{"id":"3c69ac28-d424-4581-9834-49f07e956ec5","arxiv_id":"2506.06434","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":24,"one_line_summary":"Quiescent galaxy stellar mass functions in z~1 clusters show a ~2 sigma radial dependence, with a core population shaped like an enhanced field quiescent population and an outer population shaped by mass-independent environmental quenching.","lead":"Using 17 massive galaxy clusters seen at z=0.8-1.5, the authors split quiescent galaxies by cluster-centric radius and find that the core is dominated by galaxies quenched early by mass, while the outer regions show environmental quenching of infalling field galaxies. The paper matters because it offers a way to reconcile previous contradictory measurements of cluster galaxy mass functions at cosmic noon.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The '~2σ evidence' for a radially dependent quiescent SMF comes from a model the paper's own BIC analysis disfavours (ΔBIC=-0.92 for Base+α_R,Q+M*_R,Q); without a BIC-preferred radial model, the 'distinct origins' conclusion lacks a secure observational foundation.","rationale":"After reading the manuscript carefully, I agree with the reader that the paper is a valuable analysis of a unique dataset (17 GOGREEN/GCLASS clusters) and that the overall verdict should remain CONDITIONAL. The strongest parts are the radial quenched fraction trend (Figures 6-7), the comparison to field SMFs, and the mock validation in Section 3.3. However, I identify the model-selection inconsistency in Section 4.3 as the single most load-bearing concern: the headline '~2σ evidence' for a radial dependence is computed from the combined Base+α_R,Q+M*_R,Q model, which the paper's own BIC analysis disfavours (ΔBIC=-0.92). The BIC-preferred single-parameter models each favour a different parameter, and could yield different (or no) radial trends. This is an internal inconsistency in the argument, rather than an external assumption like Eddington bias, which the paper explicitly flags in Section 5.4. The toy-model decomposition (Section 5.2, Eq. 29) is also load-bearing for the interpretation: its coefficients (A_Q, B_SF) are chosen visually rather than fitted, and the 0.24 dex core M* enhancement is the same parameter whose radial dependence is subject to the BIC objection. The paper should be credited for its public data, mock recovery tests, and an honest listing of limitations; but the 'distinct origins' conclusion should be presented as conditional on reanalysis with the BIC-preferred models. I therefore do not change the reader's CONDITIONAL verdict, but I shift the emphasis of the required revision toward the model-selection issue.","tokens_in":33417,"tokens_out":9343,"duration_ms":85846,"concrete_test":"Reanalyze the data using only the BIC-preferred single-parameter radial extensions: fit (1) Base+α_R,Q and (2) Base+M*_R,Q separately. For each, compute the posterior predictive quiescent SMF shape at R/R200=0 and R/R200=1 and the significance (in σ) of the radial difference, along with the ΔBIC relative to the Base model. If neither single-parameter model yields a ≥2σ radial difference consistent with the claimed core–edge dichotomy (higher core M* for model 2; steeper edge α for model 1), the headline radial dependence is not supported by the preferred models. Also re-derive the quoted 0.24 dex M* enhancement from model (2) with full posterior uncertainties.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires a radial dependence of the quiescent SMF shape, specifically a higher core M* and a steeper non-core α. The evidence is presented via the model 'Base + α_R,Q + M*_R,Q' (Section 4.3, Figure 12), which yields a 2.1σ core-edge separation. However, the paper's own BIC analysis (Table 3) shows this combined model has ΔBIC = -0.92 relative to the Base model: BIC prefers a model in which neither α nor M* varies with radius. Adding either α_R,Q or M*_R,Q alone moderately improves BIC (ΔBIC = +3.27 and +5.97 respectively), but adding both is penalized because the parameters are degenerate and the likelihood gain does not justify the extra parameter. The paper acknowledges this degeneracy and states it 'elects to explore' the combined model anyway. Thus the headline 2.1σ significance and the 0.24 dex M* shift are derived from a model that the analysis's own criterion rejects. Moreover, no multiple-testing correction is applied across the 17 parameter combinations in Table 3, so the effective significance of any single radial trend is lower than the nominal 2.1σ. If the radial shape dependence is not robust, the distinction between core (early mass-quenching) and non-core (environmental quenching) origins collapses; a single global quiescent SMF with enhanced normalization might suffice.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyses the stellar mass functions (SMFs) of quiescent and star-forming galaxies in 17 GOGREEN/GCLASS clusters at 0.8<z<1.5. A Bayesian model simultaneously fits Schechter-function parameters alpha and M* that vary linearly with cluster-centric radius and redshift, radial projected NFW profiles for each population, and a quenched-fraction parameterisation that depends on redshift and velocity dispersion. The main findings are that the star-forming SMF shows no radial or redshift dependence, while the quiescent SMF shows about 2-sigma evidence for radial dependence: the cluster core has a quenched fraction of about 70 per cent and a quiescent SMF similar in shape to the quiescent field but with a M* about 0.24 dex higher, whereas the outer virialised regions have a quenched fraction of about 40 per cent and a quiescent SMF similar to the star-forming field. The paper interprets this as the core being dominated by early mass-quenching and the non-core by mass-independent environmental quenching.","tokens_in":33941,"tokens_out":4638,"duration_ms":46179,"significance":"If the result is robust, it provides a spatial separation of quenching channels at z~1 that would reconcile the earlier null SMF-shape result of van der Burg et al. (2020) with local environmental-quenching expectations, and it demonstrates a useful joint modelling approach for cluster SMFs. The paper has concrete strengths: mock-data tests show that the 21-parameter model recovers input parameters, the BIC analysis does support individual radial terms for the quiescent alpha and M* (Delta BIC 3.27 and 5.97), the sample extends an existing survey with a newly characterised cluster, and field comparisons use externally published SMFs rather than being derived circularly from the cluster data. The main liability is that the combined radial model used for the headline 2.1-sigma claim is disfavoured by the paper's own BIC criterion, and the interpretive toy model is not fitted with uncertainties.","major_comments":[{"comment":"The central claim of a radial dependence of the quiescent SMF shape is presented through the model 'Base + alpha_R,Q + M*_R,Q' (Figure 12), which yields the 2.1-sigma core-edge separation and the 0.24 dex M* shift, yet Table 3 reports Delta BIC = -0.92 for this model relative to Base and a BIC that is worse than either single-radial-term model (Delta BIC = +3.27 and +5.97). Since the paper itself uses BIC to decide which parameters are informative, the model used for the headline claim is disfavoured by the paper's own criterion. Please either present the radial-dependence evidence with a BIC-supported model (for example Base + M*_R,Q alone), or justify the combined model with a model-selection procedure that accounts for the degeneracy, and report a look-elsewhere-corrected significance across the 17 parameter combinations in Table 3.","section":"Section 4.3, Table 3, and Figure 12"},{"comment":"The toy model underpinning the core/edge interpretation is not fitted statistically: the coefficients A_Q and B_SF are selected by eye so that the curves reproduce the observed SMFs, the normalisation is fixed by matching the integral, and no uncertainties are propagated. Consequently the statements that the core is 'well described primarily by early mass-quenching' and the edge 'better described through mass-independent environmental-quenching' are not quantitatively supported by this model. A likelihood fit or an explicit robustness analysis that varies the field SMFs within their uncertainties would be needed to make this interpretation load-bearing.","section":"Section 5.2, Eq. (29), and Figure 15"},{"comment":"The text acknowledges that Eddington bias is not included and that mass uncertainties could change the fit values of alpha and M*. The 0.24 dex increase in core M* is the key evidence for the early mass-quenching channel, so a quantitative estimate of the Eddington-bias effect on the radial M* gradient is needed; without it the central core/non-core distinction remains vulnerable to the plausible scenario that mass uncertainties correlate with environment.","section":"Section 5.4"}],"minor_comments":[{"comment":"The sentence 'In this paper we fit the SMFs of 170.8< z <1.5 galaxy clusters' contains a typesetting error; it should read '17 galaxy clusters at 0.8<z<1.5'.","section":"Section 1"},{"comment":"The caption describes the radial model as showing 'a strong dependence with environment', while the body text reports a 2.1-sigma effect; the wording should match the statistical significance claimed.","section":"Figure 12 caption"},{"comment":"Adding a column with the number of parameters k for each model would make the BIC comparisons and the penalty term easier to interpret.","section":"Table 3"},{"comment":"The text should state explicitly which of the two individually supported radial terms (alpha_R,Q or M*_R,Q) is preferred by the data and how the radial trend differs between those two models, rather than presenting only the combined model.","section":"Section 4.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of MNRAS and the data work is careful. The main editorial concern is the internal inconsistency between the BIC-based model selection and the model used for the headline claim; this should be resolved before publication. No citation-pattern concerns."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is the radial decomposition: a flexible Bayesian SMF fit to GOGREEN/GCLASS that splits the quiescent population into a core inside roughly 0.3 R200 and a non-core outside it. That core/edge distinction is not present in the earlier vdB20 global analysis, and it is a sensible way to square the vdB20 null shape result with Gully et al.'s deeper, lower-mass data. The star-forming SMF being invariant with radius and redshift is a clean, useful check. The mock recovery tests, the careful handling of completeness and membership weights, and the addition of SpARCS-1033 are all solid. Credit where due: the authors are also unusually upfront about Eddington bias in Section 5.4 and about the degeneracy between alpha and M* in Section 4.3.\n\nThe soft spot is the one the stress-test flags, and it is real. The 2.1 sigma core-edge separation comes from the Base + alpha_R,Q + M*_R,Q model, and Table 3 shows deltaBIC = -0.92 for that model: BIC mildly prefers the base model with no radial shape dependence. The two single-parameter radial models do have positive deltaBIC (3.27 and 5.97), so there is moderate evidence that some radial shape term is needed, but because alpha and M* are degenerate the data cannot tell which, and the combined model pays the BIC penalty. No multiple-testing correction is applied over the 17 combinations in Table 3, which lowers the effective significance of any single trend. The authors know they are exploring a disfavoured model; they say so and do it for physical reasons. That is honest, but it means the headline 'distinct origins' is built on a model that the analysis's own model-selection criterion rejects. Also, Section 5.2's toy model uses coefficients chosen by eye (A_Q, B_SF), and while it is labelled an exploration, the conclusion that the core needs a 0.24 dex M* shift depends on that choice plus the neglected Eddington bias. If mass uncertainties are environment dependent, the shift could shrink.\n\nFor all that, the paper deserves a referee. The data handling, mock tests, and literature placement look sound, and the authors have not hidden the weak points. What is missing is not competence; it is a more honest framing of what is a tentative preference. I would send this to review, asking for either a demotion of the core/edge split to a tentative trend with a combined constraint, or a model comparison that marginalises over the alpha-M* degeneracy rather than exploring one disfavoured point. The target reader works on cluster quenching at z~1, and they will want to cite this for the radial decomposition regardless of whether the full interpretation survives.","headline":"The radial SMF decomposition is a genuinely useful re-analysis, but the headline 'distinct origins' rests on a 2.1 sigma trend from a model the paper's own BIC does not prefer.","tokens_in":34513,"tokens_out":2314,"would_cite":true,"duration_ms":24988,"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 quiescent galaxies in massive z~1 galaxy clusters have two distinct origins: early mass-quenching built the core population, while environmental quenching of infalling field galaxies built the outer population.","keywords":["galaxy clusters","quenching","stellar mass function","high-redshift galaxies","GOGREEN survey","GCLASS survey","quenched fraction","Schechter function"],"falsifier":"A concrete test is to rerun the stellar mass function fits with Eddington bias included, using stellar-mass uncertainties measured separately for core and non-core galaxies and for star-forming and quiescent populations; if the ~0.24 dex elevation of the core characteristic mass is absorbed or drops below 2σ significance, the claim that early mass-quenching dominates the core loses its quantitative foundation.","tokens_in":2073,"feed_emoji":"🔭","tokens_out":2221,"duration_ms":139029,"temperature":0.7,"pith_summary":"This paper asks why quiescent (non-star-forming) galaxies are so much more common inside massive galaxy clusters at 0.8<z<1.5 than in the field, and claims the answer differs between the cluster centre and its outskirts. Analysing 17 clusters with a Bayesian fit of the stellar mass functions, the authors find a ~2σ radial dependence in the quiescent population only: outside the core the quenched fraction is ~40 per cent and the quiescent mass function resembles the star-forming field, while in the core the quenched fraction reaches ~70 per cent and the mass function resembles the quiescent field with a ~0.24 dex higher characteristic mass. They read this as two distinct origins: the core was built by early mass-quenching of massive protocluster galaxies, while the outskirts were built by mass-independent environmental quenching of infalling field galaxies. If correct, this reconciles earlier null results on cluster SMF shapes with deeper observations that found steeper low-mass slopes, and it constrains the environmental channel to contribute at most roughly 40 per cent of the cluster's quiescent population.","feed_headline":"Cluster cores quench galaxies early; outer regions quench on infall","feed_subtitle":"In 17 z~1 clusters, core galaxies were quenched early by mass; edge galaxies by environment on infall.","key_machinery":"The carrying machinery is a Bayesian likelihood that fits one Schechter function per galaxy population while letting the low-mass slope $\\alpha$ and characteristic mass $M^*$ vary linearly with cluster-centric radius and redshift, through terms such as $\\alpha=\\alpha_0+\\alpha_R\\,r'+\\alpha_z\\,z'+\\alpha_{Rz}\\,z'r'$ with $r'=R/R_{200}-1$ and $z'=z-z_0$. A projected NFW profile, with a separate concentration parameter for the quiescent and star-forming populations, supplies the radial dependence of the likelihood and fixes how the quenched fraction varies with radius; the quenched fraction itself is parameterised as a function of redshift and cluster velocity dispersion, giving 21 model parameters in total. The interpretive core is a toy model that decomposes the cluster quiescent SMF into a linear combination of the field quiescent SMF (coefficient $A_Q$, representing the early mass-quenching channel) and the field star-forming SMF (coefficient $B_{SF}$, representing mass-independent environmental quenching), and this decomposition is what assigns different origins to core and non-core galaxies.","core_discovery":"The central discovery, stated on the paper's own terms, is that the shape of the quiescent stellar mass function in 0.8<z<1.5 clusters varies with cluster-centric radius while the star-forming mass function does not. Fitting Schechter functions whose slope $\\alpha$ and characteristic mass $M^*$ vary smoothly with radius and redshift yields roughly $2\\sigma$ evidence for a radial trend in the quiescent population: at the centre the quiescent SMF is similar in shape to the quiescent field but with $M^*$ elevated by about 0.24 dex, while at the cluster edge the quiescent SMF is similar in shape to the star-forming field. The Bayesian Information Criterion rates the addition of a radial term in either $\\alpha$ or $M^*$ as moderately significant, though the two are degenerate when added together. The authors interpret this as evidence that the core is populated mainly by galaxies quenched early through a mass-dependent mechanism, while the outer virialised regions are populated by mass-independent environmental quenching of the infalling field population. A toy model that writes the cluster quiescent SMF as a linear combination of the field quiescent and field star-forming SMFs supports the interpretation: the edge is well reproduced with an unenhanced mass-quenching component and $B_{SF}\\approx0.39$, while the core requires a strongly enhanced mass-quenching component ($A_Q\\gg1$) with up to $B_{SF}\\approx0.35$.","pith_inferences":["Extending the paper's logic: the radial pattern in quiescent SMF shape is effectively a fossil record of cluster assembly time, so calibrating it against simulations of protocluster accretion histories would let future wide-field surveys read assembly history directly from mass functions.","The toy-model coefficients $A_Q$ and $B_{SF}$ are tuned by eye rather than fitted; a full Bayesian fit of the decomposition against the unbinned data, with field-SMF uncertainties propagated, would replace the 'consistent with up to 40 per cent' statement with a measured posterior on the environmental-quenching fraction.","A testable extension: the early mass-quenching interpretation predicts that merger remnants or AGN signatures should be more common among massive quiescent galaxies in the core than among their edge counterparts at the same mass and redshift, which deep imaging of these and similar clusters could check.","Because the edge result rests mainly on the low-mass slope rather than on the $M^*$ offset, the environmental-quenching channel for the outskirts would likely survive even if the core offset were fully explained by mass-measurement errors; this asymmetry is implicit in the paper but never stated."],"forward_implications":["The earlier null result from the same sample, which found no difference between cluster and field SMF shapes, is explained as an effect of averaging over the whole cluster; the core/edge split reveals structure the ensemble average hid.","The core's quiescent population cannot be produced by direct quenching of the present-day field population: its SMF shape requires a mass-quenching channel that was already enhanced in the protocluster phase, with merger-driven feedback as a plausible accelerating mechanism.","Mass-independent environmental quenching, which converts star-forming galaxies into quiescent ones, is consistent with up to roughly 40 per cent of the quiescent population across the cluster, including a comparable contribution within the core itself.","The star-forming SMF is invariant with radius and redshift and matches the field, so infalling star-forming galaxies are not measurably reshaped by the cluster environment before they quench.","The quenched fraction shows little or no dependence on cluster velocity dispersion and only a weak, insignificant redshift trend, so cluster-to-cluster variation in quenching must be driven by properties the model does not yet include."],"supporting_citations":[{"why":"Supplies the 'early mass-quenching' model, the membership-correction and completeness methods, and the null cluster-field SMF result this paper's radial split directly challenges.","marker":"van der Burg et al. (2020) (vdB20)"},{"why":"Provides the field SMF fits and field quenched fractions that serve as the comparison baseline and as the two components of the toy-model decomposition.","marker":"McLeod et al. (2021)"},{"why":"Frames the core-quenching versus starvation scenarios from the same cluster sample that motivate the radial analysis.","marker":"Baxter et al. (2023)"},{"why":"Releases the GOGREEN/GCLASS data catalogues, cluster redshifts, velocity dispersions, and membership information on which the sample is built.","marker":"Balogh et al. (2021)"},{"why":"Provides the MAMPOSSt virial radii ($R_{200}$) used to normalise cluster-centric radii across the sample.","marker":"Biviano et al. (2021)"},{"why":"Formulates the mass-independent environmental-quenching model that the toy model adopts as its star-forming-to-quiescent conversion channel.","marker":"Peng et al. (2010)"},{"why":"Measures older ages for cluster quiescent galaxies, cited as direct evidence for early mass-quenching of the core population.","marker":"Webb et al. (2020)"},{"why":"Finds with deeper data a steeper quiescent low-mass slope in four of the clusters, the result this paper's core/edge separation reconciles with the earlier null.","marker":"Gully et al. (2025)"}],"fun_headline_variants":["Core quenches by mass; edge quenches by infall in z~1 clusters","z~1 clusters: core mass-quenching, edge environmental-quenching","Distinct origins: core quenched early, outer quenched on infall","Cluster cores quench via mass; outskirts quench via environment"],"cache_read_input_tokens":36352,"weakest_assumption_plain":"The two-origin interpretation rests on the assumption that the measured excess of massive quiescent galaxies in the cluster core is real and not an artefact of stellar-mass measurement errors that differ between the dense core and the outer regions; the paper itself notes that neglecting such errors would weaken the conclusion if the errors are environment-dependent.","fun_headline_variants_meta":{"raw":{"variants":["Core quenches by mass; edge quenches by infall in z~1 clusters","z~1 clusters: core mass-quenching, edge environmental-quenching","Distinct origins: core quenched early, outer quenched on infall","Cluster cores quench via mass; outskirts quench via environment"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000827,"raw_usage":{"total_tokens":3768,"prompt_tokens":1251,"completion_tokens":2517,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":867,"completion_tokens_details":{"reasoning_tokens":2434}},"tokens_in":867,"tokens_out":2517,"duration_ms":19226,"temperature":1.0,"reasoning_tokens":2434,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T05:56:39.516730+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete test is to rerun the stellar mass function fits with Eddington bias included, using stellar-mass uncertainties measured separately for core and non-core galaxies and for star-forming and quiescent populations; if the ~0.24 dex elevation of the core characteristic mass is absorbed or drops below 2σ significance, the claim that early mass-quenching dominates the core loses its quantitative foundation.","supporting_citations":[],"review_version":1}