{"id":"1d1a3510-ef7e-42c3-a5ee-d9cf3c8ade79","arxiv_id":"2412.03533","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Impacts from disk particles can drain up to about 30-50% of a young planetesimal's spin angular momentum, but not enough alone to allow pebble cloud collapse into a single planetesimal.","lead":"This paper calculates how much spin a young planetesimal loses when particles from the surrounding disk crash into it. The effect is real but removes only 30-50% of the spin, so it helps but cannot alone make pebble clouds collapse into single planetesimals.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Torque magnitude depends on azimuthal ejecta functions (Eqs. 20–23) extrapolated from km/s iSALE and ~100 m/s sand data to 10–65 m/s and to impact angles below the 30° calibration limit; until this extrapolation is tested, the 30–50% drain estimate is not robust.","rationale":"I read the paper as a careful application of crater-scaling laws to a new regime; the central mechanism (ejecta from low-velocity impacts carries away spin angular momentum) is physically plausible and follows the Dobrovolskis-Burns logic, and the authors include explicit caveats about the unvalidated high-obliquity, low-velocity regime. The strongest quantitative claim—30–50% spin removal—depends on ϵL values computed with Eqs. 20–23. Those functions are the least secure link: they are fits to km/s iSALE and ~100 m/s sand data, applied here at 10–65 m/s and to θimp<30°, and the authors note that the resulting spin-up reversal for strengthless bodies may not be physical. Since many impacts on a sphere are oblique, the net torque is sensitive to this extrapolation. A clamping test at θimp=30° would show whether the headline result survives outside the uncalibrated region. I agree with the reader's identified weakest assumption and therefore would not change the CONDITIONAL verdict. I also note that the claimed code archive (an arXiv link) is not a functional repository, but that is a reproducibility issue rather than a scientific flaw.","tokens_in":32879,"tokens_out":8969,"duration_ms":92937,"concrete_test":"Recompute the fiducial case (Ra=10 km, ρa=235 kg m^-3, Ya=100 Pa, uhw=30 m/s) and the Figure 8c family with all impacts having θimp<30° re-evaluated at the calibration limit θimp=30° (clamping the uncalibrated grazing impacts to the boundary of Raducan et al.'s fit range). If the resulting ϵL (and hence the total drain) changes by less than ~20%, the 30–50% claim is insensitive to the extrapolation and the concern does not land; if ϵL drops by more than a factor of 2 or changes sign, the headline number is not supported without new low-velocity (10–65 m/s) oblique-impact experiments or simulations at θimp<30°.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline 30–50% spin-down (Section 7) rests on the angular momentum transfer efficiency ϵL from Eq. 36, which is computed using the azimuthally asymmetric ejecta functions μ(ζ,θ), C1(ζ,θ), Rcr(ζ,θ), and k(ζ,θ) in Eqs. 20–23. These functions are fitted to km/s-scale iSALE simulations (Raducan et al. 2022) and ~100 m/s sand-crater measurements (Quillen and Doran 2024), yet they are applied here at 10–65 m/s and to grazing impacts with θimp below 30°, outside the stated calibration range (Sec. 5.4). The spin-up reversal seen for strengthless bodies and high density ratios is explicitly acknowledged by the authors as possibly unphysical (Sec. 5.4). Because a substantial fraction of impacts on a sphere occur at oblique angles, the sign and magnitude of the net torque are controlled by exactly this unvalidated part of the model. If the butterfly pattern or the crater-radius asymmetry is weaker at low velocity, ϵL could be much smaller than the near-unity values used to derive the 30–50% claim; if the asymmetry reverses, the mechanism could even spin up the planetesimal. This is an extrapolation risk, not an internal inconsistency, but it is load-bearing because the paper's quantitative conclusion is a direct product of these functions.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies crater ejecta scaling laws to compute the mass loss and angular momentum transfer experienced by a planetesimal bombarded by low-velocity headwind particles (10–65 m/s) in a protoplanetary disk. The authors integrate oblique-impact ejecta distributions over a distribution of impacts on a spherical planetesimal, including gravitational focusing, and define an angular momentum transfer efficiency ϵL and a mass-loss ratio ϵM. They report that the Dobrovolskis and Burns (1984) spin-down mechanism operates in this low-velocity regime, with ϵL ∼ 1 for a wide range of parameters, and that under assumed pebble mass-flux conditions the cumulative angular momentum drain is 30–50% of the initial spin. The paper closes with caveats about highly oblique impacts and notes that the adopted oblique-impact scaling functions are inconsistent with one source figure.","tokens_in":33296,"tokens_out":2587,"duration_ms":30210,"significance":"If the result holds, the paper provides a quantitative pathway by which low-velocity impacts can remove a substantial fraction of a planetesimal's spin angular momentum before pebble-cloud collapse, connecting impact ejecta physics to the angular-momentum problem in streaming-instability planetesimal formation. The work has clear strengths: the transfer efficiency is computed rather than fitted to a target spin-down, the gravitational-focusing treatment is explicit, the code is archived with the submission, and the paper identifies its own most fragile assumption in Section 5.4. These strengths make the paper a useful contribution even though the headline 30–50% number inherits uncertainty from the ejecta-scaling extrapolation.","major_comments":[{"comment":"The central quantitative claim — ϵL ∼ 1 and the resulting 30–50% spin-down — is controlled by the azimuthally asymmetric ejecta functions μ(ζ,θ), C1(ζ,θ), Rcr(ζ,θ), and k(ζ,θ). These functions are calibrated to km/s-scale iSALE simulations (Raducan et al. 2022) and roughly 100 m/s sand-crater measurements (Quillen and Doran 2024), yet they are applied here at 10–65 m/s and to impact angles below the 30° calibration limit. As the authors note, the resulting spin-up reversal for strengthless bodies may not be physical. Because a large fraction of impacts on a sphere occur at highly oblique angles, the sign and magnitude of the net torque are set by exactly this unvalidated part of the model. I ask the authors to either (a) add a sensitivity test that suppresses or reverses the butterfly-pattern asymmetry and shows how ϵL changes, or (b) bound ϵL using the low-velocity experimental data available for oblique impacts, or (c) explicitly degrade the 30–50% claim to a qualitative statement pending such tests. This is an extrapolation risk rather than an internal inconsistency, but it is load-bearing for the paper's main numerical conclusion.","section":"Sec. 5.4 and Eqs. (20)–(23)"},{"comment":"The text states that Eqs. (20)–(23), taken from Raducan et al. (2022), 'are not consistent with their Figure 12 but appear consistent with the distributions shown in their other figures.' This is an explicit admission that the adopted crater-radius asymmetry function Rcr(ζ,θimp) (Eq. 22) disagrees with a stated calibration figure. Since Rcr sets the integration domain in Eq. (31), the inconsistency is not cosmetic. The authors should either reconcile the adopted functional form with the source's Figure 12, show quantitatively how much the discrepancy changes Rcr and hence ϵL, or replace Eq. (22) with a form that is actually consistent with the cited calibration. At minimum, the source of the discrepancy needs to be explained, not merely noted.","section":"Sec. 2.4.1, Eq. (22)"},{"comment":"The 30–50% total spin-down estimate combines the computed ϵL with the pebble mass flux estimate in Eq. (57), but the paper provides no uncertainty quantification on either factor. The mass-flux estimate depends on fp = 10^2, Δt = 10^4 yr, the disk model of Quillen et al. (2024), and the neglect of gas drag on incoming and outgoing particles; the ϵL values enter near unity in the fiducial regime. Because the authors do not propagate uncertainties or show how the 30–50% range follows from the parameter ranges explored in Figure 8, this headline number is presented with unjustified precision. I recommend adding a short propagation or bracketing calculation that translates plausible ranges in fp, Δt, and the oblique-scaling uncertainty into a range for the cumulative spin-down, with the caveats of Section 5.4 folded in.","section":"Sec. 7 and Eq. (57)"}],"minor_comments":[{"comment":"The caption says 'Panels (a,b) varies π4, (b,c) varies planetesimal strength Ya, (e,f) varies the distant headwind velocity uhw, and (g,h) varies the planetesimal spin rate.' The panel labels are inconsistent: the second pair should be (c,d), not (b,c). Please correct the caption.","section":"Figure 8 caption"},{"comment":"There are many typographical artifacts from the manuscript preparation, such as 'ine fficient', 'su fficiently', 'di fferent', and 'e ffect'. These should be cleaned up before publication.","section":"Throughout"},{"comment":"The paper says in Section 2.9.1 that neglecting gravitational defocusing 'does not influence our estimate of the angular momentum transfer efficiency', but in Section 5.2 it suggests that gravitational defocusing 'could be an even stronger effect on the escaping ejecta than on the incoming projectile.' Please clarify whether defocusing affects the angular momentum computation or only the linear momentum computation; the current statements are easy to read as contradictory.","section":"Sec. 2.9.1 and Sec. 5.2"},{"comment":"The code is said to be 'archived at https://arxiv.org/abs/2412.03533', which is the paper's own arXiv page rather than a code repository. Please provide a persistent code repository (e.g., Zenodo or GitHub) with a version identifier.","section":"Acknowledgements"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of Icarus and the central idea is worth publishing once the extrapolation issue is treated honestly. The authors themselves flag the two weakest links: the inconsistency of Eqs. (20)–(23) with the cited Raducan et al. figure, and the possibly unphysical spin-up reversal for strengthless bodies. These are not fatal defects, but they directly control the paper's quantitative headline, so I do not see how the 30–50% claim can stand without either new low-velocity oblique-impact validation or a sensitivity analysis that removes the load-bearing weight of the extrapolated functions. The paper's framing as a mechanism that partially, not fully, solves the angular-momentum problem is appropriate and should be preserved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take. This is a genuine extension of Dobrovolskis and Burns to a regime they never touched: planetesimals embedded in a protostellar disk, hit from a preferred direction at 10-65 m/s. The central claim—that low-velocity impacts can drain a substantial fraction of spin—is plausible and survives a skeptical read. The headline 30-50% number, though, is an order-of-magnitude estimate resting on extrapolated azimuthal ejecta functions and a rough disk mass-flux model, not a robust prediction.\n\nWhat's actually new: the paper integrates ejecta distributions with modern azimuthal sensitivity (Raducan et al. 2022, Quillen and Doran 2024), includes gravitational focusing via hyperbolic orbits, and treats a unidirectional headwind perpendicular to the spin axis. That combination matters because at low velocity the ejecta escape fraction is sensitive to rotation and impact angle. The efficiency epsilon_L is computed from crater scaling, not fitted to a desired spin-down, and the paper covers density ratio, strength, headwind speed, spin rate, and a filament velocity gradient. Credit where due: the integrals are internally consistent, and the authors are unusually candid. Section 5.4 flags the possible unphysical spin-up reversal from the butterfly pattern, Section 2.9 lists neglected physics, and they explicitly note their adopted Raducan functions are inconsistent with that paper's Figure 12.\n\nThe soft spots are real but not fatal. The load-bearing issue: Eqs. 20-23 are fitted to km/s iSALE simulations and ~100 m/s sand experiments, and are applied at 10-65 m/s and impact angles below the 30-degree calibration limit. Since a large fraction of impacts on a sphere are oblique, the net torque's sign and magnitude are controlled by exactly this unvalidated territory. The stress-test note lands here; if the azimuthal asymmetry weakens at low velocity, epsilon_L could drop well below the near-unity values used for the 30-50% claim. There is no uncertainty quantification, and the mass-flux normalization depends on fp~100 and a 10^4 yr streaming-instability epoch—reasonable but soft. Minor: the code archive is just the arXiv page, not a runnable repository, and the ejecta angle is fixed at 45 degrees.\n\nOverall: the mechanism operates, and the paper does a service by showing it is partial—30-50% drain is not enough to collapse a single planetesimal from typical cloud angular momentum at 45 AU. That conclusion is probably robust even if the exact percentage is not.\n\nRecommendation: send to peer review. A serious referee should push for sensitivity tests of the azimuthal functions at low velocity/obliquity and some uncertainty range on the final estimate. If those tests show the qualitative result survives, this is a solid contribution to planetesimal formation.","headline":"A careful extension of angular momentum drain to low-velocity embedded planetesimals; the mechanism is real and the 30-50% estimate is plausible but rests on unvalidated ejecta extrapolations.","tokens_in":33803,"tokens_out":3398,"would_cite":true,"duration_ms":35648,"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":"Low-velocity headwind impacts spin down young planetesimals by 30–50%, helping pebble-cloud collapse without completing it.","keywords":["ejecta","impact phenomena","collisional processes","planetesimals","angular momentum drain","oblique impacts","crater scaling","protoplanetary disk"],"falsifier":"Run oblique-impact experiments into weakly cohesive and strengthless granular targets at impact speeds of 10–65 m/s and impact angles down to grazing, measuring the azimuthal ejecta mass and velocity distributions; compare the measured distributions to Equations 20–23. If the low-speed, high-obliquity ejecta do not follow these functions, the predicted spin-down efficiency, the 30–50% estimate, and the spin-up reversal for strengthless bodies would all need revision.","tokens_in":32682,"feed_emoji":"☄️","tokens_out":9325,"duration_ms":87505,"temperature":0.7,"pith_summary":"This paper asks whether low-velocity impacts from particles in a protoplanetary disk can spin down a young, embedded planetesimal. The authors combine crater-scaling ejecta laws with azimuthally asymmetric ejecta distributions for oblique impacts, integrate over a realistic headwind impact population, and compute the torque on an Arrokoth-like body. They find that the asteroidal angular-momentum-drain mechanism operates at 10–65 m/s impact speeds, with an angular-momentum transfer efficiency near unity for many parameter choices. Integrating over the expected collisional mass during a streaming-instability epoch, they estimate that impacts remove 30–50% of the initial spin angular momentum. That is real but partial: it can help a pebble cloud collapse into a single planetesimal only when the cloud starts with low angular momentum.","feed_headline":"Headwind impacts shed 30–50% of a young planetesimal's spin","feed_subtitle":"A partial spin-down mechanism makes pebble-cloud collapse easier, but alone it cannot form single planetesimals.","key_machinery":"The load-bearing machinery is a numerical integration that combines four ingredients: normal-impact crater ejecta scaling laws that give ejecta velocity and mass as functions of radial distance from the impact point; analytic functions (Equations 20–23) that modulate the scaling parameters $\\mu$, $C_1$, crater radius $R_{cr}$, and the mass coefficient $k$ with impact angle $\\theta_{\\rm imp}$ and azimuthal angle $\\zeta$, capturing the asymmetry of oblique-impact ejecta curtains; gravitational focusing via hyperbolic projectile orbits, which changes impact location, velocity, and angle; and a dimensionless efficiency $\\epsilon_L$ that measures fractional spin angular momentum lost per unit projectile mass. Escaping ejecta are identified by comparing their inertial-frame speed to the escape velocity; because the surface rotates, ejecta leaving the trailing side move faster in the inertial frame than those on the leading side, so more trailing-side ejecta escape and spin angular momentum is carried away. The same integration yields the mass loss ratio $\\epsilon_M$, separating accretionary from erosional regimes.","core_discovery":"On the paper's own terms, the central discovery is that the angular-momentum-drain mechanism, in which escaping ejecta carry away spin angular momentum preferentially from a rotating body, works in the low-velocity, embedded-planetesimal regime, extending a mechanism previously studied for high-velocity asteroid impacts. Using crater scaling laws modified for oblique impacts and for ejecta azimuthal angle, the authors compute the mass loss ratio and the dimensionless angular momentum transfer efficiency $\\epsilon_L = -|\\langle\\Delta L_a\\rangle|/(L_a/M_a\\, m_{pj})$. They find $\\epsilon_L \\approx 1$ across a wide range of parameters, with stronger drain when the headwind is fast, when the projectile density exceeds the planetesimal density, and when cratering operates in the gravity regime. For the fiducial Arrokoth-like parameters, most impacts cause net accretion rather than erosion, yet the drain still operates; over the collisional era the expected total projectile mass is less than the planetesimal mass, giving a total spin-down of roughly 30–50%. The paper concludes that this partial despinning can facilitate collapse of pebble clouds with low initial angular momentum but is insufficient by itself to explain single-planetesimal formation.","pith_inferences":["Beyond the paper, a direct test of the azimuthal ejecta functions at 10–65 m/s into weakly cohesive granular targets would clarify whether the 30–50% estimate stands; the paper itself flags the low-velocity, high-obliquity regime as uncalibrated.","Beyond the paper, the same integration framework could be extended to compute the linear momentum transfer, effectively an impact-induced drag, on an embedded planetesimal; the paper notes this could tighten binaries or assist binary merging during the high-flux streaming-instability epoch.","Beyond the paper, if the predicted spin-up reversal for strengthless bodies at small radii is physical, it would mean extremely low-strength planetesimals are not efficiently despun, narrowing the conditions under which single-planetesimal collapse can occur."],"forward_implications":["Angular momentum drain works at 10–65 m/s impact speeds, not just at the km/s speeds typical of asteroid impacts, and it operates in both accretionary and erosional impact regimes.","The most favorable conditions are a fast headwind, a projectile denser than the planetesimal, and a planetesimal weak enough that cratering follows gravity-regime scaling; varying the projectile-to-planetesimal density ratio is the strongest lever on efficiency.","For an Arrokoth-like body at 45 AU over a 10^4-year streaming-instability epoch, about half the planetesimal's mass in pebbles impacts it, and the integrated effect is a 30–50% reduction in spin angular momentum.","That reduction is large but does not solve the angular momentum problem by itself: at 45 AU roughly 93% of the pebble cloud's angular momentum must be removed for collapse into a single planetesimal, so additional drain mechanisms are needed.","The pebble-cloud velocity gradient, even at its assumed upper bound, changes the drain efficiency and mass loss ratio by less than detectable amounts, so a unidirectional headwind is a good approximation."],"supporting_citations":[{"why":"Defines the angular momentum drain mechanism of escaping ejecta carrying away spin, which this paper extends from km/s asteroid impacts to low-velocity disk impacts.","marker":"(Dobrovolskis and Burns, 1984)"},{"why":"Supplies the normal-impact crater ejecta scaling laws for radius, velocity, and mass that underlie all ejecta integrations.","marker":"(Housen and Holsapple, 2011)"},{"why":"Provides the azimuthal and impact-angle-dependent functions for the scaling parameters used to make the ejecta distributions oblique.","marker":"(Raducan et al., 2022)"},{"why":"Gives roughly 100 m/s laboratory measurements of azimuthal ejecta mass from oblique impacts into sand, supporting the weak angular dependence of the ejecta mass coefficient used here.","marker":"(Quillen and Doran, 2024)"},{"why":"Supplies the disk model, headwind velocities, impact flux estimate, and planetesimal wind-erosion context that set the impact population and collisional mass.","marker":"(Quillen et al., 2024)"},{"why":"Provides the pebble-cloud angular momentum distribution and the collapse-to-single-planetesimal problem that the despinning mechanism is intended to help solve.","marker":"(Nesvorný et al., 2021)"},{"why":"Gives the Arrokoth-based fiducial planetesimal density, radius, and strength range used in the integrations.","marker":"(McKinnon et al., 2020)"}],"fun_headline_variants":["Angular momentum drain works in low-velocity regime","Pebble-cloud impacts spin down planetesimals by 50%","Ejecta torque can drain half a planetesimal's spin","Partial spin-down from impacts eases planetesimal collapse"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The results depend on ejecta functions that were fit to km/s-scale impact simulations and roughly 100 m/s sand impacts but are applied here at 10–65 m/s and at more extreme impact obliquities; if those functions fail in this regime, the computed torque, and even its sign, could be wrong.","fun_headline_variants_meta":{"raw":{"variants":["Angular momentum drain works in low-velocity regime","Pebble-cloud impacts spin down planetesimals by 50%","Ejecta torque can drain half a planetesimal's spin","Partial spin-down from impacts eases planetesimal collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00089,"raw_usage":{"total_tokens":3855,"prompt_tokens":973,"completion_tokens":2882,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":2813}},"tokens_in":589,"tokens_out":2882,"duration_ms":21046,"temperature":1.0,"reasoning_tokens":2813,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T22:17:35.099697+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run oblique-impact experiments into weakly cohesive and strengthless granular targets at impact speeds of 10–65 m/s and impact angles down to grazing, measuring the azimuthal ejecta mass and velocity distributions; compare the measured distributions to Equations 20–23. If the low-speed, high-obliquity ejecta do not follow these functions, the predicted spin-down efficiency, the 30–50% estimate, and the spin-up reversal for strengthless bodies would all need revision.","supporting_citations":[],"review_version":1}