{"id":"6a878cb8-367e-44db-94e1-edae60a97f38","arxiv_id":"1908.07250","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"The peak of the stellar initial mass function is set by tidal forces and turbulent fluctuations that fragment the envelope around a first hydrostatic core, yielding a characteristic stellar mass of about 10 times that core mass.","lead":"This paper proposes that the characteristic mass of stars, around a tenth of the Sun's mass, is set by local physics: the minimum mass of the first hydrostatic core and the fragmentation of gas in the surrounding envelope. New simulations and an analytical model show that nearby fragments limit accretion onto the central star, pegging the mass peak at about 10 times the first core mass.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The model's peak at ~10 M_L depends on VP_ext/VP_ram fits (Eq. 4) taken from already-formed cores and applied to progenitor fluctuations; the authors flag this in §3.3, and setting these terms to zero underproduces fragments, so the central prediction is contingent on an untested extrapolation.","rationale":"The paper does serious work: the run1a/run1b experiment directly shows that preventing fragmentation within 140-280 AU shifts the sink mass distribution, and the comparison of mean system sizes in Fig. 3 is a nontrivial check of the model's fragment-count statistics. The analytical model also has only seven parameters with modest effect from nf. But the quantitative claim—peak at ~10 M_L—is not parameter-free: A, M, b, and the VP fits are measured from the same simulations used for validation, and the VP terms are explicitly extrapolated from developed cores to progenitors. The sensitivity to VP terms is demonstrated by the fact that VP=0 underproduces fragments. This does not make the paper wrong; it makes the central prediction contingent on a testable but untested assumption. The reader's conditional verdict is appropriate; a targeted remeasurement of surface terms for progenitor fluctuations would settle it.","tokens_in":23506,"tokens_out":5972,"duration_ms":65367,"concrete_test":"Recompute N(r) (Eq. 14) and the unshielded-mass histogram (Eq. 16) using VP_ext and VP_ram evaluated from the surface integrals in Eq. (3) on candidate progenitor fluctuations—local density maxima in the envelope around sinks younger than 10^3 yr, before they are classified as cores or sinks—rather than the developed-core fits in Eqs. (4). If the peak of dN_unsh/dlogM stays within 5-10 M_L and the predicted <n_sys>(r) still matches Fig. 3, the extrapolation is validated; if the peak shifts outside this range, the central quantitative claim lacks support.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing vulnerability is the extrapolation of the external-pressure fits in Eqs. (4). Section 2.7 derives VP_ext ≃ Eth(0.6-0.1M) and VP_ram ≃ Ekin(1.3-0.18M) from HOP cores identified at density threshold 10^8 cm^-3 (Fig. 6), i.e. from cores that are 'already formed and developed'. Section 3.3 then uses these same fits in the virial threshold Eq. (9) to count the pre-collapse density fluctuations that become shielding fragments, with the authors noting that the fluctuations 'are more core progenitors'. This matters because the predicted peak is set by where the mean fragment count N(r) (Eq. 14) reaches unity, and Section 3.3.2 reports that setting VP_ext=VP_ram=0 gives too few fragments. The surface terms are therefore load-bearing, not a small correction: if true boundary terms for progenitor fluctuations differ (e.g. because an infalling envelope does not provide a settled virial boundary), the radius where N(r)=1 and hence the 10 M_L peak would shift. The paper's own flagged limitation is exactly the untested link in the chain.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes that the characteristic (peak) mass of the stellar IMF is set locally by the competition between accretion onto a first hydrostatic core (FHSC) and fragmentation of its surrounding collapsing envelope. Two AMR simulations of 1000 M_sun clumps with different initial conditions produce sink mass spectra peaking at about 0.1–0.2 M_sun, i.e. 5–10 times the FHSC mass M_L. Suppressing sink formation within 140 or 280 AU of existing sinks (run1a/run1b) shifts the peak upward by about a factor of two, supporting the fragmentation-starvation mechanism. An analytical model counts turbulent density fluctuations in the envelope that are self-gravitating against tidal, thermal, and pressure terms; the resulting unshielded-mass distribution peaks near 10 M_L, in apparent agreement with the simulations.","tokens_in":23874,"tokens_out":9348,"duration_ms":82882,"significance":"If the result holds, the paper provides a concrete, falsifiable mechanism for the universality of the IMF peak: the peak is tied to the FHSC mass and local envelope physics rather than to the large-scale Jeans mass. The run1a/run1b numerical experiments are a clean controlled test and are the strongest evidence in the paper. The analytical model is transparent and yields a testable relation between the peak mass and M_L, although it is not parameter-free. The authors also explicitly state several assumptions and limitations, including the extrapolation of surface-pressure fits from developed cores to core progenitors. The comparison to run2 with different initial conditions adds credibility, but the model's inputs remain partly calibrated on run1.","major_comments":[{"comment":"The fits in Eq. (4), VPext ≃ Eth(0.6−0.1M) and VPram ≃ Ekin(1.3−0.18M), become negative for M>6 and M>7.2, respectively. Yet §3.3.1 adopts M=8 for run1, so the virial threshold in Eq. (9) is evaluated with negative surface-pressure terms. Since negative VPext and VPram change the sign of the surface contribution, they artificially alter the critical density contrast ηcrit and therefore the fragment count N(r) that controls the predicted peak at about 10 M_L. The authors should state the Mach-number range over which the fits are valid, avoid extrapolating beyond it, or justify the negative values physically.","section":"§2.7, Eq. (4); §3.3.1"},{"comment":"The surface-term fits used in Eq. (9) are derived from HOP cores identified at a density threshold of 10^8 cm^-3 (Fig. 6), i.e. from cores that are already formed and developed, but in §3.3 they are applied to pre-collapse density fluctuations that the authors themselves call \"more core progenitors\" (§3.3). This extrapolation is load-bearing: §3.3.2 reports that setting VPext=VPram=0 gives too few fragments, so the radius where the mean number of fragments N(r)=1, and hence the 10 M_L peak, depends on surface terms whose value for progenitors is not measured. A direct measurement of VPext and VPram for progenitor-type fluctuations, or a sensitivity study over their plausible range, is required to support the central prediction.","section":"§2.7, §3.3"},{"comment":"The model is calibrated on run1 (A=10, M=8, b=0.6) and then compared with the run1 peak, so the agreement between the predicted 10 M_L peak and the run1 spectrum is partly a consistency check rather than an independent prediction. The run2 comparison is a useful partial out-of-sample test, but it shares the same functional forms and several parameters. The authors should quantify the sensitivity of the predicted peak to A, b, ε, and nf over the ranges in Table 1, and ideally present the model prediction for an additional simulation before measuring the inputs.","section":"§3.3.2, Table 1"},{"comment":"The claim that the mass spectrum is \"expected to be relatively universal\" because the relevant processes are local is not established by the model itself: the inputs A, b, ε, and the VP fits are environmental parameters, and the paper only demonstrates insensitivity to two particular initial-condition choices. Unless the authors provide a physical argument or additional evidence that these parameters occupy a universal range, the universality conclusion should be softened to a speculation.","section":"Abstract; §4"}],"minor_comments":[{"comment":"There are several typographical errors: \"respectivelly\" should be \"respectively\", the affiliation line contains \"Universrté\", and the Fig. 6 caption says \"radio\" where \"ratio\" is meant.","section":"§2.2, §2.7"},{"comment":"The caption should specify which black curve corresponds to which Mach number and to the VPext=0 case; currently the labels \"M = 10, VPext = 0\" and \"M = 10\" are ambiguous when read with multiple curves.","section":"Fig. 3 caption"},{"comment":"The sentence preceding Eq. (16) is confusing: it first mentions mass located at radius smaller than the fragment position, then states that fragments accrete mass outside the sphere. Clarify that Eq. (16) follows from each fragment accreting a fraction 1/nf of the mass outside its radius.","section":"§3.2.5, Eq. (16)"},{"comment":"The description of the neighbour statistics says the maximum number of neighbours over timesteps is taken; please explain why the maximum rather than the mean or instantaneous value is used, as this choice could bias the inferred multiplicity.","section":"§2.4"}],"recommendation":"major_revision","confidential_remarks":"The numerical experiments are solid and the analytical proposal is interesting, but the central quantitative prediction currently relies on calibrated inputs and on an extrapolation of surface-pressure fits whose domain of validity is not established. The paper is likely publishable after the sensitivity and domain issues are addressed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper claims that the characteristic stellar mass is set by fragmentation-induced starvation, and it brings both a new numerical experiment and a new analytical model. The run1a/run1b experiments are the cleanest part: preventing sink formation within 140 and 280 AU shifts the mass spectrum to higher masses by about a factor of two, showing directly that close fragmentation limits accretion onto the central object. That result stands on its own and supports the proposed mechanism independently of the model.\n\nThe analytical model is also new. It counts self-gravitating turbulent density fluctuations above the first hydrostatic core mass, including tidal and surface-pressure terms, and reproduces the peak near 10 M_L. The model gives a reasonable match to the neighbor statistics in both run1 and run2, and the paper is honest about the main discrepancy below 100 AU.\n\nThe soft spots are real, though. The surface-pressure fits in Eq. (4) are measured from already-formed cores, then applied to progenitor fluctuations that are less settled. The authors acknowledge this in Section 3.3, and the model with those terms set to zero underproduces fragments, so the predicted peak location depends on an untested extrapolation. Also, most of the model inputs—A, Mach number, b, and the pressure fits—are measured from run1 and then used to reproduce run1's peak. The run2 comparison helps, but it is the same machinery with slightly adjusted parameters. So the factor of ~10 is more a consistent fit than a strong prediction.\n\nThat said, the central mechanism is not fragile. The numerical experiment is clean, and the broad picture—that the first-core mass sets a floor and fragmentation in the inner few hundred AU sets the peak—is plausible and supported by the simulations. I would not reject the paper over the circularity, but I would want the authors to test the surface-term extrapolation more directly, perhaps with a run that isolates the pressure terms or with a different equation of state.\n\nThis paper is for IMF theorists and star-formation simulators. It deserves a serious referee and, with revisions that sharpen the independent test, is worth publishing. Send it to review.","headline":"A credible local mechanism for the IMF peak, with a clean numerical test and a useful but partly calibrated analytical model.","tokens_in":24363,"tokens_out":2483,"would_cite":true,"duration_ms":24925,"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":"Typical stellar mass comes from sibling fragments that crowd out infalling gas.","keywords":["initial mass function","first hydrostatic core","gravo-turbulent fluctuations","star formation","tidal forces","protostellar accretion","virial theorem","opacity limit"],"falsifier":"Measure the thermal and ram pressure at the boundary of density fluctuations before they collapse, not on developed cores; if the resulting virial surface terms differ from the fits in Eq. (4) enough to move the radius where the mean fragment count reaches one away from $r\\simeq 8\\text{--}9\\,r_{e,L}$, the predicted $10M_L$ peak shifts and the claimed agreement with the simulations fails.","tokens_in":23317,"feed_emoji":"⭐","tokens_out":9764,"duration_ms":95453,"temperature":0.7,"pith_summary":"The paper proposes that the characteristic mass of stars, the peak of the initial mass function around 0.2--0.3 solar masses, is set by a local process: a young protostar keeps accreting until sibling fragments condense out of the turbulent density fluctuations in its collapsing envelope and steal the gas beyond them. The minimum seed is the first hydrostatic core, $M_L \\sim 0.02\\text{--}0.03\\,M_\\odot$, whose mass is fixed by the dust opacity limit. Counting self-gravitating fluctuations in the envelope, with tidal forces and pressure support setting the instability threshold, the model finds that the mass remaining unshielded for the central object peaks at roughly $10M_L$. The paper's simulations show the same peak at $5\\text{--}10M_L$, and experiments that forbid fragmentation within 140--280 AU shift the mass spectrum to higher masses. If correct, the peak is nearly universal because it depends only on local envelope physics on scales below a few hundred AU, not on the mean Jeans mass or the large-scale cloud.","feed_headline":"Stars get their typical mass from sibling fragments","feed_subtitle":"Local shielding, not the cloud-scale Jeans mass, fixes why stars cluster near 0.2-0.3 solar masses","key_machinery":"The load-bearing object is the first hydrostatic (Larson) core, with mass $M_L$ set by the dust opacity limit; the mechanism is the statistical counting of shielding fragments in the accreting envelope. The model takes the envelope density as $\\rho_{\\rm env}(r)=A c_s^2/(2\\pi G r^2)$ and describes turbulent fluctuations by a log-normal PDF whose width is set by the local Mach number. Each fluctuation is tested with a virial criterion that includes self-gravity, the tidal fields of the core and envelope, thermal and turbulent support, and external thermal and ram-pressure surface terms fitted to the simulations. Counting fluctuations of mass at least $M_L$ in concentric shells gives the mean number of fragments $N(r)$; the radius where this number reaches a few defines the shielded accretion radius, and the distribution of unshielded mass peaks near $10M_L$.","core_discovery":"The central claim is that the peak of the stellar mass function is not inherited from the large-scale gas reservoir but generated during the final accretion phase. When the effective equation of state stiffens once dust becomes opaque, a first hydrostatic core of mass $M_L$ forms; because its cooling time is long, infalling gas piles up and the core grows. In the surrounding $r^{-2}$ envelope, turbulent density fluctuations that are at least as massive as $M_L$ and dense enough to overcome the tidal shear of the central core plus thermal, turbulent and ram-pressure support collapse into new fragments. These fragments accrete the gas beyond their own radii, so the final mass of the central object is set by the unshielded mass inside the radius where a few fragments have formed. Counting fluctuations analytically with a log-normal density PDF and a virial threshold gives a peak near $10M_L$; the two simulations show peaks at $0.1\\text{--}0.2\\,M_\\odot$, i.e. $5\\text{--}10M_L$, and suppressing sink formation within a few hundred AU shifts the distribution to twice larger masses.","pith_inferences":["If the peak tracks $M_L$, resolved stellar mass functions in low-metallicity dwarf galaxies should show only a mild shift, because the paper's opacity scaling of $M_L$ is weak; this is a testable prediction beyond the paper's own metallicity discussion.","Substellar objects may be the by-product of the same process: fragments below $M_L$ that cannot trigger the second collapse cool only if they stay isolated, otherwise they grow or merge into a first core.","The same local shielding argument could apply at the high-mass end, where fragmentation-induced starvation has been invoked to limit accretion onto massive protostars, making the two regimes different scales of one competition.","A direct numerical check of the model's key extrapolation would be to compute the virial surface terms on pre-collapse density fluctuations rather than on already formed cores; if those terms differ, the predicted fragment count and peak would need revision."],"forward_implications":["The peak stellar mass should sit near $10M_L$, about $0.2\\text{--}0.3\\,M_\\odot$ for solar-metallicity conditions, rather than being tied to the parent cloud's Jeans mass.","The peak position should vary little across star-forming environments, because only the local envelope physics within a few hundred AU enters the calculation.","Changing the effective equation of state so that $M_L$ changes should move the peak proportionally, as the simulations in this paper and its companions show.","Preventing fragmentation within a few hundred AU of existing protostars should shift the mass spectrum to about twice larger masses, matching the run1a and run1b experiments.","The high-mass tail of the distribution is not explained by this shielding mechanism and is left to the initial mass-reservoir statistics."],"supporting_citations":[{"why":"Supplies the earlier simulations showing the IMF peak is insensitive to large-scale initial conditions and that the high-mass tail follows a power law.","marker":"paper I"},{"why":"Establishes that the peak scales with the mass of the first hydrostatic core and proposes shielding by nearby fragments as the cause.","marker":"paper II"},{"why":"Provides detailed first-hydrostatic-core simulations fixing $M_L \\simeq 0.03\\,M_\\odot$ and the core lifetime used in the argument.","marker":"Vaytet & Haugbølle 2017"},{"why":"Identifies the density where gas becomes adiabatic, the basis for the opacity-limit estimate of $M_L$.","marker":"Masunaga & Inutsuka 1999"},{"why":"Gives the singular collapsing-envelope profile that the model uses for the mean density around the accreting core.","marker":"Shu 1977"}],"fun_headline_variants":["Local fragments, not cloud-scale gas, set star masses","Star mass peak traced to sibling fragments at birth","Turbulent clumps near first core dictate final stellar mass","How nearby fragments cap a star's growth: new model","Stellar masses emerge from local physics, not Jeans scale"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that the pressure squeezing measured on already formed cores also acts on the less dense clumps that are about to become fragments; if the two situations differ, the predicted number of sibling fragments and the peak stellar mass would move.","fun_headline_variants_meta":{"raw":{"variants":["Local fragments, not cloud-scale gas, set star masses","Star mass peak traced to sibling fragments at birth","Turbulent clumps near first core dictate final stellar mass","How nearby fragments cap a star's growth: new model","Stellar masses emerge from local physics, not Jeans scale"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000243,"raw_usage":{"total_tokens":1577,"prompt_tokens":1043,"completion_tokens":534,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":659,"completion_tokens_details":{"reasoning_tokens":453}},"tokens_in":659,"tokens_out":534,"duration_ms":5018,"temperature":1.0,"reasoning_tokens":453,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:21:48.414148+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the thermal and ram pressure at the boundary of density fluctuations before they collapse, not on developed cores; if the resulting virial surface terms differ from the fits in Eq. (4) enough to move the radius where the mean fragment count reaches one away from $r\\simeq 8\\text{--}9\\,r_{e,L}$, the predicted $10M_L$ peak shifts and the claimed agreement with the simulations fails.","supporting_citations":[{"cited_title":"1999, , 510, 822, 10.1086/306608","cited_arxiv_id":null,"evidence_quote":"Identifies the density where gas becomes adiabatic, the basis for the opacity-limit estimate of $M_L$."}],"review_version":1}