REVIEW 4 major objections 7 minor 2 cited by
A hydrodynamical CLONE of the Virgo cluster: II. Confronting observed and synthetic galaxy population twins in a dense environment
T0 review · 4 major / 7 minor · reviewed 2026-07-13 · grok-4.5
Pith's one-line read A constrained Virgo hydro simulation jointly recovers multiple galaxy trends in a dense environment, even without averaging or observational error bars.
desk verdict A careful multi-trend confrontation of one constrained Virgo hydro twin that jointly recovers the main environmental-quenching patterns down to low mass; residual high quenched fractions and low Z are flagged but not yet cleanly separated from subgrid effects. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The Virgo CLONE: a ~30 Mpc zoom-in RAMSES hydrodynamical simulation (down to ~350 pc) of a velocity-constrained Virgo counterpart, run with Horizon-AGN subgrid physics (star formation, SN and spin-dependent AGN feedback) without further cluster-specific recalibration, used as a controlled numerical twin rather than a stand-alone predictive model.
What would settle it
A re-run of the same constrained Virgo initial conditions with a Kroupa-based metal yield and mass return (or altered feedback efficiencies) that either closes the metallicity and quenched-fraction offsets while preserving the joint trends, or destroys the simultaneous match to the mass- and environment-dependent relations.
Extended reading notes
Core claim
Even without ensemble averaging or folding in observational uncertainties (aside from projection), the synthetic Virgo population jointly recovers key observational relations: reduced SFR at fixed stellar mass in the cluster, mass- and environment-dependent quenched fractions with a minimum at intermediate mass, a mass-metallicity break, dark-matter stripping of quenched galaxies, and gas-depletion-driven quenching of low-mass systems extending far beyond the virial radius.
Load-bearing premise
That the adopted feedback recipes and IMF, applied without cluster-specific recalibration to a single constrained realization, are adequate so residual mismatches can be read as physical rather than numerical.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents a multi-faceted confrontation between the galaxy population of a constrained zoom-in hydrodynamical Virgo replica (CLONE; Horizon-AGN-like subgrid physics, ~350 pc resolution) and observed relations for star-formation density, SFR/sSFR–M*, quenched fraction versus stellar and halo mass and cluster-centric distance, stellar-to-halo mass, gas-to-stellar mass, and mass–metallicity/age–metallicity trends, including limited cosmic-time evolution of quenching. Building on Paper I’s match to Virgo’s assembly history and mass/luminosity distributions (including an M87 counterpart), the authors argue that—despite slightly low metallicities and high quenched fractions—the synthetic population jointly recovers key observational trends even without ensemble averaging or folding in most observational uncertainties (aside from projection). They interpret residual offsets as largely projection/selection/systematics rather than fatal model failure, and position the replica as a useful laboratory for jellyfish galaxies and cluster-core gas dynamics, with a forthcoming Kroupa-yield run flagged to address metallicity and quenching rates.
Significance. If the multi-trend qualitative agreement holds under the stated caveats, the work is a valuable observation-driven test of current AMR hydro models on a well-observed, constrained Virgo twin rather than an ensemble-averaged random cluster. Strengths include: (i) consistent internal definitions applied uniformly down to lower stellar masses than many prior synthetic-population studies; (ii) explicit discussion of projection along a filament-aligned line of sight; (iii) joint coverage of SFR, quenching, gas content, DM stripping, and metallicity–age relations; and (iv) no cluster-specific recalibration of Horizon-AGN recipes, so residual mismatches can in principle diagnose the subgrid model. The constrained IC setup and companion ICM papers further increase the scientific return of this single realization. The result is significant for environmental quenching studies and for assessing whether current feedback/yield choices are adequate for cluster satellites, even if absolute zero-points (Z, QF) remain imperfect.
major comments (4)
- §3.5 and Fig. 4 (left): The claim that projection plus missing observational SFR uncertainties “can largely” explain the high quenched fractions (≳95% at high mass; still elevated after projection as dashed lines) is load-bearing for interpreting residuals as non-numerical, but remains qualitative. Please quantify the projected QF reduction more carefully (e.g., multiple sightlines, not only the filament-aligned one; estimate of the Euclid Collaboration: Cleland et al. ~20% effect applied to your sample), and state explicitly how much of the offset relative to Wetzel et al. (2012) / GOGREEN remains after those corrections. Soften or condition any language that residual high QF is primarily physical/selection if the residual is still large.
- §3.7 and Fig. 6: Absolute metallicities are shifted by ~0.26 dex to match Gallazzi et al. (2005), with a mix of known RAMSES under-enrichment, Salpeter yield, and observational aperture/S/N systematics invoked. Trend agreement (break near the minimum-QF mass, scatter) is the scientifically useful claim; the manuscript should more clearly separate “trend recovery” from “absolute Z reliability” and avoid implying that the zero-point is already adequate for quantitative chemical evolution or ICM–galaxy metal budget work until the planned Kroupa-yield run is shown. Atlas3D agreement without shift is helpful—make the aperture comparison more quantitative.
- Abstract, §3.5–3.6, and Conclusion: The “valuable tool” claim for jellyfish and core-gas studies rests on joint multi-trend fidelity of a single constrained realization with unaltered Horizon-AGN free parameters (ε_*=0.02, η_SN=0.2, metal yield 0.1, spin-dependent jets). That is defensible for qualitative environmental trends, but the manuscript should more explicitly bound what is and is not yet trustworthy (e.g., absolute QF and Z zero-points; low-mass end near resolution; total-gas vs HI comparison in Fig. 5). A short dedicated limitations paragraph tying residual offsets to free parameters and single-halo variance would strengthen, not weaken, the paper.
- §3.6 and Fig. 5: The finding that suburb galaxies are on average more gas-poor than core galaxies at fixed M* is interesting but counter to a simple “stronger stripping deeper in” picture and is used to support pre-processing/strangulation. Given that simulated gas is total gas in a fixed aperture (not HI) and Mun et al. (2021) is a stripping-selected sample, please strengthen the robustness checks (aperture variation already mentioned; also cold vs hot gas split if available) and clarify whether this radial trend could be affected by how quenched systems are assigned or by resolution of diffuse gas. The interpretation is plausible but currently under-supported relative to its weight in the abstract’s point #6.
minor comments (7)
- §3: Quenched definition (sSFR < 10^{-11} yr^{-1} over 100 Myr) is stated and literature variation is noted; please also show briefly (text or appendix) sensitivity of Fig. 4 to a factor-of-few threshold change or to a main-sequence-offset definition, since QF zero-points are central to the residual-offset discussion.
- Fig. 1–2 and IMF conversions: Salpeter rescaling of observational SFRs/masses is appropriate; state conversion factors once in a single methods paragraph and ensure all comparison samples (VESTIGE, Castignani, Goubert et al. sims) are treated consistently.
- Fig. 2 bottom / Virgo comparisons: VESTIGE bias toward core/RPS galaxies is noted; a short table of sample selection differences (mass limit, radius cut, SFR tracer) would help readers weight the slope comparisons.
- Appendix A / Fig. A.1: Quenched-fraction panel still sits ~25% above Wetzel et al.; cross-reference the main-text projection discussion so the appendix does not read as an unreconciled discrepancy.
- Presentation: Abstract and §1 use dense semicolon-separated claim lists (#1–#6); consider a short bullet or numbered summary box for readability. Fix minor typos (e.g., “CLONEof”, “hydrodynamicalCLONEof”, “V ogelsberger”, “by about+0.04 dex”).
- §2.2: BH formation and spin-dependent jet efficiencies are summarized; a compact table of subgrid parameters (with Horizon-AGN references) would aid reproducibility and comparison to TNG/EAGLE/Hydrangea results cited in Fig. 4.
- Backsplash galaxies (5% of sample, §3): Stating they are not treated differently is fine; one sentence on whether excluding them changes QF or gas trends would close a natural question.
Circularity Check
No load-bearing circularity: fixed Horizon-AGN subgrid physics plus constrained ICs produce galaxy trends that are compared to independent external surveys; residual offsets are acknowledged rather than fitted away.
full rationale
The paper runs a single constrained zoom-in hydrodynamical simulation of a Virgo counterpart (ICs from prior Sorce et al. peculiar-velocity constrained realizations; subgrid recipes taken unchanged from Horizon-AGN: SF efficiency 0.02, η_SN=0.2, metal yield 0.1, Bondi-Hoyle + spin-dependent AGN jets, Salpeter IMF) and confronts the resulting synthetic population against external observational relations (SDSS/Brinchmann, Woo, Gallazzi, Wetzel, VESTIGE/Boselli, Atlas3D/McDermid, Castignani, etc.). The claimed multi-trend agreement (lower SFR at fixed M*, quenched-fraction minimum at intermediate mass, mass-metallicity break, DM stripping of quenched systems, gas-depletion quenching of low-mass galaxies out to ~3 Rvir) is therefore an a-posteriori comparison, not a quantity forced by construction or by a self-citation uniqueness theorem. Self-citations establish only the large-scale environment and the prior mass/luminosity match (Sorce et al. 2021); they do not supply the SF/quenching/metallicity trends under test. Residual mismatches (slightly low Z, high quenched fractions) are explicitly flagged and deferred to a forthcoming Kroupa-yield run rather than absorbed into free parameters. No fitted input is renamed a prediction, no definitional identity is presented as derivation, and no ansatz is smuggled via citation to force the central claim. Score 1 reflects only the ordinary, non-load-bearing self-citation of the simulation setup.
Assumptions & free parameters
free parameters (4)
- star-formation efficiency =
0.02
- SN mass fraction η_SN and metal yield =
η_SN=0.2, yield=0.1
- quenched sSFR threshold =
10^{-11} yr^{-1}
- BH formation and feedback efficiencies
assumptions (5)
- domain assumption Planck 2014 cosmology (Ωm=0.307, ΩΛ=0.693, Ωb=0.048, H0=67.77, σ8=0.829)
- domain assumption Horizon-AGN subgrid recipes (cooling, SF, SN, AGN) without cluster-specific recalibration
- domain assumption Salpeter IMF for stellar particles, magnitudes, and metal yield
- ad hoc to paper A single constrained Virgo realization is sufficient for qualitative multi-trend comparison
- ad hoc to paper Projection along a filament-aligned line of sight plus absence of observational SFR uncertainties explain most of the high quenched-fraction offset
Cite this review
Pith. "Pith review of A hydrodynamical CLONE of the Virgo cluster: II. Confronting observed and synthetic galaxy population twins in a dense environment." pith.science (2026). https://pith.science/paper/55ELOX7S
@misc{pith2026260323606,
author = {Pith},
title = {Pith review of: A hydrodynamical CLONE of the Virgo cluster: II. Confronting observed and synthetic galaxy population twins in a dense environment},
year = {2026},
howpublished = {\url{https://pith.science/paper/55ELOX7S}},
note = {Machine review of arXiv:2603.23606}
}
read the original abstract
Galaxy clusters offer powerful laboratories for studying galaxy evolution in dense environments. In this context, the Constrained LOcal and Nesting Environment (CLONE) project provides a zoom-in hydrodynamical simulation of the Virgo cluster, including active galactic nucleus and supernova feedback, with a resolution down to 350~pc, designed to mirror Virgo's observed properties. Previous work showed that this replica and Virgo share the same history, mass, and luminosity distributions including the central M87. This study examines several observational relations extending to lower stellar masses than previous synthetic-population studies: star formation density, (specific) star formation rate, metallicity, and quenched fraction of galaxies as a function of stellar mass and cluster-centric distance. Despite the slightly low metallicity and the sufficiently high quenched fraction, simulated galaxies reproduce key observational trends even without averaging or accounting for observational uncertainties, aside from the consideration of projection effects: At fixed stellar mass, cluster galaxies form fewer stars than field counterparts; Most galaxies are quenched, except for intermediate-mass or isolated galaxies; Low-mass galaxies are highly quenched, thus implying a sharp metallicity drop, and low metallicity does not imply youth; Quenching occurs earlier for the most massive and the smallest galaxies than for those of intermediate mass, at least until they enter the cluster; Quenched galaxies have undergone dark matter stripping; Gas depletion drives quenching, especially in low-mass galaxies and the farther from the cluster center they are. Overall, the synthetic population jointly reproduces multiple observational trends, making it a valuable tool to probe processes from jellyfish galaxies to cluster-core gas dynamics.[Abridged]
Forward citations
Cited by 2 Pith papers
-
Cutting with precision -- Leveraging Collapse Volumes to generate the next generation of zoom-in initial conditions
Presents a new technique using collapse volumes to generate stable, low-contamination zoom-in initial conditions for SLOW constrained simulations, with tests on 20 regions showing pristine volumes beyond 6 virial radii.
-
Learning the Universe: Constrained simulations of the Coma galaxy cluster -- I. Radial X-ray and Compton-y signatures
50 constrained simulations of Coma cluster analogues reproduce the observed radial X-ray surface brightness and Compton-y profiles within the scatter expected from environment and assembly history.
Reviewed July 13, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.