{"id":"1027af2d-e084-4ef9-b080-030abc36f069","arxiv_id":"1908.01995","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"An analytic model predicts that stellar mass alone sets the star formation history of galaxies, with environment playing a minimal role in quenching at all redshifts down to z=0.","lead":"This paper uses an analytic galaxy formation model to argue that stellar mass, not environment, is the main driver of galaxy quenching from z=1.5 to the present. The claim matters because it speaks to a long-running nature versus nurture debate in galaxy evolution, but the result is largely built into the model's assumptions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The main result is pre-encoded in the model: Eq. 4 sets initial SFR from stellar mass only and Eq. 2 sets the quenching timescale from stellar mass only, so the predicted environment-independence of the SFR-M* relation is a direct consequence of the ansatz rather than an empirical finding.","rationale":"The reader's weakest_assumption identifies exactly the load-bearing point: the model's input equations enforce the output. I agree with that assessment. In good faith, the paper is transparent about the model and does compare against observations, and the agreement in Fig. 2 is genuine independent support. However, the strongest claim -- 'stellar mass is the main driver at any redshift probed' -- is not a falsifiable inference from this setup, because no environmental variable enters the SFR initialization or the quenching timescale. The only free environmental channel, satellite quenching, is modeled with a random fraction of the mass-dependent tau_c and a delay that are both halo-mass independent; thus even environmental quenching does not produce an environment-dependent SFR at fixed mass. The admitted failure to reproduce the quiescent stellar mass density evolution (Fig. 3, right) further weakens the empirical case. The proposed alternative-model test would show whether the conclusion is robust to the assumed functional form. Since this concern supports the reader's REJECT verdict rather than overturning it, no adjustment is needed.","tokens_in":13944,"tokens_out":9680,"duration_ms":103887,"concrete_test":"Re-run the analytic model with an environment-dependent satellite quenching timescale, e.g. tau_s = f * tau_c * (M_halo / 10^14 Msun)^beta, re-calibrating the free parameters to the same stellar mass function and the Wetzel et al. (2012) quiescent fractions as in Fig. 2. If a non-zero beta with comparable fit quality produces a clear dependence of the SFR-M* relation on halo mass, then the paper's central conclusion is an artifact of the assumed form of Eq. 3; if the environment-independent form is strongly preferred, the conclusion gains support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim -- that stellar mass is the main driver of quenching at all 0<z<1.5 and environment is marginal -- is not an independent result but a direct consequence of the model's ansatz. The initial SFR is assigned from Eq. 4, which depends only on M* and z, and the central quenching timescale in Eq. 2, tau_c = 10^11.7 (M*/Msun)^-1 (1+z)^-1.5, also depends only on M* and z. For satellites, Eq. 3 uses tau_s = f*tau_c with f uniform in [0.1,0.5] and a delay of 2-4 Gyr, with no dependence on host halo mass, local density, or clustercentric distance. Hence SFR(t) for every galaxy is fixed by M* at z_match/form, redshift, and, for satellites, infall time and a random draw; it cannot respond to environment. Figures 4-8 therefore show the environment-independence of SFR/SSFR at fixed M* and the strong M*-dependence at fixed halo mass by construction. The only empirical anchors that could break this are the comparisons to observed quiescent fractions (Fig. 2, which matches at z~0.1) and the split stellar mass density (Fig. 3, right). The latter explicitly underpredicts the quiescent stellar mass density at all redshifts; the authors attribute most of the tension to SSFR versus color selection, but this is the closest direct test of the quenching law and it is not cleanly passed. The paper thus offers a self-consistent model prediction, not a measurement of which process drives quenching.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents an analytic model of galaxy formation coupled to merger trees from an N-body simulation, using subhalo abundance matching to assign stellar masses and a prescribed star formation history. Star formation rates are initialized from an observed SFR-M* relation at z_match or z_form (Eq. 4) and decay exponentially with a quenching timescale that depends only on stellar mass and redshift for centrals (Eq. 2), while satellites follow a delayed-then-rapid quenching with a timescale that is a random fraction (0.1-0.5) of the central timescale (Eq. 3). The paper then reports that the SFR/SSFR-M* relations are independent of environment (halo mass or clustercentric distance) at all redshifts 0<z<1.5 for both star-forming and quiescent galaxies, and that SFR/SSFR at fixed halo mass depend strongly on stellar mass. The authors conclude that stellar mass is the main driver of galaxy quenching at all probed redshifts, with a minimal role for environment.","tokens_in":14391,"tokens_out":3163,"duration_ms":35073,"significance":"The paper is transparent in its construction: the analytic model is simple, the free parameters are stated, and the comparison to the observed quiescent fraction in Figure 2 is a useful check. The authors are also honest about the disagreement in Figure 3 and about the inability of their analysis to address mass-environment interdependence. However, the central claim is not an empirical measurement but a consequence of the model's prescriptions. Since both the initial SFR (Eq. 4) and the quenching timescale (Eq. 2) depend only on stellar mass and redshift, with satellite quenching also mass-dependent through Eq. 3, the environment-independence of the SFR-M* relation is built in by construction. The paper would be valuable as a self-consistent model prediction, but it does not constitute a test of whether mass or environment drives quenching.","major_comments":[{"comment":"The central result is pre-encoded in the model. Equation 4 assigns the initial SFR from stellar mass and redshift only, and Equation 2 sets the central quenching timescale as tau_c = 10^11.7 (M*/Msun)^-1 (1+z)^-1.5, also independent of environment. For satellites, Equation 3 uses tau_s = f*tau_c with f uniform in [0.1,0.5] and a delay of 2-4 Gyr, with no dependence on host halo mass, local density, or clustercentric distance. Consequently, the environment-independence of the SFR/SSFR-M* relations in Figures 4 and 5 is a direct algebraic consequence of these prescriptions, not an empirically inferred result. To make the claim non-circular, the paper would need to show that the observables could have differed if an environmental term had been included, or confront the model with independent data that specifically constrain environmental variations of SFR at fixed stellar mass. As written, the conclusion that 'stellar mass is the main driver of galaxy quenching' is a restatement of the model's input assumptions.","section":"Section 3, Figure 3 (right panel)"},{"comment":"The model underpredicts the stellar mass density of quiescent galaxies and overpredicts that of star-forming galaxies at all redshifts. The authors attribute much of the disagreement to the difference between a color-based selection and their SSFR cut, but this comparison is the closest direct test of the quenching law in the paper, and it is not cleanly passed. Because the quiescent fraction is the main environmental observable matched in Figure 2, the failure to reproduce the quiescent mass density weakens the robustness of any conclusion about environmental quenching efficiency. The paper should provide a quantitative assessment of the discrepancy (e.g., by stellar mass and redshift) and demonstrate that it does not affect the central conclusions, rather than asserting that it is not going to invalidate the analysis.","section":"Section 4 (last paragraph)"},{"comment":"The paper explicitly acknowledges that its analysis cannot confirm or rule out a mutual dependence between mass and environmental quenching efficiencies, and that the information in Figure 8 is not enough for a fair comparison with studies such as Balogh et al. (2016), Darvish et al. (2016), and Kawinwanichakij et al. (2017). This is a significant limitation on the breadth of the main conclusion. If mass and environmental quenching are interdependent, then the statement that 'stellar mass is the main driver of galaxy quenching at any redshift probed' is not supported by the modeling, because the model's mass-only prescription cannot produce such interdependence. The conclusions should be framed as predictions of the model rather than as a resolution of the nature/nurture debate.","section":"Section 3.1 (Figure 5, low-mass galaxies)"},{"comment":"The paper notes a trend for low-mass star-forming galaxies (log M* < 9.3) in which SFR and SSFR decrease with increasing halo mass at lower redshifts, but dismisses it as less than 0.2 dex and within the typical scatter. Given that the model has no environmental term in the star formation prescription, any such residual trend must originate from the assembly history and satellite fraction, not from direct environmental quenching. The claim that 'the environment does not play any role' is therefore too strong; the model can only place an upper limit on the magnitude of such an effect, not demonstrate its absence.","section":"Section 2.1 (Eq. 2)"}],"minor_comments":[{"comment":"The abstract and introduction contain several typographical errors ('tipically', 'generelly', 'instantantaneous', 'dispruted', 'stronlgy', 'altogehter', 'enviromental') that should be corrected before publication.","section":"Section 3, Figure 2"},{"comment":"The comparison in Figure 2 uses the SDSS quiescent fraction of Wetzel et al. (2012) at z~0.1, but the model is not compared to any observed quiescent fraction at higher redshift. Given that the conclusions extend to z=1.5, an additional comparison to higher-redshift quiescent fractions (e.g., from COSMOS/UltraVISTA as used in Figure 3) would make the model validation more convincing.","section":"Section 3.2, Figure 6"},{"comment":"The text states that the average difference in SFR for star-forming galaxies is ~1.6 dex between the least and most massive stellar mass bins at z=0, but it never provides the corresponding error bars or the significance of this difference. Adding quantitative uncertainties to the quoted values would make the comparison with observations more robust.","section":"Section 3.2, Figure 6"}],"recommendation":"reject","confidential_remarks":"The circularity concern is fundamental: the paper's main conclusion is an output of the model's mass-only prescriptions (Eqs. 2-4). This cannot be fixed by rewriting alone; the model would need to be extended to include environmental dependencies in the star formation law and then tested against independent data. Given the journal's standards, the central claim as presented does not hold as an empirical inference, so rejection is appropriate. The authors do show some self-awareness in Sections 3 and 4 about the quiescent-mass-density mismatch and the mass-environment interdependence, which is commendable, but the framing of the abstract and conclusions overstates what the model demonstrates."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe punchline: the central claim—that stellar mass drives quenching at all z<1.5 and environment is marginal—is to a large degree built into the model's input assumptions. Equation 4 assigns each galaxy an initial SFR from the observed SFR–M* relation, and Equation 2 makes the quenching timescale a power law in M* and redshift only. For satellites, the added quenching timescale is a random fraction of that same mass-dependent timescale, with no input from host halo mass or local density. So when Figures 4–8 show environment-independent SFR–M* relations and a strong M* dependence at fixed halo mass, that is mostly a projection of the ansatz, not an empirical finding. The stress-test note is right on this.\n\nWhat the paper does well: it is transparent about the model, it matches the observed stellar mass function and the quiescent fraction at z~0.1 (Fig. 2), and it openly acknowledges that it underpredicts the stellar mass density of quiescent galaxies (Fig. 3, right). The authors also do a reasonable sensitivity check on the satellite delay and quenching fraction. As a self-consistent model that extends a known framework to z<1.5, it is cleanly executed. The only genuinely new bit is the low-redshift extension, but that follows from the same prescriptions.\n\nThe soft spots: the main one is interpretive. The paper frames the result as a conclusion about the universe, but it is a model prediction. The mismatch in quiescent stellar mass density is not fully resolved; the authors attribute most of it to SSFR vs color selection, but it remains a direct test of the quenching law and it is not cleanly passed. I also find the claim that 'environmental processes must act on very short timescales' too strong given that the model's satellite quenching is artificially rapid and environment-blind by construction. That said, the paper is not incoherent—it is just over-interpreted.\n\nWho this is for: readers interested in semi-analytic quenching models will find a clear, reproducible example of how mass-dependent quenching prescriptions produce mass-dominated outcomes. Observers will not find a new constraint.\n\nRecommendation: I would not cite it as evidence about what drives quenching. But it deserves a serious referee, because the question matters and the paper is well-built; a referee can ask the authors to reframe the claims as model predictions and to strengthen the comparison with observed quiescent mass densities. Send it to review with that expectation.","headline":"The paper's main claim is largely a projection of its own mass-dependent quenching prescriptions, but the model is transparent and the work is useful as a model prediction, not as empirical evidence.","tokens_in":14895,"tokens_out":2772,"would_cite":false,"duration_ms":27674,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Stellar mass, not environment, quenches galaxies down to z=0","keywords":["galaxy quenching","stellar mass","environment","star formation rate","specific star formation rate","quiescent galaxies","analytic model","subhalo abundance matching"],"falsifier":"Measure the specific star formation rate of low-mass star-forming galaxies ($8.5 < \\log M_* < 9.25$) in the cores of massive clusters and in the field at the same redshift near $z\\sim0.5$. The model predicts no systematic offset at fixed stellar mass; a difference larger than the typical 0.2 dex scatter would contradict the environment-independent SFR–$M_*$ relation.","tokens_in":13752,"feed_emoji":"🌌","tokens_out":6037,"duration_ms":54249,"temperature":0.7,"pith_summary":"This paper tries to settle which of two suspected causes, a galaxy's stellar mass or its surroundings, drives the shutting down of star formation. Using an analytic model grafted onto a dark-matter merger tree, it finds that the star formation rate and specific star formation rate at a given stellar mass are the same in every environment, from field to massive clusters, across redshifts $0<z<1.5$. The same quantities, however, change by large factors when stellar mass changes at fixed environment. The authors conclude that stellar mass is the main driver of galaxy quenching at all redshifts probed, extending a result previously argued mainly for $z>1$, and that environmental processes must act too quickly to alter the star formation of active galaxies, only raising the chance that a galaxy is already dead.","feed_headline":"Stellar mass, not environment, quenches galaxies down to z=0","feed_subtitle":"An analytic model finds star formation is set by stellar mass, not halo mass, from z=1.5 to today.","key_machinery":"The machinery is an analytic quenching model built on a subhalo abundance matching assignment of galaxies to haloes in a high-resolution N-body merger tree. Each galaxy is placed on an observed, redshift-dependent SFR–$M_*$ relation (Eq. 4, from Tomczak et al. 2016) and then its star formation decays exponentially with a quenching timescale $\\tau_c = 10^{11.7}\\,(M_*/M_\\odot)^{-1}\\,(1+z)^{-1.5}$ Gyr for centrals; satellites follow the same decay for a 2–4 Gyr delay after infall and then quench on a faster timescale $\\tau_s$ drawn as 0.1–0.5 times $\\tau_c$. Because the quenching timescales depend on stellar mass and redshift but not on halo mass or clustercentric distance, all environment dependence enters only through whether a galaxy becomes a satellite and how long it stays active. This is what lets the paper read its output as evidence that mass, not environment, drives quenching.","core_discovery":"The central claim is that stellar mass, not environment, controls when galaxies stop forming stars. In the model, the SFR–$M_*$ relation and the specific SFR–$M_*$ relation are independent of halo mass for both star-forming and quiescent galaxies at every redshift probed, and the same holds when distance from the cluster core replaces halo mass as the environmental measure. Conversely, the SFR at fixed halo mass depends strongly on stellar mass, with about $\\sim1.6$ dex separating the lowest and highest stellar-mass bins at $z=0$. The picture the authors draw is that environmental quenching is real but effectively instantaneous and binary: it does not bend the star formation rates of still-active galaxies, it only increases the probability that a galaxy belongs to the quiescent population.","pith_inferences":["The model builds part of the answer into its inputs: because $\\tau_c$ and $\\tau_s$ depend on stellar mass and redshift but not on environment, the finding that environment is marginal is partly a consequence of the prescription, and the sharpest independent test is the observed quiescent fraction at fixed stellar mass across environments.","A concrete target: low-mass star-forming galaxies in cluster cores should have the same specific star formation rate as field galaxies at the same mass and redshift; a cluster survey in H$\\alpha$ around $z\\sim0.5$ could check this directly.","The paper states it cannot yet assess whether mass- and environment-quenching efficiencies are mutually dependent, as some recent cluster studies claim; comparing quenched fractions at fixed stellar mass with clustercentric distance would directly address that gap."],"forward_implications":["At fixed stellar mass, predicted star formation rates of both star-forming and quiescent galaxies are identical across halo masses from $z=1.5$ to $z=0$, with scatter around 0.2 dex.","At fixed halo mass or fixed distance from a cluster core, predicted SFR varies by roughly 1.6 dex between low- and high-stellar-mass bins at $z=0$, a direct signature of mass quenching.","Environmental quenching, if present, must be fast and nearly binary: it leaves the SFR of active galaxies unchanged while raising the quiescent fraction, as seen in observed quiescent fractions.","The result extends mass-quenching dominance to the present day, not only to $z>1$ as commonly claimed."],"supporting_citations":[{"why":"Supplies the analytic model that populates haloes with galaxies and evolves them along merger trees.","marker":"Contini et al. 2017a"},{"why":"Introduces the exponential tau-decay quenching prescription and its mass- and redshift-dependent timescale.","marker":"Contini et al. 2017b"},{"why":"Provides the N-body simulation and dark-matter merger tree the model runs on.","marker":"Kang et al. 2012"},{"why":"Fixes the initial redshift-dependent SFR–stellar mass relation (Eq. 4) from which every galaxy's star formation begins.","marker":"Tomczak et al. 2016"},{"why":"Supplies the exponentially decaying star formation history that the central-galaxy quenching prescription adapts.","marker":"Noeske et al. 2007"},{"why":"Provides the delayed-then-rapid satellite quenching scheme the model revises.","marker":"Wetzel et al. 2013"},{"why":"Provides observed quiescent fractions used to calibrate and compare the model at $z\\sim0.1$.","marker":"Wetzel et al. 2012"},{"why":"Provides observed stellar mass densities used to check the model's redshift evolution.","marker":"Muzzin et al. 2013"},{"why":"Establishes the mass-versus-environment quenching separation that the paper extends and compares against.","marker":"Peng et al. 2010"}],"fun_headline_variants":["Mass rules: galaxies switch off by stellar mass alone","Environment barely matters: mass drives quenching","Quenching is mass-driven, not environment-driven","Stellar mass, not halo, flips off galaxies","One parameter sets galaxy quenching: stellar mass"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conclusion rests on the assumption that every galaxy's initial star formation rate comes from an environment-independent observed SFR–stellar mass relation and then falls on a timescale set only by stellar mass and redshift; if the initial rate or the decay time actually depend on environment, the result would change.","fun_headline_variants_meta":{"raw":{"variants":["Mass rules: galaxies switch off by stellar mass alone","Environment barely matters: mass drives quenching","Quenching is mass-driven, not environment-driven","Stellar mass, not halo, flips off galaxies","One parameter sets galaxy quenching: stellar mass"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000368,"raw_usage":{"total_tokens":1998,"prompt_tokens":990,"completion_tokens":1008,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":606,"completion_tokens_details":{"reasoning_tokens":937}},"tokens_in":606,"tokens_out":1008,"duration_ms":8974,"temperature":1.0,"reasoning_tokens":937,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:56:26.813615+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the specific star formation rate of low-mass star-forming galaxies ($8.5 < \\log M_* < 9.25$) in the cores of massive clusters and in the field at the same redshift near $z\\sim0.5$. The model predicts no systematic offset at fixed stellar mass; a difference larger than the typical 0.2 dex scatter would contradict the environment-independent SFR–$M_*$ relation.","supporting_citations":[{"cited_title":"P., & Elahi, P","cited_arxiv_id":null,"evidence_quote":"Provides the N-body simulation and dark-matter merger tree the model runs on."},{"cited_title":"R., Quadri, R","cited_arxiv_id":null,"evidence_quote":"Fixes the initial redshift-dependent SFR–stellar mass relation (Eq. 4) from which every galaxy's star formation begins."},{"cited_title":"2013, ApJ, 777, 18 13","cited_arxiv_id":null,"evidence_quote":"Provides observed stellar mass densities used to check the model's redshift evolution."},{"cited_title":"J., Kovaˇ c, K., et al","cited_arxiv_id":null,"evidence_quote":"Establishes the mass-versus-environment quenching separation that the paper extends and compares against."}],"review_version":1}