{"id":"a3d9273e-a25a-435b-aa8e-9edd9dc279c4","arxiv_id":"1908.08808","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"The mass-radius relation of stellar systems is the intersection of constant-density collapse tracks with a redshift-dependent maximum halo mass from the halo mass function.","lead":"This paper tries to explain why galaxies and galaxy clusters of very different sizes line up on a single mass-radius plot instead of scattering randomly. It argues the line comes from the physics of collapsing gas clouds combined with the maximum mass that dark matter halos can reach at each cosmic time.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The CGS boundary is calibrated rather than predicted: Ns=10^-2 is the total SDSS galaxy density, not an ETG central-halo density, and Eq. (32) contains a free normalization; the claimed coincidence with the Bernardi MRR may therefore be tuning.","rationale":"The paper makes a genuinely interesting proposal: the near-universal MRR across ~10 orders of magnitude in mass may be set by the upper edge of the halo mass function (CGS) intersecting collapse tracks. Independent support comes from the qualitative agreement of Illustris, monolithic, and early-hierarchical models, and from the fact that different observational samples give similar slopes. The reader's conditional verdict is appropriate. I focused on the CGS normalization because every quantitative version of the central claim passes through it. The authors are transparent: Sec. 4.2 labels the one-halo-one-ETG hypothesis 'a strong assumption' and describes Ns as 'a very crude estimate.' That transparency does not remove the difficulty. If Ns is really the total galaxy density, then it overcounts ETG hosts by including satellites and non-ETGs; the resulting CGS is pulled to higher masses/radii than an ETG-specific boundary. Since the paper states that lower Ns shifts the locus away from the data, the successful fit may reflect tuning. The same applies to the free normalization in Eq. (32), whose value is not derived and whose gamma changes between the text and Fig. 9 caption (10 vs 50). The proposed test, replacing total galaxy counts with central-ETG halo counts, is directly implementable and would settle whether the boundary is predicted or fitted. If the test fails, the qualitative picture may survive, but the claim of an analytical demonstration of the MRR should be weakened. Therefore the verdict stays CONDITIONAL, with the condition being an independent determination of Ns_ETG.","tokens_in":36498,"tokens_out":5352,"duration_ms":61093,"concrete_test":"Recompute the CGS using an empirical central-ETG halo density rather than total-galaxy density: take a volume-limited SDSS subsample with the Bernardi et al. (2010) selection (fracDeV=1, b/a>0.6, Ms >= 10^10 M_sun), count central ETGs, divide by the survey volume, and use that Ns in n(MD,z)=Ns to re-derive log MD(z) and eqs. (12)-(14). If the resulting CGS at 9.5 <= log Ms <= 12.5 departs from eq. (2) by more than its own rms (~0.09 dex), the claimed coincidence is a calibration choice. If it still tracks eq. (2), the one-halo-one-ETG objection is quantitatively harmless.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The load-bearing quantitative step is the placement of the Cosmic Galaxy Shepherd (eqs. 12-14) and the 'analytical demonstration' in Fig. 9 that its intersections with the M-MRR manifold reproduce the Bernardi et al. (2010) MRR. That placement rests on Ns = 10^-2 halos per (Mpc/h)^3, which the authors derive from the total SDSS galaxy count (~10^6 galaxies over ~10^8 Mpc^3). The text explicitly acknowledges (Sec. 4.2) that this assumes each halo hosts exactly one galaxy and that the galaxy is an ETG, and calls this 'a strong assumption.' Real SDSS samples contain satellites and a large spiral/dwarf population, so the density of central ETG-hosting halos is plausibly 10-30 times lower. Because the CGS moves to lower masses/radii when Ns is lowered, a smaller, more realistic Ns will displace the predicted curve from the edge of the ETG data. In addition, Eq. (32) contains an unspecified normalization MN, and Sec. 5 uses gamma=10 in the text but gamma=50 in the Fig. 9 caption, so the boundary can be slid along the mass-radius plane. With two adjustable normalizations in the chain, the agreement in Fig. 9 is evidence of consistency, not a parameter-free derivation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper (the third in a series) aims to explain the observed Mass-Radius Relation (MRR) spanned by globular clusters, early-type galaxies, and galaxy clusters/groups. The authors assemble literature data (Burstein et al. 1997; Bernardi et al. 2010; WINGS) and compare them with hydrodynamical models (Illustris, plus older monolithic and early-hierarchical simulations). They propose that the observed MRR is the envelope of two complementary effects: (i) lines of constant initial collapse density (M-MRRs, described by eqs. 10 and 31) that map halo mass and redshift to stellar mass and radius, and (ii) a redshift-dependent upper mass limit for halos, which they call the Cosmic Galaxy Shepherd (CGS), derived from the halo growth function via the condition n(MD,z)=Ns (eqs. 12-14). The intersections of the M-MRR manifold with the CGS boundary are claimed to reproduce the empirical ETG MRR (eq. 2), implying that objects along the MRR are in virial equilibrium and passive evolution. The paper is explicitly exploratory: several assumptions are labeled 'strong,' 'crude,' or 'provisional,' and the authors caution that the quantitative location of the CGS depends on uncertain ingredients.","tokens_in":36836,"tokens_out":2946,"duration_ms":33850,"significance":"If the proposed picture is correct, it provides a conceptually appealing and unifying explanation for why objects spanning ~11 orders of magnitude in stellar mass lie on a single, relatively narrow MRR: the relation emerges from the convolution of collapse physics (initial density) with cosmological halo growth (maximum halo mass at each epoch). The paper's qualitative argument that the MRR is a boundary rather than a fundamental scaling relation is interesting and worth testing. The authors deserve credit for using a wide range of observational data, comparing independent simulation schemes, and providing explicit analytical fits (e.g., eqs. 12-14, 31-33) that make the proposed mechanism falsifiable. However, the quantitative claims are not yet parameter-free: the CGS location depends on a number density Ns estimated from the total SDSS galaxy count rather than the ETG-hosting halo density, on a provisional z-independent stellar-to-halo mass relation, and on a free normalization in the Press-Schechter cut-off mass. These admitted uncertainties directly affect the central figures (Figs.","major_comments":[{"comment":"The placement of the Cosmic Galaxy Shepherd uses Ns = 10^-2 halos per (Mpc/h)^3, derived from the total SDSS galaxy count (~10^6 galaxies over ~10^8 Mpc^3). The text itself calls the assumption that each halo hosts exactly one early-type galaxy 'a strong assumption.' Because the SDSS sample includes satellites and a large spiral/dwarf population, the number density of halo central galaxies that are early types is plausibly an order of magnitude or more lower. Since the condition n(MD,z)=Ns is steeply decreasing in MD, reducing Ns by a factor of 10-30 shifts the CGS to lower halo masses by roughly 0.5-1 dex, moving the predicted locus away from the Bernardi et al. (2010) points. Thus, the apparent coincidence between the CGS and the observed MRR in Fig. 8 is a tuning statement unless a separate, observationally motivated estimate of the ETG-hosting halo number density is provided.","section":"Section 4.2, eqs. (12)-(14) and Fig. 8"},{"comment":"Equation (32) defines the cut-off mass as MCO_T = MN × (1+z)^(-6/(n+3)), with MN described only as 'a suitable normalization mass scale.' In the Press-Schechter formalism the normalization is set by the power-spectrum amplitude (e.g., sigma_8) and the collapse threshold, so MN is not a free parameter; here it is left undetermined. Additionally, the text of Section 5 states that the curves are computed for gamma = 10 (MT = gamma * MCO_T), while the caption of Fig. 9 states that the dashed black lines correspond to 'total mass equal to 50 × MCO(z).' This discrepancy means that the maximum-mass boundary can be slid in the mass-radius plane by changing either MN or gamma. With two effective free normalizations in this step, the 'analytical demonstration' in Fig. 9 is not a parameter-free prediction; it is an exercise in matching the data with adjustable constants.","section":"Section 5, eqs. (32)-(33) and Fig. 9"},{"comment":"The stellar-to-halo mass ratio is a load-bearing ingredient because it converts the CGS halo masses into stellar masses via eq. (10). The adopted relation, log m = 0.062 log MD + 0.429, is fitted to Illustris galaxies at z = 0 and has no redshift dependence. The authors acknowledge this is provisional and note that other relations (e.g., Fan et al. 2010, Shankar et al. 2006, Girelli et al. 2020) disagree strongly at low halo masses and at high redshifts. Because the CGS is built from intersections over a wide redshift range, the z=0-only m(MD) may introduce a systematic bias in the predicted Ms and hence in the shape of the CGS in the MR-plane. The paper should quantify how the CGS changes when the redshift-dependent relations (A.1) or the Girelli et al. (2020) relations are used instead.","section":"Appendix A, eq. (A.5), and Section 4.1"},{"comment":"The dissipation-less collapse derivation reproduces the slope of the observed MRR (eq. 23), but the zero-point constant k_gamma in eqs. (25)-(27) is explicitly left undetermined ('cannot be determined unless the constant k is specified by fixing the initial conditions'). Consequently, the agreement in slope between eq. (26) and the Bernardi et al. (2010) relation does not by itself constrain the MRR zero point; the latter is absorbed into k_gamma. This limits the strength of the claim that eq. (2) represents the locus of ideal dissipation-less collapse, because the zero-point remains a free parameter in the derivation.","section":"Section 5, eqs. (23)-(27)"}],"minor_comments":[{"comment":"The abstract contains several typographical errors, including 'comsic' for 'cosmic' and 'evolution' in the final sentence where 'evolutions' or 'evolution processes' would be clearer. These should be corrected.","section":"Abstract"},{"comment":"The paper uses 'Section 3' twice: once for the description of simulations and once for 'The MRR of collapsing proto-galaxies.' The section numbering in the text should be made consistent.","section":"Introduction, Section 3 heading"},{"comment":"Equations (13) and (14) are cubic polynomials in log Ms with many digits; it would be helpful to state the valid mass range for each fit and to give the fitting uncertainties, as is customary for empirical scaling relations.","section":"Section 4.2, eq. (13)"},{"comment":"The caption states that the dashed black lines correspond to '50 × MCO(z),' while the body text says the calculation uses gamma = 10. This inconsistency should be resolved to avoid confusion about which parameter value was actually used.","section":"Figure 9 caption"},{"comment":"The conclusions section skips from item (iv) to item (vi); the numbering should be corrected or the missing item should be restored.","section":"Section 6, item (vi)"},{"comment":"The description of the SCDM cosmology in Chiosi & Carraro (2002) could be confusing to readers because the rest of the paper adopts the Illustris Λ-CDM parameters; the distinction between the cosmological parameters used in the older models and those used in the present analysis should be clarified.","section":"Appendix B.1"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is honest about its limitations, but the central quantitative claim (that the CGS is derived from cosmology and predicts the observed MRR) rests on several adjustable inputs: Ns, m(MD), and the normalization MN in eq. (32). As written, the paper is more an interesting suggestion than a rigorous derivation. I would encourage the authors to recast the CGS as an illustrative envelope and to either estimate Ns from a halo occupation model for ETGs or show how the results vary with Ns and MN. The paper is within scope for A&A and the qualitative picture is likely to stimulate discussion, but it needs substantial revision before it makes a solid quantitative claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the paper's central image is probably right, but the demonstration is partly a calibration. The mass-radius relation from globular clusters to clusters of galaxies looks like the envelope set by the redshift-dependent upper cutoff of the halo mass function, convolved with constant-initial-density collapse tracks. That's a real idea and worth engaging with. However, the quantitative location of this 'Cosmic Galaxy Shepherd' is adjusted to the data, not derived from first principles. The authors are honest about this in places, but the paper's framing—'analytical demonstration'—overreaches.\n\nWhat's genuinely good: the data compilation spans eleven orders of magnitude and brings together Burstein, Bernardi, and WINGS samples with Illustris in one plane. The comparison among monolithic, early-hierarchical, and fully hierarchical models showing they land on similar M-MRR loci is a useful result. The paper also flags its own weak points: Sec. 4.2 calls the one-halo-one-ETG assumption strong, labels the SDSS density estimate crude, and eq. (26) admits k_gamma is undetermined. That candor is to their credit.\n\nSoft spots, in proportion: the stress-test concern lands. Ns = 10^-2 per (Mpc/h)^3 comes from the total SDSS galaxy count, not from central ETG-hosting halos. Including satellites and spirals, the true number density of relevant halos could be ten to thirty times lower, which slides the CGS off the observed edge. The stellar-to-halo mass relation in eq. (A.5) is provisional and redshift-independent, and the paper itself notes that Girelli et al. provide redshift-dependent fits. Eq. (32) contains an unspecified normalization MN, and there is an internal inconsistency: Section 5 uses gamma=10 in eq. (33), while the Fig. 9 caption says gamma=50. None of these destroy the qualitative picture, but they mean the agreement in Fig. 9 is consistency after tuning, not a parameter-free prediction.\n\nWho gets value: researchers working on scaling relations, the halo mass function, or the claimed structural continuity from star clusters to galaxy clusters. I'd cite the data compilation and the qualitative envelope argument, but not the specific quantitative boundary without further checks.\n\nRecommendation: send to peer review. A serious referee should ask the authors to separate what is fitted from what is predicted, adopt a defensible ETG halo number density, resolve the gamma discrepancy, and either fix or remove the undetermined k_gamma. With that, this becomes a solid synthesis paper rather than a new discovery claim. My own verdict would be conditional accept after major revision.","headline":"A plausible but partially tuned explanation for the universal mass-radius locus; worth refereeing if the authors clearly separate fitting from prediction.","tokens_in":37408,"tokens_out":3593,"would_cite":true,"duration_ms":38231,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The mass-radius relation of galaxies and star clusters is the intersection of constant-density collapse tracks with a cosmic maximum-halo-mass boundary, the 'Cosmic Galaxy Shepherd'.","keywords":["mass-radius relation","early-type galaxies","galaxy clusters","globular clusters","halo growth function","Cosmic Galaxy Shepherd","virial equilibrium","passive evolution"],"falsifier":"Take a mass-complete, volume-limited sample of quiescent early-type galaxies and compare its high-mass edge on the mass-radius plane with the Cosmic Galaxy Shepherd calculated from the survey's own halo number density. If the edge falls inside the zone of avoidance at the expected number density, or if the edge does not move with redshift as the adopted cut-off mass changes, the proposed mechanism is wrong. A simpler check is to recompute $N_s$ using only early-type galaxies; if the resulting boundary no longer tracks the observed relation, the coincidence rests on the one-galaxy-per-halo assumption.","tokens_in":1861,"feed_emoji":"🌌","tokens_out":4161,"duration_ms":116661,"temperature":0.7,"pith_summary":"This paper asks why globular clusters, dwarf and early-type galaxies, and galaxy clusters and groups all crowd onto a narrow mass-radius relation instead of spreading across the plane. Its answer is that the relation is the combined product of two evolutions: the stellar component moves along collapse tracks of fixed initial density, while the dark-matter halo component follows a cosmological growth function that sets a maximum halo mass at each redshift. Translated onto the mass-radius plane, that maximum mass becomes a boundary the authors call the Cosmic Galaxy Shepherd, and the observed relation is the intersection of the collapse tracks with that boundary. Objects on the intersection are in virial equilibrium and are evolving passively, which explains both the changing slope of the relation and the existence of a zone of avoidance below it.","feed_headline":"Cosmic 'Shepherd' sets the galaxy mass-radius relation","feed_subtitle":"A single cosmic cutoff explains why objects from star clusters to galaxy clusters share one track.","key_machinery":"The central machinery is the Cosmic Galaxy Shepherd: the locus obtained by intersecting the halo growth function $n(M_D,z)$ at a fixed number density $N_s$ with the analytic collapse relation of eq. (10), which ties a stellar system's half-mass radius to its halo mass, stellar-to-halo mass ratio, Sersic profile and formation redshift. A second ingredient is the manifold of constant-initial-density tracks, the M-MRRs, summarized by eq. (31) as lines whose slope and zero point depend on the formation redshift. The observed mass-radius relation is the intersection of these two families of curves, and can also be reproduced by the Press-Schechter cut-off mass as a function of redshift; along the intersection, systems satisfy the virial condition and evolve passively.","core_discovery":"The paper claims that the observed mass-radius relation of stellar systems is not set by the virial theorem alone, nor by any single formation mechanism, but emerges where two ingredients cross. The first ingredient is a manifold of model mass-radius relations, each traced by objects of the same initial density (or formation redshift); the second is the redshift-dependent cut-off mass of the halo mass function, evaluated along a constant halo number density of $N_s = 10^{-2}$ halos per $(\\mathrm{Mpc}/h)^3$. The cut-off line, named the Cosmic Galaxy Shepherd, runs close to the empirical mass-radius relation of the SDSS early-type galaxy sample and merges smoothly into the regions populated by globular clusters and galaxy clusters. Along this line objects are in virial equilibrium and passive evolution, and the slope of the line varies from about 0.5 to about 1 with increasing mass; galaxies with active star formation and strong winds, mainly dwarfs, lie above it.","pith_inferences":["If the Shepherd is recalibrated using only the number density of early-type galaxy hosts rather than all SDSS galaxies, its position will shift; this turns the paper's claimed coincidence into a quantitative, testable prediction rather than a posterior fit.","The same boundary logic should apply to any population with a known collapse relation and a formation cutoff, such as globular cluster systems within a single galaxy or satellites within a cluster; checking whether their mass-radius envelopes show analogous shepherds would extend the idea beyond homogeneous galaxy samples.","Surveys with different volumes or depths should see the high-mass edge of the mass-radius relation move: deeper volume coverage reaches rarer, more massive halos at fixed redshift and should push the envelope to larger masses and radii.","A direct numerical convolution of the full halo mass function with the manifold of collapse tracks, rather than the linear intersection used in Section 5, would give a sharp test of the predicted slope and zero point across the whole mass range."],"forward_implications":["The slope and curvature of the observed mass-radius relation, from roughly 0.5 for low-mass systems to about 1 for the most massive ones, are inherited from the shape of the halo growth function rather than from the virial theorem itself.","The zone of avoidance below the relation is a cosmic-time statement: at each redshift, halos massive enough to populate that part of the plane are statistically too rare to be seen.","Compact massive galaxies at high redshift should lie near the Cosmic Galaxy Shepherd for their epoch, above the present-day relation, rather than contradicting it; the paper indicates this will be examined in a companion study.","Dwarf galaxies are expected to sit above the relation while star formation and galactic winds inflate them, and to settle onto it once they quench and reach mechanical and thermal equilibrium.","The same construction spans stellar masses from globular clusters to galaxy clusters, about eleven orders of magnitude, with a single physical origin."],"supporting_citations":[{"why":"Supplies the reference observational MRR of about 60,000 SDSS early-type galaxies with which the predicted boundary is compared.","marker":"Bernardi et al. (2010)"},{"why":"Provides the halo growth function $n(M_D,z)$ from which the Cosmic Galaxy Shepherd and the redshift-dependent cut-off mass are derived.","marker":"Luki´c et al. (2007)"},{"why":"Gives the analytic collapse relation (their eq. 10) linking halo mass, stellar mass, formation redshift and half-mass radius that maps halo cut-offs onto the MR-plane.","marker":"Fan et al. (2010)"},{"why":"Contributes the monolithic hydrodynamical models and initial-density tracks that define the M-MRR manifold and the earlier extension of the ETG relation.","marker":"Chiosi & Carraro (2002)"},{"why":"Supplies the halo mass-function formalism and cut-off mass whose redshift variation is used to construct the boundary and the intersection demonstration.","marker":"Press & Schechter (1974)"},{"why":"The companion study establishing the stellar-component side, virial equilibrium combined with luminosity evolution, that the present paper folds into the combined-evolution picture.","marker":"D'Onofrio et al. (2019b)"},{"why":"Provides the large-scale hydrodynamical simulations used to calibrate the stellar-to-halo mass ratio and to compare the predicted MRR with simulations.","marker":"Vogelsberger et al. (2014a,b)"},{"why":"Supplies the dissipation-less collapse scaling $R \\propto M^{(5+n)/6}$ used to derive the MRR slope from first principles.","marker":"Gott & Rees (1975)"}],"fun_headline_variants":["Cosmic shepherd line sets mass-radius relation","Halo cutoff explains mass-radius track across scales","One cosmic line rules stellar system sizes","Shepherd line: why star clusters to clusters align"],"cache_read_input_tokens":39424,"weakest_assumption_plain":"The construction assumes that every halo counted at the adopted number density $N_s=10^{-2}$ per $(\\mathrm{Mpc}/h)^3$ hosts one and only one early-type galaxy of the kind in the comparison sample; because $N_s$ is taken from total SDSS galaxy counts rather than the early-type fraction alone, a lower true number density of early-type-hosting halos would move the Cosmic Galaxy Shepherd and could turn the agreement into a tuned coincidence.","fun_headline_variants_meta":{"raw":{"variants":["Cosmic shepherd line sets mass-radius relation","Halo cutoff explains mass-radius track across scales","One cosmic line rules stellar system sizes","Shepherd line: why star clusters to clusters align"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00047,"raw_usage":{"total_tokens":2357,"prompt_tokens":978,"completion_tokens":1379,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":1322}},"tokens_in":594,"tokens_out":1379,"duration_ms":11499,"temperature":1.0,"reasoning_tokens":1322,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:41:42.542660+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a mass-complete, volume-limited sample of quiescent early-type galaxies and compare its high-mass edge on the mass-radius plane with the Cosmic Galaxy Shepherd calculated from the survey's own halo number density. If the edge falls inside the zone of avoidance at the expected number density, or if the edge does not move with redshift as the adopted cut-off mass changes, the proposed mechanism is wrong. A simpler check is to recompute $N_s$ using only early-type galaxies; if the resulting boundary no longer tracks the observed relation, the coincidence rests on the one-galaxy-per-halo assumption.","supporting_citations":[{"cited_title":"& Carraro , G","cited_arxiv_id":null,"evidence_quote":"Contributes the monolithic hydrodynamical models and initial-density tracks that define the M-MRR manifold and the earlier extension of the ETG relation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the dissipation-less collapse scaling $R \\propto M^{(5+n)/6}$ used to derive the MRR slope from first principles."}],"review_version":1}