{"id":"76c46818-2d43-48fa-b9d5-ac89ec45e328","arxiv_id":"2607.09848","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Star formation rates required by the evolving stellar mass function yield a main sequence that agrees with JWST spectroscopy and theory but differs from prior concordance relations.","lead":"This paper infers the star-forming main sequence of galaxies from the observed buildup of the stellar mass function across cosmic time, using JWST-era data. The resulting relation better matches independent spectroscopy and theory than older 'concordance' compilations, implying many prior analyses need revisiting.","discovery_kind":"extension","skeptic_critique":{"model":"grok-4.5","headline":"Merger absorption α(m⋆) calibrated only for z~0 progenitors is the load-bearing extrapolation for high-z SFRs.","rationale":"The Reader correctly isolates the Behroozi α extrapolation as the weakest assumption. The paper is otherwise carefully constructed: continuous Schechter fits, public code, independent spectroscopic and theoretical benchmarks, and a forward-evolution sanity check (App. A) all support the claim. Because the merger term is literature-standard and the paper already flags the asymmetry, the concern does not overturn the result; it simply keeps the verdict CONDITIONAL pending a high-z merger calibration. No stronger internal inconsistency or circularity is present. The concrete test above would settle whether the concern actually moves the high-z MS enough to weaken the JWST concordance.","tokens_in":35255,"tokens_out":690,"duration_ms":6493,"concrete_test":"Re-run the progenitor-tracking pipeline of §3.3 (steps 1–6) twice: once with the fiducial Behroozi α and once with α shifted by ±0.15 (the Clauwens/Wang SF-dependent range) at all z>2. Recompute the MS grid and the sSFR(z) track at m⋆=10^9.5 M⊙ (Fig. 1). If the high-z (z>3) points move outside the gray 16–84% SMF posterior band or lose overlap with the Clarke+25 Hα/UV points, the claimed spectroscopic agreement is not robust to the merger model.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that the SMF-required MS matches JWST/NIRSpec (Clarke+25) and SED (Leja+22/Simmonds+25) results rests on the continuity solution of Eq. (14). That solution subtracts the merger absorption term A = n_SF ln(10) (dz/dt) α(m⋆) (Eq. 13), where α is taken from Behroozi+13 and is stated to be valid only for backward tracking of z~0 progenitors (§3.2). The paper itself notes that forward/backward evolution is asymmetric and that α depends on the descendant SFR (Clauwens+16; Wang+23), yet applies the same α across 0<z<9. Because the high-z SMF is sparse and the inferred SFRs at z≳s 5 are obtained by integrating this term over large Δz, a systematic error in α of order the literature range (Δα~0.1–0.2) would shift the high-z MS normalization by several tenths of a dex—comparable to the claimed improvement over Popesso+23. The low-z agreement with Leja+22 is less sensitive, but the z~2–7 spectroscopic concordance that underpins the strongest claim is directly exposed.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The paper infers the star-forming main sequence (MS) from the redshift evolution of compiled star-forming and quiescent stellar mass functions (SMFs) over 0.1 < z < 9, rather than from direct SFR indicators. After fitting continuous Schechter forms (with Eddington bias) to multi-survey SMF data, the authors solve a continuity equation that includes star formation, quenching, and merger-driven mass growth and number-density absorption (Eqs. 7–14), then fit the resulting SFR–mass grid with a double power-law plus high-mass cutoff (Eq. 19). The inferred MS agrees with JWST/NIRSpec Hα/UV measurements (Clarke+25) at z ∼ 2–7 and with Prospector-based SED results (Leja+22; Simmonds+25) at z ≲ 3, while lying below the Speagle+14 and Popesso+23 concordance compilations and closer to theoretical models driven by halo accretion. Forward evolution of the SMF with literature MS relations (App. A) is used as a consistency check, and implications for SFE and the cosmic SFRD are discussed.","tokens_in":35601,"tokens_out":1312,"duration_ms":10378,"significance":"If the result holds, it supplies an independent, SMF-driven MS that is self-consistent with the observed mass buildup and that can be used as a benchmark for both observations and models. The public Python implementation of the SMF and MS fits, the joint-posterior uncertainty propagation, and the explicit forward-evolution test in Appendix A are concrete strengths that make the result reusable and falsifiable. The claimed concordance with spectroscopic JWST SFRs and the caution against uncritical use of concordance MS compilations would affect a wide range of demographic and semi-empirical applications.","major_comments":[{"comment":"§3.2–3.3, Eqs. (12)–(15): The merger absorption term A uses Behroozi+13 α(m⋆) calibrated only for backward tracking of z ∼ 0 progenitors, yet is applied over 0 < z < 9. The text itself notes that forward/backward evolution is asymmetric and that α depends on descendant SFR (Clauwens+16; Wang+23). Because high-z SFRs are obtained by integrating this term over large Δz, a systematic Δα ∼ 0.1–0.2 would shift the z ≳ 5 MS normalization by several tenths of a dex—comparable to the claimed offset from Popesso+23. A quantitative sensitivity test (varying α within the literature range, or restricting the comparison to z ≲ 3 where the calibration is safer) is needed before the z ∼ 2–7 spectroscopic agreement can be treated as robust.","section":null},{"comment":"§4.1, Table 1 note c and Fig. 2: At z ≳ 6 the total SMF is treated as purely star-forming. While the quiescent fraction is expected to be small, the quiescent SMF is already poorly constrained above z ∼ 2.5, and any residual quiescent contribution would reduce the inferred star-forming growth rates. The paper should either (i) propagate a plausible high-z quiescent fraction as a systematic or (ii) show that the z ∼ 6–9 MS is insensitive to a few-percent quiescent contamination.","section":null},{"comment":"§5.4, Eq. (20) and Fig. 5: Integrating the SMF-derived MS over the same SMF yields an SFRD whose peak at cosmic noon is lower by a factor of ∼2–3 than the Madau & Fragos (2017) luminosity-density compilation. Because the MS is constructed to reproduce the SMF evolution, this tension is not independent of the input data; it re-expresses the longstanding stellar-mass-density vs. SFRD discrepancy. The discussion should clarify what new information the SFRD comparison adds and whether the discrepancy is driven by the MS shape, the SMF faint-end slope, or the UV-to-SFR conversion assumptions.","section":null}],"minor_comments":[{"comment":"§2 and Fig. 1: The phrase “star-forming star-forming main sequence” in the section title is a typographical duplication.","section":null},{"comment":"Eq. (6) and surrounding text: The return fraction R ≈ 0.36 is assumed constant and instantaneous; a short note on the sensitivity of the MS normalization to R (or to a time-dependent return fraction) would help readers compare with models that track integrated vs. surviving mass differently.","section":null},{"comment":"Table 1 / §4.1: Systematic offsets between surveys (stellar-mass pipelines, cosmic variance) are acknowledged but not folded into the likelihood. Even a simple extra variance term or a leave-one-survey-out test would strengthen the uncertainty budget.","section":null},{"comment":"Fig. 2 caption: The statement that the solid curves under-estimate the data at the massive end because they show the intrinsic SMF is correct but easy to miss; a parenthetical reminder in the main text would help.","section":null},{"comment":"Appendix C / Fig. 7: The seven-parameter piecewise redshift evolution of the MS parameters is flexible; reporting the reduced χ² or residual scatter of the fit to the discrete MS grid would reassure readers that the functional form is not over-fitting.","section":null}],"recommendation":"major_revision","confidential_remarks":"The central scientific claim is interesting and the methodology is carefully executed, but the high-z result is load-bearing on an extrapolation of Behroozi α that the authors themselves flag as asymmetric and SFR-dependent. I would not accept without a quantitative sensitivity analysis of that term (and of the pure-SF assumption at z ≳ 6). If those tests show the z ∼ 2–7 agreement is stable, the paper becomes a strong contribution; if not, the claim should be restricted to z ≲ 3–4 where the calibration is safer. Scope is appropriate for a major astrophysics journal."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"The useful takeaway is a new, publicly coded main-sequence fit (Eq. 19 + Table 3) obtained by continuous joint Schechter modeling of star-forming + quiescent SMFs to z=9 and solving the continuity equation with explicit merger and quenching terms. It sits closer to Clarke+25 Hα/UV points and Leja/Simmonds Prospector results than to the Speagle/Popesso concordance curves, and it tracks halo accretion rates more cleanly.\n\nWhat is actually new is not the continuity idea (Bell, Leja, Davidzon already did versions of it). It is the first post-JWST continuous SF+Q fit, the careful Eddington-bias convolution, the forward-evolution sanity check in Appendix A, and the released Python implementation. The math is clean: the Eulerian continuity equation is standard, the MCMC posteriors are propagated into the MS uncertainty band, and the double-power-law + exponential form is a transparent description of the grid they recover. Citations are dense and fair; they flag the known circularity and the external spectroscopic benchmarks that break it.\n\nThe soft spot the stress-test flags is real and load-bearing at high z. Behroozi α(m⋆) is calibrated only for backward tracking of z~0 progenitors and is known to be SFR-dependent and asymmetric. They apply it across 0<z<9. A Δα of 0.1–0.2 can move the z≳5 normalization by a few tenths of a dex—comparable to the claimed improvement over Popesso. High-z SMFs are also treated as pure star-forming, and R and ξ_min are fixed. That said, the z≲3 agreement with Leja and the overall shape are less sensitive to this term, so the central low-to-mid-z claim still stands. The SFRD peak being lower than Madau is an honest residual tension they do not paper over.\n\nThis is for people who actually use the MS as an input (semi-analytics, abundance matching, SFH priors). It deserves a serious referee; the code and tables make it easy to stress-test. I would cite the fit parameters and the caution against uncritical use of the old concordance relations. Engage with it.","headline":"Solid, public SMF-to-MS re-derivation that matches JWST spectroscopy better than Speagle/Popesso; the Behroozi α extrapolation is the real soft spot but does not sink the low-z result.","tokens_in":36263,"tokens_out":597,"would_cite":true,"duration_ms":7386,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"The star-forming main sequence required by the evolving stellar mass function matches JWST spectroscopy and theory, not older concordance relations.","keywords":["star-forming main sequence","stellar mass function","galaxy evolution","JWST","star formation rate density","galaxy-halo connection","quenching","mergers"],"falsifier":"If new high-redshift mass functions or independent merger-rate measurements require a substantially different progenitor number-density evolution than the adopted alpha, the inferred high-z main sequence will shift and lose agreement with JWST spectroscopy.","tokens_in":36103,"feed_emoji":"🌌","tokens_out":832,"duration_ms":7359,"temperature":0.7,"pith_summary":"The paper asks what star formation rates galaxies must have had if the observed buildup of stellar mass over cosmic time is to be self-consistent. Using a continuous fit to star-forming and quiescent stellar mass functions from z = 0.1 to 9, and tracing galaxy growth while accounting for mergers and quenching, the authors invert the continuity equation to recover the average star-forming main sequence from 10^8 to 10^11 solar masses. That inferred sequence agrees with independent JWST spectroscopic SFR measurements from z ~ 2-7 and with non-parametric SED-fitting results at lower redshift, and it tracks the redshift evolution expected from dark-matter halo accretion. It differs systematically from the widely used Speagle and Popesso concordance compilations, which over-predict the growth of the mass function when integrated forward. The result supplies a phenomenological check on measured star formation rates and cautions that applications of older compilations need to confront these systematic offsets.","feed_headline":"Mass-function growth forces a gentler star-forming main sequence","feed_subtitle":"It matches JWST spectroscopy and halo accretion, not the older concordance fits.","key_machinery":"A continuity equation for the star-forming stellar mass function that includes star-formation-driven mass growth, merger-driven absorption parameterized by Behroozi-style number-density evolution, and quenching measured from the growth of the quiescent mass function; solved by backward-tracking progenitor abundances.","core_discovery":"The star formation rates required by the redshift evolution of the stellar mass function, once mergers and quenching are included, produce a main sequence that matches JWST/NIRSpec and non-parametric SED results while lying below and evolving more gently than the Speagle+14 and Popesso+23 concordance relations.","pith_inferences":["The residual SFRD peak mismatch implies that either dust corrections in UV-based densities or the high-mass end of the mass function still carry unaccounted systematics.","Because the procedure is mass-complete wherever the SMF is measured, it can forecast the low-mass SFMS that future wide-field surveys will test once they reach 10^8 solar masses at high z.","If the Behroozi alpha itself evolves with redshift or with star-formation rate, the present inference would need a second iteration that couples the SFMS and the merger tracks self-consistently."],"forward_implications":["Analyses that adopt Speagle or Popesso main sequences as ground truth for galaxy ages or star-formation histories will systematically mis-estimate early and intermediate SFRs.","The gentler redshift evolution of the inferred sequence brings the long-term SFMS into closer alignment with mean dark-matter halo accretion rates.","The integral of the new sequence over the mass function under-predicts the peak of the cosmic star-formation-rate density relative to UV-luminosity compilations, reopening that classic tension.","Star-formation efficiency as a function of halo mass shows a low-mass slope that flattens from near unity at z greater than 7 toward the energy-regulated value of 2/3 by z = 0."],"fun_headline_variants":["Mass functions force gentler main sequence matching JWST not concordance","SMF growth yields main sequence below Speagle and Popesso fits","Inferred MS from mass evolution softens older concordance relations","Galaxy mass histories require slower SFR main sequence than Speagle","Self-consistent SMF paths produce shallower star-forming sequence"],"cache_read_input_tokens":32896,"weakest_assumption_plain":"The merger correction uses a number-density evolution calibrated only for the progenitors of present-day galaxies and is applied unchanged out to z = 9, even though forward and backward tracks are known to differ.","fun_headline_variants_meta":{"raw":{"variants":["Mass functions force gentler main sequence matching JWST not concordance","SMF growth yields main sequence below Speagle and Popesso fits","Inferred MS from mass evolution softens older concordance relations","Galaxy mass histories require slower SFR main sequence than Speagle","Self-consistent SMF paths produce shallower star-forming sequence"]},"model":"grok-4.5","effort":"low","cost_usd":0.00633,"raw_usage":{"total_tokens":1677,"prompt_tokens":837,"num_sources_used":0,"completion_tokens":69,"cost_in_usd_ticks":63300000,"prompt_tokens_details":{"text_tokens":837,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":771,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":837,"tokens_out":69,"duration_ms":8173,"temperature":1.0,"reasoning_tokens":771,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T15:02:57.105265+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"If new high-redshift mass functions or independent merger-rate measurements require a substantially different progenitor number-density evolution than the adopted alpha, the inferred high-z main sequence will shift and lose agreement with JWST spectroscopy.","supporting_citations":[],"review_version":1}