Articles

GALAXY EVOLUTION IN THE MID-INFRARED GREEN VALLEY: A CASE OF THE A2199 SUPERCLUSTER

, , , , , , and

Published 2015 February 11 © 2015. The American Astronomical Society. All rights reserved.
, , Citation Gwang-Ho Lee et al 2015 ApJ 800 80 DOI 10.1088/0004-637X/800/2/80

0004-637X/800/2/80

ABSTRACT

We study the mid-infrared (MIR) properties of the galaxies in the A2199 supercluster at z = 0.03 to understand the star formation activity of galaxy groups and clusters in the supercluster environment. Using the Wide-field Infrared Survey Explorer data, we find no dependence of mass-normalized integrated star formation rates of galaxy groups/clusters on their virial masses. We classify the supercluster galaxies into three classes in the MIR color–luminosity diagram: MIR blue cloud (massive, quiescent, and mostly early-type), MIR star-forming sequence (mostly late-type), and MIR green valley galaxies. These MIR green valley galaxies are distinguishable from the optical green valley galaxies in the sense that they belong to the optical red sequence. We find that the fraction of each MIR class does not depend on the virial mass of each group/cluster. We compare the cumulative distributions of surface galaxy number density and cluster/group-centric distance for the three MIR classes. MIR green valley galaxies show the distribution between MIR blue cloud and MIR star-forming (SF) sequence galaxies. However, if we fix galaxy morphology, early- and late-type MIR green valley galaxies show different distributions. Our results suggest a possible evolutionary scenario of these galaxies: (1) late-type MIR SF sequence galaxies → (2) late-type MIR green valley galaxies → (3) early-type MIR green valley galaxies → (4) early-type MIR blue cloud galaxies. In this sequence, the star formation of galaxies is quenched before the galaxies enter the MIR green valley, and then morphological transformation occurs in the MIR green valley.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Color bimodality of galaxies appears in several wavelengths, including the optical (Strateva et al. 2001; Baldry et al. 2004; Choi et al. 2007), ultraviolet (Wyder et al. 2007; Brammer et al. 2009), and mid-infrared (MIR) bands (Johnson et al. 2007; Walker et al. 2013). In optical color–magnitude diagrams galaxies are often separated into two classes: a red sequence and a blue cloud. The evolution between the two populations is one of the key issues in the recent study of galaxy formation and evolution (Faber et al. 2007; Tinker et al. 2013).

Galaxies in a region between the red sequence and the blue cloud (i.e., green valley) seem to be an important population to understand this evolution. They might be in a transition phase from the blue cloud toward the red sequence; the transition would occur during a short timescale of less than 1 Gyr (Springel et al. 2005; Faber et al. 2007; Pan et al. 2013). Recently, Schawinski et al. (2014) suggested that the green valley galaxies are not in a single transitional state from the blue cloud toward the red sequence. They found that early- and late-type galaxies in the green valley of optical color–mass diagrams have significantly different ultraviolet–optical colors. They pointed out that this color difference could be explained by different quenching timescales of star formation; the quenching timescale for early-type galaxies (<250 Myr) is much shorter than that for late-type galaxies (>1 Gyr).

Galaxies can also be divided into several classes based on MIR colors such as AKARI N3(3.2 μm) − S11(10.4 μm) (Shim et al. 2011; Ko et al. 2012) and Spitzer/IRAC 3.5 μm–8 μm (Johnson et al. 2007; Walker et al. 2013). These MIR colors are useful indicators of specific star formation rates (SFRs; Duc et al. 2002; Chung et al. 2011) and of mean stellar ages (Piovan et al. 2003; Temi et al. 2005; Ko et al. 2009). For example, Hwang et al. (2012a, hereafter H12) studied the galaxies of the A2199 supercluster in the [3.4]–[12] color versus the 12 μm luminosity diagram using the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) data. They found that early- and late-type galaxies are well separated into two groups; late-type galaxies with MIR red colors form the "MIR star-forming (SF) sequence," but early-type galaxies mainly form the "MIR blue cloud."

As in the optical color–magnitude diagrams, galaxies with intermediate MIR colors seem to be transition populations (e.g., Shim et al. 2011; Ko et al. 2012, 2013). Recently, Alatalo et al. (2014) defined the infrared transition zone in the WISE [4.6]–[12] color space, and suggested that the galaxies in this zone are in the late stages of transition across the optical green valley. However, how this transition differs depending on galaxy morphology is not well understood. This is the goal of this study; we define "MIR green valley" galaxies in the MIR color–luminosity diagram, and study their properties in connection to other populations.

On the other hand, galaxy properties including star formation activity and galaxy morphology are strongly affected by environment (Park et al. 2007; Blanton & Moustakas 2009). Galaxy clusters are an ideal place to study this environmental dependence of galaxy properties, and several physical mechanisms in the cluster environment have been proposed to explain the dependence: ram pressure stripping (Gunn & Gott 1972), thermal evaporation (Cowie & Songaila 1977), strangulation (Larson et al. 1980), galaxy harassment (Moore et al. 1996), starvation (Bekki et al. 2002), and cumulative galaxy–galaxy hydrodynamic/gravitational interactions (Park & Hwang 2009).

Some cluster galaxies might have experienced the environmental effects even before they enter the cluster regions. For example, Zabludoff & Mulchaey (1998) suggested that massive early-type galaxies are formed by galaxy–galaxy mergers in a group-scale environment, and then fall into clusters later. The galaxy interactions and mergers should be very frequent and effective because of low velocity dispersion of galaxies in the group environment (Hickson et al. 1992; Sohn et al. 2013). This "pre-processing" in groups is supported by numerical simulations (Fujita 2004) and by observations (Koyama et al. 2011; Lemaux et al. 2012; Hess & Wilcots 2013; Mahajan 2013). Therefore, it is important to study the galaxies in galaxy groups to better understand the environment effects on galaxy properties.

Because superclusters of galaxies contain several clusters and groups, they can be excellent laboratories for studying the environmental dependence of galaxy properties in a full range of environments from cluster to field (e.g., Haines et al. 2011; Biviano et al. 2011). In this study, we use a multiwavelength catalog of the galaxies in the A2199 supercluster at z = 0.03 to study the MIR properties of galaxies and their environmental dependence. The main goal is to examine how the galaxies in the MIR green valley differ from other populations as a function of environment.

Section 2 describes the observational data we use. We derive the physical parameters of the groups and clusters in the A2199 supercluster in Section 3. We investigate the dependence of the mass-normalized integrated SFRs of galaxy groups and clusters on their virial masses in Section 4. We classify the galaxies in the MIR color–luminosity diagram, and study the environmental dependence of their MIR properties in Section 5. We discuss our results in Section 6, and conclude in Section 7. Throughout, we adopt flat ΛCDM cosmological parameters: H0 = 70  km s−1 Mpc−1, ΩΛ = 0.7, and Ωm = 0.3.

2. DATA

We used the galaxy catalog of the A2199 supercluster given by H12. H12 compiled the multi-wavelength data of the galaxies at mr < 17.77 from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7; Abazajian et al. 2009) and the WISE data. The WISE provides all-sky survey data in four MIR bands (3.4, 4.6, 12, and 22 μm) with a much better sensitivity than previous infrared surveyors. Thus, the data cover the entire region of the A2199 supercluster at z = 0.03 (12°  ×  12°  ≃  27 Mpc  ×  27 Mpc).

H12 selected member galaxies using the caustic method (Diaferio & Geller 1997; Diaferio 1999; Serra et al. 2011). They also compiled redshift data for the galaxies fainter than mr = 17.77 from other large spectroscopic surveys (Rines et al. 2002; Rines & Geller 2008). However, these surveys cover only the central region (R < 50 arcmin) of the A2199 supercluster. Therefore, we did not use these faint galaxies, but used the galaxies at mr < 17.77 to have a galaxy sample with uniform depth in the entire supercluster region. The final sample contains 1529 member galaxies, including 559 early-type (E/S0) and 970 late-type galaxies (S/Irr).

Galaxy morphology information is mainly from the Korea Institute for Advanced Study (KIAS) DR7 value-added galaxy catalog (VAGC; Choi et al. 2010). The morphology information is based on the automatic classification by Park & Choi (2005), which uses u − r color, g − i color gradient, and iband concentration index. The completeness and reliability of this method are >90%. H12 conducted an additional visual inspection to improve the classification results and to classify the galaxies not included in the KIAS DR7 VAGC. The detailed description of the galaxy catalog is in H12.

All the 1529 member galaxies are detected at the WISE 3.4 and 4.6 μm with signal-to-noise ratio (S/N) ⩾ 3. However, there are 1151 (75.3%) and 552 galaxies (36.1%) among the members detected at 12 and 22 μm with S/N ⩾ 3, respectively. We provide physical parameters of the member galaxies used in this paper. Table 1 lists SDSS-related parameters, and Table 2 lists WISE-related parameters.

Table 1. SDSS-related Physical Parameters of Member Galaxies

ID SDSS ObjID R.A. Decl. z Morpha log(Mstar/M) u − rb logSFRSDSS Dn4000
(DR7) (J2000) (J2000)
1 587725992501969046 16:46:57.61 41:56:22.7 0.032715 2 9.703 1.075 1.288 1.257
2 587725992503017882 16:53:52.93 39:49:46.7 0.033889 1 9.816 2.543 0.035 1.655
3 587725992503017689 16:53:56.32 39:48:45.4 0.033287 1 10.794 2.876 0.049 1.961
4 587725992502820896 16:52:29.10 40:18:43.0 0.029247 1 10.765 2.857 0.140 2.006
5 587725993036939572 16:36:03.51 45:44:47.6 0.030912 2 10.524 2.027 1.188 1.750

Notes. aMorphology classification: 1—early-types (E/S0); 2—late-types (S/Irr). bExtinction-corrected u − r.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

Table 2. WISE-related Physical Parameters of Member Galaxies

ID WISE ID [3.4]–[12] log(L12 μm/L) (S/N)12 μm MIRcla logSFRWISE (S/N)22 μm
1 J164657.63+415622.9 1.389 8.564 21.60 3 0.683 6.30
2 J165352.94+394946.7 −0.600 7.680 1.40 0 0.255 −1.70
3 J165356.33+394845.4 −1.879 8.106 8.40 1 0.404 3.30
4 J165229.09+401843.0 −1.934 8.096 11.40 1 0.311 1.20
5 J163603.52+454447.8 0.406 8.753 34.20 3 0.739 7.70

Note. aMIR classification: 0—(S/N)12 μm < 3, 1—MIR blue cloud galaxies, 2—MIR green valley galaxies, 3—MIR star-forming sequence galaxies.

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

3. GALAXY GROUPS AND CLUSTERS IN THE A2199 SUPERCLUSTER

The A2199 supercluster at z = 0.03 contains several galaxy groups and clusters. Rines et al. (2001, 2002) identified seven galaxy systems including three clusters (A2199, A2197W/E) and four X-ray bright groups (NRGs385, NRGs388, NRGs396, and the NGC 6159 group). They derived the physical parameters of these systems. Figure 6 displays the spatial distribution of three clusters and groups. The A2199 supercluster covers a wide range of environments: e.g., −0.85 < logΣ5(Mpc−2) < 2.85, if we define the environment by the galaxy surface number density (see Section 5.3 for details). The density range is wider than the range that we used in the study of A2255 (Shim et al. 2011).

A2199 is located in the center of the supercluster. It is a typical (rich and regular) cluster of galaxies in the nearby universe. The X-ray emission around A2199 is quite symmetric, suggesting that A2199 did not experience recent major merger (Rines et al. 2001). There is a cooling flow and a radio jet in A2199 (Markevitch et al. 1999; Johnstone et al. 2002; Kawano et al. 2003), which may be associated with a cD galaxy, NGC 6166 (Kelson et al. 2002).

Because we have recent redshift data for the supercluster, we recalculate some of the physical parameters of these systems. For example, we compute the mean redshifts, v[ km s−1] (biweight location of Beers et al. 1990), using the galaxies at R ⩽ 0.25 h−1 Mpc of each system. The radius for computing the mean redshifts is small enough so that the galaxy systems do not overlap each other, as used in Rines et al. (2002). We also compute the projected velocity dispersions, σp[ km s−1] (biweight scale of Beers et al. 1990), using the galaxies at R ⩽ 1.0 h−1 Mpc. We use this radius, smaller than the one used in Rines et al. (2002; i.e., 1.5 h−1 Mpc), to avoid the overlap between the systems. If we use R = 1.5 h−1 Mpc aperture, σp changes within the uncertainty.

Our estimates of v for the systems in the A2199 supercluster are in the range 8767–9507 km s−1; these agree well with previous measurements in Rines et al. (2002). On the other hand, our σp ranges from 244 km s−1 to 675 km s−1; these are on average ∼20% smaller than previous measurements in Rines et al. (2002). This difference seems to result mainly from the difference in the data. For example, σp of A2199 in this study is 675 ± 35 km s−1. This is smaller than $796^{+38}_{-33} {\rm \,km \,s^{-1}}$ in Rines et al. (2002), but is similar to $676^{+37}_{-32} {\rm \,km \,s^{-1}}$ in Rines & Diaferio (2006) who used the same SDSS data. Table 3 lists our measurements of v and σp for each system.

Table 3. Galaxy Groups/Clusters in the A2199 Supercluster

System R.A.2000 Decl.2000 v σp R200 M200 logLX
(J2000) (J2000) [ km s−1] [ km s−1] [degree] [Mpc] [1014M] [h−2 ergs−1]
A2199 16:28:38 39:33:05 9140 ± 135 675 ± 35 0.73 ± 0.04 1.65 ± 0.09 3.16 ± 0.48 44.1
A2197W 16:27:41 40:55:40 9482 ± 144 537 ± 22 0.58 ± 0.02 1.31 ± 0.05 1.63 ± 0.19 42.5
A2197E 16:29:43 40:49:12 8767 ±  88 543 ± 22 0.59 ± 0.02 1.32 ± 0.05 1.68 ± 0.20 42.4
NRGs385 16:17:43 34:58:00 9321 ± 100 426 ± 21 0.46 ± 0.02 1.04 ± 0.05 0.83 ± 0.12 42.8
NRGs388 16:23:01 37:55:21 9399 ± 146 549 ± 33 0.60 ± 0.04 1.34 ± 0.08 1.74 ± 0.30 42.3
NRGs396 16:36:50 44:13:00 9507 ± 110 244 ± 33 0.26 ± 0.04 0.60 ± 0.08 0.17 ± 0.07 42.0
NGC 6159 16:27:25 42:40:47 9397 ±  94 243 ± 31 0.26 ± 0.03 0.59 ± 0.08 0.16 ± 0.06 42.3

Notes. Right ascension and declination are mainly from Rines et al. (2001, 2002). In the cases of NRGs385 and the NGC 6159 group, their coordinates are taken from the NASA/IPAC Extragalactic Database. The uncertainties in v and σp are determined by the bootstrap re-sampling method.

Download table as:  ASCIITypeset image

To estimate the size of each system, we convert σp into R200 (approximately the virial radius) using the formula given in Carlberg et al. (1997). We then determine the virial mass (M200) of each system using the M200 − σp relation given in Rines et al. (2013): M200[1014M] = 0.093 ×p/200)2.90 ± 0.15. The derived values of R200 and M200 of A2199 are 1.65 ± 0.09 Mpc and (3.16 ± 0.48) × 1014M, respectively. These values agree well with the estimates in Rines & Diaferio (2006): R200 = 1.44 Mpc and M200 = (3.41 ± 1.10) × 1014M. The M200 for other systems in the A2199 supercluster are in the range 0.15–2 × 1014M, similar to the masses of galaxy groups at z < 0.1 (Rines & Diaferio 2010). Table 3 lists R200 and M200 of the groups and clusters in the A2199 supercluster.

4. INTEGRATED SFRs OF GALAXY SYSTEMS IN THE A2199 SUPERCLUSTER

To study the dependence of the star formation activity of clusters on cluster mass, many studies have used the mass-normalized integrated SFRs of clusters (i.e., ΣSFR/M200) as a proxy for the global SFRs (e.g., Finn et al. 2005; Geach et al. 2006; Bai et al. 2007, 2009; Chung et al. 2010; Koyama et al. 2010; Shim et al. 2011; Biviano et al. 2011; Chung et al. 2011). In this section, we study this dependence for the galaxy systems in the A2199 supercluster using the WISE data.

To compute the integrated SFR of each system, we use the SFRs of galaxies converted from total IR (8–1000 μm) luminosities (LIR). We compute the total IR luminosities of galaxies using the WISE 22 μm flux densities and the spectral energy distribution templates of Chary & Elbaz (2001). Following H12, we use the relation in Kennicutt (1998) for the conversion of the total IR luminosities into SFRs with the assumption of a Salpeter (1955) initial mass function with a power-law index x = 2.35: SFRWISE(M yr−1) = 1.72 × 10−10LIR(L). These IR-based SFRs (SFRWISE) of galaxies agree well with those derived from the optical spectra (H12; Hwang et al. 2012b; Lee et al. 2013).

To perform a fair comparison to other studies (e.g., Finn et al. 2005; Bai et al. 2007, 2009), we first use the same selection criteria adopted in other studies to compute the integrated SFRs (i.e., SFRWISE  >  2 M yr−1 and R ⩽ 0.5R200). We do not include active galactic nucleus (AGN)-host galaxies such as Seyferts, LINERs, Type I active galactic nuclei (AGNs) with broad Balmer lines (FWHM >1000 km s−1), and MIR-selected AGNs ([3.4]–[4.6] > 0.44, H12) because their LIR (based on the WISE 22 μm flux densities) can be dominated by AGNs rather than by star formation (Brand et al. 2009; Lee et al. 2012). However, we include composite galaxies because the AGN contribution to the 22 μm flux densities in these galaxies is small (∼6%; Lee 2012).

As listed in Table 4, we could obtain ΣSFR/M200 only for two galaxy systems: 4.78 ± 1.45 for A2199 and 1.97 ± 0.84 for A2197E in units of 10−14 yr−1. These values are similar to other clusters at similar redshifts (e.g., Coma and A1367; Bai et al. 2007). However, we could not determine ΣSFR/M200 for the other systems because there are no galaxies satisfying the selection criteria (i.e., SFRWISE > 2 M yr−1 and R ⩽ 0.5R200). It also should be noted that ΣSFR/M200 for A2199 and A2197E are based only on a small number of galaxies (i.e., three galaxies and one galaxy for A2199 and A2197E, respectively).

Table 4. ΣSFR and ΣSFR/M200 of Groups/Clusters in A2199 Supercluster

Galaxy R ⩽ 0.5R200   RR200
Ngal ΣSFR(> 2.0) ΣSFR/M200 Ngal ΣSFR(> 0.5) ΣSFR/M200 Ngal ΣSFR(> 0.5) ΣSFR/M200
System [M yr−1] [10−14 yr−1] [M yr−1] [10−14 yr−1] [M yr−1] [10−14 yr−1]
A2199 3 15.09 ± 3.98 4.78 ± 1.45 15 24.67 ± 4.16 7.81 ± 1.76   34 51.97 ± 5.56 16.44 ± 3.04
A2197W 0 ⋅⋅⋅ ⋅⋅⋅ 12 10.68 ± 1.41 6.56 ± 1.17   33 40.01 ± 3.44 24.57 ± 3.61
A2197E 1 3.32 ± 1.36 1.97 ± 0.84 12 13.37 ± 1.90 7.95 ± 1.47   37 48.82 ± 4.24 29.03 ± 4.24
NRGs385 0 ⋅⋅⋅ ⋅⋅⋅ 6 3.73 ± 0.63 4.49 ± 1.00   12 12.50 ± 1.77 15.03 ± 3.02
NRGs388 0 ⋅⋅⋅ ⋅⋅⋅ 4 5.87 ± 1.25 3.38 ± 0.93   15 16.66 ± 2.05 9.60 ± 2.05
NRGs396 0 ⋅⋅⋅ ⋅⋅⋅ 2 1.44 ± 0.42 8.74 ± 4.28   6 8.90 ± 2.22 53.83 ± 25.11
NGC 6159 0 ⋅⋅⋅ ⋅⋅⋅ 2 1.20 ± 0.35 7.33 ± 3.46   2 1.20 ± 0.35 7.33 ± 3.46

Download table as:  ASCIITypeset image

Thus, we change the selection criteria to increase the sample size: SFRWISE > 0.5 M yr−1 and R ⩽ 0.5R200. Most galaxies with SFRWISE > 0.5 M yr−1 are detected at 22 μm with S/N22 μm ⩾ 3 (94%), suggesting no significant bias due to incompleteness. This limit is much lower than the one used in a recent study based on similar WISE data (e.g., 4.6 M yr−1 in Chung et al. 2011).

Table 4 lists the ΣSFR and ΣSFR/M200 of each system based on samples of galaxies with SFRWISE > 0.5 M yr−1. We derive the values for two cases: R ⩽ 0.5R200 and RR200. Figure 1(a) shows that ΣSFR/M200 for R ⩽ 0.5R200 does not depend on M200. Similarly, ΣSFR/M200 for RR200 in panel (b) does not show M200 dependence. The Spearman's coefficient for the data in panel (a) is −0.18, while that in panel (b) is 0.07. In addition, the probability of obtaining the correlation by chance is over 70% in the both case. These results suggest that there is no significant correlation between ΣSFR/M200 and M200.

Figure 1.

Figure 1. ΣSFR/M200 for (a) R ⩽ 0.5R200 and (b) RR200 of galaxy groups/clusters as a function of M200 using galaxies with SFRWISE > 0.5 M yr−1.

Standard image High-resolution image

5. MIR PROPERTIES OF GALAXIES IN GALAXY SYSTEMS OF THE A2199 SUPERCLUSTER

To compare the star formation activity of galaxy systems in the A2199 supercluster in detail, we investigate the MIR properties of the galaxies in each system in this section. We first divide galaxies into three MIR classes in the MIR color–luminosity diagram (Section 5.1), and then compute the fraction of galaxies in each MIR class for the galaxy groups/clusters in the A2199 supercluster (Section 5.2). We also examine the environmental dependence of the MIR properties of galaxies in Section 5.3.

5.1. Galaxy Classification in the MIR Color–Luminosity Diagram

The MIR colors such as AKARI N3 − S11 or WISE [3.4]–[22], [4.6]–[12] are useful indicators of the specific SFRs and of the presence of intermediate-age stellar populations (e.g., Lee et al. 2009; Ko et al. 2009, 2012, 2013; Shim et al. 2011; Donoso et al. 2012, H12). WISE [3.4]–[12] colors are also a good indicator of the specific SFRs because the 3.4 and 12 μm luminosities, respectively, trace stellar masses and SFRs (Donoso et al. 2012; Hwang et al. 2012b). The other WISE colors including [4.6]–[12], [3.4]–[22], and [4.6]–[22] can also be used. However, we use [3.4]–[12] in this study because 3.4 μm probes the old stellar components better than 4.6 μm, and 12 μm sensitivity is better than 22 μm.

The left panels in Figure 2 show the distribution of galaxies in the [3.4]–[12] versus 12 μm luminosity diagram. As shown in H12, late-type SF galaxies form a linear "MIR SF sequence," while most early-type galaxies form a "MIR blue cloud" with low 12 μm luminosities. For robust classification of the galaxies in this diagram, we use a large sample of SDSS galaxies at 0.025 < z < 0.035 including the A2199 supercluster galaxies. The contours in Figures 2(a) and (b) indicate the number density distributions of early- and late-type galaxies, respectively. The linear fit to this large galaxy sample gives a relation,

Equation (1)

This relation is similar to the one derived in H12, but based on a much larger sample.

Figure 2.

Figure 2. (a) MIR color–luminosity diagram for (a) late-type and (b) early-type galaxies. Contours represent the number density distribution of galaxies at 0.025 < z < 0.035, while color symbols are member galaxies in the A2199 supercluster. Dashed lines are the "MIR SF sequence," which is a linear fit for late-type galaxies with SF nuclei, defined in Hwang et al. (2012a). (c) Histogram of Δ([3.4]–[12]), the deviation of [3.4]–[12] from the MIR SF sequence fit, for late-type galaxies. We use the Gaussian mixture modeling (GMM; Muratov & Gnedin 2010) to decompose the histogram with multi-Gaussian functions. We overplot three Gaussians and their sum, and list the mean and standard deviation of each Gaussian. (d) The probability distribution from the GMM that indicates which Gaussian galaxies should belong to (e) Histogram of [3.4]–[12] for early-type galaxies, and three Gaussians and the sum obtained from the GMM. Arrows represent the predicted [3.4]–[12] colors of early-type galaxies with different mean stellar ages from the SSP models (Piovan et al. 2003). (f) The probability of [3.4]–[12] belonging to each Gaussian.

Standard image High-resolution image

We plot the distribution of vertical offsets of late-type galaxies from the MIR SF sequence, Δ([3.4]–[12]), in panel (c). The histogram is negatively skewed. We use the Gaussian mixture modeling (GMM; Muratov & Gnedin 2010) to decompose the histogram with multi-Gaussian functions. The histogram is well described by the sum of three Gaussians rather than the sum of two Gaussians. We overplot the three Gaussians with their sum, and list the mean and standard deviation of each Gaussian.

Figure 2(d) shows the probability distribution from the GMM that indicates which Gaussian galaxies should belong to. At Δ([3.4]–[12]) < −1.38, the probability that galaxies belong to the left Gaussian (long-dashed line) is larger than 50%. At −1.38 < Δ([3.4]–[12]) < −0.41, the probability that galaxies belong to the middle Gaussian is over 50%. We adopt Δ([3.4]–[12]) = −1.38 to separate the galaxies in the left tail from the majority of late-type, SF galaxies.

On the other hand, the majority of early-type galaxies are distributed around [3.4]–[12] = −2, which is the MIR blue cloud. Some early-type galaxies have MIR colors that are much redder than the MIR blue cloud. Figure 2(e) shows the distribution of [3.4]–[12] colors for early-type galaxies at 0.025 < z < 0.035. The histogram is well described by the sum of three Gaussians. The Gaussian in the left covers a narrow color range centered at [3.4]–[12] ≃ −1.9 (i.e., MIR blue cloud). The Gaussians in the middle and right are for those in the green valley (to be defined in this section) and in the MIR SF sequence, respectively.

We overplot the predicted [3.4]–[12] colors from the single stellar population (SSP) models of Piovan et al. (2003). This model considers the MIR emission from dusty circumstellar envelopes of asymptotic giant branch stars when constructing the spectral energy distribution of early-type galaxies. Four arrows indicate mean stellar ages of 1, 2, 5, and 10 Gyr ([3.4]–[12] = −0.22, −0.51, −1.22, and −1.55). This suggests that the Gaussian in the left mainly traces the galaxies with stellar populations older than 10 Gyr, while the Gaussian in the middle traces the galaxies with intermediate-age stellar populations.

To separate the early-type galaxies with intermediate and old stellar populations, we adopt a cut of [3.4]–[12] = −1.55 that corresponds to the color of 10 Gyr stellar populations. This color cut is the one where the probability that the galaxies belong to the Gaussian in the left becomes larger than 50%, as shown in panel (f).

Based on the results in Figure 2, we divide the galaxies into three MIR classes in the MIR color–luminosity diagram as shown in Figure 3. We use two demarcation lines, Δ([3.4]–[12]) = −1.38 from the MIR SF sequence and [3.4]–[12] = −1.55. We call the galaxies at Δ([3.4]–[12]) > −1.38 from the MIR SF sequence "MIR SF sequence galaxies," and those with [3.4]–[12] < −1.55 "MIR blue cloud galaxies." Then the galaxies between the two lines are "MIR green valley galaxies."

Figure 3.

Figure 3. MIR color–luminosity diagram for galaxies with S/N12 μm ⩾ 3. Contours represent the number density distribution of galaxies at 0.025 < z < 0.035, while circles and crosses are early- and late-type member galaxies of the A2199 supercluster, respectively. We divide galaxies into three classes: MIR SF sequence galaxies, MIR green valley galaxies, and MIR blue cloud galaxies. Long- and short-dashed lines are demarcation lines for this MIR galaxy classification. The latter corresponds to the color of 10 Gyr stellar populations (Piovan et al. 2003). A dotted line indicates a linear fit for the MIR SF sequence.

Standard image High-resolution image

Among 997 member galaxies in this diagram, the number of MIR SF sequence, MIR green valley, and MIR blue cloud galaxies are 678 (68.0%), 126 (12.6%), and 193 (19.4%), respectively. The majority (93%, 628/678) of MIR SF sequence galaxies are morphologically late-types. In contrast, MIR blue cloud galaxies are mainly (90%, 174/193) early-types. Two-thirds of MIR green valley galaxies are early-type galaxies (68%, 86/126), and one-third of them are late-types (40/126).

To compare the MIR classification with that in the optical bands, we plot the galaxies in each MIR class in the optical color–mass diagram. Figure 4 shows an extinction-corrected u − r color as a function of stellar mass. To correct for dust extinction, we use E(BV) values from the stellar continuum fits of Oh et al. (2011) and apply the reddening law of Cardelli et al. (1989). Grayscale contours indicate the number density of the SDSS galaxies at 0.025 < z < 0.035. Two dashed lines define the optical green valley adopted from Schawinski et al. (2014).

Figure 4.

Figure 4. Extinction-corrected (ur) color vs. stellar mass diagram. Gray contours represent the number density distribution of galaxies at 0.025 < z < 0.035. Circles and crosses represent early- and late-type A2199 supercluster members segregated by their MIR classes: (a) MIR blue cloud galaxies, (b) MIR green valley galaxies, (c) MIR SF sequence galaxies, and (d) 12 μm undetected galaxies (S/N12 μm < 3). Dashed lines indicate the optical green valley defined in Schawinski et al. (2014).

Standard image High-resolution image

As expected, MIR SF sequence galaxies are mainly distributed in the optical blue cloud (see panel c). Some of them are in the optical green valley. MIR blue cloud galaxies are distributed in the optical red sequence (panel a). Interestingly, the majority of MIR green valley galaxies roughly follow the optical red sequence despite a scatter in color (panel b). These are consistent with the results in other studies (Ko et al. 2012; Walker et al. 2013).

The stellar mass range of MIR green valley galaxies is nearly the same as that of MIR blue cloud galaxies (log(Mstar/M) > 10). This means that the discrimination between MIR blue cloud and MIR green valley galaxies based only on optical parameters is difficult. It is interesting to note that the MIR green valley galaxies are not the same as the optical green valley galaxies. The optical green valley galaxies are included in the MIR SF sequence class. Panels (a)–(c) also show that the distributions of early- and late-type galaxies in each panel are not significantly different.

In panel (d) we plot the distribution of 12 μm undetected galaxies (i.e., S/N12 μm < 3); they are not included in the MIR classification. They are distributed either in the optical blue cloud or in the red sequence. They appear to have stellar masses lower than other MIR classes in a fixed morphology. The majority (∼80%) of the 12 μm undetected galaxies have no strong emission lines in their optical spectra; most of their spectral types are undetermined. These results suggest that they are a mixture of young galaxies with small stellar masses and of relatively massive galaxies that are quiescent; both of them have 12 μm flux densities, not high enough to be detected with WISE.

We also show the galaxies in each MIR class in the plot of SFRs derived in the optical band and stellar mass (see top panels of Figure 5). We use SFRSDSS adopted from the MPA/JHU DR7 value-added galaxy catalog (Brinchmann et al. 2004). The grayscale number density contours based on the galaxies at 0.025 < z < 0.035 show two distinctive sequences. The upper sequence (solid lines) is for SF galaxies, which is called a "main sequence" of SF galaxies (Noeske et al. 2007; Elbaz et al. 2011). The lower sequence is for those with SFRs much smaller than the main sequence galaxies.

Figure 5.

Figure 5. Top: stellar mass vs. SFRSDSS diagram. Contours represent the number density distribution of galaxies at 0.025 < z < 0.035. Circles and crosses represent early- and late-type member galaxies of the A2199 supercluster segregated by their MIR classes: (a) MIR SF sequence galaxies, (b) MIR green valley galaxies, and (c) MIR blue cloud galaxies. Solid lines are the linear fit to galaxies with SF nuclei, while dashed lines indicate 1/5 and 1/20 levels, respectively. Bottom: stellar mass vs. Dn4000 distribution of (d) MIR SF sequence, (e) MIR green valley, and (f) MIR blue cloud galaxies.

Standard image High-resolution image

As expected, MIR SF sequence galaxies are well distributed along the main sequence (panel a). In contrast, MIR blue cloud galaxies are on the sequence in the bottom. Similarly, MIR green valley galaxies are mainly found on the sequence in the bottom. The SFRs of MIR green valley galaxies appear to be slightly higher than those for MIR blue cloud galaxies in a give stellar mass. However, the difference is not significant.

We also examine the distribution of Dn4000 for the three MIR classes in bottom panels. The Dn4000 is a useful measure of the mean stellar age of galaxies (Balogh et al. 1999; Shim et al. 2011). Galaxies with Dn4000 < 1.5 contain young stellar populations with ≲ 1 Gyr (Kauffmann et al. 2003). On the other hand, Dn4000 larger than 1.5 indicates mean stellar age larger than 1 Gyr (Kranz et al. 2010). As expected, most MIR SF sequence galaxies have Dn4000 smaller than 1.5. In contrast, Dn4000 of the majority of the MIR green valley and MIR blue cloud galaxies is larger than 1.5. The mean Dn4000 of the MIR green valley galaxies (1.84 ± 0.03) is not significantly different from that of the MIR blue cloud galaxies (1.90  ±  0.01). However, the rms of Dn4000 for the MIR green valley galaxies is larger (0.10 ± 0.01) than that for the MIR blue cloud galaxies (0.07 ± 0.01), suggesting that the dispersion of mean stellar ages for the MIR green valley galaxies is larger than that for the MIR blue cloud galaxies. However, the difference in mean stellar ages between the two populations is more evident in the MIR color space.

5.2. Galaxies in Groups/Clusters and Their MIR Classes

In this section, we examine the spatial distribution of three MIR classes in the A2199 supercluster. We compute the fraction of each MIR class in galaxy groups/clusters, and examine the variation of the fraction with the mass of each system.

We first show the spatial distribution of galaxies in each MIR class over the entire region of the A2199 supercluster in Figure 6. We mark the positions of three clusters (A2199, A2197W/E) and four X-ray bright groups (NRGs385, NRGs388, NRGs396, and the NGC 6159 group) as thick circles. The size of the circles is proportional to their virial radii (R200). We also mark the positions of four X-ray faint groups (NRGs389, NRGs395, NRGs399, and NRGs400; Rines et al. 2001) as triangles. To better show the clustering of galaxies, we overplot the grayscale galaxy number density map using all the member galaxies regardless of their MIR classes in each panel.

Figure 6.

Figure 6. Spatial distribution of (a) the member galaxies undetected at 12 μm (S/N12 μm < 3), (b) MIR blue cloud galaxies, (c) MIR green valley galaxies, and (d) MIR SF sequence galaxies in the A2199 supercluster. Green and white represent the galaxy number density distribution. Large circles represent the positions of X-ray bright clusters and groups, and their sizes are proportional to virial radii (R200). Triangles represent the positions of four X-ray faint groups. Small circles and crosses represent early- and late-type galaxies, respectively.

Standard image High-resolution image

The 12 μm undetected galaxies (panel a) are mainly found in high-density regions. Early-type galaxies are highly concentrated in central regions of groups/clusters. Note that these early-type galaxies are in the optical red sequence with low stellar masses. Similarly, MIR blue cloud galaxies (panel b) and MIR green valley galaxies (panel c) are mainly found in groups and clusters. In contrast, MIR SF sequence galaxies are distributed over the entire supercluster region. Some of them are in groups and clusters, but they do not show a central concentration like MIR blue cloud and MIR green valley galaxies. Four X-ray faint groups do not show any significant clustering of galaxies.

We then calculate the fraction of each MIR class in the X-ray bright groups and clusters. We use the galaxies in the inner region (R ⩽ 0.5R200) and in the outer region (0.5R200 < RR200) separately. For comparison, we also calculate the fraction of each MIR class for the galaxies not associated with any group or cluster (i.e., R > 3 × R200) in the supercluster region; these are the galaxies in "underdense regions." We use only 12 μm detected galaxies for computing fractions. If we use all the galaxies to compute the fractions by including the 12 μm undetected galaxies, the results do not change much.

The left panels of Figure 7 show the fraction of each MIR class in the inner region as a function of M200 of each galaxy system. The fractions of MIR blue cloud galaxies (panel a) and of MIR SF sequence galaxies (panel c) vary from ∼20% to ∼60%, but do not show any dependence on M200. The fractions of MIR green valley galaxies in the galaxy systems are similar (∼25%).

Figure 7.

Figure 7. Fraction of (top) MIR blue sequence, (middle) MIR green valley, and (bottom) MIR SF sequence galaxies in the groups/clusters as a function of M200. We use galaxies in the inner region (R ⩽ 0.5R200, left) and the outer region (0.5R200 < RR200, right), separately, of groups/clusters. The symbols are the same as those in Figure 1. Hatched regions represent the fraction of each MIR class in the underdense region.

Standard image High-resolution image

Similarly, the right panels for the galaxies in the outer region of groups/clusters show no dependence of the fraction on M200. One unique point is that the fraction of each MIR class is nearly the same for all groups and clusters. Although we calculate the fraction using the galaxies at RR200 of each group and cluster (i.e., sum of the left and right panels), the M200 dependence of the fraction is not apparent.

On the other hand, panel (a) shows that the fractions of MIR blue cloud galaxies in the inner region of groups/clusters are significantly higher than in the underdense region (11.1%). However, the fractions in the outer region (panel d) are on average ∼ 18.5%, close to the fraction in the underdense region. The fraction of MIR green valley galaxies in the underdense region is 10.9%. This value is smaller than those in the inner region (panel b), but is comparable to those in the outer region (panel e). The fractions of MIR SF sequence galaxies in the inner region (panel c) are much smaller than that in the underdense region (78%). However, the fractions in the outer region are close to the value in the underdense region (panel f).

Thus Figure 7 suggests that massive galaxies with low SFRs (i.e., MIR blue cloud galaxies and MIR green valley galaxies) are mainly in the inner regions of groups/clusters regardless of their M200. Some of them are found even in the low density regions (i.e., the underdense region), but their fraction is not significant (up to ∼20%).

5.3. Environmental Dependence Based on Σ5 and R/R200

To further examine the environmental dependence of the MIR classes in the A2199 supercluster, we use two environmental indicators in this section; the surface number density (Σ5) and the cluster-(or group-)centric distance normalized by virial radius (R/R200). The surface number density, Σ5, is defined by $\Sigma _5=5(\pi D^2_{p,5})^{-1}$, where Dp, 5 is the projected distance to the fifth-nearest neighbor. The fifth-nearest neighbor to each galaxy is identified in the sample of A2199 supercluster member galaxies. We compare the cumulative distribution of each MIR class as a function of two environmental indicators. We also use the Kolmogorov–Smirnov two-sample test to quantify the difference in the distribution, if any.

Figure 8 shows the cumulative distribution of each MIR class as a function of Σ5. Panel (a) shows the case for all of the 12 μm detected galaxies. MIR SF sequence galaxies (red long-dashed line) tend to be in low-density regions, and MIR blue cloud galaxies (blue solid line) prefer high-density regions. MIR green valley galaxies (green dashed line) are between the two. This result is consistent with the expectations from the morphology- and/or SFR-density relation (e.g., Dressler 1980; Lewis et al. 2002; Park & Hwang 2009; Hwang et al. 2010).

Figure 8.

Figure 8. Cumulative distribution function of the three MIR classes as a function of the surface galaxy number density (Σ5) for MIR blue cloud galaxies (solid line), MIR green valley galaxies (short-dashed line), and MIR SF sequence galaxies (long-dashed line). The galaxy sample used in each panel is (a) all 12 μm detected galaxies, (b) early-type galaxies, (c) late-type galaxies, (d) massive galaxies greater than log(Mstar/M) > 10.5, (e) massive early-type galaxies, and (f) massive late-type galaxies. We list the number of three MIR class galaxies used in each panel.

Standard image High-resolution image

Panel (b) shows the same Σ5 distribution of each MIR class, but only for early-type galaxies. It shows that the difference between MIR blue cloud and MIR green valley galaxies is noticeably reduced. The confidence level to reject the hypothesis that the two distributions are extracted from the same parent population changes from 92% in panel (a) into 52% in panel (b). This suggests that the Σ5 values of MIR green valley early-type galaxies do not differ significantly from those of MIR blue cloud early-type galaxies. On the other hand, the cumulative distribution of MIR SF sequence early-type galaxies is still distinguishable from those of MIR blue cloud and MIR green valley early-type galaxies.

Figure 8(c) shows the cumulative distribution only for late-type galaxies. The confidence level to reject the hypothesis that the two distributions are extracted from the same parent population is 93% between MIR green valley and MIR blue cloud late-type galaxies, and is 97% between MIR green valley and MIR SF sequence late-type galaxies.

Because galaxy properties also strongly depend on stellar mass as well as on environment (Tasca et al. 2009; Hwang et al. 2010; Cucciati et al. 2010; Bolzonella et al. 2010; Peng et al. 2010; Sobral et al. 2011), we restrict our analysis to only massive galaxies with log(Mstar/M) > 10.5 to reduce the mass effect on the MIR properties. Figure 8(d) shows the cumulative distribution for the massive galaxies with 12 μm detection. Although the number of galaxies in each MIR class is reduced, the result is similar to that in panel (a). Panel (e) shows the case of massive early-type galaxies. MIR blue cloud and MIR green valley galaxies again show similar distributions. However, the cumulative distribution of MIR green valley galaxies among the massive late-type galaxies (panel f) is still between MIR blue cloud and MIR SF sequence galaxies.

We perform a similar analysis in Figure 9 using the normalized cluster/group-centric distance R/R200 instead of Σ5. Panel (a) shows the cumulative distribution of each MIR class using all of the 12 μm detected galaxies. The R/R200 distribution of MIR SF sequence galaxies is significantly different from those of MIR blue cloud and MIR green valley galaxies; the hypothesis that the two distributions are extracted from the same parent population can be rejected at a confidence level of >99%. However, the distributions of MIR blue cloud and MIR green valley galaxies are very similar. They show similar distributions even if we consider only early-type galaxies in panel (b), massive galaxies in panel (d), and massive early-type galaxies in panel (e).

Figure 9.

Figure 9. Same as Figure 8, but the cumulative distribution as a function of cluster/group-centric distance normalized by virial radius (R/R200).

Standard image High-resolution image

On the other hand, when we consider late-type galaxies, the similar distribution between MIR blue cloud and MIR green valley galaxies disappears. The confidence level to reject the hypothesis that the two distributions are extracted from the same parent population is 92% in panel (c) and 73% in panel (f). The cumulative distribution of MIR green valley late-type galaxies appears to be similar to that of MIR SF sequence late-type galaxies, and still to be between the MIR blue cloud and the MIR SF sequence galaxies (panels c and f).

Both Figures 8 and 9 show similar results. The cumulative distributions of Σ5 and R/R200 for MIR green valley galaxies are between those for MIR blue cloud galaxies and for MIR SF sequence galaxies. However, among early-type galaxies the difference between MIR blue cloud and MIR green valley galaxies disappears. In contrast, when considering only late-type galaxies, the MIR green valley galaxies are still between the MIR blue cloud and the MIR SF sequence galaxies.

6. DISCUSSION

6.1. Star Formation Activity of Galaxy Systems and Their Dependence on Virial Mass

We compare the star formation activity of X-ray bright galaxy groups and clusters in the A2199 supercluster using ΣSFR/M200 (Section 4) and the fraction of MIR SF sequence galaxies (Section 5.2). Interestingly, both ΣSFR/M200 and the fraction do not seem to depend on virial masses of galaxy systems (see Figures 1 and 7). This is seen both in inner (R ⩽ 0.5R200) and outer regions (0.5R200 < RR200).

Our results are consistent with those in other studies. For example, Goto (2005) used the SDSS data for 115 nearby clusters, and found no dependence of ΣSFR/M200 and of blue/late-type galaxy fraction on cluster virial mass (see also Bai et al. 2009, 2010; Balogh & McGee 2010). A recent study based on WISE data for 69 clusters at z < 0.1 also find no correlation between ΣSFR/M200 and M200 (Chung et al. 2011).

No correlation between the fraction of passive (or red) galaxies and the halo masses of groups/clusters could result from a simultaneous increase in the number of member galaxies and of massive quiescent elliptical galaxies with increasing M200 of the systems (De Lucia et al. 2012). No correlation between the two can also be explained by the pre-processing mechanism (Zabludoff & Mulchaey 1998), which is that giant elliptical galaxies or cD galaxies are formed by galaxy–galaxy mergers in poor groups, and then fall into clusters later. This pre-processing occurs in poor groups with M200  <  1013M. This results in no M200 dependence of the star formation activity in groups/clusters when the poor groups fall into more massive clusters with M200 > 1013M (Balogh & McGee 2010); this covers the mass range of the galaxy systems in the A2199 supercluster.

On the other hand, some studies show that there is mass dependence of the star formation activity of galaxy groups/clusters. For example, Finn et al. (2005) and Koyama et al. (2010) found an anti-correlation between ΣSFR/M200 and M200 for clusters at z < 1 (see also Bai et al. 2007 and Poggianti et al. 2008). This anti-correlation is seen in clusters with 1013M < M200 < 3 × 1014M, which covers our galaxy systems in the A2199 supercluster. Furthermore, Popesso et al. (2012) used far-infrared Herschel data for nine clusters (M200 ⩾ 3 × 1014M) and nine groups/poor clusters (1013M < M200 < 3 × 1014M) at z < 1.5, and found that ΣSFR/M200 is higher for groups/poor clusters than for massive clusters at a given redshift. They suggested that this mass-dependent star formation activity results from earlier quenching of the star formation activity in galaxies associated with more-massive halos (i.e., clusters).

Similarly, some studies show that the fraction of SF galaxies decreases with increasing cluster mass (e.g., Martínez et al. 2002; Rasmussen et al. 2012). Poggianti et al. (2006) used the clusters at low z (0.04 < z < 0.08) and at high z (0.4 < z < 0.8), and found an anti-correlation between the fraction of SF galaxies and velocity dispersion (σp) of the clusters; the anti-correlation is significant for low-z clusters with σp < 550 km s−1 (M200 ≲ 2 × 1014M). Most galaxy systems in the A2199 supercluster have σp < 550 km s−1. However, the fraction of SF galaxies (i.e., MIR SF sequence galaxy fraction) in galaxy groups/clusters of the A2199 supercluster does not show any strong dependence on M200.

It seems that the mass dependence of star formation activity of groups/clusters is still inconclusive. Some previous studies used heterogeneous data with different wavelength, sky coverage, redshift, and sensitivity (see the Appendix of Bai et al. 2007). Finn et al. (2005) and Koyama et al. (2010) suggested that the mass dependence of ΣSFR/M200 appears to be due to the strong redshift dependence (high-z clusters have higher ΣSFR/M200 and smaller M200 in their data). In addition, Poggianti et al. (2006) showed that their mass dependence of the fraction of SF galaxies becomes unclear when they used another cluster sample from Miller et al. (2005). Thus, more effort is still required to examine whether or not cluster mass is a critical factor in determining their star formation activity.

If the star formation activity in groups/clusters indeed depends on M200, no M200 dependence for the galaxy systems in the A2199 supercluster indicates that there could be other mechanisms that have affected the star formation activity of the systems. One possible mechanism is the interaction between infall groups and the supercluster. If the star formation activity of group galaxies has gradually decreased through the gravitational and/or hydrodynamical interaction with the supercluster, the current star formation activity of the groups would depend on when the groups fell into the supercluster environment. Thus, the different infall time of the groups in the supercluster could explain the different star formation activity of groups at a given M200.

Similarly, the interaction between groups/clusters can change their star formation activity (Bekki 1999; Owen et al. 2005; Miller 2005; Ferrari et al. 2005; Johnston-Hollitt et al. 2008; Hwang & Lee 2009). In fact, there are three clusters close to the central region of the A2199 supercluster: A2199, A2197E, and A2197W. In particular, A2197E and A2197W overlap each other both on the sky and in redshift space; two clusters may be currently interacting. A2199 also seems to interact with A2197W/E and with another nearby group, NRGs388. Thus, the star formation activity of A2199, A2197W/E, and NRGs388 could be enhanced by the interaction among them. Note that the merging/interaction stages are also important to determine the star formation activity of galaxy groups/clusters (Hwang & Lee 2009). Therefore, the different merging/interaction stages of the groups/clusters in the supercluster can add significant scatter to the anti-correlation between the star formation activity and M200, and can result in no correlation between the two.

6.2. Morphology Dependence of the Galaxies in the MIR Green Valley

In Section 5.3 we compare the cumulative distribution of Σ5 and R/R200 for different MIR classes. The different distributions between MIR blue cloud galaxies and MIR SF sequence galaxies are consistent with the expectations from the morphology–density and SFR–density relations; most MIR blue cloud galaxies are morphologically early-types with low SFRs, but MIR SF sequence galaxies are generally late-type galaxies with high SFRs (see Figure 5). Interestingly, the distribution of MIR green valley galaxies are between MIR blue cloud and MIR SF sequence galaxies in both environmental indicators. Because the environment strongly affects galaxy morphology and star formation activity (Park et al. 2007; Blanton & Moustakas 2009), the MIR green valley galaxies between MIR blue cloud and MIR SF sequence galaxies seem to be in a transition stage.

However, when we fix galaxy morphology, the behavior of the MIR green valley galaxies changes. The MIR green valley galaxies with early-type morphology show a distribution very similar to that of the early-type MIR blue cloud galaxies. However, the MIR green valley galaxies with late-type morphology still show a distribution between MIR blue cloud and MIR SF sequence galaxies. This trend is also apparent even if we use only massive galaxies with log(Mstar/M) > 10.5. The similarity between MIR green valley and the MIR blue cloud galaxies among the early-types is also seen in their stellar masses: log(Mstar/M) = 10.66 ± 0.04 for MIR green valley early-types and log(Mstar/M) = 10.75 ± 0.03 for MIR blue cloud early types.

To better understand the physical meaning of the observed galaxy colors, we show NUV−r and u − r colors of early-type galaxies as a function of WISE [3.4]–[12] colors in Figure 10 (see also Ko et al. 2013). We model that the MIR emission among early-types (i.e., MIR green valley early-type galaxies) can be explained by the presence of a small fraction of young stellar populations. Panel (a) shows the NUV−r versus [3.4]–[12] color–color distribution of early-type member galaxies and the SWIRE templates of Polletta et al. (2007). We overplot the two-component SSP models of Piovan et al. (2003) with various combinations of an old stellar population (12 Gyr) and young populations (0.5, 1, and 2 Gyr) with solar metallicity. The model grids indicate different fractions of a young population from 0.1% to 50%.

Figure 10.

Figure 10. (a) NUV−r vs. [3.4]–[12] and (b) u − r vs. [3.4]–[12] color–color distributions for early-type member galaxies (gray dots). Solid lines show the color prediction from the two-component SSP model (Piovan et al. 2003). These colors are obtained using several combinations of an old stellar population (12 Gyr) and young populations (0.5, 1, and 2 Gyr) with varying fraction (0.1, 1, 5, 10, 30, and 50%) of young populations. Crosses represent the SWIRE templates (Polletta et al. 2007). The vertical dashed line represents [3.4]–[12] = −1.55 that is the border line between MIR blue cloud and MIR green valley galaxies, while the horizontal dashed line is the NUV excess cut, NUV−r = 5.4, used in Ko et al. (2013).

Standard image High-resolution image

The plot shows that only >1% of the 0.5 Gyr or >5% of the 1 Gyr population can make [3.4]–[12] colors green (> − 1.55). The plot also shows that about 50% of MIR green valley early-type galaxies have no NUV excess emission (NUV−r > 5.4). Ko et al. (2013) suggested that young stellar populations with 0.5 and 1 Gyr can make the galaxies have both MIR and NUV excess, but a 2 Gyr population does not produce any NUV excess. This means that half of MIR green valley early-type galaxies stopped forming stars at least 2 Gyr ago. This >2 Gyr old recent star formation is detected only in the MIR (neither in the NUV nor in the optical) because the MIR traces star formation over a longer timescale (>2 Gyr) than the NUV and optical bands (Ko et al. 2013; Schawinski et al. 2014).

Figure 10(b) displays the u − r versus [3.4]–[12] distribution of early-type member galaxies. The model grids show that u − r colors change very little when increasing the fraction of young populations. This result is consistent with the results in Figures 4 and 5; the MIR green valley and MIR blue cloud early-type galaxies show similar optical properties.

The environmental dependence of these MIR classes and their optical properties shown in Section 5 suggest a possible evolutionary scenario of galaxies from late-type SF galaxies to early-type quiescent galaxies: (1) late-type MIR SF sequence galaxies → (2) late-type MIR green valley galaxies → (3) early-type MIR green valley galaxies → (4) early-type MIR blue cloud galaxies. Because the timescale for the suppression of star formation activity is shorter than the timescale for morphological transformation (Blanton & Moustakas 2009), we expect that the star formation activity of supercluster galaxies decreases first (1 → 2), and then the morphological transformation would be followed (2 → 3). Skibba et al. (2009) found that the color–density relation is more fundamental than the morphology–density relation. Bamford et al. (2009) also obtained a similar result, and suggested that the color transformation driven by environment should take place on shorter timescales than the morphological transformation. This agrees well with our scenario.

Masters et al. (2010) studied the environmental dependence of (optically) red (or passive) spiral galaxies that are considered as transition objects between normal blue spirals and red early-type galaxies. They showed that red spirals are mainly found in intermediate density regions and have higher stellar masses than blue spirals. These characteristics of red spirals are very similar to those of late-type MIR green valley galaxies (see Figures 58, and 9). They also suggested that strangulation or starvation (removal of gas in the outer halo and no further accretion of cold gas) is a plausible mechanism for producing red spirals. These mechanisms can explain the quenching of star formation even in low-mass groups (Kawata & Mulchaey 2008; Wetzel et al. 2013), which accounts for the presence of late-type MIR green valley galaxies located outside the cluster regions (see Figure 6). However, low-mass spirals cannot keep their spiral structures from the environmental processes (Bekki et al. 2002). This explains why late-type MIR green valley galaxies (red spirals) have higher stellar masses than late-type MIR SF sequence galaxies (normal blue spirals).

When galaxies enter the MIR green valley, they undergo morphological transformation from late-types to early-types. Early-type MIR green valley galaxies are more likely to exist in high-density regions than late-type MIR green valley galaxies (see Section 5.3). This suggests that the environment can play a role in the morphological transformation in the MIR green valley. George et al. (2013) suggested that the observed morphological dependence of galaxies on environment cannot be fully explained by strangulation and disk fading (Quilis et al. 2000). Instead, galaxy mergers and close tidal interactions should be required for the morphological transformation. Park & Hwang (2009) also found that cumulative galaxy–galaxy gravitational/hydrodynamic interactions are the main driver for the morphological transformation in the cluster environment. Bekki & Couch (2011) showed that the transformation in group environment can be driven by repetitive slow encounters with group members. These processes seem suitable for MIR green valley galaxies that are mainly in groups and clusters.

To summarize, the star formation of galaxies is quenched before the galaxies enter the MIR green valley, which is mainly driven by strangulation or starvation. Then, the morphological transformation from late-types to early-types occurs in the MIR green valley, especially for massive late-type galaxies with log(Mstar/M) > 10. The main environmental mechanisms for the morphology transformation are galaxy–galaxy mergers and interactions. After the transformation, early-type MIR green valley galaxies keep the memory of their last star formation for several Gyr.

7. SUMMARY AND CONCLUSIONS

Using the multi-wavelength data covering the entire region of the A2199 supercluster, we study the star formation activity of galaxy groups/clusters in the supercluster and the MIR properties of the supercluster galaxies. We determine ΣSFR/M200 of groups and clusters using WISE 22 μm data, and find no dependence of ΣSFR/M200 on M200.

We classify galaxies into three MIR groups in the WISE [3.4]–[12] color versus 12 μm luminosity diagram: MIR blue cloud galaxies, MIR SF sequence galaxies, and MIR green valley galaxies. MIR blue cloud galaxies are mostly (90%) early-types, while MIR SF sequence galaxies are predominantly (93%) late-types. MIR green valley galaxies consist of early- (68%) and late-type galaxies (32%). SFRs for MIR blue cloud and MIR green valley galaxies are similar, but are much smaller than those for MIR SF sequence galaxies. It is important that the MIR green valley galaxies are distinguishable from the optical green valley galaxies, in the sense that they belong to the optical red sequence. Thus, star formation quenching of galaxies occurs before the galaxies enter the MIR green valley.

We calculate the fraction of each MIR class in groups/clusters. We could not find any dependence of the fraction of each MIR class on M200, consistent with the trend of ΣSFR/M200. These results suggest that group/cluster mass does not play an important role in controlling the global star formation activity of the systems. On the other hand, there could be other mechanisms affecting the global star formation activity of the systems. These include the interaction between infalling groups and the supercluster, and between groups/clusters.

We compare the cumulative distribution of Σ5 and R/R200 for the three MIR classes. MIR green valley galaxies show the distribution between MIR blue cloud and MIR SF sequence galaxies. When only considering early-type galaxies, the difference between MIR blue cloud and MIR green valley galaxies disappears. However, late-type galaxies in the MIR green valley still show the distribution between MIR blue cloud and MIR green valley galaxies. These results suggest a possible evolutionary scenario: (1) late-type MIR SF sequence galaxies → (2) late-type MIR green valley galaxies → (3) early-type MIR green valley galaxies → (4) early-type MIR blue cloud galaxies. Thus, the MIR green valley is the site where the morphology transformation of galaxies mainly appears to occur.

We thank Margaret Geller for many helpful discussions. G.H.L. acknowledges the support by the National Research Foundation of Korea (NRF) Grant funded by the Korean Government (NRF-2012-Fostering Core Leaders of the Future Basic Science Program). M.G.L was supported by the NRF grant funded by the Korea Government (MEST; No. 2012R1A4A1028713). H.S.H. acknowledges the Smithsonian Institution for the support of his post-doctoral fellowship. J.S. was supported by the Global Ph.D. Fellowship Program through an NRF funded by the MEST (No.2011- 0007215). H.S. is supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science, ICT & Future Planning (2014R1A1A1038088). A.D. acknowledges partial support from the INFN grant In-Dark and from the grant Progetti di Ateneo/CSP TO Call2 2012 0011 "Marco Polo" of the University of Torino. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, Cambridge University, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.

Please wait… references are loading.
10.1088/0004-637X/800/2/80