Massive Quiescent Galaxies Overview

• Massive quiescent galaxies are high stellar mass systems (M_* >10^10–10^11 M_☉) with extremely low specific star formation rates, serving as key benchmarks in galaxy evolution.
• They form rapidly via intense, dust-obscured starbursts and quench quickly (within 50–300 Myr), often influenced by efficient AGN feedback and environmental factors.
• Observations and simulations show that MQGs exhibit compact, spheroidal morphologies with evolving kinematics and clustering properties, challenging current models of galaxy assembly.

Massive quiescent galaxies (MQGs) are galaxies with high stellar masses (typically $M_* \gtrsim 10^{10}$--$10^{11}\,M_\odot$) whose star formation has been shut down to low or negligible rates, as quantified by specific star formation rates (sSFR) below thresholds such as $\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$ at $z>5$ or $\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$ at lower redshift. Their presence from the epoch of reionization through cosmic noon and to the present underpins key questions in galaxy evolution, including the regulation of baryonic inflow, the efficiency of feedback mechanisms, and the early assembly of the red sequence. MQGs have become central benchmarks for models seeking to reproduce the emergence of the most massive galaxies and the diversity of galaxy quenching channels across cosmic time (Chittenden et al., 28 Apr 2025, Feldmann et al., 2016, Sherman et al., 2020, Straatman et al., 2013).

1. Definitions, Demographics, and Observational Criteria

The technical definition of MQGs depends on both stellar mass and star formation activity, with selection thresholds motivated by redshift-dependent stellar mass functions and the cosmic main sequence of star-forming galaxies.

• Mass Criterion:

High-redshift studies often use $M_* \geq 10^{10}\,M_\odot$ or $10^{10.5}\,M_\odot$ at $z>3$ (Baker et al., 4 Jun 2025, McConachie et al., 28 Oct 2025). At cosmic noon ($1.5$M_*\geq 10^{11}\,M_\odot$ (Sherman et al., 2020). For low-$10^{11}\,M_\odot$0 samples, $10^{11}\,M_\odot$1 aligns with the transition above the Schechter knee (Xu et al., 2020).

• Quiescence Criterion:

Common sSFR thresholds include $10^{11}\,M_\odot$2 or $10^{11}\,M_\odot$3 at $10^{11}\,M_\odot$4 (Binh et al., 16 Dec 2025), and $10^{11}\,M_\odot$5 at $10^{11}\,M_\odot$6 (Sherman et al., 2020). Some works use main-sequence offsets, e.g., $10^{11}\,M_\odot$7 dex below the mean (Sherman et al., 2020), or color–color criteria (rest-frame UVJ or NUV–r–J cuts) (Straatman et al., 2013, Xu et al., 2020).

• Sample Sizes and Density Evolution:

Large photometric redshift samples identify hundreds to thousands of MQGs at $10^{11}\,M_\odot$8 (Xu et al., 2020, Baker et al., 4 Jun 2025), and over 700 at $10^{11}\,M_\odot$9–$\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$0 in wide-area JWST surveys (Baker et al., 4 Jun 2025). The number density declines steeply with redshift: for $\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$1, $\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$2--$\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$3 at $\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$4, dropping by nearly two orders of magnitude by $\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$5 (Baker et al., 4 Jun 2025, Straatman et al., 2013).

• Spectroscopic Confirmation:

MQGs are confirmed with high S/N rest-frame optical/NIR spectroscopy, identifying strong Balmer absorption (H$\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$6, H$\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$7, H$\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$8), a strong 4000 Å break, and lack of nebular emission ([O II], H$\mathrm{sSFR} \lesssim 1\,\mathrm{Gyr}^{-1}$9, [O III]), with SFRs constrained to $z>5$0–$z>5$1 (Glazebrook et al., 2017, Forrest et al., 2019, Carnall et al., 2023).

2. Formation Histories, Assembly Paths, and Quenching Timescales

MQGs across cosmic time exhibit a two-phase assembly and quenching history: rapid early star formation is followed by a quenching phase, the timescale and mechanism of which vary with environment, mass, and redshift.

• Star Formation Histories:

MQGs at $z>5$2 show evidence for extremely rapid assembly, often requiring SFRs $z>5$3–$z>5$4 over $z>5$5–$z>5$6 Myr, typically in dust-enshrouded starbursts (Glazebrook et al., 2017, Tanaka et al., 2023, Carnall et al., 2023). The bulk of stellar mass forms at $z>5$7–$z>5$8, within $z>5$9 Gyr of the Big Bang (Chittenden et al., 28 Apr 2025, McConachie et al., 28 Oct 2025, Straatman et al., 2013).

• Quenching Timescales:

Observed and simulated MQGs have quenching timescales (SFR drop by $\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$0) of $\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$1–$\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$2 Myr at high redshift (Chittenden et al., 28 Apr 2025, Forrest et al., 2019, Carnall et al., 2023), much shorter than the typical gas depletion times of $\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$3–$\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$4 Myr in star-forming progenitors.

• Mechanisms:
  • Smooth Mass Accretion: In reionization-era models (e.g. Thesan-1), haloes in dense cosmic web nodes assemble via rapid, smooth filamentary inflow, with negligible major mergers ($\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$5) (Chittenden et al., 28 Apr 2025).
  • Feedback: AGN-driven feedback, particularly from rapidly growing SMBHs ($\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$6–$\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$7), injects energy (via thermal/kinetic modes) and expels cold gas, quenching star formation on Myr timescales (Chittenden et al., 28 Apr 2025, Rong et al., 2017, Carnall et al., 2023).
  • Cosmological Starvation: At slightly later times and lower redshifts, when halo specific accretion rates $\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$8–$\mathrm{sSFR} < 10^{-11}\,\mathrm{yr}^{-1}$9, the inflow of fresh gas diminishes below consumption rates, naturally quenching SF without explicit feedback (Feldmann et al., 2016).
• Environmental and Merging Effects:
  • Dense Overdensities: At $M_* \geq 10^{10}\,M_\odot$0, about half of MQGs reside in protocluster or overdense environments, where mergers and deep potentials favor both rapid assembly and quenching (McConachie et al., 28 Oct 2025, Tanaka et al., 2023), with ex-situ growth (major mergers) contributing significantly to their stellar masses.
  • Minor Mergers: At $M_* \geq 10^{10}\,M_\odot$1, dry (dissipationless) minor mergers dominate mass and size growth (Zahid et al., 2019, Patel et al., 2017). These events increase effective radii and randomize stellar orbits, especially for the most massive systems (Ji et al., 2024).

3. Structural, Kinematic, and Chemical Properties

MQGs are compact, spheroidal, and dynamically distinct from star-forming galaxies of similar mass, but exhibit significant evolution in size, concentration, and angular momentum from $M_* \geq 10^{10}\,M_\odot$2 to $M_* \geq 10^{10}\,M_\odot$3.

• Sizes and Morphologies:
  • $M_* \geq 10^{10}\,M_\odot$4 MQGs typically have rest-frame optical effective radii $M_* \geq 10^{10}\,M_\odot$5–$M_* \geq 10^{10}\,M_\odot$6 kpc, $M_* \geq 10^{10}\,M_\odot$7–$M_* \geq 10^{10}\,M_\odot$8 times smaller than local quiescent galaxies of similar mass (Patel et al., 2017, Carnall et al., 2023).
  • Sérsic indices $M_* \geq 10^{10}\,M_\odot$9–$10^{10.5}\,M_\odot$0, with axis ratios $10^{10.5}\,M_\odot$1–$10^{10.5}\,M_\odot$2 indicating round, spheroidal morphologies (Patel et al., 2017, Xu et al., 2020).
  • Size growth to $10^{10.5}\,M_\odot$3 is modest for the most massive MQGs (power-law $10^{10.5}\,M_\odot$4 or $10^{10.5}\,M_\odot$5) and is driven by minor mergers (Patel et al., 2017).
• Stellar Kinematics:
  • Quiescent galaxies at $10^{10.5}\,M_\odot$6 show a mass-dependent evolutionary trend: those with $10^{10.5}\,M_\odot$7 retain significant rotational support (projected $10^{10.5}\,M_\odot$8 decreasing to $10^{10.5}\,M_\odot$9 with age); higher-mass quiescent galaxies lose angular momentum due to incoherent mergers and become dispersion-supported (slow rotators) (Ji et al., 2024).
  • At $z>3$0, direct measurements of velocity dispersion (e.g., $z>3$1 km s$z>3$2 for a $z>3$3 MQG) confirm that dynamical and stellar masses are consistent within uncertainties (Tanaka et al., 2023).
• Chemical Abundances:
  • High-$z>3$4 MQGs are $z>3$5-enhanced ([O/Fe]$z>3$60.2 dex), with broad stellar metallicity distributions ([Z/H]$z>3$7 to $z>3$8), matching or exceeding the enhancement seen in local massive ellipticals (Lucia et al., 11 Nov 2025).

4. Environment and Clustering Across Cosmic Time

MQGs occupy a wide range of environments, with their spatial and environmental properties varying as a function of redshift, mass, and assembly history.

• High-Redshift Diversity:
  • $z>3$9 MQGs are found in cosmic overdensities (protoclusters), but also in filamentary and even void-like regions (Binh et al., 16 Dec 2025, Lucia et al., 11 Nov 2025).
  • Approximately 50% of high-mass MQGs reside in peaks with $1.5McConachie et al., 28 Oct 2025).
  • Environmental conformity (enhanced quiescence among neighbors) observed at low-$1.5Binh et al., 16 Dec 2025).
  • Dense groups of MQGs consistent with emerging red sequences are observed at $1.5cluster cores (Tanaka et al., 2023).
• Descendants and Stochastic Assembly:
  • By $1.5$M_*\geq 10^{11}\,M_\odot$0–$M_*\geq 10^{11}\,M_\odot$1, with $M_*\geq 10^{11}\,M_\odot$2 remaining permanently quenched; rejuvenation (SF reignition) is merger-driven and more common in overdense regions (Lucia et al., 11 Nov 2025).
  • Halo mass assembly at $M_*\geq 10^{11}\,M_\odot$3 for MQGs is primarily due to dry minor mergers, with empirically calibrated relations between stellar mass, velocity dispersion, and halo mass (Zahid et al., 2019).

5. Physical Quenching Channels and Simulations

The baryonic and dynamical processes responsible for massive galaxy quenching have been investigated in detail with both hydrodynamical and semi-analytic simulations, revealing key physical insights as well as enduring limitations.

• Early-epoch (z>5) MQGs:
  • Radiative transfer (RT) hydrodynamics (e.g. Thesan-1) finds rapid, smooth accretion in the densest cosmic supernodes, fast black hole growth, and AGN-driven outflows as the dominant quenching channel, with minimal major merger involvement (Chittenden et al., 28 Apr 2025).
  • Kinetic AGN feedback injects $M_*\geq 10^{11}\,M_\odot$4–$M_*\geq 10^{11}\,M_\odot$5 erg, evacuating cold gas and driving quenching on $M_*\geq 10^{11}\,M_\odot$6–$M_*\geq 10^{11}\,M_\odot$7 Myr, tightly correlated with SMBH growth episodes.
• Cosmic Noon (1.5<z<3.0):
  • The “cosmological starvation” model posits that when host halo $M_*\geq 10^{11}\,M_\odot$8 drops below $M_*\geq 10^{11}\,M_\odot$9 Gyr$10^{11}\,M_\odot$00, SFR declines rapidly—without requiring explicit AGN feedback—matching observed quiescent fractions ($10^{11}\,M_\odot$01–$10^{11}\,M_\odot$02) for $10^{11}\,M_\odot$03–$10^{11}\,M_\odot$04 (Feldmann et al., 2016).
• Simulations and Model Performance:
  • Hydrodynamical simulations (IllustrisTNG, SIMBA, EAGLE, Magneticum): generally reproduce the abundance and structure of MQGs at $10^{11}\,M_\odot$05–$10^{11}\,M_\odot$06 but under-predict high-$10^{11}\,M_\odot$07 ($10^{11}\,M_\odot$08) MQG number densities by factors 5–100 (Baker et al., 4 Jun 2025, Lustig et al., 2022).
  • Semi-analytic models: SHARK can match observed MQG number densities at $10^{11}\,M_\odot$09 after introducing Gaussian mass/SFR scatter, but most other SAMs (GAEA, GALFORM) fail at both the abundance and mass scale (Baker et al., 4 Jun 2025, Lustig et al., 2022).
  • AGN feedback (radio/kinetic mode) is required in models to reproduce rapid and deep quenching, with merger-driven black hole growth a key trigger at high redshift (Rong et al., 2017, Forrest et al., 2019).
  • Simulations that include coupled RT and explicit modeling of cosmic reionization, such as Thesan-1, uniquely produce MQGs in the densest environments at $10^{11}\,M_\odot$10, absent in otherwise identical runs lacking RT (Chittenden et al., 28 Apr 2025).

6. Open Challenges and Future Prospects

Despite remarkable progress, significant tensions and uncertainties remain in the theoretical and observational study of MQGs.

• Selection and Purity:
  • Rest-frame UVJ diagram, commonly used for MQG selection, suffers from $10^{11}\,M_\odot$1130\% incompleteness (misses young, recently quenched galaxies) and up to $10^{11}\,M_\odot$1260\% contamination by dusty star-forming interlopers at $10^{11}\,M_\odot$13 in typical deep-field photometry (Lustig et al., 2022).
  • sSFR-based selection recovers a substantially higher abundance, especially among lower-mass ($10^{11}\,M_\odot$14) systems (Baker et al., 4 Jun 2025).
• Simulation–Observation Tensions:
  • Simulations systematically underproduce MQGs at $10^{11}\,M_\odot$15 and require more efficient, earlier AGN feedback, or additional mechanisms (cosmic ray heating, shock-induced morphological transformations) to reach observed abundances and quenching rates (Forrest et al., 2019, Straatman et al., 2013, Baker et al., 4 Jun 2025).
  • MWGs detected at $10^{11}\,M_\odot$16–$10^{11}\,M_\odot$17 have higher stellar mass fractions and shorter assembly timescales than simulated analogs, implying that theoretical models underestimate baryon conversion efficiency and/or cannot trigger rapid enough quenching (Glazebrook et al., 2017, Carnall et al., 2023).
• High-redshift Cluster Assembly:
  • The first spectroscopically confirmed MQG-centric proto-clusters at $10^{11}\,M_\odot$18 provide constraints on simultaneous and synchronized quenching across $10^{11}\,M_\odot$19 Mpc, inconsistent with current large-volume simulations (e.g., Illustris-TNG300) (Tanaka et al., 2023).
• Kinematic and Structural Evolution:
  • MQGs at $10^{11}\,M_\odot$20 undergo mass-dependent dynamical transformation from fast to slow rotators due to cumulative dry merging, but the initial origins of rotational support and its cosmic evolution remain incompletely mapped (Ji et al., 2024).

Planned and ongoing JWST NIRCam and NIRSpec surveys, in combination with ALMA deep fields, will further clarify the demographics, physical state, and evolutionary fates of the MQG population. Simulations incorporating next-generation feedback models, high dynamic range, and explicit radiative transfer will be essential to match the observed abundance, internal properties, and environment dependence of the most massive quenched galaxies at high redshift.

---

References:

See (Chittenden et al., 28 Apr 2025, Feldmann et al., 2016, Sherman et al., 2020, Binh et al., 16 Dec 2025, Tanaka et al., 2023, McConachie et al., 28 Oct 2025, Lucia et al., 11 Nov 2025, Baker et al., 4 Jun 2025, Rong et al., 2017, Zhang et al., 2019, Carnall et al., 2023, Glazebrook et al., 2017, Straatman et al., 2013, Patel et al., 2017, Lustig et al., 2022, Ji et al., 2024, Zahid et al., 2019, Xu et al., 2020, Forrest et al., 2019).