{"id":"6ce5f797-8111-4075-8293-42328ddec3bf","arxiv_id":"1908.00552","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A purely analytic galaxy formation model links halo growth to stellar mass via a time-independent efficiency and reproduces key galaxy statistics, with admitted failures at the highest redshifts.","lead":"This paper introduces a compact set of analytic equations that predict a galaxy's stellar mass from the growth of its dark matter halo, using a simple efficiency curve for star formation. If the framework holds, researchers can quickly test how feedback from stars and black holes shapes galaxy populations across cosmic time.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Neither efficiency model reproduces the GSMF at both low and high redshift; the abstract overstates high-redshift agreement.","rationale":"The paper presents a compact and useful analytic framework connecting halo growth to stellar mass, and the z=0 calibration and low/intermediate-z behavior are genuinely informative. The most load-bearing concern is the one the reader identified: the time-independent efficiency assumption, combined with calibration only to z~0, is exactly what produces the high-z underpredictions shown in Fig. 11. The paper itself acknowledges this limitation in Section 4.1, so the abstract's wording overstates what is actually demonstrated. I do not see an additional internal inconsistency that would warrant rejection; the double-Schechter discussion is largely a restatement of the assumed double-power-law efficiency, and the lack of code is a reproducibility inconvenience, not a correctness flaw. Because the reader's CONDITIONAL verdict already captures this overstatement and the need for toned-down claims, my stress-test does not change the verdict. A quantitative per-redshift goodness-of-fit check would settle the matter cleanly and should be reported in any revision.","tokens_in":27614,"tokens_out":11825,"duration_ms":117725,"concrete_test":"Recompute the Model I and Model II GSMFs from Eqs. (2), (5), (13) and (16), and evaluate reduced chi-squared per redshift bin against the Behroozi et al. (2019) standardized data used in Fig. 11, splitting at z=4 (and additionally 1 2 for z>4, and Model II yields chi-squared_nu > 2 for 14. The abstract should be toned down to 'up to z≈4' for Model I, or the claims should be qualified to describe Model I's low/intermediate-z success and Model II's high-z success separately.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper develops a fully analytic model of galaxy formation in which the growth of stellar mass is tied to the growth of dark matter haloes through an effective star formation efficiency. The efficiency is either a time-independent function of halo mass (Model I, Eq. 4) or a function of virial temperature (Model II, Eq. 8), with four parameters calibrated to the z~0 galaxy stellar mass function (GSMF) in Table 2. From analytic halo growth histories and a Press-Schechter mass function, the authors derive stellar mass growth, the GSMF, the cosmic star formation rate density, and the specific star formation rate, and compare these with observations and with EAGLE simulation variants. The paper claims that the model reproduces the shape and evolution of the cosmic SFR density, the sSFR, and the GSMF both at the present time and at high redshift, and that the double Schechter function arises naturally from two efficiency regimes in the stellar-to-halo mass relation.","tokens_in":27857,"tokens_out":6439,"duration_ms":69192,"significance":"The analytic machinery is elegant and potentially useful: the derivations in Appendix A are transparent, the model is invertible and fast, and the out-of-sample agreement with the cosmic SFR density and the sSFR is a genuine success. The comparison with EAGLE variants in Figs. 5, B1, and B2 also gives a welcome sanity check. However, the paper's headline claim about reproducing the high-redshift GSMF is not supported by either model as presented, and the double-Schechter 'explanation' is largely a restatement of the assumed double power-law parametrization. If those claims are properly qualified, the framework is a valuable contribution to analytic galaxy formation modelling.","major_comments":[{"comment":"The abstract's claim that the model reproduces the GSMF 'both at the present time and at high redshifts' is contradicted by the paper's own results. Section 4.1 states that Model I 'reproduces very well the evolution of the GSMF up to redshift z≈4, but significantly under predicts the abundance of distant galaxies,' while Model II 'provides a good fit both at low and high redshift, but the evolution is too rapid at intermediate redshift (z=1 to z=4).' No single efficiency prescription in the paper works across the full claimed range. The abstract and the similar concluding sentence in Section 5 should be revised to present the two models as complementary (Model I to z≈4; Model II at high z but with an intermediate-redshift tension), or the claim of reproducing the high-redshift GSMF should be removed.","section":"Abstract; Section 4.1; Fig. 11"},{"comment":"The 'explanation' of the double Schechter function as two efficiency regimes is largely circular. The double power-law form of epsilon_* in Eq. (4), and the resulting slope transition in Eq. (7), are fitted to the z~0 GSMF (Table 2), so the appearance of a bump at the knee of the GSMF is built into the parametrization rather than independently predicted. The paper should state explicitly that the double-Schechter-like shape is a consequence of the assumed functional form, and should frame the genuine prediction as the redshift evolution of that shape, which is tested in Fig. 11 up to z≈4. As written, the 'origin' claim overstates what the model independently establishes.","section":"Section 2.1, Eqs. (4), (7), (15); Section 4"},{"comment":"The high-redshift comparison in Fig. 11 is only qualitative: no chi-squared or other quantitative statistic is given for the z>0 panels. Given that the abstract claims a reproduction of the high-redshift GSMF, the authors should either quantify the agreement per redshift bin, or explicitly identify the panels where the model deviates by more than the observational uncertainties. This would also clarify whether the phrase 'reproduces very well' is supported across the full mass range or only at certain masses.","section":"Section 4.1; Fig. 11"}],"minor_comments":[{"comment":"The text repeatedly writes 'viral temperature' where 'virial temperature' is meant; please correct these typographical errors.","section":"Section 4.1.1"},{"comment":"The caption contains the garbled string '10th^a˘A¸S90th'; this should be '10th to 90th percentiles'.","section":"Figure 5 caption"},{"comment":"The symbol epsilon is used for the logarithmic slope of the stellar-to-halo mass relation in Eq. (1) and for the star formation efficiency in Eq. (3). These are distinct quantities, and using the same base symbol is confusing; consider denoting the slope by, for example, s(M_h,t).","section":"Eq. (1) and Eq. (3)"},{"comment":"The variable z is defined as (M_h/M_crit)^(alpha+beta), which conflicts with the standard use of z for redshift throughout the paper. The authors note the clash, but a different symbol, such as y, would be clearer.","section":"Eq. (4)"},{"comment":"The reduced chi-squared values are quoted to one decimal place; giving the raw chi-squared and the number of data points would allow readers to assess the fit quality more precisely.","section":"Table 2"},{"comment":"The text contains the typo 'In. this section'; it should read 'In this section'.","section":"Section 4.1"}],"recommendation":"major_revision","confidential_remarks":"The central derivation appears internally consistent, and the out-of-sample SFRD and sSFR comparisons are real strengths. The main issue is that the abstract overstates the high-redshift GSMF agreement in a way that conflicts with the paper's own Fig. 11, and the double-Schechter explanation is partly circular. Both are fixable with careful rewriting and more quantitative high-redshift comparisons, so I see no reason for rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is a useful analytic repackaging of abundance matching, not a new physical discovery. The authors derive closed-form expressions for halo growth and stellar mass assembly from a double-power-law efficiency, and they show that a time-independent epsilon*(Mh) with four parameters fit to the z~0 GSMF also gets the cosmic SFRD and sSFR right without fitting them. Those out-of-sample checks are real and worth taking seriously.\n\nWhat is actually new: the analytic Press-Schechter machinery, with the hypergeometric solution for M*(Mh), and an invertible set of ODEs that lets you vary feedback parameters quickly. The comparison to EAGLE variant runs (no SN, no AGN) is a nice sanity check. The framework will be convenient for people who want fast experiments on feedback physics before running simulations.\n\nWhere the paper is soft: the abstract's claim about reproducing the GSMF \"both at the present time and at high redshifts\" is contradicted by Fig. 11 and the text of Section 4.1. Model I reproduces z<4 but underpredicts galaxies at z>4; Model II does better at high-z but evolves too fast at z=1-4. The stress-test note is accurate; the claim should be qualified. The double-Schechter \"explanation\" is also partly circular: the inflection in the abundance is a direct consequence of the slope transition in Eq. (7), and the same four parameters were fit to the z~0 GSMF, so saying the model explains the knee is weaker than it sounds. And while the paper describes the model as fully analytic, no implementation is provided, so \"easily extended and inverted\" is not immediately reproducible by a reader. It is also worth noting that the final SHMR is essentially a double-power-law abundance-matching relation, so the conceptual distance from Moster et al. and Behroozi et al. is small.\n\nNone of this kills the paper. The math is internally consistent, the calibration is transparent, and the uncalibrated SFRD and sSFR agreements are genuine evidence that a time-independent mass-dependent efficiency captures the average behaviour. The central argument, that much of galaxy-halo co-evolution can be captured by a simple efficiency curve plus halo growth, holds up for z<4 at least.\n\nWho this is for: anyone doing quick phenomenological experiments on feedback, or teaching galaxy formation with a transparent analytic model. It deserves a serious referee, but the referee should ask for a toned-down abstract, a clearer separation of fitted versus predicted statements, and ideally a code release.\n\nRecommendation: engage with it and send to review with moderate revision. I would not sink it, but I would push on the overclaim.","headline":"A handy analytic wrapper around abundance matching with genuine out-of-sample checks, but the high-redshift claim in the abstract needs reining in and the double-Schechter explanation is largely built into the fitting function.","tokens_in":28432,"tokens_out":2023,"would_cite":true,"duration_ms":20806,"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":"This paper argues that a single time-independent star-formation efficiency curve, peaking at Milky-Way-scale haloes and falling off on both sides, reproduces the observed evolution of galaxies from the present day to redshift four, and…","keywords":["galaxy formation","star formation efficiency","stellar mass function","double Schechter function","feedback regulation","dark matter haloes","analytic model","cosmic star formation history"],"falsifier":"Measure the galaxy stellar mass function or the stellar-to-halo mass relation at $z=4{-}8$ with small uncertainties: if the abundance of high-redshift galaxies lies systematically above the time-independent efficiency model while the virial-temperature model also overshoots at intermediate redshifts, the central claim fails. More directly, abundance matching that shows the peak efficiency mass $M_{\\rm crit}$ shifting with redshift by more than the model's fixed value would violate Eq. (4).","tokens_in":27394,"feed_emoji":"🌌","tokens_out":5217,"duration_ms":56637,"temperature":0.7,"pith_summary":"This paper argues that most of the galaxy formation story, from the rise and fall of cosmic star formation to the characteristic knee of the galaxy stellar mass function, follows from a single time-independent function: the efficiency with which a dark matter halo of a given mass turns its infalling gas into stars. The efficiency is a double power law that peaks near Milky-Way-scale haloes (about $10^{12}\\,M_\\odot$) and falls off at low masses, where stellar feedback expels gas, and at high masses, where black-hole feedback suppresses cooling. From this one curve plus an analytic model of halo growth, the authors reproduce the cosmic star formation rate density, the specific star formation rate, and the stellar mass function from the present day to about $z=4$. They also show that the double Schechter form of the low-redshift stellar mass function is not an added ingredient but a direct consequence of the two efficiency regimes. A variant where efficiency depends on virial temperature instead of mass does better at $z>4$ but evolves too rapidly between $z=1$ and $z=4$.","feed_headline":"A fixed halo-mass efficiency law reproduces galaxy history","feed_subtitle":"Two feedback regimes, stellar winds and black holes, create the double Schechter mass function.","key_machinery":"The central object is the effective star formation efficiency $\\epsilon_*(M_h)$, modelled as a double power law (Eq. 4) with four free parameters: normalisation $\\epsilon_N$, peak halo mass $M_{\\rm crit}$, low-mass slope $\\alpha$, and high-mass slope $\\beta$. This efficiency is combined with an analytic Press--Schechter description of halo growth and abundance, including a Taylor-expansion solution of the Friedmann equations and a power-law approximation to the density-field variance $S\\propto M_h^{-\\gamma}$. The stellar mass function is then obtained through the Jacobian relation $\\varphi(M_*)=\\epsilon^{-1}\\varphi_h(M_h)$, where $\\epsilon$ is the logarithmic slope of the stellar-to-halo mass relation; the discontinuity in this slope between $1+\\alpha$ and $1-\\beta$ creates the inflection point that becomes the double Schechter bump. A separate factor $f_{\\rm SFR}$ splits stellar growth into in-situ star formation and accreted stars, allowing the star formation rate and cosmic star formation rate density to be computed.","core_discovery":"Galaxy formation can be reduced to a product of cosmology and a single astrophysical object: the effective star formation efficiency $\\epsilon_*(M_h)$, defined as the fraction of infalling baryons converted into stars. Taking this efficiency to be a fixed, time-independent double power law of halo mass, with four parameters calibrated to the present-day stellar mass function, the model reproduces the cosmic star formation rate density, the specific star formation rate of galaxies, and the galaxy stellar mass function at both low and high redshift. The paper's key explanatory result is that the double Schechter function emerges from the mapping between halo mass and stellar mass: the logarithmic slope of the stellar-to-halo mass relation changes sharply at the efficiency peak, producing an inflection point, or 'bump', at the knee of the stellar mass function. Physically, this bump is the pile-up of galaxies around the halo mass where star formation efficiency peaks, between a stellar-feedback-regulated regime at low masses and a black-hole-regulated regime at high masses.","pith_inferences":["The model's failure at $z>4$ under the fixed-efficiency assumption, noted in the paper, is itself a measurement: it brackets the redshift where the effective star formation efficiency must begin to evolve with cosmic time.","Because all ingredients are analytic and differentiable, one could invert the model directly against the observed stellar mass function to recover the efficiency function without running hydrodynamical simulations, which would make the method a practical tool for survey data.","A natural extension is to let $M_{\\rm crit}$ or the slopes $\\alpha,\\beta$ vary smoothly with time while keeping the double power-law form; the paper's two limiting cases bracket the allowed behaviour of such an extension.","The virial-temperature model's overly rapid evolution at intermediate redshifts could be cured by a mild redshift dependence of $T_{\\rm crit}$ itself, for instance through metallicity or dust effects on cooling, a modification the paper leaves unexplored."],"forward_implications":["If the fixed, mass-dependent efficiency is correct, then the redshift evolution of the galaxy stellar mass function from $z=0$ to $z\\approx 4$ is driven almost entirely by the growth and abundance of dark matter haloes, not by evolving baryonic physics.","The double Schechter function of low-redshift galaxy surveys would have a mechanistic explanation in terms of two feedback stages, making it a prediction rather than an empirical fitting function.","Supernova feedback, through the low-mass slope $\\alpha$, sets the rise and peak of the cosmic star formation history; removing it turns the cosmic star formation rate density into a power law with no peak.","Black-hole feedback, through the high-mass slope $\\beta$, is primarily responsible for the sharp knee of the stellar mass function; removing it leaves the mass function shallow and featureless at the high-mass end.","The virial-temperature version of the model identifies a physical critical temperature $T_{\\rm crit}$ at which stellar outflows stall and black-hole feedback switches on, offering a concrete threshold that can be sought in hydrodynamic simulations."],"supporting_citations":[{"why":"Supplies the analytic halo mass function and collapse threshold from which halo growth rates and abundances are derived.","marker":"Press & Schechter 1974"},{"why":"Provides the double power-law parametrisation of star formation efficiency used in Eq. (4).","marker":"Moster et al. 2010"},{"why":"Provides the analytic halo accretion-rate formula that links halo growth to the stellar mass growth rate.","marker":"Correa et al. 2015"},{"why":"Supplies the Taylor-expansion solution of the Friedmann equations underlying the analytic cosmology and growth factors.","marker":"Salcido et al. 2018"},{"why":"One of the $z\\sim 0$ galaxy stellar mass function datasets used to calibrate the four efficiency parameters.","marker":"Baldry et al. 2012"},{"why":"Second $z\\sim 0$ galaxy stellar mass function dataset used in the calibration fit.","marker":"Moustakas et al. 2013"},{"why":"Compiles the observed cosmic star formation rate density to which the model's predictions are compared.","marker":"Behroozi et al. 2013"},{"why":"Describes the reference hydrodynamical simulation whose behaviour the analytic models are shown to approximate.","marker":"Schaye et al. 2015"},{"why":"Underlies the entropy and buoyancy argument for the critical halo mass and provides the no-supernova simulation variant used in comparison.","marker":"Bower et al. 2017"}],"fun_headline_variants":["Two feedback regimes explain galaxy mass function","Analytic model maps halo growth to star formation","Stellar and black hole feedback set galaxy masses","One efficiency law reproduces cosmic galaxy growth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a halo's effective star formation efficiency is a fixed function of halo mass alone, calibrated to the $z\\sim 0$ stellar mass function and held constant with cosmic time; if that efficiency evolves, the high-redshift predictions break down, and Figure 11 shows the fixed model under-predicting galaxy abundances at $z>4$.","fun_headline_variants_meta":{"raw":{"variants":["Two feedback regimes explain galaxy mass function","Analytic model maps halo growth to star formation","Stellar and black hole feedback set galaxy masses","One efficiency law reproduces cosmic galaxy growth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000324,"raw_usage":{"total_tokens":1838,"prompt_tokens":986,"completion_tokens":852,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":602,"completion_tokens_details":{"reasoning_tokens":796}},"tokens_in":602,"tokens_out":852,"duration_ms":9055,"temperature":1.0,"reasoning_tokens":796,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:47:55.163660+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the galaxy stellar mass function or the stellar-to-halo mass relation at $z=4{-}8$ with small uncertainties: if the abundance of high-redshift galaxies lies systematically above the time-independent efficiency model while the virial-temperature model also overshoots at intermediate redshifts, the central claim fails. More directly, abundance matching that shows the peak efficiency mass $M_{\\rm crit}$ shifting with redshift by more than the model's fixed value would violate Eq. (4).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Taylor-expansion solution of the Friedmann equations underlying the analytic cosmology and growth factors."}],"review_version":1}