We study the environments of a sample of 61 extremely rare z~1.6 Ultra-Massive Passively Evolving Galaxies (UMPEGs: stellar masses M_stars >10^11.5 M_sun) which -- based on clustering analysis presented in Cheema et al. (2020) -- appear to be associated with very massive (M_halo ~ 10^14.1 h^-1 M_sun) dark matter halos that are likely to be the progenitors of z~0 massive (Coma- and Virgo-like) galaxy clusters. We find that UMPEGs on average have fewer than one satellite galaxy with mass ratio M_sat : M_UMPEG >~ 1:5 (i.e., M_sat >~ 10^10.8 M_sun) within 0.5 Mpc; the large mass gap that we observe between the typical UMPEG and its most massive satellite implies that the z~1.6 UMPEGs assembled through major mergers. Using observed satellite counts with merger timescales from the literature, we estimate the growth rate due to mergers with mass ratio of >~ 1:4 to be ~13% Gyr^-1 (with a ~2x systematic uncertainty). This relatively low growth rate is unlikely to significantly affect the shape of the massive end of the stellar mass function, whose evolution must instead be driven by the quenching of new cohorts of ultra-massive star-forming galaxies. However, this growth rate is high enough that, if sustained to z~0, the typical z~1.6 M_UMPEG=10^11.6 M_sun UMPEG can grow into a M_stars~10^12 M_sun brightest cluster galaxy (BCG) of a present-day massive galaxy cluster. Our observations favour a scenario in which our UMPEGs are main-branch progenitors of some of the present-day BCGs that have first assembled through major mergers at high redshifts and grown further through (likely minor) merging at later times.

We use a 27.6 deg$^2$ survey to measure the clustering of $gzK_s$-selected quiescent galaxies at $z\sim1.6$, focusing on ultra-massive quiescent galaxies. We find that $z\sim1.6$ Ultra-Massive Passively Evolving Galaxies (UMPEGs), which have $K_s(AB)<19.75$ (stellar masses of $M_{stars}$ $>\sim 10^{11.4}M_{\odot}$ and mean $<$$M_{stars}$$>$ = $10^{11.5}M_{\odot}$), cluster more strongly than any other known galaxy population at high redshift. Comparing their correlation length, $r_0 = 29.77 \pm 2.75$ $ h^{-1}$Mpc, with the clustering of dark matter halos in the Millennium XXL N-body simulation suggests that these $z\sim1.6$ UMPEGs reside in dark matter halos of mass $M_{h}\sim10^{14.1}h^{-1}M_{\odot}$. Such very massive $z\sim1.6$ halos are associated with the ancestors of $z\sim0$ massive galaxy clusters such as the Virgo and Coma clusters. Given their extreme stellar masses and lack of companions with comparable mass, we surmise that these UMPEGs could be the already-quenched central massive galaxies of their (proto)clusters. We conclude that with only a modest amount of further growth in their stellar mass, $z\sim1.6$ UMPEGs could be the progenitors of some of the massive central galaxies of present-day massive galaxy clusters observed to be already very massive and quiescent near the peak epoch of the cosmic star formation.

We introduce the largest to date survey of massive quiescent galaxies at redshift z~1.6. With these data, which cover 27.6 deg^2, we can find significant numbers of very rare objects such as ultra-massive quiescent galaxies that populate the extreme massive end of the galaxy mass function, or dense environments that are likely to become present-day massive galaxy clusters. In this paper, the first in a series, we apply our gzK adaptation of the BzK technique to select our z~1.6 galaxy catalog and then study the quiescent galaxy stellar mass function with good statistics over M_stars ~ 10^10.2 - 10^11.7 M_sun --- a factor of 30 in mass --- including 60 ultra-massive z~1.6 quiescent galaxies with M_stars > 10^11.5 M_sun. We find that the stellar mass function of quiescent galaxies at z~1.6 is well represented by the Schechter function over this large mass range. This suggests that the mass quenching mechanism observed at lower redshifts must have already been well established by this epoch, and that it is likely due to a single physical mechanism over a wide range of mass. This close adherence to the Schechter shape also suggests that neither merging nor gravitational lensing significantly affect the observed quenched population. Finally, comparing measurements of M* parameters for quiescent and star-forming populations (ours and from the literature), we find hints of an offset (M*_SF > M*_PE), that could suggest that the efficiency of the quenching process evolves with time.

We study the angular clustering of z~2 galaxies using ~40,000 star-forming (SF) and ~5,000 passively-evolving (PE) galaxies selected from ~2.5 deg$^2$ of deep ($K_{lim}$=23-24 AB) CFHT imaging. For both populations the clustering is stronger for galaxies brighter in rest-frame optical and the trend is particularly strong for PE galaxies, indicating that passive galaxies with larger stellar masses reside in more massive halos. In contrast, at rest-frame UV we find that while the clustering of SF galaxies increases with increasing luminosity, it decreases for PE galaxies; a possible explanation lies in quenching of star formation in the most massive halos. Furthermore, we find two components in the correlation functions for both SF and PE galaxies, attributable to one- and two-halo terms. The presence of one-halo terms for both PE and SF galaxies suggests that environmental effects were producing passive galaxies in virtualized environments already by z~2. Finally, we find notable clustering differences between the four widely-separated fields in our study; the popular COSMOS field is the most discrepant (as is also the case for number counts and luminosity functions), highlighting the need for very large areas and multiple sightlines in galaxy evolution statistical studies.