The Gemini/HST Galaxy Cluster Project: Environment Effects on the Stellar Populations in the Lynx Clusters at z = 1.27

Galaxy evolution may be studied by investigating the resulting stellar populations and the galaxy structure at z ≈ 0, commonly referred to as galaxy archaeology; see recent studies by Cappellari et al. (2013a), Conroy et al. (2014), McDermid et al. (2015), and Parikh et al. (2019). Alternatively, one may study galaxies over a substantial redshift interval and attempt to establish an evolutionary timeline consistent with the properties, e.g., van Dokkum & van der Marel (2007, and references therein), Saglia et al. (2010), Jørgensen et al. (2017), and Beifiori et al. (2017). In the present paper we focus on the latter technique as it is applied to bulge-dominated galaxies.

Multiple studies have found that since z = 1 − 2 the sizes of field galaxies increase by a factor of 3–5 (Toft et al. 2009, 2012; Newman et al. 2012; Cassata et al. 2013; van der Wel et al. 2014). However, it appears that the evolution is accelerated in cluster environments and may happen at earlier epochs than in the field. Therefore the (remaining) size evolution for cluster galaxies from z ≈ 1.5 to the present is less than found for field galaxies (Papovich et al. 2012; Zirm et al. 2012; Jørgensen & Chiboucas 2013; Strazzullo et al. 2013; Delaye et al. 2014; Jørgensen et al. 2014). The recent results for three z = 1.3–1.6 clusters published by Beifiori et al. (2017) also indicate that size evolution for cluster galaxies is smaller than found in the field. Some of the apparent size growth for passive galaxies may be due to larger more recently quenched galaxies entering the samples. Belli et al. (2015) find that for galaxies with masses above 10^10.7 M_☉ about half the size growth between z = 2 and 1.25 originates from such recently quenched galaxies, while the other half represents the size growth of already quenched galaxies.

Stellar population evolution studies beyond z ≈ 1 have primarily focused on ages through studies of luminosity changes. Beifiori et al. (2017) used new data for 19 galaxies in z = 1.3–1.6 clusters obtained with the Very Large Telescope/KMOS to extend the redshift coverage of the results regarding the evolution of the mass-to-light (M/L) ratios of bulge-dominated passive galaxies. The authors used their new results together with the available literature results covering up to z = 1.3 (van Dokkum & Franx 1996; Jørgensen et al. 1999, 2006, 2014; Kelson et al. 2000; Wuyts et al. 2004; Holden et al. 2005, 2010; Barr et al. 2006; van Dokkum & van der Marel 2007; Saglia et al. 2010; Jørgensen & Chiboucas 2013) and low-redshift reference data for the Coma cluster (Jørgensen 1999; Jørgensen et al. 2006) to further solidify the evidence supporting passive evolution and a formation redshift z_form ≈ 2. The formation redshift should be understood as the epoch of the last major star formation episode. At z ≈ 1 the massive (Mass > 10¹¹ M_☉) bulge-dominated galaxies in clusters appear to be in place and mostly passively evolving. Lower mass galaxies may still be added to the red sequence and from then on passively evolve (e.g., Sánchez-Blázquez et al. 2009; Choi et al. 2014), but see also Cerulo et al. (2016) for results supporting that the red sequence well below L^⋆ is fully populated in rich clusters already at $z\approx 1.5$. Ultimately, the properties of galaxies mapped over a large fraction of the age of the universe, may constrain the models for building the galaxies. It is difficult to understand within the prevailing hierarchical model favored by the ΛCDM (cold dark matter) cosmology, the existence of such massive passive galaxies with relatively old stellar populations at z ≈ 1, while less massive galaxies appear to harbor younger stellar populations, e.g., Jørgensen et al. (2017, and references therein), see Kauffmann et al. (2003) for a discussion of this tension between the observational results and the hierarchical models of galaxy formation. However, more recent cosmological simulations like Illustris (Genel et al. 2014; Vogelsberger et al. 2014; Wellons et al. 2015) and UniverseMachine (Behroozi et al. 2019) find that massive quiescent galaxies can be in place by z ≳ 2.

Metallicities and abundance ratios for bulge-dominated galaxies up to z ≈ 1 have been studied from spectroscopy both through our project, the Gemini/Hubble Space Telescope (HST) Galaxy Cluster Project (GCP) (Jørgensen et al. 2005, 2017; Jørgensen & Chiboucas 2013), and the ESO Distant Cluster Survey (Sánchez-Blázquez et al. 2009). However, few studies have attempted to map the evolution of metallicities and abundance beyond z ≈ 1, and samples with a deep spectroscopy of z > 1 (passive) bulge-dominated galaxies are still very small, e.g., Jørgensen et al. (2014) and Beifiori et al. (2017). Onodera et al. (2015) stacked the spectra of 24 massive galaxies at $z\approx 1.6$ and established the mean age, metallicity, and abundance ratio [α/Fe] from this composite spectrum. The results suggest that the metallicities and abundance ratios of massive galaxies are set already at z ≈ 1.6, and the stellar populations may evolve passively from there to agreement with low-redshift passive galaxies. Kriek et al. (2016) studied a single massive and passive galaxy at z = 2.1 and found an abundance ratio [Mg/Fe] = 0.45, which is significantly higher than in the majority of low-redshift passive galaxies, and consistent with a very short star formation timescale.

Several authors have used lower resolution HST grism spectra to study z = 0.5–2 galaxies. Such recent studies find ages consistent with quiescent galaxies forming a large fraction of their stars at z ≳ 2 (e.g., Fumagalli et al. 2016; Ferreras et al. 2019; Estrada-Carpenter et al. 2019). Metallicities of these galaxies are found to be solar already at z ≈ 2 (Ferreras et al. 2019; Estrada-Carpenter et al. 2019).

In addition to possibly affecting the speed of the size and mass evolution, the cluster environment is also known to affect the star formation histories of the galaxies. Of particular importance for studies of z > 1 galaxies, it appears that even at z = 1.2–1.5 centers of very massive clusters are devoid of star-forming galaxies (Grützbauch et al. 2012; Quadri et al. 2012; Koyama et al. 2013) possibly because the cluster environment has contributed to the quenching of the star formation. However, our results for Lynx W (Jørgensen et al. 2014) show that this z = 1.27 cluster appears to be just in the process of quenching star formation. In addition, $z\approx 2$ protoclusters still contain actively star-forming galaxies, e.g., Wang et al. (2016). Thus, the main transformation must happen at z = 1 − 2 and depends on galaxy mass (Quadri et al. 2012; Koyama et al. 2013) and possibly also on the cluster environment (Tanaka et al. 2013).

In this paper we focus on the structure and stellar populations of the bulge-dominated galaxies in the two clusters Lynx E and W at z = 1.27. Lynx E and W are part of the Lynx supercluster. The clusters were first spectroscopically confirmed by Rosati et al. (1999) and Stanford et al. (1997), while Nakata et al. (2005) carried out an extensive photometry survey of the supercluster identifying not only the two main clusters, but also several smaller groups as part of the supercluster. van Dokkum et al. (2001) used HST imaging to investigate the Lynx W galaxy population in detail. These authors concluded that the cluster contains ≈50% bulge-dominated galaxies and that the two most luminous galaxies in the cluster may have been in recent interactions, or still in ongoing interactions. One of these is the triple-core galaxy in the center of the cluster. In their investigation of a larger part of the Lynx supercluster, Mei et al. (2012) focused on the morphologies and the color–magnitude relation. These authors also concluded that the fraction of bulge-dominated galaxies in the clusters is ≈50%. They identify some galaxies as being quenched before being transformed morphologically to bulge-dominated galaxies. They find that the bulge-dominated galaxies are smaller (at a given mass) than found in lower redshift clusters. On the other hand, Saracco et al. (2014), who also studied sizes and luminosities of Lynx W elliptical galaxies, concluded that the photometric data were consistent with passive evolution and no additional size or mass growth from z = 1.27 to the present.

The Lynx clusters are part of our z = 1 − 2 extension of the GCP. Using our deep optical spectroscopy together with HST imaging, in this paper we investigate to what extent structure evolution is required between z = 1.27 and the present, and whether the stellar populations of the galaxies are consistent with passive evolution from z = 1.27 to the present. In addition, because the two clusters are at identical redshifts but have very different cluster masses, we have the opportunity to study cluster environment effects by direct comparison of the galaxy populations in the two clusters.

The observational data are briefly described in Section 2. We make use of our new low-redshift reference sample (Jørgensen et al. 2018b), see Section 2.3. Appendices A, B, C, and D provide additional details regarding the data. In Section 3, we discuss the properties of the two clusters, in particular the masses compared with masses of other GCP clusters and our low-redshift reference clusters. Section 4 provides an overview of the methods and models used in the analysis. We then define the subsamples of galaxies in Section 5 and evaluate the selection effects for the available sample. The apparent differences between Lynx E and W are presented in Section 6. We investigate the possible structure evolution as well as the evolution of the Fundamental Plane (FP) and the M/L ratios in Section 7. In Section 8, we focus on the stellar populations both as a function of galaxy velocity dispersion and as a function of cluster environment. In particular, we construct composite spectra of samples matching the cluster environments and fit stellar population models to these. In Section 9, we discuss the results in the context of simple models for the evolution with redshift. The conclusions are summarized in Section 10.

Throughout this paper we adopt a ΛCDM cosmology with H₀ = 70 km s⁻¹ Mpc⁻¹, Ω_M = 0.3, and Ω_Λ = 0.7. Magnitudes are quoted as in the AB system (Oke & Gunn 1983), except where noted.

Sections 2.1 and 2.2 briefly summarize the observational data for Lynx E and W. The reader is referred to Appendices A and B, as well as Jørgensen et al. (2014) for details. Section 2.3 summarizes the data for the low-redshift reference sample, see also Appendix D and Jørgensen et al. (2018b), as well as other GCP data used in the analysis.

2.1. Imaging of Lynx E and W

In our analysis we use available HST/ACS imaging in the filters F775W and F850LP, obtained in 2004 through the HST program ID 9919. The three pointings obtained for Lynx E and W cover all galaxies in our spectroscopic sample, see Figure 1. We have reprocessed the HST/ACS imaging since our analysis of Lynx W in Jørgensen et al. (2014). The same methods were used as described in Chiboucas et al. (2009). The images were drizzled (Fruchter & Hook 2002), and then processed with SExtractor (Bertin & Arnouts 1996). For sample selection purposes, we adopted the SExtractor magnitudes mag_auto in F850LP as the total magnitudes z_tot,850, understanding that mag_auto misses a small fraction of the flux (see, e.g., Jørgensen et al. 2018a for detailed simulations of the effect). Aperture colors (i₇₇₅ − z₈₅₀) were derived within an aperture with diameter of 05. Based on point-spread-function (PSF) modeling as described in Chiboucas et al. (2009), we find the median FWHM = 01045 and 01050, for the drizzled stacked images in F775W and F850LP, respectively. Thus, the fixed aperture size may lead to a small systematic error on the colors. We quantified this from simulations matching the observed properties for z₈₅₀ = 21–24.6 mag galaxies, and the PSF and the noise properties of the data. The simulations give an offset in the (i₇₇₅ − z₈₅₀) aperture colors of −0.008 from the true colors. This difference is not significant for our purpose, as the colors are used primarily for selection of targets for the spectroscopic observations. All galaxies in the fields brighter than z_tot,850 = 24.6 mag were then fit with the profile-fitting program GALFIT (Peng et al. 2002). This limit corresponds to a signal-to-noise ratio (S/N) of ≈25 for the available data, resulting in radii and total magnitudes with formal fitting uncertainties of 10% or less. The limit also ensures that our main spectroscopic sample is included in the processing. We fit both r^1/4 profiles and Sérsic (1968) profiles and derive effective radii, total magnitudes, and surface brightnesses. This processing was done for the observations in F850LP only. Throughout the paper we use circularized effective radii derived from the semimajor and semiminor axes as r_e = (a_e b_e)^1/2.

Zoom In Zoom Out Reset image size

The photometry is calibrated to the AB system using the synthetic zero-points for the filters, 25.654 for F775W and 24.862 for F850LP (Sirianni et al. 2005). We adopt galactic extinction in the direction of the cluster as provided by the NASA/IPAC Extragalactic Database using the calibration from Schlafly & Finkbeiner (2011). The extinction in the two filters are ${A}_{775}=0.044$ and A₈₅₀ = 0.034. The photometric parameters derived using SExtractor and GALFIT are for the spectroscopic sample listed in Table 7 in Appendix A. The re-derived photometry is consistent with our measurements in Jørgensen et al. (2014), see Appendix A. For consistency in the current analysis we exclusively use the re-derived photometry. The full photometric catalog for the three fields will be published in a future paper (I. Jørgensen et al. 2020, in preparation).

We calibrate the photometry to rest-frame B-band (Vega magnitudes) using stellar population models from Bruzual & Charlot (2003). We use the improved U, B, and V filter functions from Maíz Apellániz (2006) as described in Jørgensen et al. (2018a). The calibration is given in Appendix A.

2.2. Spectroscopy of Lynx E and W

2.3. The Low-redshift Reference Sample and the z = 0.2–0.9 GCP Clusters

In order to determine the cluster redshifts and velocity dispersions, we first examine the distribution of the galaxies in velocity space. Figure 2(a) shows the relative radial velocity distribution of the spectroscopic sample within ±10,000 km s⁻¹ of the median redshift of the potential cluster members. The Lynx clusters are well isolated in velocity space, with a well-defined peak of galaxies within ±2000 km s⁻¹ of the median redshift. The galaxies present at relative radial velocities of −4400 kms⁻¹ and 6000 km s⁻¹, respectively, are at cluster center distances larger than 1.5 R₅₀₀ and therefore unlikely to be gravitationally bound to the clusters. Based on this, we limit the possible cluster members in redshift space to z = 1.253–1.283.

Zoom In Zoom Out Reset image size

In order to assign the galaxies to either Lynx W or E we need to adopt cluster centers. We consider (1) the BCG, (2) the X-ray centers, or (3) the centroid of positions of member galaxies. The BCGs are the triple-core galaxy in Lynx W and ID 4942 in Lynx E. These and the X-ray centers from Rosati et al. (1999) are marked in Figure 1. Due to limitations of the MOS mode used for the observations, our spectroscopic sample is too sparse in the centers of the clusters to give reliable centroids of the two clusters. Thus, to derive centroids from member galaxy positions we include both our spectroscopic sample and galaxies within three sigma of the red sequence and z₈₅₀ ≤ 24.1 mag. We then proceed to assign each possible member to either Lynx W or Lynx E, depending on its closest angular distance to the centers of the two clusters. It turns out that using the centroids as cluster centers results in the same cluster assignments as if using the BCGs. If we use X-ray centers instead of the BCGs, ID 2653 is assigned to Lynx E instead of Lynx W. All other cluster assignments are unchanged. We proceed using the BCGs as the cluster centers, while in the following commenting on to what extent our results depend on this choice.

We determine the velocity dispersions using the bi-weight method (Beers et al. 1990). The resulting median redshifts and velocity dispersion are listed in Table 1. The uncertainties are derived using a bootstrap method as detailed in Beers et al. Determinations are given for both BCGs as cluster centers and X-ray centers. Figure 2(b) highlights the radial velocity distributions for the two clusters relative to the median redshift of members of both clusters, see Table 10 in Appendix B for membership information of the individual galaxies in the spectroscopy sample. One could argue that it would be better to assign membership based on the closest distance to the cluster centers in units of the cluster radii R₅₀₀. (R₅₀₀ is the radius within which the mean over-density of the cluster is 500 times the critical density at the cluster redshift.) If we take that approach, three galaxies (IDs 2653, 2655, and 2757) would be assigned to Lynx E, for both choices of cluster centers. However, this has no significant effect on the derived cluster redshifts or velocity dispersions. Specifically, for this alternative membership assignment we find velocity dispersions of 646 km s⁻¹ and 675 km s⁻¹ for Lynx E and Lynx W, respectively. We note that Ettori et al. (2009) found a size difference of the clusters of only a factor 1.1, while Pascut & Ponman (2015) found a factor 1.9 difference. Because of the considerable uncertainty in the R₅₀₀ measurements from the literature we choose to proceed with memberships assigned using the closest angular distance to the cluster centers. Formally, the uncertainties on the cluster velocity dispersions are 70–90 km s⁻¹, see Table 1. The possible uncertainty in the cluster assignments contributes systematic uncertainties of $30\,\mathrm{km}\,{{\rm{s}}}^{-1}$ for Lynx W and 80 km s⁻¹ for Lynx E.

Table 1.  Cluster Properties

Cluster	BCGs as Centers			X-Ray Centers
	Redshift	${\sigma }_{\mathrm{cluster}}$	Nmember	Redshift	σcluster	Nmember	L500	M500	R500
		(km s−1)			(km s−1)		(1044erg s−1)	(1014 M☉)	(Mpc)
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
Lynx E	1.2642	${568}_{-90}^{+68}$	15	1.2646	${580}_{-77}^{+65}$	16	3.73	1.65	0.518
Lynx W	1.2706	${693}_{-87}^{+67}$	24	1.2706	${706}_{-79}^{+83}$	23	0.63	0.55	0.359

Note. Column 1: galaxy cluster. Column 2: cluster redshift, adopting BCGs as cluster centers. The uncertainties on the redshifts are 0.001. Column 3: cluster velocity dispersion, adopting BCGs as cluster centers. Column 4: number of member galaxies for which spectroscopy is available, adopting BCGs as cluster centers. Columns 5–7: cluster redshift, velocity dispersion, and number of member galaxies, when adopting X-ray centers. Column 8: X-ray luminosity in the 0.1–2.4 keV band within the radius R₅₀₀. Column 9: cluster mass derived from X-ray data within the radius R₅₀₀. Column 10: radius within which the mean overdensity of the cluster is 500 times the critical density at the cluster redshift.

Download table as:  ASCIITypeset image

In Figure 3, we show the cluster masses versus redshifts for the clusters, together with our GCP cluster sample at z = 0.2–1.0, the low-redshift reference clusters, the z = 1.2–1.6 GCP extension, and for reference the catalog of z < 1 clusters from Piffaretti et al. (2011). The X-ray data for the GCP clusters have been calibrated to consistency with Piffaretti et al. (2011), see Jørgensen et al. (2018a). The same calibrations are used to bring the literature data for the z = 1.2–1.6 GCP clusters, including Lynx E and W, to the same system. The Lynx X-ray data are from Ettori et al. (2004, 2009), Stott et al. (2010, Lynx E only), and Pascut & Ponman (2015). Table 1 lists the average values of the recalibrated X-ray luminosities, radii, and masses for Lynx E and W.

Zoom In Zoom Out Reset image size

The mass, M₅₀₀, of Lynx W is significantly lower than that of Lynx E, and of the z = 0.2–1.0 GCP clusters, while Lynx E has a mass similar to the lowest mass z = 0.2–1.0 GCP clusters. However, due to the expected growth of cluster masses with time (see sample models in Figure 3), both clusters are viable progenitors for clusters of masses similar to those of Coma and Perseus at z ≈ 0.

In Figure 4, we show the cluster velocity dispersions versus the X-ray luminosity and the combination of X-ray masses and radii. The figure includes available data for the GCP clusters and the four low-redshift reference clusters (Perseus, Coma, A2029, and A2142, see Jørgensen et al. 2018b). Lynx W falls on the low luminosity and mass extension of the relations for the more massive GCP clusters, while Lynx E falls significantly below the relations as its velocity dispersion is much smaller than expected from its X-ray properties. It is not clear whether this discrepancy is due to a misestimate of the X-ray properties, given the rather low S/N of the X-ray observations, or selection effects in our spectroscopic sample artificially leading to a non-representative low cluster velocity for the cluster. We inspected the Chandra X-ray imaging of the clusters. There are no X-ray point sources very close to the center of Lynx E, thus there is no reason to expect that the cluster's X-ray luminosity has been overestimated due to contamination by point sources. Jee et al. (2006) and Stanford et al. (2001) give velocity dispersions for Lynx E and W of ${740}_{-134}^{+113}\,\mathrm{km}\,{{\rm{s}}}^{-1}$ and 650 ± 170 km s⁻¹, respectively. The result from Jee et al. is based on weak lensing analysis. Within the uncertainties these results agree with ours. If we treat all spectroscopic members as members of one cluster, we find a cluster velocity dispersion of ${773}_{-90}^{+60}\,\mathrm{km}\,{{\rm{s}}}^{-1}$, which can be considered an upper limit on the cluster velocity dispersion supported by our data. Figure 4 shows this as the point labeled "Lynx E limit."

Zoom In Zoom Out Reset image size

To investigate the issue of the cluster masses further, we assess the cluster richness from the number of bulge-dominated galaxies on the red sequence within R₅₀₀ for each cluster. The ratio of the X-ray masses M₅₀₀ is approximately 3:1 with Lynx E being three times as massive as Lynx W. Figure 5 shows the color–magnitude relation for the full field of Lynx E and W covered by the HST imaging. The luminosities of the galaxies in the two clusters including only galaxies within R₅₀₀, brighter than the limit for our spectroscopic sample, and within 3σ of the red sequence gives a total luminosity of such galaxies a factor 2.5 larger for Lynx E than for Lynx W. The choice of including only galaxies out to R₅₀₀ is not critical. If we instead include galaxies to the Lynx E R₂₀₀ = 1.52 R₅₀₀ distance in both clusters, the ratio is 2. Ultimately, a larger spectroscopic sample of cluster members and deeper X-ray observations would be needed to provide more secure determinations of the cluster velocity dispersions and the cluster X-ray masses. For our purpose it is sufficient to know that, assuming that the luminosity ratio traces the mass ratio, then Lynx E is at least twice, possibly three times, as massive as Lynx W. However, see also the analysis by Jee et al. (2006), which shows weak lensing masses of similar sizes for the two clusters, and X-ray mass estimates supporting that Lynx E is three to four times more massive than Lynx W.

Zoom In Zoom Out Reset image size

Our analysis focuses on the bulge-dominated galaxies in the two clusters, with the main emphasis on the possible structure evolution and testing whether passive evolution of the stellar populations from z = 1.27 to the present is sufficient to reach properties consistent with our low-redshift reference sample. We also address possible differences between the galaxy populations in the two clusters.

To parameterize the properties of the galaxies we use (1) the FP (Djorgovski & Davis 1987; Dressler et al. 1987; Jørgensen et al. 1996), and the relations between masses, sizes, and velocity dispersions; and (2) the absorption lines as a function of galaxy velocity dispersion. The zero-points of the relations are derived separately for Lynx E and W. Finally, we investigate high S/N composite spectra of galaxy populations to further establish possible differences between stellar populations as a function of cluster environment.

In the analysis involving the photometry and structural parameters, we show results based on parameters from the fits with Sérsic profiles. None of our results change significantly, when we instead used parameters from the fits with r^1/4 profiles, and none of the conclusions depend on the choice of profile fitting.

We adopt the same method for establishing the scaling relations and associated uncertainties on slopes and zero-points as we used in Jørgensen et al. (2014, 2017), see also Jørgensen et al. (1996). The fits are determined by minimizing the sum of the absolute residuals. Unless otherwise noted, the minimization is done perpendicular to the relations and the zero-points are median zero-points. The uncertainties are derived using a bootstrap method.

In our analysis we use dynamical masses of the galaxies. We adopt ${{\rm{M}}\mathrm{ass}}_{\mathrm{dyn}}=\beta {r}_{e}{\sigma }^{2}\,{G}^{-1}$, with β = 5 (Bender et al. 1992). The results from Cappellari et al. (2006) show that the approximation gives a reasonable mass estimate. However, see also discussion in Cappellari et al. (2013b) regarding the possible merits of using a value of β dependent on the Sérsic index n_ser. Our results do not depend significantly on whether we use $\beta =5$, or adopt the expression for β (n_ser) from Cappellari et al. (2006). Appendix E contains a comparison of relations and zero-point offsets when using the two different methods of mass determination.

In the analysis we use single stellar population (SSP) models from Maraston & Strömbäck (2011) for a Salpeter (1955) initial mass function (IMF). These authors provide spectral energy distributions (SEDs) for the models. As described in Jørgensen et al. (2014), we establish model values for the indices CN3883 and Hζ_A from the SEDs. In addition we use M/L ratios for very similar models from Maraston (2005). We note that for the rest-frame wavelength interval and ages of interest for our analysis, the SEDs from Maraston & Strömbäck (2011) are almost identical to those from Vazdekis et al. (2010) for solar abundance ratio models. The main difference between the two sets of models is the inclusion of contributions from thermally pulsating asymptotic giant branch (TP-AGB) stars in the Maraston & Strömbäck models, while these stars are not included in the models from Vazkedis et al. However, the TP-AGB stars do not contribute significantly to the flux in the rest-frame wavelength interval of interest here, even for the ≈1 Gyr old stellar populations relevant for our analysis.

Jørgensen et al. (2014) list linear model relations between measurable parameters and the logarithm of the age and the metallicity [M/H], valid for the Maraston & Strömbäck (2011) models. As in that paper, we use these relations to aid our analysis. In particular, we derive the expected changes of the measurable parameters with redshift under the simple assumption of passive evolution of the stellar populations. In these models, it is assumed that after an initial period of star formation the galaxies evolve passively without any additional star formation. The models are usually parameterized by a formation redshift z_form, which corresponds to the approximate epoch of the last major star formation episode. In particular, if the Lynx galaxies and the low-redshift reference sample share a common formation epoch, at a lookback time of t_form, we can then derive that from their age difference in log space, ${\rm{\Delta }}\mathrm{log}\ \mathrm{age}$,

$$\begin{eqnarray}&&{t}_{\mathrm{form}}=\displaystyle \frac{\left({t}_{\mathrm{lookback},\mathrm{Lynx}}{10}^{{\rm{\Delta }}\mathrm{log}\mathrm{age}}-{t}_{\mathrm{lookback},\mathrm{low}-{\rm{z}}}\right)}{({10}^{{\rm{\Delta }}\mathrm{log}\mathrm{age}}-1)}\end{eqnarray} \tag{ 1 }$$

where ${t}_{\mathrm{lookback},\mathrm{Lynx}}$ and ${t}_{\mathrm{lookback},\mathrm{low}-{\rm{z}}}$ are the lookback times for the redshifts of the Lynx sample and the low-redshift reference sample, respectively. With the aid of models, the age difference may be determined from the difference in M/L ratios or line strengths. Using our adopted cosmology, t_form can be converted to the formation redshift z_form.

Finally, we use the parameterization of structural changes due to mergers established in Bezanson et al. (2009). Specifically, for minor mergers

$$\begin{eqnarray}&&{\rm{\Delta }}\mathrm{log}{r}_{e}=2\,{\rm{\Delta }}\mathrm{log}\,{\rm{M}}\mathrm{ass}=-4\,{\rm{\Delta }}\mathrm{log}\,\sigma \end{eqnarray} \tag{ 2 }$$

where r_e is the half-light radius, and σ is the velocity dispersion. For major mergers

$$\begin{eqnarray}&&{\rm{\Delta }}\mathrm{log}\,{r}_{e}={\rm{\Delta }}\mathrm{log}\,{\rm{M}}\mathrm{ass}\end{eqnarray} \tag{ 3 }$$

while the velocity dispersion is unchanged. As done by Bezanson et al. (2009), we assume that the mergers do not involve star formation, i.e., they are dry mergers.

In the absence of mergers, we implicitly assume that the galaxies we observe in the Lynx clusters can be considered progenitors to the galaxies in the low-redshift reference sample. However, the low-redshift reference sample contains galaxies with stellar populations younger than ≈9 Gyr and are therefore too young to have been passive at z = 1.27 (Jørgensen et al. 2018b). This progenitor bias is discussed in detail by van Dokkum & Franx (2001). We also consider progenitor bias for our investigation of the structure evolution originating from newly quenched galaxies being added to the sample of passive galaxies since z = 1.27 and that such galaxies may be larger than the older passive galaxies, see Belli et al. (2015). We return to both effects of progenitor bias in the discussion (Section 9).

In order to ensure consistency between final sample selection for our Lynx analysis in the present paper and the samples used in our analysis of the z = 0.2–0.9 GCP clusters we proceed as follows. Following Jørgensen et al. (2014), we divide the cluster members into subsamples according to available spectroscopic parameters, S/N, Sérsic index, and the strength of the [O ii] emission. As noted above we assign the galaxies to either Lynx E and Lynx W, depending on their proximity to the adopted cluster centers. The samples are listed in Table 2. Our main sample consists of subsamples 4 and 5, which are the bulge-dominated galaxies with EW[O ii] > 5 Å and ≤5 Å, respectively, and spectroscopic S/N ≥ 10 Å⁻¹ in the rest frame of the galaxies. The galaxies with EW[O ii] ≤ 5 Å are considered passive galaxies and in the following referred to as such. The limit in S/N is equivalent to uncertainties on the velocity dispersions of ≤0.125 dex (see Appendix B). The galaxies have total magnitudes of z₈₅₀ ≤ 24.2. Two galaxies, IDs 316 and 4508, listed as sample 6 have significant emission in the high excitation neon lines [Ne v]3426 Å and [Ne iii]3869 Å, most likely originating from AGNs (Schmidt et al. 1998; Mignoli et al. 2013). ID 4508 also coincides with a Chandra X-ray point source. These two galaxies are excluded from the analysis.

Table 2.  Subsamples

No.	Criteria	Lynx E		Lynx W
		N	IDs	N	IDs
1	Only redshift and EW[O ii] from spectra			4	129 197 2655 2834
2	nser < 1.5, $\mathrm{log}\sigma$ meas.			4	145 503 691 707
3	nser ≥ 1.5, $\mathrm{log}\sigma$ meas., S/N < 10	1	4921	3	254 545 1036
4	nser ≥ 1.5, $\mathrm{log}\sigma$ meas., S/N ≥ 10, EW[O ii] > 5 Å	3	2098 4093 4828	5	209 293 365 2653 2953
5	nser ≥ 1.5, $\mathrm{log}\sigma$ meas., S/N ≥ 10, EW[O ii] ≤ 5 Å	10	2309 2416 4095 4363 4593 4860 4926 4942 5009 5037	7	148 392 454 886 1134 2757 2849
6	nser ≥ 1.5, $\mathrm{log}\sigma$ meas., S/N ≥ 10, AGN	1	4508	1	316

Note. Samples 4 and 5 require a log mass ≥ 10.3. The lower mass n_ser ≥ 1.5 galaxy ID 4921 is listed under sample 3. ID 1036 in sample 3 does not have the $\mathrm{log}\sigma$ measured and has no significant [O ii] emission.

Download table as:  ASCIITypeset image

Except for allowing galaxies with lower S/N spectra as part of the analysis, the selection criteria for sample 5 are the same as used for our analysis of the z = 0.2–0.9 GCP clusters (Jørgensen & Chiboucas 2013; Jørgensen et al. 2017). The bulge-dominated galaxies with EW[O ii] > 5 Å (sample 4) are included in the present analysis, as we assume that their fairly weak emission originates from the now decreasing star formation rather than ongoing low level star formation. We want to investigate if these galaxies therefore can passively evolve into galaxies similar to the low-redshift reference sample galaxies. In the following, unless explicitly stated otherwise the figures show only the Lynx galaxies from samples 4 and 5.

For bulge-dominated galaxies within 0.12 mag (3σ) of the red sequence, the spectroscopic sample is 41% complete to z₈₅₀ = 24.2, the magnitude of the faintest galaxies included in the analysis. Comparing the Lynx sample of the bulge-dominated galaxies with spectroscopy to the full sample along the red sequence, we find that a Kolmogorov–Smirnov test gives a ≈30% probability that the two samples are drawn from the same parent distribution in z₈₅₀. In addition, if we offset the Lynx sample of bulge-dominated galaxies with spectroscopy by ${\rm{\Delta }}\mathrm{log}L=-0.75$ to account for the average luminosity evolution, the luminosity distribution of the offset sample and that of the low-redshift reference sample are in agreement with being drawn from the same parent sample.

To evaluate if our sample is biased in sizes or surface brightnesses, Figure 6 shows these two parameters, $\mathrm{log}{r}_{{\rm{e}}}$ and $\mathrm{log}\langle I{\rangle }_{{\rm{e}}}$, versus each other for the bulge-dominated galaxies in our spectroscopic Lynx samples 4 and 5. The figure also includes bulge-dominated galaxies with z₈₅₀ ≤ 24.2 mag and within 0.12 mag (3σ) of the red sequence for Lynx, but without spectroscopy. The Lynx data have been offset by ${\rm{\Delta }}\mathrm{log}\langle I{\rangle }_{{\rm{e}}}\,=-0.75$ to take into account the luminosity evolution. The low-redshift reference sample is overlaid. The distributions of the Lynx samples closely resemble the distribution of the low-redshift reference sample, without any obvious bias in the coverage of both parameters. The spectroscopic sample is limited at effective radii ${r}_{{\rm{e}}}\geqslant 1\,\mathrm{kpc}$, which is similar to the limits of literature samples (e.g., Saglia et al. 2010; Beifiori et al. 2017).

Zoom In Zoom Out Reset image size

In summary, we find no obvious selection effects in luminosities, effective radii, or surface brightnesses, that may bias our results.

Examining the gray-scale figure of the two clusters, Figure 1, it is striking that the core of Lynx E is dominated by passive bulge-dominated galaxies. The core of Lynx W appears to be roughly a 50:50 mix of passive bulge-dominated galaxies and disk galaxies plus emission line bulge-dominated galaxies. We investigate this closer using a phase space diagram of line-of-sight velocities relative to the cluster velocity dispersions, ${v}_{| | }/{\sigma }_{\mathrm{cl}}$, versus the cluster center distances in units of R₅₀₀, see Figure 7. On the figure, the cyan dashed lines denote caustics with constant ${v}_{| | }/{\sigma }_{\mathrm{cl}}\times {R}_{\mathrm{cl}}/{R}_{500}$ values of 0.2, 0.64, and, 1.35. Following Noble et al. (2015), galaxies between the caustics of 0.64 and 1.35 are recently accreted, while those outside 1.35 are still infalling. This is in general agreement with the majority of the low-redshift reference cluster galaxies being inside the 0.64 caustic and therefore accreted early in the formation of the clusters. Both the Lynx W and Lynx E emission line bulge-dominated galaxies are primarily recently accreted or still infalling, while the passive galaxies inside R₅₀₀ are also inside the 0.64 caustic and thus can be assumed to have been accreted earlier. The Lynx E sample contains two post-starburst galaxies (marked with yellow outlines in Figure 7). These are also outside the 0.64 caustic. The Lynx E cluster core is significantly richer in bulge-dominated galaxies on the red sequence than seen for the Lynx W core. To illustrate this we show gray-scale figures of the cluster cores with all bulge-dominated galaxies on the red sequence marked, see Figure 8.

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

As we proceed to the quantitative analysis of the data in the following sections, we highlight in the figures the passive bulge-dominated galaxies located in the cores of the two clusters (marked with black outlines in Figure 7). In Section 8.2, we directly investigate the stellar populations of the composite spectra of these galaxies compared to the rest of the passive bulge-dominated galaxies in the clusters and to the bulge-dominated galaxies with significant [O ii] emission.

All scaling relations relevant for our analysis of the possible structure evolution and the evolution of the FP are summarized in Tables 3 and 4, and shown in Figures 9–11. For all relations, we list in the tables the zero-points for the low-redshift reference sample and Lynx samples relative to the Coma relations from Jørgensen & Chiboucas (2013). We also establish relations based on our new larger low-redshift reference sample described in Section 2.3, and give zero-points for the reference sample and the Lynx samples relative to these relations. All slopes for the relations based on our new larger low-redshift reference sample are consistent within 1σ with the slopes from Jørgensen & Chiboucas (2013). On the figures we show the new relations, except for the FP where we choose to show the relation edge-on and face-on using the same coefficients as in Jørgensen & Chiboucas (2013). This choice was made to make it easier to compare our results visually with other published results since many of these use the coefficients from Jørgensen & Chiboucas (2013). The tables list the results for Lynx cluster assignments based on adopting the BCGs as cluster centers. There are no significant differences in the results if using the X-ray centers for the cluster assignments. Specifically, in all cases the zero-point changes are less than half the uncertainties on the zero-points (derived as $\mathrm{rms}\,{N}_{\mathrm{gal}}^{-1/2}$).

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Table 3.  Structure Relations and the FP

Relation	Low Redshift			Lynx E			Lynx W
	γ	Ngal	rms	γ	Ngal	rms	γ	Ngal	rms
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$\mathrm{log}{r}_{e}=(0.57\pm 0.06)\mathrm{log}\,\mathrm{Mass}+\gamma$ a	−5.836	228	0.21	−5.881	13	0.19	−5.815	12	0.17
$\mathrm{log}{r}_{e}=(0.53\pm 0.04)\mathrm{log}\,\mathrm{Mass}+\gamma$ b	−5.423	228	0.21	−5.470	13	0.19	−5.407	12	0.17
$\mathrm{log}\sigma =(0.26\pm 0.03)\mathrm{log}\,\mathrm{Mass}+\gamma$ a	−0.602	228	0.10	−0.578	13	0.09	−0.615	12	0.09
$\mathrm{log}\sigma =(0.28\pm 0.02)\mathrm{log}\,\mathrm{Mass}+\gamma$ b	−0.826	228	0.10	−0.805	13	0.09	−0.840	12	0.09
$\mathrm{log}\,{{\rm{r}}}_{e}=(1.30\pm 0.08)\mathrm{log}\sigma -(0.82\pm 0.03)\mathrm{log}\langle I{\rangle }_{e}+\gamma$ a	−0.429	228	0.09	0.154	13	0.21	0.174	12	0.21
$\mathrm{log}\,{{\rm{r}}}_{e}=(1.23\pm 0.07)\mathrm{log}\sigma -(0.84\pm 0.02)\mathrm{log}\langle I{\rangle }_{e}+\gamma$ b	−0.219	228	0.09	0.374	13	0.21	0.402	12	0.22
$\mathrm{log}\,M/L=(0.24\pm 0.03)\mathrm{log}\,\mathrm{Mass}+\gamma$ a	−1.818	228	0.12	−2.537	13	0.25	−2.563	12	0.26
$\mathrm{log}\,M/L=(0.27\pm 0.02)\mathrm{log}\,\mathrm{Mass}+\gamma$ b	−2.162	228	0.12	−2.894	13	0.24	−2.922	12	0.25
$\mathrm{log}\,M/L=(1.07\pm 0.12)\mathrm{log}\sigma +\gamma$ a	−1.592	228	0.11	−2.302	13	0.24	−2.308	12	0.27
$\mathrm{log}\,{\text{}}M/L=(0.96\pm 0.10)\mathrm{log}\sigma +\gamma$ b	−1.356	228	0.11	−2.059	13	0.24	−2.076	12	0.27

Notes. Column 1: scaling relation. Column 2: zero-point for the low-redshift sample (average of the zero-points for Perseus and Coma). Column 3: number of galaxies included from the low-redshift sample. Column 4: rms in the Y-direction of the scaling relation for the low-redshift sample. Columns 5, 6, and 7: zero-point, number of galaxies, rms in the Y-direction for the Lynx W sample. Columns 8, 9, and 10: zero-point, number of galaxies, rms in the Y-direction for the Lynx E sample.

^aSlopes adopted from Jørgensen & Chiboucas (2013), zero-points use effective radii (and surface brightnesses) from fits with Sérsic profiles, r_e = (a_e b_e)^1/2. ^bSlopes from best fit to Coma and Perseus samples. Effective radii (and surface brightnesses) from fits with Sérsic profiles, r_e = (a_e b_e)^1/2.

Download table as:  ASCIITypeset image

Table 4.  The FP and Relations for the M/L Ratios

Cluster	Relationa	${N}_{\mathrm{gal}}$	rms
Comaa	$\mathrm{log}\,{{\rm{r}}}_{e}=(1.30\pm 0.08)\mathrm{log}\sigma -(0.82\pm 0.03)\mathrm{log}\langle I{\rangle }_{e}-0.443$	105	0.08
Coma, Perseusb	$\mathrm{log}\,{{\rm{r}}}_{e}=(1.23\pm 0.07)\mathrm{log}\sigma -(0.84\pm 0.02)\mathrm{log}\langle I{\rangle }_{e}-0.219$	228	0.09
Lynx E+W	$\mathrm{log}\,{{\rm{r}}}_{e}=(0.55\pm 0.39)\mathrm{log}\sigma -(0.61\pm 0.11)\mathrm{log}\langle I{\rangle }_{e}+1.120$	25	0.14
RXJ0152.7–1357, RXJ1226.9+3332, Lynx E+Wc	$\mathrm{log}\,{{\rm{r}}}_{e}=(0.65\pm 0.21)\mathrm{log}\sigma -(0.65\pm 0.09)\mathrm{log}\langle I{\rangle }_{e}+\{1.000,1.000,1.060\}$	74	0.09, 0.11, 0.15
Comaa	$\mathrm{log}\,M/L=(0.24\pm 0.03)\mathrm{log}\,\mathrm{Mass}-1.754$	105	0.09
Coma, Perseusb	$\mathrm{log}\,M/L=(0.27\pm 0.02)\mathrm{log}\,\mathrm{Mass}-2.162$	228	0.12
Lynx E+W	$\mathrm{log}\,M/L=(0.70\pm 0.16)\mathrm{log}\,\mathrm{Mass}-7.580$	25	0.23
RXJ0152.7–1357, RXJ1226.9+3332, Lynx E+Wc	$\mathrm{log}\,M/L=(0.58\pm 0.09)\mathrm{log}\,\mathrm{Mass}-\{6.163,6.136,6.222\}$	74	0.13, 0.17, 0.22
Comaa	$\mathrm{log}\,M/L=(1.07\pm 0.12)\mathrm{log}\sigma -1.560$	105	0.11
Coma, Perseusb	$\mathrm{log}\,M/L=(0.96\pm 0.10)\mathrm{log}\sigma -1.356$	228	0.11
Lynx E+W	$\mathrm{log}\,M/L=(2.38\pm 0.89)\mathrm{log}\sigma -5.281$	25	0.30
RXJ0152.7–1357, RXJ1226.9+3332, Lynx E+Wc	$\mathrm{log}\,M/L=(2.25\pm 0.26)\mathrm{log}\sigma -\{4.775,4.787,4.989\}$	74	0.12, 0.19, 0.29

Notes.

^aRelations for Coma adopted from Jørgensen & Chiboucas (2013). ^bVelocity dispersions from Jørgensen et al. (2018b), photometry as described in Section 2.3. ^cRXJ0152.7–1357, RXJ1226.9+3332, and Lynx E+W fit with parallel relations. The zero-points and rms for the three samples are listed in the same order as the clusters.

Download table as:  ASCIITypeset image

7.1. Sizes and Velocity Dispersions

7.2. The FP

We can gain a more detailed understanding of the stellar populations of the Lynx galaxies by using the line indices, Hζ_A, CN3883, and D4000, together with the velocity dispersions and the M/L ratios. We first establish scaling relations between the indices and the velocity dispersions and compare the data to SSP models (Section 8.1). Then in Section 8.2 we assemble the spectra into three composite spectra and fit these with SSP models.

8.1. Line Strengths and M/L Ratios

Figure 12 shows the absorption line indices, Hζ_A, CN3883, D4000, and the M/L ratios versus velocity dispersion and versus each other. Table 5 lists the relations shown in panels (a)–(c). On the other panels we overlay SSP models based on SEDs from Maraston & Strömbäck (2011) and M/L information for similar SSP models from Maraston (2005). All models use solar abundance ratios.

Zoom In Zoom Out Reset image size

Table 5.  Scaling Relations for Line Indices

Relation	Low Redshift			Lynx E			Lynx W
	γ	Ngal	rms	γ	Ngal	rms	γ	Ngal	rms
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$\mathrm{log}\ {\rm{H}}{\zeta }_{{\rm{A}}}=(-0.82\pm 0.39)\mathrm{log}\sigma +\gamma$	1.919	213	0.16	2.379	13	0.21	2.338	9	0.19
$\mathrm{CN}3883=(0.18\pm 0.01)\mathrm{log}\sigma +\gamma$	−0.177	217	0.03	−0.290	13	0.08	−0.269	9	0.03
${\rm{D}}4000=(0.42\pm 0.04)\mathrm{log}\sigma +\gamma$	1.127	209	0.09	0.906	13	0.20	1.085	9	0.24

Note. Column 1: scaling relation. Column 2: zero-point for the low-redshift sample (average of the zero-points for Perseus and Coma). Column 3: number of galaxies included from the low-redshift sample. Column 4: rms in the Y-direction of the scaling relation for the low-redshift sample. Columns 5, 6, and 7: zero-point, number of galaxies, rms in the Y-direction for the Lynx W sample. Columns 8, 9, and 10: zero-point, number of galaxies, rms in the Y-direction for the Lynx E sample.

Download table as:  ASCIITypeset image

The Hζ_A-velocity dispersion relation for the Lynx sample is offset from the relation for the low-redshift reference sample, in agreement with expectations for stellar populations with ages of 1–3 Gyr. Similarly as also found in Section 7.2, the M/L ratios for the Lynx sample support ages of 1–3 Gyr. In particular, the majority of the Lynx galaxies agree with the 1–3 Gyr model locations in the Hζ_A–M/L diagram (Figure 12(f)). The RXJ0152.7–1357 and RXJ1226.9+3332 samples follow similar trends with slightly weaker Hζ_A indices as expected as their redshifts are lower than the Lynx redshift.

The CN3883-velocity dispersion relation for the Lynx sample is significantly offset from the low-redshift reference sample, while the RXJ0152.7–1357 and RXJ1226.9+3332 samples show no significant offsets, see Figure 12(b). Our results for RXJ0152.7–1357 and RXJ1226.9+3332 were originally published in Jørgensen et al. (2005) and Jørgensen & Chiboucas (2013). Additional results for the z = 0.2–0.9 GCP clusters can be found in Jørgensen et al. (2017). A small number of the Lynx galaxies, primarily the passive galaxies in the centers of the two clusters (shown with black outlines on the points in Figure 12) fall closer to the low-redshift relation than the median relation for the bulk of the Lynx galaxies. These galaxies still have strong Hζ_A and low M/L ratios in agreement with ages of 1–3 Gyr (Figures 12(a) and (f)). Thus, their unusually strong CN3883 may reflect very high [CN/Fe] or more broadly high total metallicity [M/H]. Due to the short wavelength coverage of our spectra, we cannot evaluate whether only [CN/Fe] or also the total metallicity [M/H] is high relative to the remainder of the sample.

The D4000–velocity dispersion relation, Figure 12(c), for the Lynx sample has a very large scatter. Eight of the 22 Lynx galaxies have D4000 as strong or stronger as the bulk of the galaxies in the low-redshift reference sample. The remainder of the Lynx galaxies have weaker D4000 for a given velocity dispersion. While D4000 is usually used as an indicator of stellar population age (e.g., Gallazzi et al. 2014), we note that the blue passband for D4000 covers the CN3883 absorption feature. Thus, D4000 and CN3883 can be expected to be tightly correlated, as shown in Figure 12(e). The oldest Maraston & Strömbäck models predict stronger D4000, for a given CN3883, than seen at low redshift, but nevertheless the models illustrate that variations in age and metallicity are completely degenerate in the two indices over the parameter ranges populated by the galaxies in our samples. See also our discussion in Jørgensen & Chiboucas (2013) regarding the ability of the Maraston & Strömbäck models to correctly model the CN features.

Turning to Figures 12(g) and (h), it is clear that the strong CN3883 and D4000 values for part of the Lynx sample and the majority of the RXJ0152.7–1357 and RXJ1226.9+3332 galaxies are unassociated with a change in the M/L ratio, and that they are also not well modeled with the Maraston & Strömbäck models. A possible explanation may be that these models assume solar abundance ratios, while the quiescent massive galaxies are known to have above solar α-element, carbon, and nitrogen abundances, see Conroy et al. (2014) and Jørgensen et al. (2017) and references therein. Detailed stellar population models giving both M/L ratios and SEDs, from which we can derive model CN3883 and D4000 values, are not yet available for such element-enhanced models.

8.2. Composite Spectra

To investigate further in particular the strong CN3883 indices, we derive composite spectra by stacking the available Lynx spectra according to their environment as well as the presence of significant [O ii] emission. Based on the phase space diagram, Figure 7 discussed in Section 6, we establish three composite spectra as follows. The cores composite is made up of all passive galaxies inside R₅₀₀ in either Lynx E or Lynx W. The outskirts composite is made up of all passive galaxies outside the cluster cores. Finally, the emission composite is made up of all bulge-dominated galaxies with [O ii] > 5 Å, excluding the two AGNs. As noted in Section 6 the galaxies included in both the outskirts composite and the emission composite populate the phase space diagram in areas that are either recently accreted by the clusters or still infalling. We repeat the determination of the composite spectra and the following analysis also for assignments using the X-ray cluster centers. In this case, ID 545 is assigned to the outskirts composite instead of the cores composite, and ID 148 is assigned to the cores composite instead of the outskirts composite. As the effects of this alternative assignment are minimal, the figures only show the results for assignments using the BCGs as cluster centers, while we comment on results from both in the text. We also establish a high S/N composite spectrum of SDSS spectra of 60 galaxies in the Coma cluster spanning the same range in velocity dispersion ($\mathrm{log}\sigma =2.1\mbox{--}2.4$) as the Lynx sample. Figure 13 shows all of these composite spectra, as well as the ratios between the Lynx composites and the Coma cluster composite.

Zoom In Zoom Out Reset image size

Except for the presence of the [O ii] emission, the emission composite appears similar to the outskirts composite. The equivalent width of [O ii] in the emission composite is 14.4 ± 0.2 Å, which we convert to a star formation rate (SFR) using the average luminosity of the galaxies in the composite (${L}_{B}={10}^{11.08}{L}_{B,\odot }$) and the calibrations from Kennicutt (1992) and Gallagher et al. (1989), see also Jørgensen et al. (2014). We find SFR = 3.39 ± 0.05 M_☉ yr⁻¹. The emission composite exclude the two galaxies with obvious neon emission presumably originating from AGNs in the galaxies, see Section 5. However, we cannot rule out that lower level AGN emission contributes to the strength of the [O ii] emission in the emission composite. Therefore the SFR should be taken as an upper limit. Given the limit on the SFR, the emission in the emission composite may originate from weak residual star formation, which has little influence on bulk of the stellar population, and/or from weak AGNs.

The main difference between the cores composite and the outskirts composite is the strength of the absorption within the CN3883 passband. Relative to the Coma composite, representing the low-redshift reference, the outskirts composite has very weak CN3883, as also seen from the measurements of the individual line indices. The cores composite on the other hand has stronger CN3883 approaching the absorption strength in the Coma composite. Thus, the composite spectra confirm our result from the individual line index measurements.

To quantify the origin of the difference and also explore any differences between the emission composite and the outskirts composite, we fit all three composite spectra with linear combinations of a set of SSP models from Maraston & Strömbäck (2011). The fitting was done with the same kinematics fitting software used in our determinations of the velocity dispersions, see Gebhardt et al. (2000, 2003) for a description of the software. We use SSP models with a Salpeter (1955) IMF. We use only models with ages between 1 and 3 Gyr. The absence of strong emission lines limit the ages to 1 Gyr or older and is in agreement with our results based on the M/L ratios of the galaxies. The upper limit on the ages is chosen to accommodate formation redshifts of up to ≈3.5 and is also in agreement with our results from the strength of Hζ_A and the M/L ratios (see Figure 12). In reality the composite spectra are fully described with only four models. Figure 14 shows the data and the best fits, while the linear compositions of the fits are summarized in Table 6. The uncertainties are determined using Monte Carlo simulations in which we add 1σ additional noise to the composite spectra. In each case, 100 realizations of such spectra are used, and the uncertainties are quoted as the rms scatter of the results from these realizations.

Zoom In Zoom Out Reset image size

Table 6.  SED Fits to Composite Spectra

		[M/H] = −0.3		[M/H] = 0.3
Composite	Cluster Centers	1 Gyr	3 Gyr	1 Gyr	3 Gyr	χ2
(1)	(2)	(3)	(4)	(5)	(6)	(7)
Cores	BCGs	0.38 ± 0.05	0.00 ± 0.00	0.11 ± 0.08	0.51 ± 0.03	2.6
Outskirts	BCGs	0.33 ± 0.07	0.05 ± 0.08	0.62 ± 0.14	0.00 ± 0.02	2.6
Emission	⋯	0.50 ± 0.06	0.27 ± 0.05	0.23 ± 0.11	0.00 ± 0.04	2.5
Cores	X-ray	0.28 ± 0.04	0.06 ± 0.04	0.12 ± 0.08	0.54 ± 0.02	3.1
Outskirts	X-ray	0.40 ± 0.06	0.16 ± 0.13	0.43 ± 0.13	0.00 ± 0.02	2.2

Note. The fits are listed as fractions of luminosity originating from each SSP model. Cluster centers based on BCGs or X-ray centers as noted for the cores composite and outskirts composite.

Download table as:  ASCIITypeset image

The cores composite is dominated by stellar populations of super-solar metallicity ([M/H] = 0.3), mostly with ages of 3 Gyr, combined with a ≈40% contribution by luminosity from a stellar population of subsolar metallicity ([M/H] = −0.3) and an age of 1 Gyr. The outskirts composite can be modeled almost completely by 1 Gyr old stellar populations of both super-solar and subsolar metallicity. The emission composite is modeled primarily by the models with subsolar metallicity, with only about a quarter of the flux originating from a 1 Gyr old super-solar stellar population. In Section 9 we will address to what extent these mixes of stellar populations may evolve to match our low-redshift reference sample spectra in the time available. Table 6 also lists the results for the cores composite and outskirts composite if using the X-ray centers for assignments to the composites. In this case a larger fraction of the outskirts composite best fit has subsolar metallicity. However, the main result is unchanged: the outskirts composite is dominated by stellar populations only 1 Gyr old, while the cores composite contains a ≈50% contribution from 3 Gyr old stellar populations with super-solar metallicity.

There are two main themes in our discussion of the results: (1) to what extent do the data support or allow structure evolution (sizes and velocity dispersions), and (2) can the stellar populations evolve passively from z = 1.27 to the present and be consistent with the observed stellar populations in the low-redshift reference sample.

The cores of Lynx E and W are very different, Lynx E resembling the cores of very massive clusters at similar redshifts, e.g., RDCS J1252.9–2927 at z = 1.24 (Nantais et al. 2013), and XMMU J2235.2–2557 at z = 1.39 (Grützbauch et al. 2012), dominated by passive bulge-dominated galaxies, while Lynx W lacks a well-defined core of such galaxies. However, the passive bulge-dominated galaxies within R₅₀₀ of the Lynx W center appear similar to those in the Lynx E core. Thus, in the following we focus on the evolution of the galaxy sample as a whole, rather than the difference between the two clusters.

9.1. Structure Evolution

9.2. Stellar Population Evolution Measured by the Scaling Relations

In Figure 15 we show the zero-point differences between the z > 0.2 GCP clusters, including the two Lynx clusters, and the low-redshift reference sample, for the main scaling relations. The zero-point differences for the z = 0.2–0.9 GCP clusters are fully consistent with our results in Jørgensen & Chiboucas (2013) and Jørgensen et al. (2014, 2017), but for consistency have been recalculated using the scaling relations listed in Tables 3 and 5 in the present paper. The zero-point differences for the Lynx data are shown for Lynx E and W separately, though there is no significant difference between the results for the two clusters. Predictions for passive evolution based on the models from Maraston (2005) and Maraston & Strömbäck (2011) are overlaid in Figure 15. The changes in M/L ratios and the Hζ_A indices are both consistent with formation redshifts of z_form ≈ 1.9. When comparing our results to the literature, we focus on the highest mass galaxies M_dyn ≥ 10^11.5 M_☉, for which the M/L ratios give ages of about 3 Gyr and z_form ≈ 3.7. This is in agreement with recent results from Beifiori et al. (2017) who based on galaxies in three z = 1.4–1.6 clusters found ages of 1.2–2.3 Gyr and formation redshifts of z_form = 2.3–3.

Zoom In Zoom Out Reset image size

Turning to the CN3883 index, Figure 15(c), it is clear that the index is not well modeled by the passive evolution model based on the Maraston & Strömbäck (2011) SSP SEDs. It is not clear if this is a shortcoming of the models, or due to more complex mixtures of the stellar populations being present and not well modeled by the SSPs. While in principle the Lynx cluster data on average are consistent with a formation redshift of z_form = 4, as we saw in Section 8.1, the sample contains many galaxies with CN3883 indices as strong as in the low-redshift reference sample. There are no other studies of this index at $z\gt 1$. Kriek et al. (2016) studied a high mass z = 2.1 galaxy, COSMOS-11494, based on Keck/MOSFIRE spectroscopy and found a super-solar abundance ratio of [Mg/Fe] = 0.45. A comparison of the spectrum from Kriek et al. to our composite spectra in the wavelength interval of the CN3883 index shows that COSMOS-11494 has a CN3883 feature stronger than the Lynx outskirts composite spectrum, but not as strong as the cores composite. Thus, the CN3883 feature supports the super-solar abundance ratio of this galaxy, and also hints that strong CN3883 features may be present at very early epochs in some massive galaxies. Onodera et al. (2015) measured the CN₂ index from their stacked spectrum of passive z ≈ 1.6 galaxies. Converting the very weak CN₂ index measured by these authors to CN3883 using our empirical conversion based on z = 0.5–0.9 galaxies (Jørgensen & Chiboucas 2013) gives CN3883 ≈ 0.13. This is similar to our results for the Lynx galaxies outside the cluster cores. The galaxies in the Onodera et al. sample may reside in groups or clusters, evaluating from these authors' redshift information for the galaxies. However, it is not possible for us to evaluate how massive these clusters or groups may be.

9.3. Stellar Population Evolution Probed by the Composite Spectra

To further evaluate to what extent passive evolution is sufficient to turn the stellar populations of the bulge-dominated Lynx galaxies into stellar populations similar to those in the low-redshift reference sample, we test whether the SED fits to the Lynx composite spectra can evolve into models compatible with the Coma composite spectrum, in the time available between z = 1.27 and the present. We first convert the luminosity fractions from our best fits, Table 6, into mass fractions using the M/L ratios for the models from Maraston (2005). Then under the assumption of mass conservation we derive the model spectra of the best fits aged by 8 Gyr. Finally, we attempt to fit the Coma composite spectrum with each of the three such evolved best-fit spectra. The result is shown in Figure 16. From panels (a) and (b) it is clear that after 8 Gyr of passive evolution, we can no longer detect any significant differences between what originated from the cores composite and the outskirts composite. This is a consequence of the stellar population models being quite insensitive to small relative age differences once the ages are above about 8–10 Gyr. On the other hand, the evolved emission-composite model does not fit the Coma composite well, showing significant residuals in the vicinity of the CN3883 absorption feature, see Figure 16(c). Thus, for the galaxies included in the emission composite passive evolution is not sufficient for the galaxies to become similar to the Coma galaxies. The main problem is that the metallicity is too low. For this composite to result in a spectrum consistent with the Coma composite, the star formation needs to continue in order to recycle the metals into new stars such that the metallicity increases. However, such a scenario may be difficult to accommodate considering how low the SFR is already (3.4 M_☉ yr⁻¹) and that the galaxies are in the process of falling into the clusters and presumably completing the quenching process before they cross the cores of the clusters. For a galaxy at R₂₀₀ = 1.52 R₅₀₀ in these clusters to reach the cluster core typically takes 0.9 Gyr. For simplicity we will assume that the SFR stays at 3.4 M_☉ yr⁻¹ producing another $3\cdot {10}^{9}{M}_{\odot }$ of new stars. However, this is only 1%–10% of the mass of the already existing stars and appears insufficient to change the luminosity weighted metallicity enough to reach those of the Coma cluster galaxies. A more complex star formation history involving an increased SFR followed by (rapid) quenching would be needed in order to reach the luminosity weighted metallicity of the passive Coma cluster galaxies without the resulting luminosity weighted ages becoming too low. Any prolonged star formation in these galaxies would also be in conflict with their possible evolution into the largest passive galaxies in the clusters as proposed in Section 9.1.

Zoom In Zoom Out Reset image size

We have used deep ground-based optical spectroscopy from Gemini North and HST/ACS imaging to investigate the structure and stellar populations of galaxies in Lynx E and W at redshift z = 1.27, both part of the Lynx Supercluster. Our main conclusions are as follows:

• 1.  
  The Lynx E cluster core is similar to the cores of RDCS J1252.9–2927 and XMMU J2235.2–2557 dominated by bulge-dominated passive galaxies and void of star-forming galaxies. Lynx W lacks a well-defined core of bulge-dominated passive galaxies and appears to contain a significant fraction of emission line galaxies. However, a closer inspection of the phase space diagram shows that the emission line galaxies in general are either recently accreted or still infalling.
• 2.  
  At a given dynamical mass, the galaxies in the Lynx clusters show only a very small difference in size and velocity dispersion when compared to our low-redshift reference sample. However, to produce BCGs similar to those seen in massive low-redshift clusters, the central galaxies will have to grow by at least a factor of 5, possibly through major merging with available nearby galaxies of similar sizes and masses. Alternatively, large and massive galaxies with current star formation (e.g., emission line galaxies as present in the Lynx W sample) may become quenched and enter the passive population by a later epoch. Such newly quenched galaxies may also merge with passive galaxies in the cluster centers, though in some cases their masses and sizes are already similar to those of low-redshift BCGs.
• 3.  
  The bulge-dominated galaxies in the Lynx clusters populate an FP similar to that seen for lower redshift galaxies. However, the slope is steeper; in the M/L–mass plane it is as steep as for the z = 0.9 cluster galaxies. While the data are consistent with passive evolution with a median formation redshift of ${z}_{\mathrm{form}}={1.89}_{-0.10}^{+0.14}$, the steeper slope means that the low mass galaxies may be as young as 1 Gyr, while the highest mass galaxies have stellar population ages of about 3 Gyr. The ages are also supported by the strong Hζ_A absorption line indices.
• 4.  
  The passive galaxies in the cores of the two Lynx clusters have significantly stronger CN3883 indices than found for the passive galaxies in the outskirts of the clusters. Using composite spectra, the cores composite is well modeled by approximately a 60:40 combination of a super-solar metallicity and 3 Gyr old stellar population and a subsolar metallicity and 1 Gyr old stellar population. The outskirts composite on the other hand can be modeled almost completely by 1 Gyr old stellar populations mixed 60:40 in super-solar metallicity and subsolar metallicity. Thus, the galaxies must have gone though different star formation histories. Evolving the best-fit models passively by 8 Gyr shows that both models can then equally well fit the spectra of the Coma cluster galaxies.
• 5.  
  The composite of bulge-dominated galaxies with emission lines, the emission composite, is best fit with approximately a 75:25 combination subsolar and super-solar metallicity models. Evolving such a model by 8 Gyr shows that it cannot satisfactorily fit the Coma cluster galaxies. Thus, ongoing star formation is needed in these Lynx galaxies for them to evolve to consistency with Coma cluster stellar populations. However, the current level of star formation as measured from the [O ii] emission seems to be insufficient to do so before the star formation may be quenched as the galaxies fall into the cluster core. In addition, weak AGN emission may contribute to the [O ii] strength, making the contribution from current star formation even lower.

Further insight into the stellar populations of the Lynx bulge-dominated galaxies would benefit greatly from deep near-infrared spectroscopy enabling measurements of strengths of the iron and magnesium absorption lines. This would enable us to study the total metallicity as well as abundance ratios in greater detail and possibly understand to what extent these galaxies can evolve passively into galaxies similar to those in our low-redshift reference sample.

Karl Gebhardt is thanked for making his kinematics software available. Kristi Webb and Charity Woodrum are thanked for their contributions to the calibration of the low-redshift reference sample photometry. Ricardo Demarco is thanked for comments on a draft version of this paper. The paper makes use of simulations software produced by Riley Peterson (https://github.com/rileypeterson/Gemini-Python-Scripts). The anonymous referee is thanked for comments and suggestions that helped improve this paper.

This work is based on observations obtained at the Gemini Observatory (processed using the Gemini IRAF package), which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), the National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnología e Innovación Productiva (Argentina), Ministério da Ciência, Tecnologia e Inovação (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea). The data presented in this paper originate from the Gemini programs GN-2014B-Q-22 and GN-2014B-DD-4.

In part, based on observations made with the NASA/ESA HST, obtained from the data archive at the Space Telescope Science Institute.

The paper makes use of photometry and spectroscopy from SDSS. Funding for SDSS-IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High-Performance Computing at the University of Utah. The SDSS website is www.sdss.org.

Data from the Chandra X-ray Observatory archive were used to aid the interpretation presented in this paper.

I.J., C.O., and R.C. acknowledge support from grant HST-AR-13255.01 from STScI. STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555.

The Cosmic Dawn Center (DAWN) is funded by the Danish National Research Foundation. S.T. acknowledges support from the ERC Consolidator Grant funding scheme (project Context, grant No. 648179).

Table 7 lists the photometric parameters for the spectroscopic sample as derived from the HST/ACS observations in F850LP and F775W. The F850LP images were processed to derive two-dimensional surface photometry using GALFIT (Peng et al. 2002). F775W was used for the color determinations, only. The effective radii in Table 7 are derived from the semimajor and -minor axes as r_e = (a_e b_e)^1/2. The following faint galaxies in the spectroscopic sample were not refit in the new processing of the HST imaging: IDs 129, 197, 780, 1133, and 2911. Of these IDs 129 and 197 are faint disk galaxies in Lynx W, while none of the other galaxies are members of the Lynx clusters. In Table 7, we list the photometry from Jørgensen et al. (2014) for these galaxies. All magnitudes in Table 7 have been corrected for galactic extinction.

Table 7.  Lynx Photometric Parameters from HST/ACS data

ID	IDJ2014	R.A. (J2000)	Decl. (J2000)	mtot,SEx	(i775 − z850)	${m}_{\mathrm{tot},\mathrm{dev}}$	$\mathrm{log}{r}_{{\rm{e}},\mathrm{dev}}$	${m}_{\mathrm{tot},\mathrm{ser}}$	$\mathrm{log}{r}_{{\rm{e}},\mathrm{ser}}$	nser	PA
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)	(11)	(12)	(13)
104	1698	132.15605	44.89076	24.45	0.552	24.36	−0.906	24.36	−0.906	1.5	79.9	0.58
129	1533	132.17006	44.91985	24.64	0.679	23.98	−0.401	24.51	−0.753	0.6	−88.4	0.63
145	1644	132.15846	44.90064	22.36	0.704	21.30	0.536	22.24	−0.130	0.4	−36.8	0.34
148	1748	132.15444	44.89276	23.07	0.933	22.76	−0.556	22.87	−0.643	3.0	−86.4	0.39
197	1809	132.15401	44.89896	24.35	0.931	23.96	−0.256	24.33	−0.562	1.0	−2.9	0.50
209	1763	132.14989	44.89334	21.68	0.971	21.12	0.200	21.36	0.021	2.9	−13.9	0.26
254	2015	132.14189	44.88400	24.42	0.886	24.20	−1.121	24.14	−1.086	4.8	−52.5	0.13
293	1888	132.15065	44.90477	22.14	0.919	21.98	−0.333	21.63	−0.063	5.8	44.2	0.24
316	2063	132.14194	44.89228	23.18	0.781	22.77	−0.393	23.04	−0.593	2.1	−70.4	0.28
347	2138	132.13880	44.89084	24.12	0.494	23.56	−0.284	24.12	−0.670	0.8	−24.0	0.52

Note. Column 1: galaxy ID. Column 2: ID from Jørgensen et al. (2014). Columns 3 and 4: positions consistent with USNO (Monet et al. 1998), with an rms scatter of ≈05. Column 5: total F850LP magnitude from SExtractor. Column 6: aperture color within an aperture with radius 05. Column 7: total magnitude from fit with r^1/4 profile. Column 8: logarithm of the effective radius in arcsec from fit with r^1/4 profile. Columns 9 and 10: total magnitude and logarithm of the effective radius, from fit with a Sérsic profile. Column 11: Sérsic index. Column 12: position angle of major axis measured from north through east; Column 13: ellipticity. This table is available in its entirety as a machine readable table. A portion is shown here for guidance on its content.

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

In Figure 17 and Table 8, we compare the new determinations for the Lynx W galaxies with those used in Jørgensen et al. (2014). There are no significant offsets between the two sets of parameters.

Zoom In Zoom Out Reset image size

Table 8.  Comparison with Jørgensen et al. (2014)

Parameter	Δ	rms
mtot,ser	0.001	0.10
mtot,dev	0.002	0.11
$\mathrm{log}{r}_{{\rm{e}},\mathrm{ser}}$	−0.0003	0.078
$\mathrm{log}{r}_{{\rm{e}},\mathrm{dev}}$	−0.0001	0.053
PA	−0.004	3.2
	0.0001	0.018
nser	0.005	0.54

Note. All comparisons are for the filter F850LP and include 47 galaxies. Δ list the median differences "this paper"–"Jørgensen et al. (2014)."

Download table as:  ASCIITypeset image

As in Jørgensen et al. (2014), we derive the rest-frame B-band from the observed z₈₅₀ magnitudes and colors, using calibrations established based on Bruzual & Charlot (2003) stellar population models spanning the observed color range. For the calibration, we adopt the U, B, and V filter functions from Maíz Apellániz (2006), rather than the older filter functions distributed with the Bruzual & Charlot models. For the Lynx cluster redshift of z = 1.27 the calibration using these new filter functions is

$$\begin{eqnarray}&&{B}_{\mathrm{rest}}={z}_{850}+1.416-1.038({i}_{775}-{z}_{850}).\end{eqnarray} \tag{ 4 }$$

The large color term is due to the fact that at the cluster redshift the F850LP filter spans the 4000 Å break. For Lynx member galaxies on the red sequence, this calibration will give B_rest values 0.1 mag fainter than the calibration used in Jørgensen et al. (2014). As noted in Jørgensen et al. (2018a), the effect of the new calibration for the z < 1 GCP clusters is typically 0.05 mag, with the new calibration leading to fainter rest-frame B-band magnitudes. All GCP data have been recalibrated using the Maíz Apellániz (2006) filter functions. The distance modulus for our adopted cosmology is DM(z) = 44.74 at z = 1.27. The absolute B-band magnitude, M_B, is then derived as

$$\begin{eqnarray}&&{M}_{{\rm{B}}}={B}_{\mathrm{rest}}-{DM}(z).\end{eqnarray} \tag{ 5 }$$

Techniques for how to calibrate to a "fixed-frame" photometric system are described in detail by Blanton et al. (2003).

Here we summarize processing of the spectroscopic data, and provide additional information useful for evaluating the robustness of our results. The reader is also referred to the information in Jørgensen et al. (2014), as we do not repeat the tests described in that paper.

The spectroscopic observations were obtained in MOS mode with GMOS-N during the period of UT 2014 November 27 to 2015 January 12 under the Gemini program IDs GN-2014B-Q-22 and GN-2014B-DD-4. All observations were carried out in queue mode. Table 9 summarizes the instrumentation, exposure times, and other information relevant for the observations. The observations were processed using the same methods as those used for Lynx W, and described in detail in Jørgensen et al. (2014), see also Jørgensen & Chiboucas (2013). The processing includes standard steps for bias subtraction, flat-fielding, and sky subtraction for nod-and-shuffle observations. As for Lynx W, we use a spatially dependent correction for the charge diffusion effect. All frames were then combined and one-dimensional spectra extracted, and calibrated to a relative flux scale.

Table 9.  Spectroscopic Observations

Instrument	GMOS-N
CCDs	3 × E2V DD 2048 × 4608
r.o.n.a	(3.17, 3.22, 3.46) e−
gaina	(2.31, 2.27, 2.17) e−/ADU
Pixel scale	00727 pixel−1
Field of view	55 × 55
Grating	R400_G5305
Spectroscopic filter	OG515_G0306
Wavelength rangeb	5500–10500 Å
Slit width x slit length	1'' × 275
Total exposure time	55,800 s
Number of frames	31
Image qualityc	057
Instrumental resolutiond	3.002 Å, 100 km s−1
Aperturee	1'' × 07, 053
S/Nf	23.3

Notes.

^aValues for the three detectors in the array. ^bThe exact wavelength range varies from slitlet to slitlet. ^cImage quality measured as the average FWHM at 8000 Å of the blue stars included in the masks. ^dMedian instrumental resolution derived as sigma in Gaussian fits to the sky lines of the stacked spectra. The second entry is the equivalent resolution in km s⁻¹ at 4000 Å in the rest frame of the cluster. ^eThe first entry is the rectangular extraction aperture (slit width × extraction length). The second entry is the radius in an equivalent circular aperture, r_ap = 1.025 (length × width/π)^1/2, see Jørgensen et al. (1995). ^fMedian S/N per Ångstrom for the 16 cluster members, in the rest frame of the cluster, measured in the wavelength interval 3750–4100 Å.

Download table as:  ASCIITypeset image

The calibrated spectra were fit with stellar templates as described in Jørgensen et al. (2014). This results in determination of the redshifts and the velocity dispersions. The fits were performed with the kinematics fitting software made available by Karl Gebhardt, see Gebhardt et al. (2000, 2003) for a description of the software. Spectra of member galaxies were fit in the wavelength range of 3750–4125 Å. Because the software determines the fits in pixel space, it is straightforward to mask wavelength ranges not to be included in the fit. We use this to flag areas of strong residuals from the sky subtraction as well as emission lines. The instrumental resolution is determined from stacked sky spectra processed in the same way as the science spectra. As described in Jørgensen & Chiboucas (2013), we use three template stars with spectral types K0III, G1V, and B8V. Figure 18 shows the normalized spectra, fits, and residuals in the wavelength region covered by the fits. The purpose of the set of template stars is to span the stellar populations in the galaxies. From Figure 18 we conclude that the use of these three stars accomplishes this. Further, we choose to use the same template stars as were used in our previous work to ensure consistency with our previous publications. The software determines the line-of-sight-velocity-distribution (LOSVD) from the science spectra and then the velocity dispersion is derived from the LOSVD through both a Gauss–Hermite polynomial fit and a Gaussian fit. For the 15 member galaxies with velocity dispersion determined, the offset in log space between the two measurements is −0.034 with the Gauss–Hermite polynomial fits giving the smaller velocity dispersions. For consistency with our previous data, we use the velocity dispersions derived from the Gauss–Hermite polynomial fits to the LOSVD, rather than the Gaussian fits. Aperture correction of the velocity dispersions were performed using the technique from Jørgensen et al. (1995).

Zoom In Zoom Out Reset image size

The uncertainties, ${\sigma }_{\mathrm{log}\sigma }$, on the velocity dispersions are derived by the fitting software using a bootstrap method. In Figure 19 we show the resulting uncertainties versus the S/N of the spectra. The uncertainties scale approximately with S/N as ${\sigma }_{\mathrm{log}\sigma }\approx 1.3\times {\rm{S}}/{{\rm{N}}}^{-1}$. This relation is shown in the figure. The relation is not used to determine uncertainties in our measurements, but may be useful for other researchers planning similar observations.

Zoom In Zoom Out Reset image size

The spectra cover wavelengths to 4400 Å in the rest frame. However, the sky subtraction errors limit the range for which the spectra can be used for determination of absorption line indices to ≈4100 Å in the rest frame (9250 Å observed). Thus, we can derive the following indices CN3883, CaHK, D4000, and Hζ_A. The indices CN3883 and CaHK are defined in Davidge & Clark (1994). For the high-order Balmer line index ${\rm{H}}{\zeta }_{{\rm{A}}}$ we adopt the definition from Nantais et al. (2013). For D4000 we use a shorter red passband than that used in the standard definition. The resulting D4000 indices are calibrated to consistency with the usual definition from Gorgas et al. 1999), see Jørgensen et al. (2014).

The line indices have been aperture corrected and corrected for velocity dispersion, as described in Jørgensen et al. (2005), see Jørgensen & Chiboucas (2013) for a discussion of the method applied to intermediate redshift galaxies. As for the lower order Balmer lines, we assume that Hζ_A has no aperture correction.

For galaxies with detectable emission from [O ii] we determined the equivalent width of the [O ii]λλ3726, 3729 doublet. With an instrumental resolution of σ ≈ 3 Å (FWHM ≈ 7 Å), the doublet is not resolved in our spectra and we refer to it simply as the "[O ii] line."

Table 10 lists the results from the template fitting and the derived line strengths.

Table 10.  Spectroscopic Parameters

ID	Redshift	Membera	$\mathrm{log}\sigma$	$\mathrm{log}{\sigma }_{\mathrm{cor}}$ b	${\sigma }_{\mathrm{log}\sigma }$	Template Fractions			χ2	S/Nc	CN3883	σCN3883	CaHK	σCaHK	D4000	σD4000	HζA	${\sigma }_{{\rm{H}}{\zeta }_{A}}$	EW [O ii]	${\sigma }_{\mathrm{EW}[\mathrm{OII}]}$
						B8V	G1V	K0III
2076	1.1409	0	⋯	⋯	⋯	⋯	⋯	⋯	⋯	5.8	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	56.3	28.3
2098	1.2714	1	2.304	2.328	0.063	0.44	0.42	0.14	3.4	14.0	0.126	0.015	24.58	0.91	1.863	0.011	2.52	0.40	20.2	3.2
2309	1.2682	1	2.576	2.600	0.073	0.00	0.51	0.49	3.4	21.7	0.297	0.012	19.89	0.61	2.098	0.009	2.10	0.32	⋯	⋯
2333	1.0138	0	⋯	⋯	⋯	⋯	⋯	⋯	⋯	24.1	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	3.2	0.3
2478	1.0731	0	2.186	2.210	0.077	1.00	0.00	0.00	2.0	27.6	⋯	⋯	6.53	0.45	1.388	0.004	2.87	0.17	39.7	2.5
2626	0.6777	0	2.410	2.432	0.226	0.74	0.26	0.00	1.6	10.8	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯
2688	1.2236	0	2.225	2.249	0.138	0.82	0.18	0.00	1.9	9.9	0.085	0.019	20.75	1.34	1.292	0.010	6.44	0.34	45.3	17.9
2757	1.2740	1	2.057	2.081	0.052	0.08	0.84	0.08	2.3	25.2	0.116	0.009	11.12	0.59	1.545	0.006	1.54	0.28	⋯	⋯
3047	1.0370	0	1.972	1.995	0.077	0.40	0.60	0.00	4.7	26.1	⋯	⋯	17.53	0.41	⋯	⋯	2.95	0.20	⋯	⋯
4093	1.2635	1	2.294	2.318	0.129	0.66	0.00	0.34	3.8	13.4	0.139	0.017	16.89	1.12	1.889	0.013	6.31	0.35	18.1	1.9
4095	1.2631	1	2.337	2.361	0.028	0.27	0.54	0.20	10.7	77.8	0.122	0.003	18.43	0.20	1.886	0.002	3.14	0.08	⋯	⋯
4123	1.1371	0	2.203	2.227	0.052	0.81	0.19	0.00	2.3	17.2	0.080	0.012	17.01	0.83	1.511	0.008	4.61	0.28	42.9	9.4
4128	1.1999	0	2.139	2.163	0.095	0.00	0.80	0.20	3.6	16.1	0.061	0.013	21.45	0.76	1.665	0.008	2.76	0.28	5.2	0.4
4216	1.2293	0	2.408	2.432	0.086	1.00	0.00	0.00	2.1	16.8	⋯	⋯	6.29	1.00	1.230	0.006	3.47	0.28	35.9	3.6
4298	1.1991	0	2.227	2.251	0.109	0.50	0.27	0.23	2.3	11.1	0.381	0.022	5.25	1.28	1.986	0.014	3.38	0.46	10.9	0.5
4363	1.2662	1	2.237	2.261	0.106	0.46	0.22	0.32	3.5	25.0	0.140	0.010	13.47	0.56	1.874	0.007	3.39	0.25	2.9	0.2
4508	1.2644	1	2.064	2.088	0.102	0.37	0.62	0.01	3.3	26.2	0.063	0.009	14.12	0.58	1.796	0.007	3.66	0.20	28.4	1.7
4593	1.2692	1	2.391	2.415	0.050	0.62	0.38	0.00	2.7	29.1	0.079	0.008	16.69	0.47	1.864	0.006	6.21	0.16	4.5	0.5
4628	1.1718	0	⋯	⋯	⋯	⋯	⋯	⋯	⋯	6.6	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	35.9	1.7
4638	1.5870	0	⋯	⋯	⋯	⋯	⋯	⋯	⋯	2.6	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	162.1	77.8
4722	⋯	3	⋯	⋯	⋯	⋯	⋯	⋯	⋯	1.6	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯
4828	1.2686	1	2.244	2.268	0.089	0.28	0.72	0.00	1.4	25.9	0.120	0.009	19.34	0.55	1.847	0.007	3.24	0.23	10.5	0.7
4860	1.2587	1	2.166	2.190	0.088	0.34	0.66	0.00	1.6	20.4	0.102	0.012	15.63	0.81	2.003	0.009	2.50	0.29	⋯	⋯
4906	⋯	3	⋯	⋯	⋯	⋯	⋯	⋯	⋯	2.2	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯	⋯
4921	1.2590	1	2.011	2.035	0.144	0.39	0.40	0.20	2.5	13.5	⋯	⋯	23.54	1.13	1.487	0.012	4.65	0.40	⋯	⋯
4926	1.2681	1	2.053	2.077	0.110	0.09	0.68	0.22	1.7	15.2	0.224	0.016	18.89	0.90	1.770	0.011	3.06	0.42	⋯	⋯
4942	1.2594	1	2.445	2.469	0.038	0.12	0.51	0.37	5.0	45.0	0.236	0.006	22.03	0.32	1.992	0.004	2.26	0.15	⋯	⋯
5009	1.2622	1	2.206	2.230	0.058	0.22	0.44	0.34	1.8	14.6	0.323	0.017	22.13	0.94	2.079	0.013	3.09	0.44	⋯	⋯
5037	1.2607	1	2.120	2.144	0.103	0.30	0.29	0.41	4.6	24.8	0.240	0.010	18.85	0.56	2.161	0.008	4.08	0.27	⋯	⋯

Notes. The typical uncertainties on the redshifts are 0.0001. Uncertainties on the template fractions are not derived by default by the fitting software, but our independent Monte Carlo simulations indicate uncertainties of 0.02–0.05. The absorption line indices have been corrected for galaxy velocity dispersion and aperture corrected. D4000 is measured as ${\rm{D}}{4000}_{\mathrm{short}}$, see Jørgensen et al. (2014). This table is also available as a machine readable table.

^aAdopted membership: 1—galaxy is a member of Lynx E or W; 0—galaxy is not a member of either Lynx E or W; 3—redshift cannot be determined. ^bVelocity dispersions corrected to a standard size aperture equivalent to a circular aperture with a diameter of 34 at the distance of the Coma cluster. ^cS/N per Ångstrom in the rest frame of the galaxy. The wavelength interval was chosen based on the redshift of the galaxy as follows: redshift 1.00–1.35: 3750–4100 Å; and redshift <1.00: 4100–4600 Å. For ID 4722 and 4906 a redshift of 1.27 was assumed for the S/N calculation.

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

The new spectra as well as stamp-sized images of the galaxies from the HST/ACS imaging of the cluster members are shown in Figure 20. The stamps cover the equivalent of 75 kpc × 75 kpc at the distance of the cluster. The spectra used to create Figure 20 are available in the online journal.

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

D.1. Calibration of SDSS Photometry

In Jørgensen et al. (2018b) we presented consistently calibrated spectroscopic parameters for a large low-redshift reference sample of bulge-dominated galaxies in both the Perseus and Coma cluster. The data were calibrated to consistency with our Legacy Data for the two clusters used in our previous GCP papers (Barr et al. 2005; Jørgensen et al. 2005, 2014, 2017; Jørgensen & Chiboucas 2013). In order to take advantage of this reference sample, we establish the photometric parameters for the sample, and ensure that those are also consistent with the Legacy Data. For the purpose of the use in the present paper we rely on the SDSS g'-band photometry of the galaxies. The full detail and the calibrated data will be included in a future paper (I. Jørgensen et al. 2019, in preparation). Here we summarize the main points and compare the calibrated data to the photometric data from Simard et al. (2011), see Section D.2.

The SDSS data release 14 (DR14) contains photometric parameters for the galaxies from r^1/4 profile fitting (devmag, devrad), exponential (expmag, exprad) profile fitting, and a total magnitude cmodelmag from a best-fit linear combination of these two fits with information about the fraction fracdev of the flux originating from the r^1/4 fit in the combination. We use the SDSS photometric parameters to derive parameters that mimic Sérsic (1968) parameters. These parameters are in the following referred to as pseudo-Sérsic parameters to not confuse them with parameters from actual fits to Sérsic profiles. In the main body of the paper, we simply refer to the parameters as Sérsic parameters. We proceed as follows: we adopt cmodelmag as the best-fit total magnitude, ${m}_{\mathrm{tot},\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$. We then derive the matching half-light radius, ${r}_{{\rm{e}},\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$, from the information about the linear combination of the r^1/4 profile and the exponential profile.

Using the general equation the enclosed luminosity of a Sérsic profile within a radius (Graham & Driver 2005) and the information about the linear combination of the r^1/4 and exponential profiles, we derive the following:

$$\begin{eqnarray}&&\begin{array}{l}{f}_{\mathrm{dev}}\left(\displaystyle \frac{1}{{\rm{\Gamma }}(8)}\gamma \left(8,{b}_{4}{\left(\displaystyle \frac{R}{{r}_{\mathrm{dev}}}\right)}^{0.25}\right)-0.5\right)\\ \quad +\,(1-{f}_{\mathrm{dev}})\displaystyle \frac{{L}_{\exp }}{{L}_{\mathrm{dev}}}\left(\displaystyle \frac{1}{{\rm{\Gamma }}(2)}\gamma \left(2,{b}_{1}\left(\displaystyle \frac{R}{{r}_{\exp }}\right)\right)-0.5\right)=0\end{array}\end{eqnarray} \tag{ 6 }$$

where f_dev is fracdev from SDSS, and the ratio of luminosities ${L}_{\exp }\,{L}_{\mathrm{dev}}^{-1}$ is calculated from the corresponding magnitudes devmag and expmag. The radii r_dev and r_exp are devrad and exprad, respectively, both circularized using the axis ratios from SDSS. Γ is the gamma function, γ is the incomplete gamma function, and b₁ and b₄ are specific values of the parameter b_n for which ${\rm{\Gamma }}(2n)=2\gamma (2n,{b}_{n})$ and have values of b₁ = 1.6783 and b₄ = 7.6692. Solving Equation (6) for R gives ${r}_{\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$ matching ${m}_{\mathrm{tot},\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$. From ${r}_{\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$ and ${m}_{\mathrm{tot},\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$ we can derive the mean surface brightness within ${r}_{\mathrm{pseudo}\mbox{--}\mathrm{Sersic}}$.

We correct the SDSS photometry for galactic extinction using the extinction info included in DR14. The correction is based on the dust calibration from Schlafly & Finkbeiner (2011). Finally we calibrate the SDSS r^1/4 and pseudo-Sérsic parameters to consistency with the Legacy Data by comparing the SDSS r^1/4 parameters with to the Legacy Data B-band data for the Coma cluster. The Legacy Data originally used the Burstein & Heiles (1982) extinction correction. Here we instead adopt the calibration from Schlafly & Finkbeiner (2011), which gives 0.02 mag fainter extinction corrected magnitudes.

Figure 21 shows the comparison of the Legacy Data with the SDSS r^1/4 parameters. In panels (a) and (b) the difference in total magnitudes, ${B}_{\mathrm{Legacy}}-{g}_{\mathrm{dev},\mathrm{SDSS}}^{{\prime} }$, is shown versus color and the product of the three filters' fracdev, f_g × f_r × f_i. The magnitude offset does not depend on the color. About half of the tail of points scattering to low values of ${B}_{\mathrm{Legacy}}-{g}_{\mathrm{dev},\mathrm{SDSS}}^{{\prime} }$ are from galaxies with f_g × f_r × f_i < 0.4. We take this as in indication that the growth curve method used for the Legacy Data (Jørgensen et al. 1995) and the fitting implemented in the SDSS pipeline handle galaxies with faint disks differently. The median magnitude difference is 0.398 ± 0.019 mag. In Figure 21(c), we show the difference in radii versus the difference in surface brightnesses. The two photometric bands are close enough in wavelength space that we do not expect any significant difference in radii due to the color gradients of the galaxies. As expected the differences are tightly correlated. The best fit is shown as the solid line. The intercept gives the offset between the two magnitudes, 0.444 ± 0.005 mag. For both determinations of the magnitude difference, points deviating more than 3σ have been excluded. To calibrate the SDSS $g^{\prime}$-band photometry to consistency with the Legacy Data B-band photometry, we adopt the average of those two offsets, 0.421 mag. Both r^1/4 SDSS and pseudo-Sérsic magnitudes are calibrated to the B-band using this offset. Due to the correlation between the errors on the surface brightnesses (or total magnitudes) and the effective radii, the intercept in Figure 21(c) is less affected by the tail of points at lower ${B}_{\mathrm{Legacy}}-g{{\prime} }_{\mathrm{dev},\mathrm{SDSS}}$, than the median magnitude difference itself. Based on the difference between the adopted offset and the intercept value, we estimate that this tail in the distribution affects the calibration by no more than 0.02 mag, which is insignificant for our purpose.

Zoom In Zoom Out Reset image size

The B-band photometry is calibrated to rest-frame B-band using the k-correction info from Chilingarian et al. (2010). On average the correction for Coma cluster galaxies is 0.075 mag. In our use of the Legacy Data we adopted the k-correction from Pence (1976), which is 0.111 mag for the Coma cluster. Thus, together with the difference in the adopted galactic extinction the recalibrated new reference sample is 0.056 mag fainter than the Legacy Data used in previous GCP papers.

D.2. Comparison with Literature Data

To ensure that our approach to converting the SDSS DR14 parameters to pseudo-Sérsic parameters is valid, we compare the Coma cluster data with the Sérsic parameters from the two-dimensional fitting of the SDSS imaging performed by Simard et al. (2011). The Perseus cluster is not included in the Simard et al. catalog. Figure 22 summarizes the comparisons. The agreement is very good with the total magnitudes as well as the effective radii being fully consistent with the two-dimensional fitting parameters from Simard et al. (2011). As shown in Figure 22(c) we take the opportunity to evaluate if our technique in Jørgensen et al. (2018b) for selecting the bulge-dominated galaxies using the product of the SDSS fracdev parameters, f_g × f_r × f_i, is consistent with using the Sérsic index, n_ser. The vertical and horizontal dashed lines on the figure show the limits in the two parameters used in our previous sample selections, while the red points highlight the final selection of bulge-dominated Coma cluster galaxies (Jørgensen et al. 2018b). In general the two methods are equivalent. The few galaxies marked in blue at high ${f}_{g}\times {f}_{r}\times {f}_{i}$ and high n_ser are galaxies that upon visual inspection were identified as spiral galaxies (see Jørgensen et al. 2018b).

Zoom In Zoom Out Reset image size

We here test to what extent our results are sensitive to our choice of deriving the dynamical masses as ${{\rm{m}}\mathrm{ass}}_{\mathrm{dyn}}=\beta {r}_{e}{\sigma }^{2}\,{G}^{-1}$ with β = 5, rather than using β values dependent on the Sérsic (1968) index as proposed by Cappellari et al. (2006)

$$\begin{eqnarray}&&\beta ({n}_{{\rm{s}}{\rm{er}}})=8.87-0.831\,{n}_{{\rm{s}}{\rm{er}}}+0.0241\,{n}_{{\rm{s}}{\rm{er}}}^{2}.\end{eqnarray} \tag{ 7 }$$

For this test we use the Sérsic indices from Simard et al. (2011) for the Coma cluster sample. This limits the low-redshift sample to the 80 passive bulge-dominated galaxies with ${\mathrm{log}{\rm{m}}\mathrm{ass}}_{\mathrm{dyn}}\geqslant 10.3$ and available Sérsic index. We first fit relevant relations to this smaller low-redshift sample, using β = 5. Then we repeat the fits using masses derived using β from Equation (7). For the latter fits we also derive zero-points for the Lynx E and W. The fits and zero-points are summarized in Table 11, while Figures 23 and 24 show the relations for masses using β (n_ser).

Zoom In Zoom Out Reset image size

Zoom In Zoom Out Reset image size

Table 11.  Structure Relations

Relation	Low Redshift			Lynx E			Lynx W
	γ	Ngal	rms	γ	Ngal	rms	γ	Ngal	rms
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$\mathrm{log}{r}_{e}=(0.72\pm 0.18)\mathrm{log}\,\mathrm{Mass}+\gamma$ a	−7.416	80	0.18	⋯	⋯	⋯	⋯	⋯	⋯
$\mathrm{log}{r}_{e}=(0.65\pm 0.17)\mathrm{log}\,\mathrm{Mass}+\gamma$ b	−6.712	80	0.16	−6.819	13	0.22	−6.745	12	0.20
$\mathrm{log}\sigma =(0.33\pm 0.05)\mathrm{log}\,\mathrm{Mass}+\gamma$ a	−1.343	80	0.08	⋯	⋯	⋯	⋯	⋯	⋯
$\mathrm{log}\sigma =(0.30\pm 0.04)\mathrm{log}\,\mathrm{Mass}+\gamma$ b	−1.071	80	0.09	−1.024	13	0.09	−1.087	12	0.09
$\mathrm{log}\,M/L=(0.32\pm 0.06)\mathrm{log}\,\mathrm{Mass}+\gamma$ a	−2.612	80	0.08	⋯	⋯	⋯	⋯	⋯	⋯
$\mathrm{log}\,M/L=(0.35\pm 0.06)\mathrm{log}\,\mathrm{Mass}+\gamma$ b	−2.888	80	0.07	−3.696	13	0.23	−3.720	12	0.25
$\mathrm{log}\,M/L=(0.91\pm 0.40)\mathrm{log}\sigma +\gamma$ a	−1.193	80	0.10	⋯	⋯	⋯	⋯	⋯	⋯
$\mathrm{log}\,M/L=(0.88\pm 0.32)\mathrm{log}\sigma +\gamma$ b	−1.045	80	0.11	−1.820	13	0.26	−1.820	12	0.28

Notes. Column 1: scaling relation, slopes from best-fit Coma sample with available Sérsic indices. Column 2: zero-point for the low-redshift sample (Coma). Column 3: number of galaxies included from the low-redshift sample. Column 4: rms in the Y-direction of the scaling relation for the low-redshift sample. Columns 5, 6, and 7: zero-point, number of galaxies, rms in the Y-direction for the Lynx W sample. Columns 8, 9, and 10: zero-point, number of galaxies, rms in the Y-direction for the Lynx E sample.

^aMasses derived using β = 5. ^bMasses derived using β (n_ser) from Equation (7).

Download table as:  ASCIITypeset image

A comparison to the fits to the low-redshift sample using fixed β with the fits using β (n_ser) shows that all slope differences are well within the 1σ uncertainties, see Table 11. Comparing the slopes from the fits using β (n_ser) (Table 11) to the relations derived using the full Perseus and Coma samples and β = 5 (Table 3 in the main text) shows slightly larger differences, but still well within the 1σ uncertainties on the slopes.

We then compare the zero-point offsets for Lynx E and W as derived in the main text, using β = 5, see Table 3, with the zero-point offsets using β (n_ser), see Table 11. In both cases, there are no significant differences between Lynx E and W. The median zero-point differences for the full Lynx sample using $\beta ({n}_{\mathrm{ser}})$ are ${\rm{\Delta }}\mathrm{log}{r}_{e}=-0.07\pm 0.05$ and ${\rm{\Delta }}\mathrm{log}\sigma =0.016\,\pm 0.021$, to be compared with ${\rm{\Delta }}\mathrm{log}{r}_{e}=-0.02\pm 0.04$ and ${\rm{\Delta }}\mathrm{log}\sigma =0.003\pm 0.019$ for β = 5, see Section 7.1. Similarly for the median offset in the M/L ratio at a fixed mass we find ${\rm{\Delta }}\mathrm{log}\,M/L=-0.82\pm 0.05$, to be compared with ${\rm{\Delta }}\mathrm{log}\,M/L=-0.75\pm 0.05$ for β = 5. None of the results for the two methods are different at more than the 1σ level.

In summary, the results in this section confirm that our results as presented in the main text do not significantly depend on our choice to derive the dynamical masses as ${\mathrm{Mass}}_{\mathrm{dyn}}=\beta {r}_{e}{\sigma }^{2}\,{G}^{-1}$ with β = 5.