The Evolution of the Quenching of Star Formation in Cluster Galaxies since z ∼ 1

I. Pintos-Castro; H. K. C. Yee; A. Muzzin; L. Old; G. Wilson

doi:10.3847/1538-4357/ab14ee

1. Introduction

Observationally, galaxies can be classified into two broad types: red passive galaxies dominated by early-type morphologies and blue star-forming (SF) galaxies dominated by late-type morphologies (e.g., Strateva et al. 2001; Baldry et al. 2004). The probability of a galaxy belonging to one of these types is mainly driven by its stellar mass and its surrounding environment (e.g., Dressler 1980; Balogh et al. 2004; Kauffmann et al. 2004; Baldry et al. 2006): on average, more massive galaxies are more likely to become quiescent, and galaxies in denser regions are more likely to become quenched. Another important parameter is the redshift: while local clusters are populated by red and dead massive galaxies, the fraction of blue, SF galaxies is higher in more distant clusters (e.g., Butcher & Oemler 1984; Li et al. 2009; Raichoor & Andreon 2012; Hennig et al. 2017), and there is even evidence for high star formation activity close to the cluster core of the highest-redshift galaxy clusters (e.g., Hilton et al. 2010; Santos et al. 2015; Webb et al. 2015; Wang et al. 2016).

Several studies have focused on disentangling the separate contributions of processes driven by intrinsic galaxy properties and the environmental influence. It is the "nature versus nurture" scenario of the debate (Dressler 1980). Over the past decade, a number of investigations have classified mechanisms of suppression of star formation activity into two major categories known as "mass quenching," related to internal processes that scale with stellar mass, and "environmental quenching," related to external changes as galaxies interact with their surroundings (Peng et al. 2010). Whether these two quenching processes act independently from each other is still under discussion. While observational studies show evidence that the two effects are separable in the local universe or for relatively high-mass galaxies at high redshifts (Peng et al. 2010; Muzzin et al. 2012; Quadri et al. 2012; Kovač et al. 2014), there is growing observational evidence that environmental quenching evolves and depends on stellar mass (Balogh et al. 2016; Darvish et al. 2016; Kawinwanichakij et al.2017; Lin et al. 2017).

Some of the apparent discrepancies when studying "environment" may arise from its many different definitions. One approach is to analyze galaxy properties in terms of the large-scale structure in which galaxies are located (e.g., superclusters, clusters, groups, field, voids). In particular, a commonly used proxy is the projected distance to the center of the parent halo to find at which radius galaxy properties are effectively altered (e.g., Wetzel et al. 2014; Allen et al. 2016; Vulcani et al. 2016; Bufanda et al. 2017). Another useful approach, which requires spectroscopic data to build up a picture of the phase-space profiles of clusters, is to identify whether galaxies are infalling or virialized to determine how long they have been under the influence of the host halo (Muzzin et al. 2014; Jaffé et al. 2015, 2018; Noble et al. 2016; Weinzirl et al. 2017). A common approach is to study galaxy properties in terms of the local galaxy density around each individual galaxy: in this case no previous knowledge of clusters/groups positions is needed (e.g., Peng et al. 2010; Li et al. 2012; Davies et al. 2016; Kawinwanichakij et al. 2017; Darvish et al. 2018). There is no universal method to measure the environment, and its suitability depends on the physical scale being probed (see, e.g., Muldrew et al. 2012, for a detailed discussion on galaxy environment measures).

Physical explanations for the suppression of star formation in galaxies related to mass include feedback from an active galactic nucleus (e.g., Fabian 2012) or energetic feedback from supernova explosions and stellar winds (e.g., Dekel & Silk 1986; Veilleux et al. 2005). Regarding the environmental influence, many physical processes have been proposed that might transform the morphology and star formation properties of galaxies in dense environments, acting at different timescales and physical distances from the center of the host halo (see, e.g., Boselli & Gavazzi 2006, for a review). Such mechanisms can be broadly classified into (1) the interaction of a galaxy with the gaseous component of the cluster (e.g., ram pressure stripping; Gunn & Gott 1972), (2) interactions between the galaxy and the cluster gravitational potential (e.g., tidal stripping; Merritt 1984), and (3) smaller-scale galaxy–galaxy interactions (e.g., harassment; Moore et al. 1996). The relative importance of quenching mechanisms also depends on the epoch. For example, at high redshift, in the absence of cosmological accretion, once a galaxy falls into a larger halo (Dekel & Birnboim 2006), vigorous star formation may drive outflows that will exhaust the gas supply in short timescales ("overconsumption"; see McGee et al. 2014; Balogh et al. 2016).

To have a clearer picture of the quenching mechanisms in play, one needs to track how galaxies change their star formation properties as they fall into the cluster, isolating the role of the environment from galaxies' internal processes. In this paper, we combine the SpARCS (Spitzer Adaptation of the Red-Sequence Cluster Survey; Muzzin et al. 2009; Wilson et al. 2009) sample of clusters with the first public data release from the HSC-SSP (Hiper Suprime-Cam Subaru Strategic Program; Aihara et al. 2018) in the XMM-LSS and ELAIS-N1 fields. We focus on examining the fraction of SF cluster galaxies, out to as far as 10 virial radii, with a sample sufficiently large enough to perform background corrections, to control for the galaxies' intrinsic properties (such as stellar mass), over a large range of redshift. This allows us to track how star formation is quenched as galaxies infall into clusters, providing a clearer picture of how galaxies in clusters evolve.

The paper is structured as follows. Section 2 introduces the data set, along with the cluster galaxy selection. Section 3 describes how galaxy properties have been determined, including possible systematics. In Section 4 we present our results, focusing on the dependence of the fraction of SF galaxies on stellar mass and distance to the cluster center at different redshift bins. We also quantify the environmental and mass quenching efficiencies and their dependences on stellar mass and cluster-centric radius. In Section 5 we discuss our results, and in Section 6 a summary of this work is given. Unless otherwise stated, throughout this paper we adopt a ΛCDM cosmology with H₀ = 70 km s⁻¹ Mpc⁻¹, Ω_m = 0.3, and Ω_Λ = 0.7. All magnitudes are on the AB system.

2. Data and Sample Selection

2.1. Optical and Infrared Data

We use the first public data release from the HSC-SSP⁵ (Aihara et al. 2018) project for the optical broadband photometry: g, r, i, z', and Y bands. This catalog covers a large region of the XMM-LSS SpARCS field and the majority of the ELAIS-N1 SpARCS field, as shown in Figure 1. The total overlapping area is about 16.5 deg². The HSC-SSP survey consists of three layers of depth (wide, deep, and ultradeep), with our fields of interest being part of the deep layer. The 5σ depths estimated by the HSC-SSP team for the deep layer are 26.8, 26.6, 26.5, 25.6, and 24.8 mag for the g, r, i, z, and Y bands, respectively.

**Figure 1.** Coverage of the available data in the ELAIS-N1 (top) and *XMM*-LSS (bottom) fields used in the present work. Blue regions represent the HSC-SSP deep layer from the first data release, red areas show the SpARCS/SWIRE coverage, green areas are the SERVS data, and purple regions come from UKIDSS-DSX. The clusters analyzed in this study (black circles; see Section 2.5) are in the overlapping regions of HSC and SpARCS.
Download figure:
Standard image High-resolution image

For near-IR (NIR) data we benefit from the UKIDSS project (Lawrence et al. 2007). From the final data release of DXS (Deep Extragalactic Survey, DR10), we have J and K data covering a significant portion of both fields.

Additional NIR photometry comes from the SWIRE Legacy survey (Lonsdale et al. 2003), which is the NIR-band photometry for SpARCS. This survey has imaged nearly 50 deg², in six high-latitude fields, with Spitzer at 3.6, 4.5, 5.8, and 8 μm using IRAC (Fazio et al. 2004) and at 24 μm using MIPS (Rieke et al. 2004). Additional deeper observations in IRAC 3.6 and 4.5 μm bands are available from the more recent Spitzer Extragalactic Representative Volume Survey (SERVS; Mauduit et al. 2012). IRAC data are taken from the band-merged catalog of the XMM-LSS and ELAIS-N1 fields. For these merged catalogs we estimated the 90% completeness magnitude limit values of 20.7, 20.8, 19.6, and 19.5 mag for 3.6, 4.5, 5.8, and 8 μm bands, respectively. The details of catalog construction are discussed in Section 2.4.

2.2. Spectroscopic Redshifts

The HSC-SSP DR1 catalog already includes the significant amount of spectroscopic data publicly available in these fields. For the ELAIS-N1 field data come from Sloan Digital Sky Survey (SDSS) DR12 (∼3100 spectroscopic sources). For the XMM-LSS field, data from various surveys are available: SDSS DR12, PRIMUS DR1 (Coil et al. 2011), VIPERS PDR1 (Garilli et al. 2014), and VVDS (Le Fèvre et al. 2013). The overlapping regions of these spectroscopic surveys with the HSC-SSP catalog of the XMM-LSS field contain a total of 49,930 spectroscopic sources. We use these data to verify our cluster and photometric redshifts.

2.3. Photometric Redshifts

Photometric redshifts are taken from the HSC-SSP DR1. The release contains different photometric redshifts computed using six independent codes. We use the values obtained through the DEmP algorithm (Hsieh & Yee 2014). This code combines the nearest neighbor technique in multiple color–magnitude spaces with a polynomial fitting technique. The redshift value for each object is estimated using the 40 nearest neighbors in the nine-dimensional space (five magnitudes axes and g − r, r − i, i − z', z' − Y color axes using a point-spread function [PSF] matched aperture photometry) with a linear function. The redshifts obtained with this methodology show a good correlation with their corresponding spectroscopic values (σ_(1+z) = 0.049), especially within our range of interest, between 0.2 and 1.1. Further details on the photometric redshift computation of HSC-SSP data are described in Tanaka et al. (2018).

2.4. Selection of Galaxies with Clean Photometry

For selecting sources with reliable photometry, we follow the HSC-SSP team recommendations. We select primary objects, which means objects that have been deblended and are in the inner patch and tract.⁶ We also apply a set of flags to make sure that the objects do not suffer from problematic pixels (flags_pixel_interpolated_center, flags_pixel_edge, flags_pixel_cr_center), their centroids are correct (centroid_sdss_flags), their photometric measurements are good (parent_flux_convolved_flags, cmodel_flux_flags), and they survived the junk suppression step⁷ (detected_notjunk) for the five photometric bands.

The star/galaxy separation was performed using the classification_extendedness parameter, which estimates the difference between photometry from PSF and CModel.⁸ For compact sources, the CModel measurement approaches PSF photometry asymptotically and the difference in magnitude becomes small. In the HSC catalog, the extendedness parameter takes the value 0 when the difference is sufficiently small and 1 otherwise. Thus, we remove from our sample only those objects that are classified as point sources in all five HSC bands. This stringent criterion for classifying an object as a star is due to faint small galaxies that are often classified as point sources. Furthermore, at faint magnitudes, there are many more galaxies than faint stars in high Galactic latitude fields, and the overwhelming number of objects are galaxies. The number of stars rejected this way represents 6.5% of our galaxy catalog. In Figure 2 black and gray histograms show the distributions before and after removing the stars, respectively. In particular, the peak seen in the 0.7 < z < 0.9 panel, placed at rare blue V − J colors and red U − V colors, shows the importance of removing such point sources from our catalog. Taking advantage of the ancillary J and K data, we perform a sanity check of the star classification based on the extendedness parameter. We find that 90%, 77%, 87%, and 59%⁹ of the objects classified as stars at z ∼ 0.4, 0.6, 0.8, and 1, respectively, show blue J − K colors (J − K < 1).

**Figure 2.** *UVJ* color–color diagrams for different cluster redshift bins. Red lines separate *UVJ*-quiescent (top left region) and *UVJ*-SF (bottom and right regions) galaxies. The number in the lower right corner of each central panel indicates the stellar mass completeness limit used in each redshift bin. Top and right panels show V − J and U − V histograms, respectively. Gray histograms show the final sample of galaxies employed, while black histograms show the sample before removing the point-like sources (see Section 2.4).
Download figure:
Standard image High-resolution image

Once we have our clean optical photometric catalog (griz'Y), we cross-match it with the infrared catalogs. We use the nearest neighbor technique using a maximum separation radius of 1'' for UKIDSS and 1 farcs 5 for IRAC/SWIRE+SERVS. The absolute separation distribution peaks at 0 farcs 06 and 0 farcs 2 for UKIDSS and IRAC/SWIRE+SERVS, respectively.

2.5. Cluster Sample and Membership

The clusters analyzed in this study are from two of the six SpARCS (Muzzin et al. 2009; Wilson et al. 2009) fields. Clusters are originally found using a slightly modified version of the cluster red sequence algorithm developed by Gladders & Yee (2000, 2005). The code is described in more detail in Muzzin et al. (2008). The only modification of the algorithm from Muzzin et al. (2008) is that of the optical band in the color, z' − 3.6 μm instead of R − 3.6 μm. From the full sample of cluster candidates we choose those with richness N_red > 6.0¹⁰ (equivalent to ${B}_{\mathrm{gc}}\,\gtrsim \,330$ ¹¹ ), in order to reduce the number of false clusters. The approximate cluster mass limit is 9 × 10¹³ M_⊙, computed using Equation (9) from Muzzin et al. (2007), which relates B_gc and M₂₀₀ (the mass within a radius that encloses a density 200 times the critical density of the universe).

The original SpARCS sample has 188 and 167 cluster candidates in the ELAIS-N1 and XMM-LSS fields up to redshift 1.2 and N_red > 6.0, respectively. This is chosen as our parent sample, which we further refine in several ways. The first is that the redshifts were defined using only the red sequence, but now we have full photometric redshift information. We use this information to refine the cluster redshifts as follows:

1.
We compute stacked red sequences for clusters in redshift bins of 0.1 from ${z}_{\mathrm{SpARCS}}=0.3$ up to 1.1, where z_SpARCS is the original SpARCS catalog cluster photometric redshift, by considering galaxies at a maximum distance from the cluster center of R_MAX = 1.5 Mpc and a redshift range defined by σ_Δz = 0.05, where we use both photometric and spectroscopic redshifts. Note that we use observed magnitudes, so the color–magnitude diagram employed is different for different redshift bins. We fit a line to the galaxies within the red color region (1.1–1.8, 1.3–1.9, and 1.3–1.9 in g − r vs. r for z_SpARCS ∼ 0.35, 0.45, and 0.55, respectively; 1.1–1.9, 1.4–2.1, 1.6–2.3, and 1.7–2.4 in r − z' vs. z' for z_SpARCS ∼ 0.65, 0.75, 0.85, 0.95, respectively; and, 1.9–3.0 in r − Y vs. Y for z_SpARCS ∼ 1.05) and define the 98% confidence intervals of the fit as the red sequence region.
2.
For each cluster, we utilize galaxies interior to the maximum cluster-centric radius of R_MAX = 0.5 Mpc and σ_Δz = 0.2 within the fitted red sequence region to define an initial cluster redshift, z_CL, which is the peak of the redshift probability density function resulting from the sum of the z_PDF of all these red sequence galaxies. Here we reduce the distance to the cluster center to avoid contamination due to nearby clusters in the projected plane, and we use a wider redshift range to take into account the uncertainty in the original SpARCS redshift estimation.
3.
For clusters with catastrophic photometric redshifts, defined as either a difference between the redshift of the SpARCS brightest cluster galaxy (BCG) and z_SpARCS of larger than 0.15 × (1 + z) or when fewer than three red sequence galaxies are found in the previous step, we define a provisional z_CL as the peak of the redshift PDF resulting from the sum of the z_PDF of all galaxies, not only the red sequence ones. We perform step 2 again, changing z_SpARCS to the provisional z_CL, and estimate the initial z_CL when more than two red sequence galaxies are found.
4.
We repeat steps 1 and 2, changing z_SpARCS to the initial z_CL to determine the final z_CL. The red sequence region obtained after the second iteration is shown in Figure 3 as the gray shaded area between the black curves. From the original SpARCS catalog, the persistent catastrophic z_SpARCS clusters (i.e., when fewer than three red sequence galaxies are found after the second iteration) are flagged with value 9, and those clusters outside the HSC-SSP field of view are flagged with value −99. Clusters flagged with either of these values are not included in the analysis of this paper.

**Figure 3.** Color–magnitude diagrams for different cluster redshift bins in the two fields: ELAIS-N1 (top two rows) and *XMM*-LSS (bottom two rows). Red and blue colors represent the 2D histograms of *UVJ*-quiescent and *UVJ*-SF galaxies, respectively. The number in the lower left corner of each panel indicates the number of clusters in each bin, which are stacked and used to calculate the red sequence. The shaded region between the two black curves defines the red sequence region, i.e., the 98% confidence interval of the fit (used to determine the final cluster redshift; see points 1 and 4 in Section 2.5). The dotted lines indicate the magnitude completeness limit.
Download figure:
Standard image High-resolution image

**Figure 4.** Photometric cluster redshift, estimated as detailed in Section 2.5, vs. the spectroscopic cluster redshift obtained from a subsample of spectroscopic galaxies within 15,000 km s⁻¹ of BCGs that have spectroscopic redshifts (blue circles) and from the *XMM*-LSS survey (orange squares; Adami et al. 2011; Willis et al. 2013). The size of the symbols is scaled with the cluster richness. The total number of clusters shown is 80, as we do not plot the 11 clusters from the BCG-selected sample that are already included in the X-ray sample.
Download figure:
Standard image High-resolution image

After these cuts, the final sample of clusters with which we work comprises 127 in ELAIS-N1 and 82 in XMM-LSS. We then select potential cluster galaxies by selecting galaxies that are within the cluster photometric redshift slices. We adopt spectroscopic redshifts when available and DEmP photometric values for the rest of the galaxies. The final galaxy catalog includes a total of 1,264,407 galaxies, with 60% of them in the ELAIS-N1 field and the remainder in the XMM-LSS field. For each cluster, we use galaxies that are included in a slice defined by the z_CL ± 0.05 × (1 + z_CL) range and up to 20 Mpc in cluster-centric radius. We note that since the slices associated with the clusters can overlap, a galaxy in overlapping areas can be included in multiple slices.

We compare our photometrically determined cluster z_CL with available spectroscopic measurements from the XMM-LSS survey from Adami et al. (2011) and Willis et al. (2013) for 20 clusters, and we obtain a good correlation with a slope of 1.04 and rms = 0.015. Moreover, we take advantage of the large spectroscopic sample to double-check our z_CL values. We identify clusters whose BCGs have spectroscopic redshifts and compute the peak of the distribution of the spectroscopic redshifts, including that of the BCG, within R₂₀₀ and σ_z = 0.05. Considering the 70 BCG-selected clusters with three or more spectroscopic redshifts, we find that the rms of the z_CL versus z_SPEC relation is 0.02. We note that 11 out of these 70 clusters overlap with the XMM-LSS cluster sample. Figure 4 shows these results. These comparisons show that the photometric cluster redshift estimates are very accurate.

3. Fraction of Star-forming Galaxies

3.1. Stellar Masses and UVJ Rest-frame Colors

With our 11-band ( ${grizYJK},$ 3.6, 4.5, 5.8, 8.0 μm¹² ) catalog we estimate stellar masses and rest-frame colors. For both computations, we use SED-fitting software and the photometric/spectroscopic redshifts of the galaxies themselves.

To determine stellar masses, we use the LePhare code (Arnouts et al. 1999; Ilbert et al. 2006) to fit the stellar part of the spectrum with Bruzual & Charlot (2003) population synthesis models, with star formation histories exponentially declining as e^−t/τ. The complete template library was built considering solar metallicity, nine different values of τ (from 0.1 to 30 Gyr) with 57 steps in age, and the extinction law of Calzetti et al. (2000) with values of E(B − V) ranging from 0 to 1. The templates also account for the contribution of the emission lines to the flux. We assume the initial mass function of Chabrier (2003).

The rest-frame U − V and V − J colors are computed with the EAZY (Brammer et al. 2008) program. We use a combination of a reduced set of seven templates, similar to the Brammer et al. (2008) original set: five are the reduced set from Grazian et al. (2006) templates but with emission lines to match the Ilbert et al. (2009) prescription, the sixth template accounts for dusty galaxies that do not appear in semianalytic models, and the final template is an evolved Maraston (2005) model to be red enough for massive old galaxies at z < 1. Our lower redshift limit of z = 0.3 is set to ensure a good estimation of rest-frame magnitudes, as the bluest filter we have available is the g band.

These rest-frame colors allow us to classify galaxies as UVJ-SF and UVJ-quiescent, as shown by the UVJ color–color diagrams of Figure 2 and the 2D colored distributions of Figure 3. The cuts used to define the SF and quiescent regions are slightly modified from the calibration made by Williams et al. (2009) to fit the bimodality in our data. The criteria are U − V > 1.45, V − J < 1.5, and U − V > 0.88 × (V − J) + 0.59. Figure 3 shows a good agreement between the classification made using UVJ colors and the red sequence computed in Section 2.5 for the selection of cluster galaxies, where most of the quiescent population falls within the red sequence region and some of the SF galaxies populate its faint end.

3.2. Stellar Mass Completeness

A careful study of the evolution of the SF fraction needs to take into account different completeness limits at different redshifts. Therefore, we have examined the stellar mass differential distributions, and we define the stellar mass threshold for each redshift range as the peak and, when studying evolution, compare samples in equivalent stellar mass bins. The values obtained for the stellar mass limit are log M_⋆/M_⊙ = 8.6, 8.6, 9.5, and 9.7 for z ∼ 0.4, 0.6, 0.8, and 1, respectively.

3.3. Background Subtraction

Our estimation of the number of SF galaxies over the total number of cluster galaxies is contaminated by the inclusion of field galaxies owing to our selection method based on photometric redshifts. We make use of the large coverage of our catalog to correct for this background contamination. We generate background field control samples for each cluster considering those galaxies in an annular area of inner R/R₂₀₀ radius 6 and outer R/R₂₀₀ radius 9. The selection of these radii is a compromise between being sufficiently large that no cluster galaxies are included and being not so large that the correction for the area outside our spatial coverage and the contamination due to nearby clusters are small. The results we show are stable against changes in the size of the area used for defining the background control sample, which we expect, given that the surface galaxy density as a function of radius is almost flat beyond ∼2.5 Mpc. The subtraction of the background signal is performed when computing the fraction of SF galaxies, as described in Section 3.4.

3.4. The SF Fraction with Bayesian Inference

The fraction of SF galaxies is, by definition, in the [0, 1] range, but the presence of background sources could yield unphysical values of negative fractions or values >1 (a subsample is larger than the whole sample). This occurs when the formula

$\begin{eqnarray}&&{f}_{\mathrm{SF}}(\mathrm{Cluster})=\displaystyle \frac{{n}_{\mathrm{SF}}(\mathrm{Cluster}+\mathrm{Back})-\,{n}_{\mathrm{SF}}(\mathrm{Back})}{{n}_{\mathrm{TOT}}(\mathrm{Cluster}+\mathrm{Back})-\,{n}_{\mathrm{TOT}}(\mathrm{Back})}\end{eqnarray} \tag{ 1 }$

is used to compute the SF fraction in clusters. To avoid nonsense fractions, we use Bayesian inference to calculate the posterior probability that clusters have a fraction of SF galaxies f_SF, given the same quantities that appear in Equation (1). We followed the methodology described in Andreon et al. (2006) and the mathematical details expounded in D'Agostini (2004). However, instead of taking the median value from the posterior distribution as in Andreon et al. (2006), our estimate of the cluster f_SF is the peak of the posterior distribution, and its uncertainty is defined as the FWHM interval.

We compute the f_SF for each cluster as a function of the cluster-centric distance in four cluster redshift bins (0.3 ≤ z < 0.5, 0.5 ≤ z < 0.7, 0.7 ≤ z < 0.9, 0.9 ≤ z < 1.1) and five stellar mass bins (8.6 ≤ log M_⋆/M_⊙ < 9.0, 9.0 ≤ log M_⋆/M_⊙ < 9.5, 9.5 ≤ log M_⋆/M_⊙ < 10.0, 10.0 ≤ log M_⋆/M_⊙ < 10.5, 10.5 ≤ log M_⋆/M_⊙ < 11.2). The R/R₂₀₀ bin width is 0.2, except for the highest stellar mass bin, where the width is increased to 0.4. R₂₀₀ is computed using Equation (2) from Muzzin et al. (2008), which relates B_gc and R₂₀₀. The results are presented in Figure 5, where the f_SF values for individual clusters are shown as light-colored points. To estimate the composite f_SF in each radial bin, we stack the galaxies from all clusters in each bin and apply the same Bayesian inference approach as for individual clusters. The result obtained for the composite f_SF is completely compatible with the mean of the f_SF of individual clusters. The composite values are plotted as dark circles. The diamond symbols and the light horizontal lines represent the near field, from R/R₂₀₀ = 4 to 6, which we use as a comparison "field environment" and refer to as "field" hereafter. These f_SF field values have also been background corrected. The squares denote the value of the "background field" region employed for the background correction from R/R₂₀₀ = 6 to 9. We note here that our f_SF composite values at R/R₂₀₀ ≳ 2 lie along the value of the field. We note that this behavior is not a bias of the method because the points have net counts over that of the background field values.

**Figure 5.** Fraction of SF galaxies as a function of the distance to the center of the cluster (R/R₂₀₀). Panels from left to right increase in redshift and from top to bottom in stellar mass, as indicated by top titles and labels on the right. SF galaxies are selected by their *UVJ* rest-frame colors. As described in Section 3.4, the fraction is estimated for individual clusters (smaller light-colored circles) using Bayesian inference to account for the background contamination. The composite fraction is computed by stacking galaxies from all clusters and applying the same Bayesian formula (bigger dark circles), and the errors are the FWHM of the posterior distribution. In the right secondary panels, diamond and square symbols indicate the mean f_SF of the "field environment" (also the horizontal line in main panels), defined by the region 4 ≤ R/R₂₀₀ < 6, and the control "background field," defined by the region 6 ≤ R/R₂₀₀ < 9. The width of R/R₂₀₀ bins is 0.2, except for the highest stellar mass bin, where it is increased to 0.4.
Download figure:
Standard image High-resolution image

3.5. Elimination of BCGs

As we are analyzing the effect of the cluster environment on the quenching of the star formation in its member galaxies, we note that the BCGs are special objects whose formation and evolution are complex: they are involved in a combination of galaxy merging, cluster cooling flows, AGN feedback, and star formation processes (e.g., Rawle et al. 2012; Webb et al. 2015; McDonald et al. 2016). Thus, these processes will affect our results, particularly in a sensitive region as is the innermost R/R₂₀₀ of the highest stellar mass bin. We select BCGs as the brightest galaxy in the same magnitude used for estimating the red sequence (i.e., r, z', and Y bands for z = 0.3–0.6, 0.6–1.0, and 1.0–1.1, respectively) and remove them from our cluster galaxy sample.

4. Results

In this section we present the results on the dependence of the SF fraction on galaxy stellar mass, distance to the cluster center, and redshift. Our aim is to trace how the evolution of the fraction of galaxies that quench their star formation is connected to its cluster environment by using cluster-centric radius. For this purpose, we need to distinguish the effect of the cluster environment from the secular evolution of galaxies. Our approach is to measure the SF fraction as a function of the environment at fixed stellar mass and vice versa.

4.1. The Fraction of Star-forming Galaxies: Environment and Mass Dependency

The large amount of data available allow us to study the f_SF in four epochs, in R/R₂₀₀ bins of 0.2 from the cluster core to the field, and over a large range of galaxy stellar masses at the same time, while maintaining small errors, as shown in Figure 5. Note that for each redshift we only show stellar mass bins down to their completeness limit. The panels in this figure give us an idea of the role of the stellar mass and the environment in quenching star formation activity and their evolution.

Focusing on a single epoch, i.e., looking at one column in Figure 5, the composite f_SF decreases drastically from 1 to below 0.5 from low to high stellar masses. This drop is more dramatic for the lower-redshift bins: from log M_⋆/M_⊙ ≃ 9.7 to 10.8 the field SF population is reduced by a factor of ∼4. Over the redshift range from 0.7 to 0.9 the field f_SF is reduced by a factor of ∼3 over the same mass range. The decline of the f_SF from low to high stellar masses is evident not only in the reference value of the field but also in the cluster core, where the SF population decreases by a factor of ∼4 from log M_⋆/M_⊙ ≃ 9.7 to 10.8 in all three lower-redshift bins. This trend is seen more clearly in Figure 6, where the R/R₂₀₀ dimension has been binned into three regions and the f_SF is presented as a function of the stellar mass using narrower stellar mass bins. The regions used are the cluster core up to R/R₂₀₀ = 0.4, the cluster outskirts combining 0.4 ≤ R/R₂₀₀ < 1, and the field as estimated between 4 ≤ R/R₂₀₀ < 6. We deliberately exclude the plot of the 1 ≤ R/R₂₀₀ < 4 region because it overlaps with the field curve, as one can infer from Figure 5.

**Figure 6.** Fraction of SF galaxies as a function of stellar mass in three different environmental regions: 0 ≤ R/R₂₀₀ ≤ 0.4 (cluster core; dark-colored circles), 0.4 ≤ R/R₂₀₀ ≤ 1 (cluster outskirts; medium-colored crosses), and 4 ≤ R/R₂₀₀ ≤ 6 (field; light-colored diamonds). We include the field value at the highest-redshift bin as a reference in the other panels (dashed gray curve). Arrows indicate the upper limit of the SF fraction at stellar mass bins where the SF population is complete but the passive population is not. Panels from left to right increase in redshift as indicated by top titles. The dot-dashed horizontal lines in the right panel are the locus of the fractions used to define ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ . As in the previous figure, the fraction is computed summing up galaxies from every cluster and using Bayesian inference to account for the background contamination, and the errors are the FWHM of the posterior distribution. The width of stellar mass bins is 0.3 dex.
Download figure:
Standard image High-resolution image

With these data, we can quantify two characteristic stellar masses: the mass at which the quenching process starts, ${M}_{\mathrm{INI}}^{* }$ , and the mass at which almost all galaxies are quenched, ${M}_{\mathrm{END}}^{* }$ , using a simple procedure. We estimate ${M}_{\mathrm{INI}}^{* }$ at where f_SF becomes lower than 0.8 and ${M}_{\mathrm{END}}^{* }$ at where f_SF goes below 0.2 (as marked in the right panel of Figure 6). The result is presented in Figure 7. Note that all measured values are above our stellar mass completeness limit. Within our redshift range of study, we observe a clear trend in both measures: these characteristic stellar masses evolve to lower values toward lower redshifts.

**Figure 7.** ${M}_{\mathrm{INI}}^{* }$ (blue bottom lines) and ${M}_{\mathrm{END}}^{* }$ (pink top lines) values, as defined in Section 4.1, as a function of redshift for three different environments: 0 ≤ R/R₂₀₀ ≤ 0.4 (cluster core; dark-colored circles), 0.4 ≤ R/R₂₀₀ ≤ 1 (cluster outskirts; medium-colored crosses), and 4 ≤ R/R₂₀₀ ≤ 6 (field; light-colored diamonds). The two arrows indicate that the measured value is just below our mass completeness limit.
Download figure:
Standard image High-resolution image

Figure 6 reveals information not only on the quenching effects caused by stellar mass but also about the environmental quenching process that takes place in the cluster environment. In these four panels, the curves representing the cluster core show a shortage of the SF population compared to the field at all stellar masses. We observe little difference in the f_SF versus stellar mass relation between the outskirts of the cluster and the field in Figure 6 for the two higher-redshift bins. However, the relation for the cluster outskirts becomes slightly lower for the 0.5 ≤ z < 0.7 bin and significantly lower for the lowest-redshift bin at all stellar masses. This indicates that little or no environmental quenching effect as measured by f_SF is evident beyond R/R₂₀₀ ∼ 1 for clusters at z ≥ 0.7.

4.2. Evolution of the SF Fraction: Environment and Mass Dependency

Focusing on a single stellar mass bin in Figure 5, i.e., analyzing rows, we can study the evolution of the SF fraction with redshift. In Figure 8, in a fashion similar to previously done in Figure 6, we bin the R/R₂₀₀ dimension in three regions to show the fraction of SF galaxies as a function of epoch. In these panels we observe the strong evolution of the fraction of SF galaxies with redshift, first seen by Butcher & Oemler (1984) and subsequently in many other studies (e.g., Balogh et al. 1999; Poggianti et al. 2006; Saintonge et al. 2008; Brodwin et al. 2013; Wagner et al. 2017; Jian et al. 2018). Furthermore, our large sample of clusters covering a relatively large range of redshift demonstrates clearly the dependence of this evolution on both galaxy stellar mass and the cluster-centric radius. In the far right panel of Figure 8, we see how more massive cluster galaxies (log M_⋆/M_⊙ ≥ 10.5) are mostly quenched by z ∼ 0.6 (f_SF ≤ 0.1). We find that the change of f_SF is fractionally larger for higher stellar mass galaxies, in both the field and the cluster environment.

The large area coverage of our data allows us to analyze large radii and define the field and the cluster in a consistent way, being able to confidently compare the f_SF evolution curves of the different environments. The four panels in Figure 8 show how the cluster environmental impact shifts the f_SF to lower values compared to the field for all redshift and stellar mass bins, but such a shift also depends on redshift and stellar mass. The shift between the f_SF of the cluster core and the field is fractionally larger for lower-redshift bins and higher-mass galaxies. Thus, the environment causes a steepening of the f_SF versus redshift curves.

4.3. Possible Systematics and Their Effects on the ${f}_{\mathrm{SF}}$ Results

4.3.1. Photo-z Uncertainties for Different Populations

A quantity that can be different for SF and quiescent populations is the photometric redshift dispersion σ_Δz, as the form of the SED of quiescent galaxies makes it easier to estimate the correct redshift through photometric data fitting. We employ the available spectroscopic data to examine the dispersion of the z_SPEC versus z_PHOT relation for these two populations. We find that, although the SF sample contains more outliers (20% versus 14% of the galaxies are outside the 0.05 × (1 + z) interval), the σ_Δz for the SF and quiescent populations are not significantly different and, down to our mass completeness limit, are less than 0.05. This means that we are not biasing our result when selecting cluster galaxies in redshift slices of σ_Δz = 0.05.

4.3.2. Cluster Sample: Missing Clusters?

As mentioned in Section 2.5, our initial sample of clusters is limited to those with N_red > 6.0, and then it is cleaned by excluding clusters that fall outside of the HSC area or do not have enough red sequence galaxies to determine z_CL by means of a photometric redshift distribution. This procedure helps us to reduce the contamination of false clusters, but it might bring the question of the lack of clusters with very blue populations, as by construction our sample is biased toward clusters with well-defined red sequences. Nevertheless, the richness limit of N_red = 6.0 is rather low, and we require only a small number of red sequence galaxies to estimate the cluster redshift, just >3. Even for the case of a cluster with a very high fraction of blue galaxies, we still expect it to have a number of red galaxies, and since we find no literature reporting the detection of such blue clusters, we can safely argue that the number of missing clusters is very low. In addition, when comparing X-ray-selected clusters from the XMM-LSS distant cluster survey and clusters from the red sequence method of the SpARCS cluster survey, Willis et al. (2018) find that 10/10 spectroscopically confirmed XMM-LSS clusters at 0.8 < z < 1.2 are matched to a SpARCS cluster at the same redshift range. Foltz et al. (2015) observed no measurable differences between the properties of quiescent cluster members in clusters selected by different methods at z ∼ 1, and they found a remarkable agreement in color–magnitude relation zero-point throughout their cluster samples, indicating that the red sequence method does not preferentially select older, more evolved systems. Furthermore, we note that the f_SF in the field is much lower than 1 for high-mass galaxies at all redshifts (Figure 8), so even clusters that show no additional quenching above what would be expected for a field population would still be detected using the red sequence method.

4.3.3. Cluster Sample: Progenitor Bin

Lin et al. (2017) probed a similar redshift range using a sample of clusters constructed with another red sequence algorithm and concluded that the top N clusters from redshift bins occupying similar comoving volumes are a good approach for studying the evolution of clusters and their galaxy populations. Our data allow us to examine the same redshift bins as in Lin et al. (2017; z = 0.3–0.6, 0.6–0.77, 0.77–0.9, 0.9–1.02) and select the top N = 30 halos to probe the solidity of our results. We obtain compatible results: the f_SF reference value (horizontal lines in Figure 5) remained practically unaltered, with variations of around 3%. The larger differences appear at the highest-mass bin of the 0.6–0.77 redshift bin, where f_SF values at R/R₂₀₀ ∼ 1.0 are higher and also the field f_SF value is around 20% higher, which could be explained by the slightly narrower and higher-redshift bin, since f_SF values increase with redshift. In general, the overall shape of the SF fraction as a function of R/R₂₀₀ is almost indistinguishable from that of using the whole sample.

4.3.4. Cluster Center

As we are analyzing the effect of the environment as measured by the distance to the center of the cluster, normalized by its R₂₀₀ value, the definition of that center could affect our results. Thus, we have also tested the dependence of our results on the identification of the center of the cluster. We define the cluster center as the center of the overdensity of galaxies, as estimated by the SpARCS red sequence algorithm, but we can also calculate the distance to the BCG and use this as a proxy for the R/R₂₀₀ parameter. We can safely conclude that there are no systematics due to the exact definition of the center of the cluster, as using the BCGs does not change our results on the radial dependence of the SF fraction. We observe comparable shapes in all the projections of the f_SF, R/R₂₀₀, M_⋆, and z parameter space, showing no environmental quenching beyond $R/{R}_{200}\sim 1.5$ in the redshift range we cover.

4.3.5. UVJ Colors

The UVJ diagram exhibits a strong dependence on the stellar mass: smaller (bluer) U − V values are populated by lower stellar mass galaxies, while the highest (reddest) U − V values are populated by the most massive galaxies. Thus, our results, especially when binning in stellar mass, could strongly depend on the UVJ boundaries selected (see Section 3.1). We have observed that slight variations in the UVJ limits for classifying SF/quiescent galaxies change the reference value of the SF fraction but not its shape, i.e., the horizontal line in Figure 5 moves up/down as we shrink/spread the limits with respect to the passive population. The highest stellar mass bin is particularly sensitive to the V − J threshold, and using threshold values larger than 1.6 makes the f_SF zero even for larger cluster-centric distances. Nevertheless, we note that the radial gradient of the SF fraction is always apparent, independently of the UVJ limits used.

5. Discussion

Considering "mass quenching" as part of the secular evolution of a galaxy, we identify the field f_SF as the result of the processes related to the "nature" of a galaxy. On top of that fiducial f_SF value we measure the effect of the cluster environment, which can be interpreted as an acceleration of the quenching process, the "nurture" component. In this scenario, we analyze the dependence of the f_SF on galaxies' stellar mass and position within the cluster halo and how it evolves.

5.1. The Dependence of the ${f}_{\mathrm{SF}}$ on Galaxy Stellar Mass

Focusing on the f_SF curves representing the field environment in Figure 6, we observe how the mass quenching pushes down the f_SF starting with the high-mass galaxies at earlier epochs and affecting increasingly lower stellar mass galaxies as we move to more recent epochs. We introduce the parameters ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ (Section 4.1) as a direct way to quantify the dependence of f_SF on stellar mass. We are witnessing the mass-dependent evolution of the quenched population: more massive galaxies turn off their star formation, becoming red and passive first, as has been observed by many others (e.g., Cowie et al. 1996; Marchesini et al. 2009; Muzzin et al. 2013; Sobral et al. 2014; Tomczak et al. 2014). Our approach allows us to determine the f_SF dependence on stellar mass for different cluster-centric radii and for different redshift bins. In the same way as described by Fang et al. (2018) through their UVJ diagram grid, we clearly see the buildup of the quiescent population in both the cluster and the field in Figure 5. By choosing a stellar mass and redshift bin and searching for the most similar plot at higher mass, we will always identify it with that of an earlier epoch (e.g., panel 9.5 ≤ log M_⋆/M_⊙ < 10 at 0.5 ≤ z < 0.7 is most similar to panel 10 ≤ log M_⋆/M_⊙ < 10.5 at 0.9 ≤ z < 1.1).

The rise of ${M}_{\mathrm{END}}^{* }$ (and ${M}_{\mathrm{INI}}^{* }$ ) with redshift also confirms that more massive galaxies evolve more quickly (Figure 7); while log M_⋆/M_⊙ > 11.1 galaxies are already mostly quenched by z ∼ 1, most of the log M_⋆/M_⊙ ≳ 10.6 galaxies are also quiescent by z ∼ 0.4. The ${M}_{\mathrm{END}}^{* }$ value we have defined is equivalent to the upper stellar mass limit of the blue cloud recently measured by Haines et al. (2017) using the d4000–M_⋆ plane for field galaxies. The authors find that the upper stellar mass limit of the blue population is steadily retreating toward lower redshifts, from log M_⋆/M_⊙ = 11.2 at z ∼ 0.9 to log M_⋆/M_⊙ = 10.9 at z ∼ 0.58. This is in good agreement with our field ${M}_{\mathrm{END}}^{* }$ values. Remarkably, the trend we observe for field ${M}_{\mathrm{INI}}^{* }$ is also compatible with the local (z < 0.055) study of Geha et al. (2012), since these authors find that quenched galaxies in the field do not exist below log M_⋆/M_⊙ = 9.

Our quantitative measures of field ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ for different redshifts could be useful to test galaxy formation models and constrain the physical mechanisms responsible for mass quenching at 0.3 ≤ z < 1.1. These kinds of estimates demonstrate the capability of photometric redshift studies to produce useful information by providing sample sizes that are much larger than those from spectroscopic surveys over a wide redshift range.

5.2. The Dependence of the ${f}_{\mathrm{SF}}$ on Environment

5.2.1. Radial Gradient of ${f}_{\mathrm{SF}}$

On top of the "natural" buildup of the quiescent population, we observe an excess of quiescent galaxies toward the cluster core that gives evidence to a quenching process driven by the cluster environment. Many previous works have shown the radial dependence of the red/blue fraction in groups and clusters of galaxies. Using a sample of ∼1000 clusters at 0.45 ≤ z ≤ 0.9 drawn from the first Red-Sequence Cluster Survey (RCS1), Loh et al. (2008) found that the red fraction of R/R₂₀₀ < 0.25 galaxies was much greater than at larger radii and that there is a moderate evolution in the red fraction at radii smaller than 0.5R₂₀₀ over the small redshift range. Raichoor & Andreon (2012), using an X-ray-selected cluster sample of 25 clusters up to z ∼ 1, showed that the dependence of the fraction of blue galaxies with cluster-centric distances held for different galaxy mass bins. More recently, Lin et al. (2014) and Jian et al. (2017) studied galaxy groups in the Pan-STARRS1 Medium-Deep Survey and found that the quiescent fraction slightly decreased as the radius increased and that the slope of f_Q (quiescent fraction) as a function of the cluster-centric radius was steeper for less massive galaxies (log M_⋆/M_⊙ ≲ 10.1). Our work confirms and extends these findings to lower stellar mass bins.

We observe a confinement of the cluster effect to within R ≃ R₂₀₀. From panels in Figure 5, we can infer that the effect is mostly interior to R/R₂₀₀ ∼ 0.5, except for the lowest-redshift bin or the highest stellar mass bin, where the effect of the environment extends to R/R₂₀₀ ∼ 1–1.5. This is also observable in Figure 6, where the cluster f_SF versus log M_⋆/M_⊙ curves show that the strongest impact happens at the very cluster core, while the outskirts curves show that the cluster environment affects larger R/R₂₀₀ primarily toward later epochs. Allen et al. (2016) found a similar extent of the cluster effect: cluster SF galaxies within 0.5R_vir have lower mass-normalized average sizes and a higher fraction of Sérsic indices with n > 1 than field SF galaxies. However, they measure consistent mass-normalized average star formation rate (SFR) for SF galaxies in the cluster outskirts and core, but this is elevated by a factor of two in the field, suggesting that environmental effects extend to 2R_vir. The authors argue that galaxies in the cluster outskirts have begun to quench, but that is not yet reflected in their sizes/morphologies, and they remain SF for 2–4 Gyr after their first cluster infall, in accordance with Wetzel et al. (2013). The f_SF parameter used in this study is sensitive effectively only to completely quenched galaxies. Our result of the radial gradient and our data are equally consistent with two scenarios. The radial gradient of f_SF, flat from the field to R ≲ R₂₀₀ and then strongly decreasing toward the cluster core, is consistent with a delayed environmental quenching scenario in which galaxies remain unaltered after first entering the cluster outskirts and then rapidly quench when approaching the core. However, the observed radial gradient is also consistent with a scenario in which galaxies after entering the outskirts remain SF, but at a lower mass-normalized rate, and completed quenching does not occur until significantly well past the virial radius.

From Figure 5 alone, we cannot definitively deduce whether the extent of the effect of the environment, as indicated by the R/R₂₀₀ at which the f_SF drops, evolves or not. The larger extent of the cluster effect at our lowest redshift could be a sign of environmental quenching evolution, but it also could be caused by the cluster sample selection. In reanalyzing the data using the Lin et al. (2017) cluster sample setup (see Section 4.3.3), the differences between the curves in Figure 6 representing the field and the cluster outskirts still remain for clusters in the two lower-redshift bins (z < 0.77). Further analysis of Figure 5 also gives a hint of a potential dependence of the cluster sphere of influence with galaxies' stellar mass. In the next paper of this series we will investigate such evidence and its implications on the environmental quenching effect.

5.2.2. The Environmental Dependence of ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$

As shown in Figure 7, the stellar masses ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ are lower in the cluster core than in the field at all redshifts. This shift indicates that galaxies in the cluster core have their quenching accelerated over what is expected from mass quenching in the field, i.e., cluster galaxies of certain stellar mass complete their quenching earlier than field galaxies of the same stellar mass. Such a shift is at ~0.3–0.4 in log M_⋆/M_⊙, except for ${M}_{\mathrm{INI}}^{* }$ at the highest-redshift bin, where it is reduced to 0.14, and for the ${M}_{\mathrm{END}}^{* }$ at the lower-redshift bins, where it is increased to 0.8–0.9. This effect produces f_SF versus M_⋆ curves in Figure 6 that are bunched up at the highest-redshift bin and more spread out at lower redshifts. While we observe evolution with redshift for both characteristic masses in all environments, the difference between them, which is a proxy for the slope of the f_SF versus M_⋆ curves, is almost constant for the field and the outskirts of the cluster. For the cluster core, it appears that the f_SF versus M_⋆ curves could be flattening toward lower redshifts, but we are limited by the stellar mass completeness of our two high-redshift bins to arrive at a more definitive result.

5.2.3. The Effect of the Local Galaxy Density

The global environment is sensitive to the cluster-centric radius parameter, which gives the position with respect to the gravitational potential of the parent halo. Thus, R/R₂₀₀ is useful to evaluate the role of mechanisms such as ram pressure stripping, strangulation, and global tidal effects from the cluster dark matter halo mass. In order to have a complete picture of the quenching processes, one also needs to parameterize the local galaxy density to measure the more immediate environment. Such a parameter will be sensitive to mechanisms that are dependent on the existence of neighboring galaxies, e.g., harassment, galaxy–galaxy interactions, and mergers.

While our sample size is large enough to study the effect of the stellar mass and the position within the cluster separately, a similar careful analysis of the local galaxy density, at fixed cluster-centric radii, particularly in the cluster core, requires a larger sample of clusters. Li et al. (2012) studied the evolution of the red fraction in a sample of 905 galaxy groups and clusters with 0.15 ≤ z < 0.52 from the first Red-Sequence Cluster Survey. They analyze the f_red dependence on four parameters—galaxy stellar mass, total group stellar mass, normalized group-centric radius, and local galaxy density—and find that the relative effect of local density on f_red is roughly constant as a function of group-centric radius, suggesting that one can separate the influence between local galaxy density and the global effects of the group environment. Their results also indicate that both local galaxy density and cluster/group-centric radius are significantly correlated with galaxy population. Therefore, even though a detailed study of the effect of local galaxy density is essential to a clearer picture of galaxy evolution, we expect that the conclusions derived from the analysis of the cluster-centric distances are correct.

5.3. The Evolution of the f_SF

The picture that galaxies evolve with redshift has been largely demonstrated since Butcher & Oemler (1984) showed that the fraction of blue galaxies in clusters is higher at earlier epochs (see Section 4.2). Our large sample and the use of a large cluster-centric reference f_SF (i.e., "field") allow us to break down the dependence of the evolution of f_SF in detail by measuring f_SF versus z for different stellar mass and cluster-centric radius bins simultaneously.

Raichoor & Andreon (2012) already showed the benefit of using narrow stellar mass bins when studying the evolution of the blue fraction in clusters. In particular, in their highest stellar mass bin (log M_⋆/M_⊙ > 10.5, the same as in Butcher & Oemler 1984), these authors observed no evolution of the blue fraction up to z = 0.43, and only when including one cluster at z = 1.05 did they measure an increase of the blue population. Our data allow us to study in detail their redshift gap and locate the rise of the f_SF of the most massive galaxies at z > 0.7. Although we cannot make a direct comparison, since their sample of clusters is X-ray selected and their inner cluster bin is up to R/R₂₀₀ = 0.5, our results indicate an increment of the f_SF with redshift and toward lower stellar masses shallower than the extrapolation of their model (left panels in their Figure 13).

Figure 8 shows that the galaxy evolution observed in clusters, reflected by the increase of the blue population toward higher redshifts, occurs as an extension of the galaxy evolution of the field galaxies. This is not surprising considering that clusters of galaxies are continually accreting galaxies from the field and that a non-negligible population of field galaxies is blended into the cluster population. This evolutionary trend is also easily seen from Figure 6. We have shown that the total fraction of the SF population is driven by lower-mass galaxies; however, the stellar mass at which galaxies begin to quench is not constant but evolves to higher values with higher redshift, as demonstrated by ${M}_{\mathrm{INI}}^{* }$ . This means that at a fixed stellar mass range, as when calculating f_SF of a mass-selected sample, we will obtain higher f_SF values for higher redshifts. The f_SF versus z relation may appear to be flat if the stellar mass range in which it is computed is selected such that the upper limit is $\lesssim {M}_{\mathrm{INI}}^{* }$ of the lowest-redshift bin or the lower limit is $\gtrsim {M}_{\mathrm{END}}^{* }$ of the highest-redshift bin. Consequently, if a study is not careful with the stellar mass limit, one may overestimate the effect (e.g., using a blue luminosity limit instead of a stellar mass limit) or not see any evolutionary effect (e.g., if the stellar mass limit is too high for the redshift range studied).

5.4. Environmental and Mass Quenching Efficiencies

5.4.1. Measuring ${\epsilon }_{\mathrm{env}}$ and ${\epsilon }_{\mathrm{mass}}$

Following the approach of previous works, we estimate the environmental quenching efficiency, _env, to explore the degree to which cluster galaxies have suppressed their star formation relative to similar galaxies in the field (van den Bosch et al. 2008; Peng et al. 2010; Phillips et al. 2014). The formula that we employ is ${\epsilon }_{\mathrm{env}}(r,m)=({f}_{Q}(r,m)\,-\,{f}_{Q}({r}_{\mathrm{field}},m))/{f}_{\mathrm{SF}}({r}_{\mathrm{field}},m)$ , which accounts for the fraction of galaxies with stellar mass m at a certain distance r from the cluster in excess of those in the field. As in previous sections, the field environment is defined as between R/R₂₀₀ = 4 and 6. In the same way, the stellar mass quenching efficiency, _mass, is defined as the fraction of galaxies that are quenched compared to the SF population at low stellar mass: ${\epsilon }_{\mathrm{mass}}(r,m)=({f}_{Q}(r,m)-\,{f}_{Q}(r,{m}_{0}))/{f}_{\mathrm{SF}}(r,{m}_{0})$ , where the reference mass m₀ is the lowest stellar mass at which (almost) all galaxies at cluster-centric radius r are forming stars. In practice, we use the stellar mass completeness limit for each redshift bin. While these definitions are usually employed when studying clusters, _env and _mass are usually defined by means of the local galaxy density when analyzing field galaxies. For _env, works in the literature usually compare the top local galaxy density quartile to the bottom quartile; e.g., Kawinwanichakij et al. (2017) use δ₇₅ and δ₂₅ instead of r_cluster and r_field. Similarly, _mass is defined in field studies avoiding galaxies within the densest environments, e.g., using only galaxies in the three bottom local galaxy density quartiles, δ < δ₇₅. These different definitions complicate direct comparisons between different works and the resulting interpretations.

In Figures 9 and 10 we plot the 2D profile of the two quenching efficiencies as a function of galaxy stellar mass and cluster-centric radius for different redshift bins. We note that _env and _mass are degenerate as the fraction of SF galaxies, f_SF(r, m), approaches 0, as both parameters would approach the value of 1. When (almost) all galaxies are quenched, we cannot distinguish which of the two quenching mechanisms, mass or environment, has been more efficient.

**Figure 9.** Mass quenching efficiency as a function of the stellar mass and the cluster-centric distance for four redshift bins, as indicated by the top labels. Purple curves show the contours at fixed _mass = 0.25 and 0.5. Horizontal dashed lines indicate the stellar mass completeness limit of each redshift bin.
Download figure:
Standard image High-resolution image

epsilon — **Figure 9.** Mass quenching efficiency as a function of the stellar mass and the cluster-centric distance for four redshift bins, as indicated by the top labels. Purple curves show the contours at fixed _mass = 0.25 and 0.5. Horizontal dashed lines indicate the stellar mass completeness limit of each redshift bin.
Download figure:
Standard image High-resolution image

**Figure 10.** Environmental quenching efficiency as a function of the stellar mass and the cluster-centric distance for four redshift bins, as indicated by the top labels. Horizontal dashed lines indicate the stellar mass completeness limit of each redshift bin.
Download figure:
Standard image High-resolution image

5.4.2. Are Mass and Environmental Quenching Efficiencies Independent?

In the case of _mass being independent of the environment we expect in Figure 9 horizontal stripes with a higher mass quenching efficiency for more massive galaxies and lower efficiency toward lower stellar masses. We do observe a gradient with stellar mass, but the stripes are not completely flat. This effect is clearer in Figure 11, where the mass quenching efficiency is plotted as a function of R/R₂₀₀ for various stellar mass bins. We can see how _mass appears as flat lines on top of each other as we move toward higher stellar mass bins for environments at R ≳ R₂₀₀, indicating that _mass is independent of R/R₂₀₀. Peng et al. (2010), using galaxy samples from SDSS and zCOSMOS, concluded that _mass is completely independent of the environment, where environment is defined as quartiles in local density. However, as we approach the cluster core, a dependence of ${\epsilon }_{\mathrm{mass}}$ on R/R₂₀₀ becomes evident: _mass is clearly more efficient toward the cluster center. This fact suggests that in very dense areas, such as the cores of clusters, galaxy formation history and evolution may be different from that of the field. We may expect galaxies in the cluster core, on average, to be older than those in the field, since some of them are likely formed "in situ" (e.g., Maulbetsch et al. 2007). These galaxies may be formed earlier, and it is reasonable that they follow a different mass quenching history from that of the "universal clock" of field galaxy mass quenching. Peng et al. (2010) likely did not find this population of galaxies formed in very dense environments because their galaxy density did not sample sufficiently dense environments, or the effect was masked by the large uncertainty in estimating local galaxy density.

**Figure 11.** Mass quenching efficiencies as a function of cluster-centric distance for four redshift bins, as indicated by the top labels, and different stellar mass bins.
Download figure:
Standard image High-resolution image

Looking at Figures 9 and 10, we see how _mass clearly dominates over _env outside R ∼ R₂₀₀. But inside R₂₀₀, _env grows faster than _mass, such that for the intermediate to low stellar mass range (M_⋆ ≲ 10^10.3–10^9.7 M_⊙ depending on the redshift) very close to the cluster center they appear to be equally effective. We also observe how both effects become less efficient with decreasing stellar masses, in agreement with, e.g., Jian et al. (2017). Such a clear dependence of _env on stellar mass has been reported by Kawinwanichakij et al. (2017) at z > 1. These authors found that more massive galaxies experience stronger environmental quenching, and they concluded that the environmental quenching effects are not separable from stellar mass at z ≳ 1. However, at 0.5 < z < 1, they measured a _env nearly independent of stellar mass. Van der Burg et al. (2018) found also that there is no clear stellar mass dependence of _env in any radial bin for a sample of Planck-selected clusters at 0.5 < z < 0.7. However, the stellar mass dependence that we find in this redshift bin is mainly observed below their lowest stellar mass bin, i.e., ≲10^9.7 M_⊙. We also note that their sample of clusters consists of the most massive halos, very likely more massive than the average of the halos we are probing. Papovich et al. (2018), using the same ZFOURGE data set as Kawinwanichakij et al. (2017), studied the stellar mass function of the quiescent population and deduced that the observed shape requires an environmental quenching that depends on stellar mass even for the lowest-redshift bin of 0.5 < z < 1. Our measurement of the environmental quenching is consistent with their deduction and extends to lower redshifts and the very dense environment of the cluster cores.

Our results suggest that, within the cores of clusters, mass and environmental quenching efficiencies, as defined, are interdependent and not separable. However, the environmental influence in the core of clusters is unlikely to arise from simply the much higher local galaxy density, as the gravitational potential of the parent halo must also exert a significant influence. In future work, we will study this in greater detail by isolating the effects of local galaxy density at different cluster-centric radii, as attempted by Li et al. (2012), which may provide key insights into the possible different quenching mechanisms at play.

5.4.3. The Evolution of the Environmental Quenching Efficiency

From the _env displayed in Figure 10, we observe how the influence of the cluster diminishes with increasing redshift. At z > 0.7, _env becomes significant only at R/R₂₀₀ ≲ 0.6. In Figure 12 we show the evolution of the environmental quenching efficiency in the cluster core (R/R₂₀₀ < 0.4) and its stellar mass dependence. By estimating _env in three stellar mass bins, we see how the environment is more effective at quenching more massive galaxies over the whole redshift range, as clearly revealed by the color _env map of Figure 10. Figure 12 shows a strong evolution of the environmental quenching in our redshift range of study, for all stellar masses. We also observe that such evolution of _env is similar for different stellar masses. By definition, _env accounts for the fraction of galaxies that have been quenched by the cluster in excess of those in the field across all time. Hence, the evolution with redshift we measure for _env could just reflect the increased time that galaxies have been in clusters by z ∼ 0.4 compared to z ∼ 1. However, to attribute the observed evolution to this, the fading timescale needs to be longer than the infall timescale. Therefore, the changes we see in _env with redshift likely reflect a true evolution of the environmental quenching efficiency.

**Figure 12.** Environmental quenching efficiency as a function of redshift for different stellar mass bins (red, teal, and indigo regions) in the cores ( $R/{R}_{200}\lt \,0.4$ ) of our cluster sample and various other samples from the literature. The hatched area indicates that the stellar mass is below our stellar mass completeness limit. Black and gray regions show the environmental quenching efficiency of galaxies from Peng et al. (2010) for their highest galaxy density quartile for galaxies with 9.0 < log M_⋆/M_⊙ < 11.0 from SDSS and with 10.2 < log M_⋆/M_⊙ < 11.0 from zCOSMOS, respectively. Red circles show the measurements from Nantais et al. (2017) for galaxies in clusters with log M_⋆/M_⊙ < 10.3, while the red rectangle region shows the result of van der Burg et al. (2013) using the same sample up to z = 1.34. Turquoise circles show the intermediate stellar mass bin (9.8 < log M_⋆/M_⊙ < 10.2) of galaxies in the highest-density quartile from Kawinwanichakij et al. (2017).
Download figure:
Standard image High-resolution image

**Figure 12.** Environmental quenching efficiency as a function of redshift for different stellar mass bins (red, teal, and indigo regions) in the cores ( $R/{R}_{200}\lt \,0.4$ ) of our cluster sample and various other samples from the literature. The hatched area indicates that the stellar mass is below our stellar mass completeness limit. Black and gray regions show the environmental quenching efficiency of galaxies from Peng et al. (2010) for their highest galaxy density quartile for galaxies with 9.0 < log M_⋆/M_⊙ < 11.0 from SDSS and with 10.2 < log M_⋆/M_⊙ < 11.0 from zCOSMOS, respectively. Red circles show the measurements from Nantais et al. (2017) for galaxies in clusters with log M_⋆/M_⊙ < 10.3, while the red rectangle region shows the result of van der Burg et al. (2013) using the same sample up to z = 1.34. Turquoise circles show the intermediate stellar mass bin (9.8 < log M_⋆/M_⊙ < 10.2) of galaxies in the highest-density quartile from Kawinwanichakij et al. (2017).
Download figure:
Standard image High-resolution image

Due to systematic effects such as stellar mass completeness limits, cluster sample selection, and different definitions of environment (e.g., local galaxy density, different R/R₂₀₀ distances), it is difficult to compare our _env curves with those in the literature directly. For field galaxies, Peng et al. (2010) and Kawinwanichakij et al. (2017) show lower _env values and a smoother evolution as illustrated in Figure 12. These studies define environment through local galaxy density and do not focus on galaxy clusters; hence, our ${f}_{\mathrm{SF}}({r}_{\mathrm{cluster}},m)$ estimate is very likely measured in a much denser environment than the highest-density quartile that Kawinwanichakij et al. (2017) used to calculate _env. Conversely, our ${f}_{\mathrm{SF}}({r}_{\mathrm{field}},m)$ reference value is measured as the average between R/R₂₀₀ = 4 and 6, which is likely closer to the 50th percentile in local density than the lowest-density quartile used by Kawinwanichakij et al. (2017). This difference would at least partly compensate for the ${f}_{\mathrm{SF}}({r}_{\mathrm{cluster}},m)$ measured in a denser environment; therefore, if we could use an ${f}_{\mathrm{SF}}({r}_{\mathrm{field}},m)$ comparable to that of Kawinwanichakij et al. (2017), we would expect our data to show an even higher _env. Note that we only include the Kawinwanichakij et al. (2017) intermediate-mass bin in Figure 12. The steeper slopes of our curves compared to the trends of Peng et al. (2010) and Kawinwanichakij et al. (2017) could be interpreted as the acceleration of the quenching produced by the cluster with respect to the field/group environment.

Nantais et al. (2017) analyzed a sample of spectroscopically confirmed SpARCS galaxy clusters at 0.86 < z < 1.65 with log M_⋆/M_⊙ > 10.3, showing a substantial evolution of the environmental quenching efficiency within their redshift range. Their trend is compatible with the strong evolution that we observe (Figure 12), though the _env we measure for our highest stellar mass bin in the overlapping redshift bin at z ∼ 1 is considerably lower. However, our measurements are in better agreement with van der Burg et al. (2013) using the same data as Nantais et al. (2017) at 0.86 < z < 1.34. The large discrepancy is very likely due to how the environmental quenching efficiency is computed, e.g., Nantais et al. (2017) used different stellar mass completeness limits than van der Burg et al. (2013). Furthermore, we define the cluster core within R/R₂₀₀ ≤ 0.4 and use the same data set to account for the field, while Nantais et al. (2017) and van der Burg et al. (2013) consider all galaxies within 1 Mpc, and their field galaxies are from the UltraVISTA/COSMOS catalog.

5.4.4. The Evolution of the Mass Quenching Efficiency

In Figure 13 we focus on the evolution of the mass quenching efficiency in three different environments: the cluster core (R/R₂₀₀ < 0.4), the outskirts of the cluster (0.4 ≤ R/R₂₀₀ < 1), and the field (4 ≤ R/R₂₀₀ < 6). The mass dependence is also analyzed by using two stellar mass bins. As with the environmental quenching, the efficiency of the mass quenching diminishes with higher redshift. We observe how for the two stellar mass bins the strongest evolution in our redshift range happens between z ∼ 0.6 and 1. However, the slope of the lower stellar mass bin is steeper, suggesting that the evolution of _mass is faster for lower stellar masses. The R/R₂₀₀ bins show again that the mass quenching within the cluster environment (R ≲ R₂₀₀) is different from that of the field. While it still depends on the mass of the galaxy, i.e., mass quenching is more efficient for more massive galaxies, galaxies in the cluster core are, on average, in a more advanced stage of mass quenching. The mass quenching dependence on environment within the cluster core and stellar mass is clearly present in our entire redshift range. However, for each mass bin, the slopes of the three R/R₂₀₀ bins are very similar, suggesting that the environment does not affect the evolution of the mass quenching efficiency, consistent with the scenario that there is a population of galaxies formed "in situ" in the cluster cores in earlier time, as described in Section 5.4.2.

Our results for the outskirts and field environments (blue and green hatched curves in Figure 13) are in reasonable agreement with previous works on field galaxies. Our two stellar mass field bins are compatible at lower redshift with results from Peng et al. (2010) and at higher redshift with the high- and intermediate-mass bins of Kawinwanichakij et al. (2017), in both amplitude and evolutionary rate with redshift. The differences can likely be attributed to the slightly different stellar mass limits, or the difference in stellar mass computation.

The interpretation of redshift evolution is complex owing to selection effects, and while the threshold of N_red > 6.0 takes care of the sensitivity to lower-mass clusters with decreasing redshift, by using a fixed mass limit we are not analyzing a progenitor/descendant sample. However, when using the setup of Lin et al. (2017) described in Section 4.3.3, we obtain very similar evolutionary trends to those in Figure 13 (mass-limited sample). This suggests that the effect of a mass-limited cluster sample does not mask the trend with redshift.

5.4.5. Clues to the Quenching Process

Our results show a sharp decline of f_SF at R ≲ R₂₀₀ after a remarkably steady field value, which extends up to 6R/R₂₀₀. The constant f_SF value from the outskirts to as far as 6R/R₂₀₀ indicates that little or no environmental quenching is acting outside the cluster, or if there is an environmental quenching process acting out to 6R/R₂₀₀, it is independent of the distance to the cluster center. The latter could be the case if, e.g., pre-processing in groups or high-density regions is occurring. At lower redshifts, several works have found that the fraction of SF galaxies in clusters is suppressed with respect to the field value even at 3R/R₂₀₀ (e.g., Chung et al. 2011; Haines et al. 2015). The comparison with these works is complex since, for example, the sample selection is different, as we select SF galaxies by the UVJ diagram while they use an SFR_IR threshold. Nevertheless, in the bottom panels of Figure 5 we observe a trend of increasing R/R₂₀₀ in the effect of the environment from high to low redshift, which could be compatible with the environmental effect found at 3R/R₂₀₀ for z < 0.3 clusters. To come to a more definitive conclusion, we need a larger sample that allows us to study the fraction of SF galaxies as a function of the local galaxy density at fixed cluster-centric radii. If f_SF versus local galaxy density shows a correlation that is independent of R/R₂₀₀, at least between R₂₀₀ and 6R₂₀₀, pre-processing would be a suitable mechanism for explaining the results.

More interestingly, within the cluster core, we show (1) that the f_SF versus M_⋆ curve of the cluster is shifted with respect to the field curve, (2) that _mass increases at smaller R/R₂₀₀, and (3) that _env depends on the galaxy stellar mass and evolves with redshift to higher values at later epochs. From condition 1, we deduce that the cluster environment accelerates the quenching, and condition 2 indicates that part of the advanced quenching level in the cluster core may be the outcome of an earlier population formed "in situ" in the primordial overdensities where clusters are located. Additionally, condition 3 may indicate that the efficiency of the stellar-mass-dependent physical process(es) responsible for environmental quenching changes with time, being more efficient at lower redshifts. Strangulation and starvation are environmental mechanisms that can produce such a mass-dependent efficiency as galaxies become quenched after consuming all their cold gas reservoir, once their hot gas supply was cut off when falling into the cluster host halo. For example, in the low-redshift model of Davies et al. (2016), starvation is the favorite model for quenching intermediate-mass galaxies. At higher redshift (z ≳ 1), our results are consistent with the model of Balogh et al. (2016), where the truncation of new gas supply limits the SF lifetime of a galaxy as it exhausts its reservoir through star formation and outflows (the "overconsumption" model; McGee et al. 2014). In such a model, the higher SFRs of more massive galaxies lead to shorter quenching times, leading to an _env that rises with stellar mass. Our results fit in this model and may suggest that, at intermediate redshifts (0.3 < z < 1.1), starvation and strangulation become more efficient at quenching intermediate to more massive galaxies.

6. Summary and Conclusions

We have analyzed the fraction of SF galaxies, f_SF, in a sample of 209 SpARCS clusters at 0.3 $\leqslant$ z_CL ≤ 1.1. The large area and depth of the HSC-SSP data, along with our wide redshift range, have allowed us to study the dependence of the f_SF on galaxy stellar mass and distance to the cluster center and the evolution of these parameters with redshift. Our major results can be summarized as follows:

1.
On a grid of stellar masses and redshifts, we demonstrate that the f_SF decreases from R/R₂₀₀ ∼ 1 toward the cluster core beyond that measured in field galaxies of similar stellar masses and redshifts. More massive cluster galaxies show a decrease at earlier redshifts, while lower-mass cluster galaxies begin such a decrease at lower redshifts.
2.
Using the same data set to measure the evolution of f_SF with redshift in the field and the cluster, we find that the environmental effect produces not only a shift to lower values but also a faster evolution rate. The relative decrease in f_SF between the cluster core and the field is strongly dependent on the galaxy stellar mass and cluster redshift.
3.
We estimate ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ , the fiducial stellar masses at which SF galaxies begin to quench and at which almost all galaxies (80%) are already quenched, respectively, in bins of cluster-centric radii, as a function of redshift. These fiducial masses show a clear evolution in redshift, with lower values toward lower redshifts, indicating that more massive galaxies begin and complete their quenching process earlier. Furthermore, we observe a definitive dependence of ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ on the environment, as defined by the distance to the cluster core, with both ${M}_{\mathrm{INI}}^{* }$ and ${M}_{\mathrm{END}}^{* }$ being lower in the cluster core compared to the field. This strong environmental dependence demonstrates that galaxies of a given mass in clusters begin and complete their quenching of star formation much earlier than their counterparts in the field.
4.
We calculate environmental and mass quenching efficiencies as a function of stellar mass and R/R₂₀₀. We find that outside the cluster core, mass and environmental quenching efficiencies are separable. However, in the cluster core, _mass increases by a significant amount relative to that of the field, suggesting that galaxies in the cluster core are in a more advanced stage of mass quenching. This difference may be the result of an initial population of galaxies formed "in situ" in the cores of the clusters, well ahead of the general population of field galaxies. Furthermore, the environmental quenching in the cluster core is less efficient for lower-mass galaxies. We conclude that, within the cluster core, mass and environmental quenching efficiencies are not separable. We also observe that low-mass galaxies in general have lower _env, suggesting that the effect of environmental quenching has a dependence on galaxy mass even outside the cluster core.
5.
We find that both mass and environmental quenching processes evolve with redshift, becoming more efficient at later epochs. The evolution of the environmental quenching efficiency is similar for different stellar masses but more rapid in the cluster core than for the field/group environments. On the other hand, while the mass quenching efficiency also evolves in a similar manner for galaxies in different environments over our redshift range, the evolution of _mass is more rapid for low-mass galaxies.

In this work, we have used cluster-centric radius as an environmental parameter and have not isolated the effects of local galaxy density simultaneously. Such an investigation can provide key insights into the quenching process and mechanisms. A data set with a larger sample of clusters is required to carry out such an analysis. Such a data set will also allow us to separate the sample into different cluster richnesses, allowing the examination of the dependence of these results on the cluster mass. A larger redshift baseline with deeper data will provide crucial information on the evolution of the more massive cluster galaxies at higher redshifts and the beginning of the quenching of the less massive galaxies in different epochs.

In memory of Prof. Juan-Luis Pintos: sociologist and beloved father. We thank the anonymous referee for the careful reading of the manuscript and providing many useful comments. We also thank Remco van der Burg for valuable discussion and constructive feedback on the manuscript.

This work is based on data collected at the Subaru Telescope and retrieved from the HSC data archive system, which is operated by the Subaru Telescope and Astronomy Data Center at the National Astronomical Observatory of Japan. This work is also based on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory of the California Institute of Technology under NASA contract 1407. This work also uses data obtained at UKIRT, which is owned by the University of Hawaii (UH) and operated by the UH Institute for Astronomy; operations are enabled through the cooperation of the East Asian Observatory. This work is also based in part on observations obtained with MegaPrime/MegaCam, a joint project of CFHT and CEA/DAPNIA, at the Canada–France–Hawaii Telescope (CFHT), which is operated by the National Research Council (NRC) of Canada, the Institut National des Science de lUnivers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii.

The research presented in this paper is supported by grants from the Discovery Grant Program and the Canada Research Chair Program of the National Science and Engineering Council of Canada (NSERC), and the Faculty of Arts and Science of University of Toronto to H.K.C.Y. Further, L.O. acknowledges a European Space Agency (ESA) Research Fellowship. The research of G.W. is supported by the National Science Foundation through grant AST-1517863, by HST program Nos. GO-13677/14327.01 and GO-15294, and by grant No. 80NSSC17K0019 issued through the NASA Astrophysics Data Analysis Program (ADAP). Support for program Nos. GO-13677/14327.01 and GO-15294 was provided by NASA through a grant from the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555.

This research made use of NumPy, SciPy, Astropy, and Matplotlib.

The Evolution of the Quenching of Star Formation in Cluster Galaxies since z ∼ 1

Article metrics

Permissions

Share this article

Author e-mails

Author affiliations

ORCID iDs

Dates

Abstract