{"id":"c07509e8-8f04-4e15-8802-736ad244033e","arxiv_id":"2502.06942","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"An abundance-matched SFR-halo accretion relation plus halo quenching reproduces low-redshift galaxy abundances and quenched fractions, with extra low-mass quenching required for the faint end.","lead":"A new semi-empirical galaxy model, DECODE, turns dark matter halo growth into star formation using an abundance-matched relation between star formation rate and halo accretion rate, then switches star formation off in massive haloes. The model matches observed galaxy mass functions and the fraction of dead galaxies at low redshift, and suggests an extra quenching process is needed for the smallest haloes.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Result 3 (low-mass quenching is needed) is not robust to the choice of input SFR function: Sect. 2.2 uses Mancuso+16, while Appendix A states Sargent+12 and Fujimoto+23 give opposite low-mass behaviors, so the ad hoc 10^11 Msun cut in Sect. 4.2.1 may be fixing an artifact of that choice.","rationale":"The TNG test in Sect. 3 is genuine supporting evidence for a monotonic SFR-HAR relation and for the abundance-matching algorithm, but it cannot validate the observed SFR function; the observational anchor remains the input. The paper itself flags the alternative-SFR-function sensitivity in Appendix A, and that flagged limitation is decisive for the third conclusion. This is exactly the reader's weakest assumption, so I agree with the reader's CONDITIONAL verdict and recommend no change to it. A rerun with all three SFR functions would settle whether the low-mass quenching requirement is robust.","tokens_in":24277,"tokens_out":7263,"duration_ms":68826,"concrete_test":"Re-run the complete DECODE pipeline (Sects. 2.2-2.4) with identical HMF, merger trees, scatter, and halo-quenching recipe, changing only the input SFR function among Mancuso+16, Sargent+12, and Fujimoto+23. For each run compute the z=0 SMF, quenched fraction, and SMHM relation, and quantify agreement with Bernardi+17, Leja+20, and Weaver+23. Then toggle the 10^11 Msun low-mass quenching cut off and on and record the change in goodness of fit for each input. If the cut is required, or its required sign changes, only for the Mancuso input, Result 3 should be reframed as a property of the chosen SFR function rather than a physical necessity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's headline claim that 'additional quenching processes in the least massive haloes are needed' (Abstract Result 3, Sect. 4.2.1) is load-bearing for the paper's physical conclusion, but it is established only for one input SFR function. In Sect. 2.2 the model adopts Mancuso et al. (2016) as the census of star-forming galaxies and removes haloes above Mh,lim = 10^12 Msun from the abundance matching. Appendix A then reports that switching to Fujimoto et al. (2023) yields a high excess of low-mass galaxies and lower quenched fractions, while Sargent et al. (2012) predicts 'opposite results.' Because Eq. (1) maps the input SFR function directly to the SFR-HAR relation used to grow galaxies, the predicted low-mass SMF excess, and therefore the hand-added 10^11 Msun low-mass quenching cut intended to cure it, is a function of the assumed SFR function rather than an independent physical inference. The high-mass halo-quenching result may still be robust since it is controlled by the bright end and the Mh,lim truncation, but Result 3 as stated is not supported unless the same requirement emerges for all SFR-function choices. The 0.4 dex quenching dispersion is also a fitted parameter, further softening the 'sufficient' language.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents DECODE, a semi-empirical model that assigns galaxy star formation rates along dark matter halo accretion histories via an abundance-matched SFR-HAR relation, integrates these SFRs to build stellar masses, and tests a halo-quenching scenario against observed SMFs, quenched fractions, and SMHM relations from z=0 to z~2. The authors validate the monotonic SFR-HAR assumption against TNG100, use the Mancuso et al. (2016) SFR function as the observational input, and conclude that halo quenching alone can reproduce the high-mass galaxy population while additional low-mass quenching is needed to match the steep low-mass end of the SMHM relation and the flat low-mass end of the SMF.","tokens_in":24616,"tokens_out":5392,"duration_ms":49240,"significance":"If the central claims hold, the paper provides a useful, transparent framework with minimal free parameters, and the TNG-based check gives independent support for the monotonic SFR-HAR hypothesis and for the abundance-matching procedure itself. The model's ability to predict quenched fractions and SMHM relations as outputs rather than inputs is a genuine strength, as is the direct integration of SFRs along merger trees. However, the headline low-mass-quenching conclusion is currently tied to one specific input SFR function, and the quenched-fraction comparison depends on a dispersion parameter that is tuned to the data; these caveats reduce the significance of the paper's third conclusion until the robustness is quantified.","major_comments":[{"comment":"Result 3, that 'additional quenching processes in the least massive haloes are needed,' is not robust to the choice of input SFR function. The abundance matching in Eq. (1) uses Mancuso et al. (2016) as the census of star-forming galaxies, and Sect. 2.2 states that switching to alternative SFR functions yields 'qualitatively similar results.' Appendix A, however, reports that Fujimoto et al. (2023) produces a high excess of low-mass galaxies and lower quenched fractions, while Sargent et al. (2012) predicts 'opposite results.' Since the low-mass excess that motivates the ad hoc 10^11 Msun quenching cut in Sect. 4.2.1 is generated by the SFR-HAR relation derived from this input, the inference is input-dependent. I request a quantitative re-analysis with at least one alternative SFR function, or a clear and quantitative justification for preferring Mancuso et al. (2016), before Result 3 is presented as a general finding. The Conclusions bullet that 'some of the latest current data on the SMF do not require this additional ingredient' should also be reconciled with the abstract's categorical statement of Result 3.","section":"Sect. 2.2 and Appendix A"},{"comment":"The quenched-fraction comparison is not fully independent because the 0.4 dex dispersion around the quenching threshold is explicitly fitted to the data ('we found a value of 0.4 dex to best suit the outputs to the data'). Since the model's claim that halo quenching is 'sufficient' to reproduce the statistics of quenched galaxies is assessed through this comparison, the fitted nature of this parameter should be stated more prominently, and the sensitivity of the high-mass conclusions to sigma_quench should be quantified rather than only described qualitatively as flattening or sharpening the quenched-fraction curves.","section":"Sect. 2.4"},{"comment":"The abundance matching removes all haloes above Mh,lim = 10^12 Msun from the HAR function, so the SFR-HAR relation at high HAR is not constrained by observations but is effectively imposed by the quenching hypothesis under test. This means the predicted high-mass SMF and quenched fraction partly encode the assumption being tested. I ask the authors to clarify how this truncation affects the claimed agreement at the high-mass end and to state explicitly what part of the high-mass prediction is a test of halo quenching rather than a consequence of the input truncation.","section":"Sect. 2.2"}],"minor_comments":[{"comment":"The notation used for the scatter is defined only in the text; please define mu = d log SFR/d log HAR directly after Eq. (1) and clarify how mu is evaluated when the SFR-HAR relation bends at low redshift, as in Fig. 6.","section":"Eq. (1)"},{"comment":"The axis labels in the upper and lower panels of Fig. 2 appear as 'Mh [M yr 1]' and should be typeset as log10(Mh/[M_sun yr^-1]) and log10(SFR/[M_sun yr^-1]) for clarity.","section":"Fig. 2"},{"comment":"The phrase 'by including the energy release, for example from strong SN feedback' overstates what is implemented; the model simply shuts down star formation in haloes below 10^11 Msun, and this should be described as an ad hoc test rather than a feedback model.","section":"Sect. 4.2.1"},{"comment":"There is a typo in the sentence 'a strong and instantaneous halo quenching could be not a realist scenario'; this should read 'may not be a realistic scenario.'","section":"Sect. 5"},{"comment":"The footnote defining 'surviving' as 'those centrals that have become satellites' is confusing; please clarify whether these are central galaxies that have survived as distinct objects after infall or something else.","section":"Sect. 2.3.2"},{"comment":"The TNG test in Sect. 3 is a self-consistency check of the abundance-matching technique using simulation inputs; the text should state explicitly that this test does not independently validate the observational SFR function or the halo-quenching recipe.","section":"Sect. 3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a useful contribution with a transparent methodology, but the robustness of Result 3 to the input SFR function is the key issue. The requested re-analysis or a substantial reframing of the conclusion is feasible within the paper's scope, so I do not recommend rejection at this stage."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. First, the core ingredient is genuinely new: they build a monotonic SFR-HAR relation by abundance matching the observed SFR function against the halo accretion-rate function, then integrate along individual halo histories. The TNG test is a real check, not just an appeal to authority, and it supports both the monotonicity and the abundance-matching procedure itself. Second, the halo-quenching story for the high-mass end survives the paper's own robustness tests, but the low-mass quenching claim does not; it depends on which SFR function you feed in, and the authors say as much in Appendix A.\n\nWhat the paper does well: the framework is transparent, uses few parameters, and predicts the SMF, SMHM relation, and quenched fractions without inputting the z=0 SMF. The distinction from TOPSEM is meaningful, and the merger-tree integration plus satellite treatment is careful. The high-mass result, that halo quenching with a threshold around 10^12 Msun flattens the SMHM and steepens the SMF bright end, is robust and consistent with other semi-empirical models.\n\nSoft spots, in proportion. The 0.4 dex quenching dispersion is fitted to the quenched-fraction data, so those comparisons are partly tuned. The hand-set 10^11 Msun low-mass cut is introduced to fix an excess that Appendix A shows reverses if you swap in Sargent+12 or Fujimoto+23; the low-mass quenching requirement is therefore a function of the Mancuso+16 SFR function, not an independent physical inference. The abstract's third result overstates this, though the Conclusions hedge correctly by noting some current SMF data do not require the extra ingredient. Removing haloes above 10^12 Msun from the abundance matching is also a modeling choice worth stress-testing, since the scatter around the quenching threshold could in principle put star-forming galaxies there.\n\nWho is this for: people building semi-empirical models or interpreting JWST/Euclid quenched fractions. The method is a useful, fast tool. The paper deserves a serious referee; it is not desk-rejectable, but the 'sufficient' framing should be softened and the low-mass claim recast as conditional on the chosen SFR function.","headline":"New abundance-matching-based SFR-HAR framework with a real TNG validation; the high-mass halo-quenching story holds up, but the paper's low-mass quenching claim is input-dependent and should be recast as conditional.","tokens_in":25171,"tokens_out":1905,"would_cite":true,"duration_ms":18021,"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":"A single monotonic star-formation–halo relation plus halo quenching explains galaxy statistics up to $z\\sim2$.","keywords":["star formation rate-halo accretion rate relation","abundance matching","galaxy quenching","halo quenching","stellar mass function","stellar mass-halo mass relation","semi-empirical model","galaxy formation"],"falsifier":"A decisive check would be a direct, per-galaxy measurement of star formation rate and host halo accretion rate at $z\\sim1$: if the rank-ordered relation is non-monotonic, bends at high accretion rates, or shows asymmetric scatter well beyond 0.4 dex, the abundance-matching anchor collapses. A cheaper calculation is to re-run the pipeline with an alternative low-mass SFR function census of the kind used in the paper's Appendix A; the paper itself reports that this reverses the low-mass conclusion, so the choice of census decides whether low-mass quenching is actually required.","tokens_in":24081,"feed_emoji":"🌌","tokens_out":8410,"duration_ms":69248,"temperature":0.7,"pith_summary":"This paper argues that a single, monotonic relation between a galaxy's star formation rate and the rate at which its dark matter halo accretes mass, calibrated by matching observed star-forming number counts to simulated halo counts, can account for the bulk of the local star-forming galaxy population. It then asks whether the classic halo-quenching picture, where star formation shuts off once a halo crosses roughly $10^{12}\\,M_\\odot$, explains the quenched population. The model integrates these assigned star formation rates along dark-matter merger trees to predict stellar mass functions, quenched fractions, and the stellar mass–halo mass relation. The authors find that monotonic star formation rate–halo accretion rate matching plus halo quenching reproduces the high-mass end and the quenched statistics up to $z\\sim2$, but that the low-mass end of the stellar mass function and the steep low-mass stellar mass–halo mass slope require additional quenching in the least massive haloes. The payoff is a minimal, transparent description of galaxy assembly that marks precisely where new physics must enter.","feed_headline":"A monotonic SFR–halo relation fits galaxy populations to z~2.","feed_subtitle":"DECODE shows halo quenching handles massive galaxies; the lowest-mass galaxies need extra quenching.","key_machinery":"The load-bearing object is the star formation rate–halo accretion rate (SFR–HAR) relation: a monotonic mapping constructed by abundance matching the observed SFR number density with the simulated halo accretion rate number density via a one-equation formalism, with Gaussian scatter in SFR at fixed HAR as the only adjustable parameter. This relation converts dark-matter accretion histories into galaxy star formation histories. The same mapping is validated for self-consistency: feeding a large cosmological hydrodynamical simulation's own SFR and HAR functions into the abundance-matching recipe recovers the simulation's mean SFR–HAR relation at all redshifts studied, and a comparison of co-evolving SFR and HAR tracks shows no appreciable time delay. Quenching enters as a truncation: when the host halo crosses the threshold mass, the assigned SFR is dropped to zero, and the built-up stellar mass thereafter grows only by mergers.","core_discovery":"The central discovery is that the galaxy–halo connection can be derived, not fitted: rank-order the observed star formation rate function against the halo accretion rate function at each redshift, with a single 0.4 dex scatter, and the resulting relation assigns a star formation rate to every halo in a cosmological merger tree. Stellar mass is then the time integral of that rate, truncated instantaneously when the host halo crosses a quenching threshold near $10^{12}\\,M_\\odot$ with 0.4 dex of dispersion. This recipe predicts the local and $z\\sim1\\text{--}2$ stellar mass functions, quenched fractions, and stellar mass–halo mass relations in broad agreement with data, supporting three specific claims: (1) a monotonic star formation rate–halo accretion rate relation is sufficient for the number densities of the bulk of star-forming galaxies; (2) halo quenching alone accounts for the quenched population and for the flat high-mass end of the stellar mass–halo mass relation and the steep bright end of the stellar mass function; and (3) matching the steep low-mass end of the stellar mass–halo mass relation and the flat faint end of the stellar mass function requires an extra quenching process in haloes below roughly $10^{11}\\,M_\\odot$.","pith_inferences":["A natural next test would replace the rank-ordered abundance matching with a direct, per-galaxy measurement of SFR and host halo accretion rate at a few redshifts, to check whether the 0.4 dex symmetric scatter and monotonicity hold outside the simulation used for calibration.","If low-mass quenching is confirmed, the same machinery could quantify how much of the faint-end slope is produced by supernova-driven outflows versus satellite stripping by rerunning the model with a mass-dependent quenching threshold.","Since a delayed-quenching variant with 1–2 Gyr exponential decline leaves the outputs essentially unchanged, the model's discriminating power appears to reside in the threshold mass and its scatter rather than the quenching timescale, an interpretation worth testing with more gradual quenching recipes.","Starting the assembly at higher redshift using early-Universe luminosity functions would test whether the overabundance of high-redshift quenched galaxies seen in recent surveys is a missing second channel for quenching or a problem with the assumed SFR census."],"forward_implications":["If correct, the bulk of the observed galaxy population at $z\\lesssim2$ can be predicted from halo accretion histories alone, without a large calibrated set of baryonic feedback parameters.","The high-mass end of the stellar mass function and the flat high-mass stellar mass–halo mass slope are explained by halo quenching, so the model locates the need for additional feedback mainly in low-mass haloes and at $z\\gtrsim2$.","Because stellar masses are integrals of assigned SFRs, the model produces self-consistent predictions for the stellar mass–halo mass relation and its scatter, including decreasing scatter with halo mass, which can be tested against weak-lensing and dynamical measurements.","The same pipeline can be rerun with updated empirical SFR censuses from new surveys; Appendix A already shows that alternative inputs reverse the low-mass conclusion, making the low-mass quenching requirement a directly testable data-driven claim.","The predicted quenched fractions at $z\\sim2$ fall below the latest high-redshift data, suggesting the model can be used to single out the redshift range where a second quenching channel becomes necessary."],"supporting_citations":[{"why":"Supplies the abundance-matching formalism (their Eq. 37) that converts the SFR and HAR functions into the SFR-HAR relation without free fitting.","marker":"Aversa et al. (2015)"},{"why":"Provides the observed star formation rate function used as the primary input to the abundance matching.","marker":"Mancuso et al. (2016)"},{"why":"Defines the shock-heating halo quenching threshold near 10^12 Msun and the cold-stream physics that motivates the quenching recipe.","marker":"Dekel & Birnboim (2006)"},{"why":"Establishes the shock-heating mechanism that sets the quenching threshold applied in the model.","marker":"Birnboim & Dekel (2003)"},{"why":"Provides the SatGen analytic merger trees from which each halo's accretion rate history is drawn.","marker":"Jiang et al. (2021)"},{"why":"Supplies the algorithm underlying the SatGen merger-tree construction.","marker":"Parkinson et al. (2008)"},{"why":"Gives the halo mass function used to sample the parent halo catalogue.","marker":"Tinker et al. (2008)"},{"why":"Provides the TNG100 simulation used to test whether the assumed monotonic SFR-HAR relation and the abundance-matching recovery hold.","marker":"Nelson et al. (2019)"},{"why":"Defines the previous DECODE model whose satellite merging recipe and fudge factors are reused here.","marker":"Fu et al. (2022)"},{"why":"Presents the TOPSEM semi-empirical model whose abundance-matching philosophy DECODE builds on.","marker":"Boco et al. (2023)"}],"fun_headline_variants":["SFR–halo relation derived, not fitted, in DECODE","Halo quenching only part: low-mass haloes need extra","Abundance matching yields monotonic SFR–halo link","One SFR–halo relation explains galaxy masses to z~2","DECODE: quenching recipe for massive plus low-mass haloes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything downstream rests on treating the adopted star formation function as a complete census of star-forming galaxies and on removing all haloes above $10^{12}\\,M_\\odot$ from the abundance matching, so the inferred SFR-HAR relation is calibrated only below the quenching threshold and then extrapolated to every halo's earlier accretion history.","fun_headline_variants_meta":{"raw":{"variants":["SFR–halo relation derived, not fitted, in DECODE","Halo quenching only part: low-mass haloes need extra","Abundance matching yields monotonic SFR–halo link","One SFR–halo relation explains galaxy masses to z~2","DECODE: quenching recipe for massive plus low-mass haloes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000223,"raw_usage":{"total_tokens":1549,"prompt_tokens":1132,"completion_tokens":417,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":748,"completion_tokens_details":{"reasoning_tokens":327}},"tokens_in":748,"tokens_out":417,"duration_ms":4735,"temperature":1.0,"reasoning_tokens":327,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T14:16:35.169384+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive check would be a direct, per-galaxy measurement of star formation rate and host halo accretion rate at $z\\sim1$: if the rank-ordered relation is non-monotonic, bends at high accretion rates, or shows asymmetric scatter well beyond 0.4 dex, the abundance-matching anchor collapses. A cheaper calculation is to re-run the pipeline with an alternative low-mass SFR function census of the kind used in the paper's Appendix A; the paper itself reports that this reverses the low-mass conclusion, so the choice of census decides whether low-mass quenching is actually required.","supporting_citations":[{"cited_title":"2016, , 823, 128","cited_arxiv_id":null,"evidence_quote":"Provides the observed star formation rate function used as the primary input to the abundance matching."},{"cited_title":"2008, , 383, 557","cited_arxiv_id":null,"evidence_quote":"Supplies the algorithm underlying the SatGen merger-tree construction."},{"cited_title":"TOPSEM, TwO Parameters Semi Empirical Model: Galaxy Evolution and Bulge/Disk Dicothomy from Two-Stage Halo Accretion","cited_arxiv_id":"2307.13036","evidence_quote":"Presents the TOPSEM semi-empirical model whose abundance-matching philosophy DECODE builds on."}],"review_version":1}