{"id":"6e407259-e554-48f6-b60e-b272ee10f45c","arxiv_id":"1908.02513","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Boulders between 1 and 100 m can transport water ice about 10% inward of the classical snowline before fully sublimating, depending on disk heating and dust-mantle assumptions.","lead":"Meter-to-hundred-meter icy boulders drifting inward through a protoplanetary disk can cross the water snowline and carry water ice to a region about 10% closer to the star than the classical snowline before they evaporate. The paper gives a quantitative estimate of this snowline smearing and shows a simple surface-sublimation formula matches a full thermal model for clean, mantle-free bodies.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing assumption is that boulders survive collisional erosion while drifting; the paper's own Sect. 4.1 estimate gives erosion timescales of about 10 yr, far shorter than the modeled 10^2-10^3 yr sublimation phase.","rationale":"The paper's stated goal is to quantify how far inward boulder-sized bodies can carry water ice before complete disintegration. That quantification has two prerequisites: boulders exist in disks, and they survive the 10^2-10^3 yr drift phase long enough to cross the snowline. Existence is handled as an explicit postulate, but survival is built into the no-collision ablation model. The reader identifies this as the weakest assumption, and I agree. The manuscript itself performs the relevant collision calculation in Sect. 4.1 and finds that the collision-free assumption fails on roughly 10 yr timescales under the Windmark et al. (2012) erosion fit. This is not an external challenge to the paper; it is the paper's own limiting caveat, and it sits at the load-bearing point of the central claim. A secondary issue is the overstated disk-independence of the '10%' result: the non-irradiated disk case gives only 3%, so the abstract's blanket wording needs qualification. But that is a reporting defect, not a defeat of the mechanism. Collisional erosion, by contrast, can remove the boulders before they deposit water inward of the snowline. The analytic ablation formula and the favorable comparison with the full cometary nucleus model are real supporting evidence for the sublimation step, and the parameter study is useful. No public code or data is provided, so independent reproduction would require reimplementation. The appropriate disposition is unchanged from the reader's conditional recommendation: the mechanism is plausible and the model is competent, but the headline result should be conditioned on the collision-erosion caveat that the authors themselves identify.","tokens_in":28269,"tokens_out":4937,"duration_ms":56217,"concrete_test":"Implement the paper's nominal 10 m case (Tables 1 and 2) with the cometary nucleus model, adding the Windmark et al. (2012) erosion rate quoted in Sect. 4.1 to the time-dependent radius evolution and radial drift, using the same assumed impactor size distribution as in Appendices B and C. Record the location where the radius first falls to 10 cm and the cumulative water mass released inside the snowline. Compare with the no-erosion run: if the inward boundary moves from roughly 0.9 r_snowline to within a few percent of the snowline, or if the inside-snowline water yield drops by more than a factor of two, the collision-free premise is indeed load-bearing. Also rerun with the lower Krijt et al. (2015) erosion rate to bracket the uncertainty.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that 1-100 m boulders cross the snowline and deposit water ice down to about 0.9 of the classical snowline radius depends on the drifting body losing mass only by sublimation. The paper's Section 4.1 shows that this premise is not secure in the collisional environment it assumes. Using the Windmark et al. (2012) erosion fit, the nominal 10 m target encounters roughly 4e-4 of its own mass per year in small impactors and loses about 8e-2% of its mass per year. The authors state that this makes the collision-free assumption valid only on timescales of about 10 yr, while the modeled drift and sublimation phase lasts 10^2-10^3 yr. Over that interval, collisional erosion can remove a substantial fraction of the boulder's mass before it reaches the snowline crossing region, so the water may be released outside the snowline rather than transported inward. The paper is transparent about this limitation, noting that the results are 'only strictly valid if either the surface density of solids is reduced ... or erosion is less efficient in the relevant mass regime.' The concern is not a peripheral robustness issue; it attacks the core premise of the mechanism. The analytic sublimation formula and the favorable comparison with the cometary nucleus model are genuine supporting evidence for the ablation step, but they do not address whether the boulder survives long enough for that ablation to matter.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript couples an α-disk evolution model to the radial drift of spherical, dusty-water-ice bodies and to two treatments of sublimation—an analytic Hertz-Knudsen-Langmuir surface-ablation formula and the full Marboeuf et al. (2012) cometary nucleus model—to ask how far inward of the classical H2O snowline drifting boulders can carry water ice. The authors find that bodies with radii between roughly 1 m and 100 m drift inward faster than the snowline moves and, in the absence of a dust mantle, shrink to 10 cm at distances about 10% inside the snowline, with the ratio remaining stable as the disk evolves. Disk mass, lifetime, and photo-evaporation affect this inward penetration at the percent level, but the non-irradiated disk case gives only about 3% inward penetration. The analytic sublimation formula is shown to reproduce the full thermal model for homogeneous bodies up to about 100 m in radius, and the authors separately assess dust-mantle, water-vapor, frictional-heating, and collision effects, identifying collisions with small impactors as the main limitation on the mechanism.","tokens_in":28552,"tokens_out":6946,"duration_ms":78221,"significance":"If the central result is accepted, the paper provides a concrete, quantitative mechanism for smearing the water snowline and for delivering water to terrestrial-planet-forming regions, with implications for disk chemistry, planetesimal formation, and exoplanet compositions. A genuine strength is that the 10% inward penetration is not a fitted parameter: it emerges from the coupled drift-sublimation simulations and is stable across several disk parameters. The validation of the analytic surface-sublimation expression against the full cometary nucleus model is also nontrivial, because the latter includes heat conduction and vapor diffusion that the former omits. The authors are transparent about their assumptions, including the collision caveat. The significance is nevertheless contingent on the survival of boulders against collisional erosion over the 10^2–10^3 yr drift phase: the paper's own Section 4.1 estimates a nominal erosion timescale of about 10 yr, so the headline claim currently applies to a restricted regime unless erosion is explicitly included in the trajectory calculation.","major_comments":[{"comment":"The collision analysis in Section 4.1 is load-bearing for the central claim. The authors state that the nominal 10 m target encounters about 4e-4 of its own mass per year in small impactors and, using the Windmark et al. (2012) erosion fit, loses about 8e-2% of its mass per year, making the collision-free assumption valid only on timescales of about 10 yr. This is one to two orders of magnitude shorter than the 10^2–10^3 yr sublimation/drift phase that produces the 10% inward penetration. Since erosion directly reduces the body's radius and mass, it can shorten the distance actually traveled before disintegration; this is not a peripheral parameter but a condition on the validity of the headline result. I recommend that the authors either add an erosion term to the radius-evolution equation (e.g., Eq. 14) and recompute the disintegration locations, or explicitly present the 10% result as conditional on the low-solid-surface-density or reduced-erosion regimes named in Section 4.1. The abstract and Conclusions currently state the 10% result without carrying this caveat forward.","section":"Sec. 4.1"},{"comment":"The abstract and Conclusion item 1 state that the 10% inward penetration holds 'almost independently' of disk properties, but the non-irradiated disk in Fig. 10 reaches only about 3% inside the snowline, as acknowledged in Section 3.3.2. A factor-of-three difference between 10% and 3% is not a 'percent-level' effect relative to the headline value. The claim should be qualified to irradiated disks with the nominal pressure-gradient profile, or the reported range of inward penetration should be given explicitly (e.g., 3–10% depending on irradiation and pressure gradient).","section":"Sec. 3.3.2 / Fig. 10"},{"comment":"The '10%' inner boundary is defined by the adopted 'complete disintegration' radius of 10 cm, after which the body is assumed to have a very short lifetime. Since the total inward drift is integrated from the initial radius down to this cutoff, the reported inward penetration depends on the choice of R_stop = 10 cm. A convergence check with a smaller cutoff radius (e.g., 1 cm) or a stated estimate of the residual drift time below 10 cm would make the headline distance more robust. I expect the effect to be small because sub-meter bodies drift more slowly, but the sensitivity is not currently quantified.","section":"Sec. 3.3.1 / Fig. 8"}],"minor_comments":[{"comment":"The heading contains a typo: 'Comparision' should be 'Comparison'.","section":"Sec. 3.1 heading"},{"comment":"The caption for Fig. 7(a) should explain that the log-ratio is artificially set to 12 once the snowline starts moving outward; currently this is stated only in the main text.","section":"Fig. 7(a) caption"},{"comment":"The sentence describing the no-vapor case as an 'upper boundary for the dynamical snowline location' is easy to misread. Since the presence of vapor slows sublimation and lets the body penetrate further inward, the no-vapor result is an outer (more conservative) limit on the dynamical snowline; please rephrase for clarity.","section":"Sec. 4.4"},{"comment":"The initial placement of bodies at 10% outside the snowline is justified as a relaxation condition, but the manuscript does not explicitly state that the pre-sublimation drift from farther out has no effect on the resulting disintegration position because sublimation is inactive below the 150 K threshold. A sentence stating this explicitly would remove a possible concern about the initial condition.","section":"Sec. 2.4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is transparent about the collision caveat, which is a credit. However, the abstract and conclusions currently present the 10% result without the condition that the collision analysis places on it, and the non-irradiated disk case gives a significantly smaller value. The manuscript is not beyond repair: the mechanism may well be valid in reduced-solid-density or low-erosion regimes, but the headline claim should be reframed or the erosion term should be included in the model. I would not recommend rejection at this stage."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the short version: this is a competent, readable parameter study of a genuinely neglected size regime, and the headline result—that 1–100 m boulders crossing the water snowline by radial drift can carry water ice about 10% closer to the star than the classical snowline—is plausibly real. But the paper's own collision analysis in Sect. 4.1 gives erosion timescales around 10 yr, far shorter than the 100–1000 yr sublimation/drift phase. That means the central mechanism only operates if the solid surface density is low or the erosive effect is weaker than the Windmark et al. fit. The authors are transparent about this, but the abstract and conclusion oversell the 10% as almost disk-independent, when the non-irradiated disk case gives only 3%.\n\nWhat is actually new: previous pebble-drift models stopped at around decimeter sizes, and Schorghofer's asteroid study was static. This is the first time the 1–100 m drift+ablation regime gets a quantitative treatment. The analytic Hertz-Knudsen-Langmuir surface sublimation formula is validated against the full cometary nucleus model with internal heat conduction and vapor diffusion, and the agreement in the no-mantle, pre-heated limit is solid. The parameter study covers disk mass, lifetime, irradiation, and mantle thickness; the disk-variation results are mostly robust at the percent level, and the snowline-normalized scaling is a nice way to decouple the result from disk evolution. There are no fitted parameters in the central claim; it comes out of the time-dependent simulation. The collision section deserves credit for showing its work and stating the survival problem explicitly.\n\nThe soft spots are real but not fatal to the paper as a parameter study. First, the collision timescale problem is load-bearing: if boulders don't survive to drift across the snowline, the 10% water delivery doesn't occur. The paper can't fix that without a full treatment of the solid size distribution and collisional evolution. Second, the abstract's \"almost independently\" overstates the disk dependence. Third, no code or data is released, which makes re-use harder. None of these make the modeling careless—the analysis is internally consistent within its stated assumptions.\n\nWho should read it: anyone modeling disk composition, pebble/planetesimal formation, or water delivery to terrestrial planet formation regions. It deserves serious peer review, with a request to correct the summary generalization and to state the collision-survival condition front and center rather than in a caveat.","headline":"A solid, honest parameter study of a neglected size regime, but the headline 10% snowline smearing is conditional on boulders surviving erosion, and the abstract oversells the disk-independence.","tokens_in":29120,"tokens_out":3246,"would_cite":true,"duration_ms":35099,"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":"Meter-sized boulders can carry water ice about ten percent closer to the star than the classical snowline.","keywords":["protoplanetary disks","water snowline","radial drift","boulder-sized bodies","sublimation","dust mantle","Hertz-Knudsen-Langmuir","planet formation"],"falsifier":"The claim would be settled by measuring the impact-erosion lifetime of meter-sized icy bodies in the relevant velocity regime: if such bodies lose more than a few percent of their mass per year, they would disintegrate before crossing the snowline and the predicted ten-percent inward water-ice zone would not form; if erosion is much slower, the zone is expected. A complementary observation is a spatially resolved water-vapor or ice map just inside the snowline of a boulder-rich disk, looking for the predicted inward extension rather than a sharp cutoff.","tokens_in":28058,"feed_emoji":"❄️","tokens_out":11683,"duration_ms":117878,"temperature":0.7,"pith_summary":"This paper tries to establish that boulder-sized solids, roughly 1 m to 100 m across, are not simply lost from a protoplanetary disk: they drift inward faster than the water snowline recedes, cross it, and shed water ice as they ablate, so that ice can be present up to about ten percent closer to the star than the classical snowline location. A reader should care because disk composition, planetesimal formation, and the water budgets of forming planets are usually computed with a sharp snowline set by temperature and pressure alone; if this result holds, that boundary is smeared inward by a dynamical process that depends on the presence and size distribution of boulders. The paper also shows that for an unmantled body made of dust and water ice, a simple analytic surface-sublimation formula matches a full cometary-nucleus thermal model, because bodies below about 100 m reach internal thermal equilibrium in about a year.","feed_headline":"Boulders carry water ice 10 percent past the snowline","feed_subtitle":"If correct, planet-forming disks hold ice-bearing solids well inside the classical snowline, changing where water is available.","key_machinery":"The argument runs on a race between two speeds: the radial drift speed of a solid body, computed from gas-drag stopping times in the Epstein, Stokes, and quadratic regimes, and the inward speed of the water snowline during disk evolution. On the ablation side, the Hertz-Knudsen-Langmuir formula $\\phi(T)=P_s(T)/\\sqrt{2\\pi m_{\\mathrm{H_2O}} R_g T}$ gives the free surface sublimation rate of water ice; because heat conduction equilibrates bodies below roughly 100 m within about a year, this surface formula alone tracks the full cometary-nucleus model. The feedback that carries the result is the drift-ablation loop: sublimation shrinks the radius, a smaller radius drifts faster, down to the fastest meter-size bodies, and faster drift brings the body into hotter regions, accelerating ablation until the body reaches the 10 cm disintegration limit.","core_discovery":"The central claim is that the dynamical snowline for water, the place where drifting solids actually lose their ice, is not the classical temperature-pressure snowline for bodies between about 1 m and 100 m. These bodies drift starward faster than the snowline moves inward, so they cross it; once inside, sublimation shrinks them and the smaller body drifts even faster, with the fastest drift near meter size, so they carry water inward until they shrink to roughly 10 cm and disintegrate. In the paper's terminology, the region polluted with water ice extends to ten percent closer to the star than the snowline location, and this result is nearly independent of disk mass, disk lifetime, and irradiation. A non-eroding dust mantle would let bodies reach about half the classical snowline distance, but the paper argues such mantles are stripped by collisions with pebble-sized objects. For homogeneous dust-water-ice bodies without a mantle, the paper finds the interior is isothermal on timescales of about a year, so the analytic Hertz-Knudsen-Langmuir sublimation expression reproduces the full numerical thermal model.","pith_inferences":["If boulders are as abundant as the paper postulates, the water vapor they release just inside the snowline could locally raise the vapor pressure and slow further sublimation, potentially pushing the polluted region deeper inward than ten percent; the paper only tests this with an imposed vapor profile.","A full coagulation-erosion model could remove boulders before they cross the snowline, which would turn the ten-percent zone into an upper limit; the paper identifies collisional erosion as its main uncertainty.","The same drift-and-ablation mechanism should act on other volatiles with their own snowlines, such as CO2 or CO, and because each volatile has a different sublimation temperature, the smearing width should differ by species in a way that may be observable in resolved disk chemistry.","The near time-independence of the result suggests a practical rule for planet formation codes: deposit water at roughly 0.9 times the snowline radius whenever boulders are present, rather than at the snowline itself."],"forward_implications":["The water ice available to planet formation begins about ten percent closer to the star than the classical snowline whenever meter-to-100 m bodies exist, so sharp-snowline disk models underestimate the inner water reservoir.","The volatile flux carried by drifting boulders can be estimated from the boulder size distribution and composition and should be added to disk chemistry and planet formation models.","For unmantled homogeneous dust-ice bodies up to 100 m, the analytic surface-sublimation formula can replace the full thermal model, reducing computational cost in future simulations.","Because the disintegration location is nearly constant when measured in units of the snowline distance, the predicted inward pollution zone scales across disks with different masses and lifetimes.","If a non-eroding dust mantle is present, icy bodies can reach about half the classical snowline distance, but pebble collisions should remove such mantles, keeping the unmantled result closer to reality."],"supporting_citations":[{"why":"supplies the basic gas-drag radial drift mechanism for solids in the disk.","marker":"Weidenschilling 1977"},{"why":"provides the drift solutions underlying the radial drift formula used throughout.","marker":"Adachi et al. 1976"},{"why":"gives the specific radial drift formula and stopping-time expressions adopted for the three drag regimes.","marker":"Chambers 2008"},{"why":"provides the analytic midplane temperature expressions used to locate the snowline and drive sublimation.","marker":"Nakamoto & Nakagawa 1994"},{"why":"supplies the disk temperature model with viscous and stellar irradiation heating.","marker":"Hueso & Guillot 2005"},{"why":"contributes the cometary nucleus model that is the full thermal benchmark against which the analytic ablation formula is tested.","marker":"Marboeuf et al. 2012"},{"why":"shows a prior application of similar sublimation reasoning to ice survival on small bodies.","marker":"Schorghofer 2008"},{"why":"establishes that small bodies reach gas-temperature equilibrium and quantifies the negligible frictional heating.","marker":"D'Angelo & Podolak 2015"},{"why":"supplies the laboratory-based erosion efficiencies used to bound the collision-free assumption.","marker":"Windmark et al. 2012"}],"fun_headline_variants":["Drifting boulders push water ice 10% past snowline","Meter-sized rocks carry water beyond the snowline","Boulder drift extends water ice region by 10%","Water ice rides boulders past the classical snowline","Boulders deliver H2O ice 10% closer to star"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a drifting boulder survives intact while it crosses the snowline; the paper's own collision-rate calculation finds that with standard laboratory erosion rates a 10 m body loses roughly $8\\times10^{-2}$ percent of its mass per year, making the collision-free assumption valid only for about ten years, far shorter than the hundreds to thousands of years of sublimation modeled.","fun_headline_variants_meta":{"raw":{"variants":["Drifting boulders push water ice 10% past snowline","Meter-sized rocks carry water beyond the snowline","Boulder drift extends water ice region by 10%","Water ice rides boulders past the classical snowline","Boulders deliver H2O ice 10% closer to star"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000217,"raw_usage":{"total_tokens":1503,"prompt_tokens":1078,"completion_tokens":425,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":694,"completion_tokens_details":{"reasoning_tokens":342}},"tokens_in":694,"tokens_out":425,"duration_ms":5542,"temperature":1.0,"reasoning_tokens":342,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:42:00.311587+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The claim would be settled by measuring the impact-erosion lifetime of meter-sized icy bodies in the relevant velocity regime: if such bodies lose more than a few percent of their mass per year, they would disintegrate before crossing the snowline and the predicted ten-percent inward water-ice zone would not form; if erosion is much slower, the zone is expected. A complementary observation is a spatially resolved water-vapor or ice map just inside the snowline of a boulder-rich disk, looking for the predicted inward extension rather than a sharp cutoff.","supporting_citations":[{"cited_title":"& Nakagawa, Y","cited_arxiv_id":null,"evidence_quote":"provides the analytic midplane temperature expressions used to locate the snowline and drive sublimation."},{"cited_title":"& Guillot, T","cited_arxiv_id":null,"evidence_quote":"supplies the disk temperature model with viscous and stellar irradiation heating."},{"cited_title":"2012, A & A, 542, A82","cited_arxiv_id":null,"evidence_quote":"contributes the cometary nucleus model that is the full thermal benchmark against which the analytic ablation formula is tested."},{"cited_title":"2008, ApJ, 682, 697","cited_arxiv_id":null,"evidence_quote":"shows a prior application of similar sublimation reasoning to ice survival on small bodies."},{"cited_title":"2012, A & A, 540, A73","cited_arxiv_id":null,"evidence_quote":"supplies the laboratory-based erosion efficiencies used to bound the collision-free assumption."}],"review_version":1}