{"id":"ef3040e3-61b0-44ca-8be7-b11324b0d0c2","arxiv_id":"1908.02608","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"A 0-d toy model combining dust growth, pebble-driven planetesimal formation, and planetesimal collisions predicts planetesimals dominate the solid mass within about a million years, after which collisions resupply small particles.","lead":"This paper uses a simple local computer model to track how dust, pebbles, and planetesimals in a protoplanetary disk evolve together over time. It finds that planetesimals can quickly become the dominant form of solid matter and that their collisions later replenish dust.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Default parameters violate the model's own locality condition: at R=10 AU, eps=0.01, the pebble conversion length exceeds the drift length by about a factor of 25, so radial transport cannot be neglected.","rationale":"The reader's weakest assumption correctly identified the neglect of radial transport as the most fragile premise. The present stress-test sharpens that concern into a quantitative failure of the model's own validity criterion: with default parameters, the pebble conversion length l = d/eps exceeds the drift length R by a factor of about 25 at R=10 AU, so pebbles would be advected inward before being converted into planetesimals under the model's own formation recipe. This is not an external objection to using a toy model; it is an internal check that the stated regime of validity is violated in the default parameter study. The concern bears directly on the strongest claim because the reported mapping from dust and disk mass to planetesimal mass, and the timescales on which planetesimals dominate, are computed for eps values where the 0-d mass conservation in Eq. (2) cannot hold. The authors are transparent about the transport limitation and Appendix A compares to Lenz et al. (2019), which includes radial transport and a particle size grid; that comparison gives some support for the early growth phase, but it does not validate the collision-dominated phase or the quantitative mapping. I also verified the reader's note that Section 4.2.3 contains a contradictory statement about the disk-mass dependence of the normalized planetesimal fraction, but the transport violation is the more load-bearing issue. Because the paper presents itself as an explicitly local toy model and discloses its main limitation, a conditional verdict remains appropriate; the condition should be sharpened to require tau_conv < tau_drift (or, equivalently, l < R) and the parameter study should be restricted or corrected accordingly.","tokens_in":25888,"tokens_out":13723,"duration_ms":154362,"concrete_test":"Implement a minimal 1-D extension of Eqs. (35a)-(35c) by adding the radial advection term -(1/R) d/dR (R v_drift Sigma_pbb) for pebbles (and, if desired, for dust), using v_drift from Eq. (18), on a grid R=1-100 AU with the default parameters of Table 1 (eps=0.01, Mdisk=0.02 Mstar, St=0.1). Compare the radially integrated planetesimal mass at t=1e5 yr and t=1e6 yr with the 0-D results in Fig. 6. If the planetesimal mass drops by a large factor or the inward pebble flux at the inner boundary becomes large, then the local approximation is invalid in the default regime and the advertised mapping cannot be applied without including radial transport. As a simpler preliminary check, the ratio tau_conv/tau_drift = l/R should be evaluated for the parameter grid; if it exceeds unity, the stated validity condition is violated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that planetesimals form fast and everywhere and that the dust and total disk mass can be mapped to hidden planetesimal mass for a given epsilon. Both depend on the local 0-d premise that radial transport is unimportant, which the authors state is valid only if pebble conversion is faster than inward drift. Using the paper's own equations (13), (18), (20), (21) and Table 1, I estimate at R=10 AU, St=0.1, eps=0.01: the drift velocity from Eq. (18) is about 640 cm/s, so the drift timescale R/|v_drift| is roughly 7.5e3 yr. The conversion timescale is l/|v_drift| with l=d/eps=250 AU, giving roughly 1.9e5 yr. Thus the conversion length is about 25 times the drift length, and the conversion timescale is about 25 times longer than the drift timescale. A pebble would reach the star before traversing the roughly 100 trap separations required for conversion. Even for eps=0.1, the conversion timescale (about 1.9e4 yr) still exceeds the drift timescale (about 7.5e3 yr) at 10 AU; only for eps on the order of 0.25 or larger does the stated condition become self-consistent. The parameter study in Figs. 4-7 uses eps=0.001-0.5 and R=10 AU, so most of the reported regime violates the model's own locality condition. The absence of the radial advection term in Eq. (2) would therefore remove the pebble supply on timescales shorter than the conversion timescale, changing the predicted planetesimal fractions, the onset times, and the dust-to-planetesimal mapping. The authors disclose the transport limitation in Section 5, but the disclosure does not quantify that the default parameter set falls outside the valid regime.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper formulates a 0-dimensional, local model of solid evolution in a protoplanetary disk, coupling dust growth to pebbles (Birnstiel et al. 2012), pebble-flux-regulated planetesimal formation (Lenz et al. 2019), and destructive planetesimal collisions. The model is written as a small set of coupled ODEs for the column densities of dust, pebbles, and planetesimals, with two scenarios for whether collisional dust can regrow. The authors integrate the model locally for a grid of radii and combine the solutions to estimate global mass fractions. They report that planetesimals form quickly and everywhere, dominate the solid mass after roughly 10^4–10^6 yr depending on parameters, and that collisions later resupply dust and pebbles. They propose that, for a given planetesimal-formation efficiency, the observed dust content and disk age can be mapped to the hidden planetesimal mass.","tokens_in":26247,"tokens_out":10549,"duration_ms":110332,"significance":"If the central claims survive scrutiny, the paper offers a simple heuristic tool for estimating unobservable planetesimal mass from observable dust mass and age, and a clear framework for testing how trapping efficiency and disk mass affect the solid inventory. Its strengths are the transparent derivation of the rate equations from cited models, the explicit two-scenario treatment, a parameter study that spans relevant disk conditions, and an honest limitations section. The main quantitative conclusions are, however, contingent on the locality assumption and on the unconstrained efficiency parameter, and the default parameter set lies outside the regime where the locality assumption is valid; the paper is therefore best viewed as a conceptual toy model at present rather than a calibrated predictor.","major_comments":[{"comment":"The model's own validity condition is that radial transport can be neglected only when pebble conversion into planetesimals is faster than inward drift. Using the paper's expressions, with the default parameters R=10 AU, St_pbb=0.1, eps=0.01, and d=5 h_g, the conversion length l=d/eps is approximately 250 AU, whereas a pebble drifting at the speed in Eq. (18) traverses only about 10 AU in the same time; equivalently, the conversion timescale is about 1.9e5 yr and the drift timescale is about 7.5e3 yr. The conversion length thus exceeds the drift length by roughly a factor of 25, so the no-transport assumption behind Eq. (2) fails for the default parameters and for most of the parameter space explored in Figs. 4–7 (eps=0.001–0.5 at R=10 AU). Because the reported onset times, maximal planetesimal fractions, and the dust-to-planetesimal mapping in Section 6 all depend on the local pebble supply, this is a load-bearing inconsistency, not merely a limitation; the authors should either restrict the study to the self-consistent regime (approximately eps ≥ 0.25 at 10 AU with the default disk parameters) or include the radial advection term and revisit the results.","section":"Introduction, §2.3, §5 (Eqs. (2), (18), (20), (21); Table 1; Figs. 4–7)"},{"comment":"The sentence \"Less massive disks have a lower relative planetesimal fraction than more massive disks\" is the opposite of the abstract's finding that planetesimal collisions are more significant in more massive disks and lead to lower relative planetesimal fractions compared to less massive disks; it also contradicts the bullet summary in Section 6 and the trend visible in Fig. 5, top-left panel. This appears to be a sign error in the wording, but because it is one of the paper's headline parameter-study results, it must be corrected and checked against the figure.","section":"§4.2.3"},{"comment":"The only direct benchmark, Appendix A, shows that the two-population growth model makes planetesimal formation set in earlier than in Lenz et al. (2019), which includes radial transport and a resolved size grid. Since the paper's central timescales (e.g., planetesimal dominance by 10^4–10^6 yr) feed directly into the proposed age-based mapping in Section 6, the authors should state explicitly whether the predicted onset times are to be read as lower limits and should quantify the offset from a transport-including model.","section":"Appendix A"}],"minor_comments":[{"comment":"The default parameter set is inconsistently reported: Table 1 lists eps=0.01 and Mdisk=0.02 M_sun, while Fig. 2 and Section 4.1 use eps=0.1 and Mdisk=0.01, and Fig. 4 also says Mdisk=0.01. The default set should be defined once and used consistently in captions and text.","section":"Table 1, Fig. 2, Fig. 4"},{"comment":"Equation (C.5) appears to have a bracket or formatting error in the printed polynomial; check the typesetting.","section":"Eq. (C.5)"},{"comment":"The limitations paragraph correctly identifies the locality assumption as a strong limitation, but the abstract and conclusions do not hedge the headline claims accordingly; consider adding a sentence that the quantitative mapping is only valid where the local approximation holds.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest about its limitations, but the headline claims in the abstract and conclusion should be scaled back until the parameter study is restricted to the regime where the local approximation is self-consistent. The authors may be able to address the main concern with a moderate revision, either by adding a radial advection term or by re-framing the study as valid only for high trap efficiencies."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a transparent 0-d toy model combining Birnstiel et al. (2012) two-population dust growth with Lenz et al. (2019) pebble-flux-regulated planetesimal formation and adding planetesimal collisions. That combination is new, and the paper earns credit for showing that collisions resupply dust and pebbles and push the system through three phases: dust growth, planetesimal dominance, then collisional fragmentation. The differential equations are clearly derived, the parameter study is honest, and the limitations section is unusually upfront.\n\nThe soft spots, in order of severity.\n\nFirst, the model's own locality condition is violated by the default parameters. At R=10 AU with eps=0.01 and St=0.1, the conversion length l=d/eps is tens of AU, while the drift timescale R/|v_drift| is about 10^4 yr and the conversion timescale is several times longer. So a pebble would reach the star before traversing enough traps. The paper discloses that radial transport is neglected and says this is fine if the latter process is sufficiently efficient, but it does not check that the default values satisfy that condition. Most of the parameter study sits in the violating regime, so the quantitative outputs—onset times, mass fractions, the dust-to-planetesimal mapping—are not self-consistent. This is not a fatal blow for a toy model, but it means the advertised mapping is conditional on an epsilon that is both unconstrained and, at the default value, incompatible with the local approximation.\n\nSecond, Section 4.2.3 has a sentence that says less massive disks have a lower relative planetesimal fraction than more massive disks. The figure and the abstract say the opposite, and the reasoning given (collisions more important in massive disks) supports the opposite. Looks like a typo, but it is confusing.\n\nThird, the code is not provided, so the figures cannot be exactly reproduced. For a 0-d model with few ODEs, that would be easy to fix.\n\nThe central claim—planetesimals form fast and come to dominate—follows from the assumed rate equations, so internally the paper is sound. But the practical predictive power is weaker than the title suggests, because epsilon is a free parameter and the local approximation fails for most of the explored range.\n\nWho this is for: people building initial conditions for late-stage disk models or interpreting dust observations in terms of hidden mass. It deserves a serious referee. I would accept it with mandatory revision: quantify the local-approximation regime, move the parameter study into it or add a radial advection term, and fix the 4.2.3 sentence.\n\nWorth bringing to reading group as an example of how to present a toy model honestly.","headline":"Transparent local toy model of dust, pebbles, and planetesimals, but the default parameters sit in a regime where the neglected radial drift is faster than the assumed conversion, so the quantitative mapping needs a major caveat.","tokens_in":26838,"tokens_out":3993,"would_cite":true,"duration_ms":40620,"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":"A local model of a protoplanetary disk shows that planetesimals can form fast and everywhere, locking more than 98% of the solid mass into bodies too faint to observe within about a million years.","keywords":["protoplanetary disks","planetesimal formation","dust growth","pebble flux","collisional fragmentation","circumstellar matter","solid mass distribution","disk age"],"falsifier":"Measure the millimeter-sized dust mass of a sample of protoplanetary disks with well-determined ages in the range $10^4$–$10^6$ yr and compare it with the model's predicted small-particle fraction: the model predicts that by roughly $10^5$–$10^6$ yr the dust-and-pebble reservoir falls to a few percent of the initial solid mass (before collisions replenish it at later times), so disks of that age that still hold dust close to their initial dust reservoir would contradict fast, universal planetesimal formation.","tokens_in":25627,"feed_emoji":"🪐","tokens_out":19012,"duration_ms":168067,"temperature":0.7,"pith_summary":"Dust in a protoplanetary disk may not stay dust. This paper argues that when pebble trapping and gravitational collapse are efficient, a local patch of the disk converts most of its solid mass into km-sized planetesimals within $10^4$–$10^6$ years, with planetesimals reaching more than 98% of the local solid column density (mass per unit disk area). The authors build a deliberately simple zero-dimensional model at a single radius, combining dust growth with pebble-flux-regulated planetesimal formation and collisional fragmentation, and they vary the distance to the star, the disk mass, and the formation efficiency. The payoff is a mapping: for a given formation efficiency, the dust and pebble mass observed at a given disk age, together with an independent estimate of the total disk mass, translates into the hidden mass stored in planetesimals. That mapping matters because telescopes measure only micrometer- and millimeter-sized grains, while planet formation theory needs the large bodies.","feed_headline":"98% of disk solids hide in planetesimals within a million years","feed_subtitle":"A simple model connects the dust and pebbles telescopes can see to the planetesimal mass hidden below detection.","key_machinery":"The load-bearing object is a closed set of coupled rate equations for the column densities of three solid species — dust, pebbles, and planetesimals — evolved at one radius while holding the total solid column density fixed, since radial transport is neglected. Three rates carry the argument: dust grows into pebbles on an exponential growth timescale from the two-population dust-growth model, with growth rate $\\dot\\Sigma_{\\rm growth}=\\Sigma_{\\rm dst}/\\tau_{\\rm growth}\\propto\\Sigma_{\\rm dst}^2$; pebbles are converted into planetesimals at the rate $\\dot\\Sigma_{\\rm form}=|v_{\\rm drift}|\\,\\Sigma_{\\rm pbb}/l$, where the conversion length $l=d/\\epsilon$ is the trap spacing divided by the trap efficiency; and planetesimals fragment collisionally at a rate $\\dot\\Sigma_{\\rm col}\\propto\\Sigma_{\\rm pls}^2$, with fragments distributed among the three species by a collisional-cascade power law of slope $\\xi=1.83$. In the scenario used for the parameter study, collisional dust is treated as too compact to grow again, so the system never reaches equilibrium: primordial dust drains monotonically, planetesimals dominate for a phase, and collisional debris accumulates at late times. The identity that enables the mapping is strict local mass conservation, $\\partial(\\Sigma_{\\rm dst}+\\Sigma_{\\rm pbb}+\\Sigma_{\\rm pls})/\\partial t=0$, which ties any observed decrease in small-particle column density directly to the hidden planetesimal column density.","core_discovery":"The paper's central claim is that planetesimals form quickly and at every radius in a protoplanetary disk: in the fiducial run the planetesimal column density exceeds 98% of the total solid column density starting around $5\\times10^4$ yr, and across the explored parameters a planetesimal-dominated phase begins between roughly $10^4$ and $10^6$ yr. Planetesimal collisions take over after about $10^6$ yr and resupply dust and pebbles, so a late disk can look dust-rich even though most of its mass was once locked in large bodies. The authors therefore propose a conditional relation between observables and hidden mass: given the formation efficiency $\\epsilon$ and an independent estimate of the total disk mass, the observed dust-and-pebble fraction at a known disk age determines how much mass sits in planetesimals, with more observed dust implying relatively less mass in large bodies. Quantitatively, the timing of the planetesimal-dominated phase and its peak mass depend strongly on the distance to the star $R$, the initial disk mass, and $\\epsilon$, and more massive disks end up with lower relative planetesimal fractions because their stronger collision activity recycles solids back into small particles.","pith_inferences":["Editorial inference: the 'over 98% within about $10^5$ yr' result is a model prediction tied to the chosen trap parameters, not an observed fact; comparing millimeter dust masses of roughly 1 Myr-old disks against the predicted small-particle fraction would directly bracket the effective trap efficiency $\\epsilon$.","Editorial inference: the model exposes a practical degeneracy — pairs of formation efficiency and disk mass can give nearly identical dust evolution, so real observations constrain a curve of hidden planetesimal mass rather than a single value unless total disk mass is measured independently.","Editorial inference: the omitted process with the largest lever is inward pebble drift; a transport-inclusive version of the model should delay planetesimal growth in the outer disk, which would turn the inferred planetesimal masses into upper limits for disks near 1 Myr.","Editorial inference: the two scenarios diverge only after about $10^6$ yr, so multi-wavelength observations separating small grains from larger pebbles in older disks could test whether collisional fragments regrow or stay inert."],"forward_implications":["Because planetesimals can absorb more than 98% of the solid column density within $10^5$–$10^6$ yr, millimeter-continuum dust mass is not a proxy for the total solid mass of a disk.","After roughly $10^6$ yr, planetesimal collisions resupply dust and pebbles, so an older disk can appear dust-rich even though its solids were earlier locked in planetesimals; dust-based mass estimates must include the collision channel.","For a fixed formation efficiency, the observed dust-and-pebble fraction at a known disk age, plus an independent total disk mass, determines the hidden planetesimal mass — more observed dust implies relatively less mass in planetesimals.","Planetesimal formation proceeds inside-out: the inner disk depletes its pebbles and forms planetesimals first, and the peak planetesimal column density migrates outward with time.","More massive disks have lower relative planetesimal fractions because their higher collision rates recycle solids back into small fragments sooner; the lowest-mass disks in the study convert solids into planetesimals most efficiently."],"supporting_citations":[{"why":"Supplies the two-population dust-growth model with its growth timescale and sets the pebble Stokes number of 0.1 used throughout.","marker":"Birnstiel et al. (2012)"},{"why":"Supplies the pebble-flux-regulated planetesimal formation recipe (conversion length and formation rate) that this paper couples to growth and collisions.","marker":"Lenz et al. (2019)"},{"why":"Supplies the collisional-cascade fragment size distribution with slope 1.83 that partitions collision fragments among dust, pebbles, and planetesimals.","marker":"Dohnanyi (1969)"},{"why":"Provides the 50-km planetesimal radius and the choice of Hill velocity as the collision relative velocity.","marker":"Morbidelli et al. (2009)"},{"why":"Provides the radial drift velocity expression that sets the pebble flux into traps and hence the planetesimal formation rate.","marker":"Nakagawa et al. (1986)"},{"why":"Provides the numerical trap distance of 5 scale heights and trap lifetime of 100 orbits used to set the conversion length.","marker":"Dittrich et al. (2013)"},{"why":"Supplies the mean-free-path collision timescale used to compute the planetesimal collision rate.","marker":"Birnstiel et al. (2016)"},{"why":"Supplies the irradiated temperature profile that sets disk scale height and growth timescales at each radius.","marker":"Chiang & Goldreich (1997)"},{"why":"Supplies the Minimum Mass Solar Nebula disk mass baseline against which the parameter study's disk masses are calibrated.","marker":"Weidenschilling (1977)"}],"fun_headline_variants":["Most disk solids hide in planetesimals within a million years","Planetesimal formation hides most disk solids within 1 Myr","Dust and pebbles reveal hidden planetesimal mass","Model ties observable dust to invisible planetesimal mass","Most solid mass in disks becomes planetesimals quickly"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Pebbles must be converted into planetesimals faster than they drift inward toward the star; if conversion is too slow, solid material cannot be treated as staying at one radius, and the predicted timescales and mass fractions change.","fun_headline_variants_meta":{"raw":{"variants":["Most disk solids hide in planetesimals within a million years","Planetesimal formation hides most disk solids within 1 Myr","Dust and pebbles reveal hidden planetesimal mass","Model ties observable dust to invisible planetesimal mass","Most solid mass in disks becomes planetesimals quickly"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00058,"raw_usage":{"total_tokens":2817,"prompt_tokens":1116,"completion_tokens":1701,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":732,"completion_tokens_details":{"reasoning_tokens":1620}},"tokens_in":732,"tokens_out":1701,"duration_ms":14198,"temperature":1.0,"reasoning_tokens":1620,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:39:39.738942+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the millimeter-sized dust mass of a sample of protoplanetary disks with well-determined ages in the range $10^4$–$10^6$ yr and compare it with the model's predicted small-particle fraction: the model predicts that by roughly $10^5$–$10^6$ yr the dust-and-pebble reservoir falls to a few percent of the initial solid mass (before collisions replenish it at later times), so disks of that age that still hold dust close to their initial dust reservoir would contradict fast, universal planetesimal formation.","supporting_citations":[{"cited_title":"F., Nesvorný, D., & Levison, H","cited_arxiv_id":null,"evidence_quote":"Provides the 50-km planetesimal radius and the choice of Hill velocity as the collision relative velocity."},{"cited_title":"2013, , 763, 117","cited_arxiv_id":null,"evidence_quote":"Provides the numerical trap distance of 5 scale heights and trap lifetime of 100 orbits used to set the conversion length."}],"review_version":1}