COLIBRE Simulations of a z=4 Quiescent Proto-cluster

• COLIBRE Simulations are advanced computational models that replicate a proto-cluster of massive, passively evolving galaxies at z=4 with rapid starburst and quenching phases.
• The simulations integrate multi-band imaging and precise spectroscopic methods to determine stellar masses, velocity dispersions, and structural parameters of early quiescent galaxies.
• Detailed SED fitting and dynamical measurements from the simulations challenge current cosmological models, indicating a need to revise theories on early galaxy assembly and quenching.

A proto-cluster of massive quiescent galaxies at z=4 represents a pivotal observational discovery in galaxy evolution, revealing the existence of a concentration of compact, passively evolving galaxies in the early Universe. Identified in the Subaru/XMM-Newton Deep Field, this system comprises at least five massive ($M_*\gtrsim10^{10}\,M_\odot$) quiescent galaxies co-located on $\sim$1 Mpc physical scales at $z\approx4$, with spectroscopic confirmation via Keck/MOSFIRE for the brightest ($M_*\sim1.1\times10^{11}\,M_\odot$) member. The population’s properties—stellar masses, velocity dispersions, sizes, derived dynamical and halo masses, over-density metrics, and nascent red sequence—directly constrain models for early starburst and quenching, hierarchical structure formation, and the limits of current cosmological simulations (Tanaka et al., 2023).

1. Survey Data, Photometric Selection, and Spectroscopic Confirmation

The observational strategy employs multi-band, deep near-IR and optical imaging (Subaru HSC, UKIRT, VISTA, Spitzer) spanning 1.3 deg² in the SXDS. Photometric redshifts are obtained using the Tanaka et al. (2015) template fitting code, based on Bruzual & Charlot (2003) models and the Chabrier IMF. This yields a precision $\sigma_{\Delta z/(1+z)}\approx0.036$ with $<5\%$ outliers against VANDELS spectroscopy. Quiescent candidates are selected with sSFR$<10^{-9.5}\,\rm yr^{-1}$.

A search for overdensities reveals a $20\sigma$ concentration of quiescent galaxies within $<2$ physical Mpc for $3.7<z_{\rm phot}<4.3$. Five of these undergo Keck/MOSFIRE K-band spectroscopy (16 hr total, $0.7''$ slit, ABBA dither), with reduction including custom flat fielding, A–B sky subtraction, wavelength/telluric calibration, and optimal extraction. The primary confirmed galaxy has $z_{\rm spec}=3.99$, and companions $z_{\rm phot/spectro}\sim4$.

2. SED Modeling, Star Formation Histories, and Quenching

Stellar mass estimation combines joint photometric and binned spectroscopic SED fits. The method integrates the star formation rate over cosmic time, incorporating mass return $R\approx0.3$ for a Chabrier IMF: $M_* = (1-R)\int_0^{t_{\rm obs}}\mathrm{SFR}(t)\,dt$ The recent SFH is parameterized as a two-phase model representing a brief starburst followed by rapid exponential quenching: $$\mathrm{SFR}(t) = \begin{cases} \mathrm{SFR}_{\rm burst} & t < t_q \ \mathrm{SFR}_{\rm burst}\exp[-(t-t_q)/\tau_q] & t \geq t_q \end{cases}$$ For the $z=3.99$ galaxy, the SED fit yields $M_* \simeq 1.1 \times 10^{11}\,M_\odot$, a starburst $\sim$500 Myr prior to observation, and a quenching timescale $\tau_q<200$ Myr. Similar histories are found for companions at $z=4.01$.

3. Velocity Dispersion, Structural Parameters, and Dynamical Masses

Absorption-line fitting (pPXF, Vazdekis et al. 2010 SSPs) measures a stellar velocity dispersion $\sigma_*=305\pm103$ km/s for the $z=3.99$ system. GALFIT on ground-based HSC imaging provides an effective radius $R_e=1.04\pm0.43$ kpc. The dynamical mass formula, using $k(n)\simeq5$ for Sersic index $n\sim$2–4,

$M_{\rm dyn}=k(n)R_e\sigma_*^2/G$

yields $M_{\rm dyn}\approx M_*$, supporting a stellar-mass-dominated potential. Such compactness and high dispersion are characteristic of quiescent "red nuggets" at high $z$.

4. Halo Masses and Large-Scale Environment

Halo occupation models and abundance matching (Shuntov et al. 2022) relate $M_*$ to $M_{\rm halo}$: $\log M_{\rm halo} \approx \alpha\log M_* + \beta$ For these galaxies, $M_{\rm halo}\sim\mathrm{few}\times10^{12}\,M_\odot$ individually, and a collapsed total $M_{\rm halo}\sim10^{13}\,M_\odot$ for the proto-cluster, corresponding to a group-scale halo.

The galaxies are embedded in a larger-scale N–S filament traced by spectroscopically confirmed VANDELS galaxies ($3.95

5. Overdensity Signatures, Red Sequence, and Proto-Cluster Classification

Within a projected radius $\approx1$ Mpc, the five quiescent galaxies achieve the highest known $z\ge4$ surface density ($\sim1.5$ Mpc$^{-2}$). H–K vs K color–magnitude diagrams show a tight color locus ($\Delta(\mathrm{H-K})<0.1$ mag), establishing a red sequence at $z=4$. The localized quiescent fraction approaches $35\%$, far exceeding the $1.5\%$ baseline field value—a strong Butcher–Oemler effect extension.

No analogs of such a dense quiescent proto-cluster are found in the 300 cMpc volume of Illustris-TNG300 at $z=4$—it contains only 11 galaxies with $M_*>10^{10}\,M_\odot$ and $\log\,\mathrm{sSFR}<-9.5$, none within 1 Mpc of another.

6. Implications for Early Cluster Formation and Quenching Physics

This system demonstrates that by $z=4$—only 1.5 Gyr after the Big Bang—massive, compact, quiescent galaxies (with evidence of a red sequence) can form and concentrate in group-scale halos. The presence of a rapid starburst and subsequent quenching ($\tau_q<200$ Myr), supported by the measured velocity dispersions and structural parameters, imposes strong constraints on feedback, burst-triggering mechanisms (possibly mergers or AGN), and early environmental processes.

The rarity of such systems in current simulations suggests that (i) present simulation volumes may be insufficient, or (ii) early rapid quenching in dense peaks is under-predicted, requiring recalibration of high-redshift feedback and star-formation physics. The high local quiescent fraction, compactness, and mass make the system a critical target for testing hierarchical assembly, AGN feedback, and the emergence of red sequences in the early Universe.

7. Future Directions and Survey Implications

The identification of a $z=4$ proto-cluster dominated by quiescent galaxies offers immediate avenues for spectroscopic and environmental follow-up to probe kinematics, chemical abundances, relic AGN, and subsequent mergers. Larger-volume cosmological simulations, coupled with deep spectroscopic surveys, are needed to quantify the statistical frequency of such concentrations and their implications for cluster formation and early galaxy quenching. This discovery sets new benchmarks for the timing, physical drivers, and environmental dependence of galaxy evolution in the first Gyr of cosmic history (Tanaka et al., 2023).