{"id":"b81d370f-8e55-48c7-b02a-51dfe00d7e48","arxiv_id":"2507.16739","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Soft-sphere contact simulations of gravitationally collapsing pebble clouds form binary planetesimal systems and resolve, for the first time, their spins and shapes.","lead":"Simulations of collapsing pebble clouds with contact physics instead of instant mergers produce planetesimals with measured spins, shapes, and tighter binary orbits than earlier models allowed. If correct, the results give a first look at the initial rotation and shape statistics that relict Kuiper Belt objects and asteroids may preserve.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"1-AU collapse likely misses gas drag and disk tides; spin/shape predictions rest on an untested orbital-radius assumption.","rationale":"The reader identified the heliocentric-orbit scaling assumption as the weakest point, and this stress-test agrees. The paper's core claim—resolving spins, shapes, and tight binaries from gravitational collapse—requires a valid mapping from the simulated 1-AU collapse to the observed Kuiper Belt. The mapping is argued but never tested: no 30-AU SSDEM simulation is presented, and the argument omits gas drag and tidal effects that differ strongly between 1 AU and 30 AU. The concern is not speculative; the paper itself acknowledges in the Discussion that gas drag has been proposed to harden binaries, and its own limitation statement says the long-term evolution is left to future work. This is an externally testable assumption, so the appropriate verdict remains CONDITIONAL: the paper's orbital and binary results are internally coherent, but the spin and shape comparisons to relict populations are not yet validated. No ad hominem or theatrical framing is needed; the issue is simply that the central comparison is built on an untested scaling step. Agreement is full because the reader's weakest_assumption is precisely this scaling step.","tokens_in":36855,"tokens_out":1283,"duration_ms":13920,"concrete_test":"Run the SSDEM collapse at a⊙ = 30 AU with identical cloud mass, particle number, particle radius, rotation rate f, and contact parameters, and compare accreted planetesimal mass ratios, binary orbits, spin periods, and shape frequencies to the 1-AU simulations. If the orbital and shape distributions differ by more than the run-to-run scatter, the 1-AU-to-30-AU scaling assumption fails and the comparisons to Kuiper Belt binaries, spins, and shapes cannot stand without accounting for gas drag and tides. Alternatively, run a 1-AU simulation with gas drag and tidal terms included and show they do not change the outcome.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The paper's novel claims depend on treating the 1-AU simulations as proxies for Kuiper Belt collapse: Section 2 asserts that heliocentric orbit does not significantly affect accreted planetesimal properties, arguing that collision outcomes scale with surface escape velocity. This ignores gas drag, which at 1 AU is vastly stronger than at 30 AU and can preferentially remove angular momentum from the collapsing cloud, and it ignores the heliocentric tidal field, whose ratio to the cloud's self-gravity changes with orbital radius. The statement is not backed by a dedicated simulation set—no 30-AU SSDEM runs or a controlled orbital-radius comparison are presented. The spin and shape distributions, including the 10-hr average period that is highlighted in the abstract and compared to relict populations, are therefore only secure if these environmental differences are genuinely negligible. Since no direct test is shown, the paper's own limitation passage leaves the central comparison to Kuiper Belt observations unsupported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the soft-sphere discrete element method (SSDEM) within PKDGRAV to simulate the gravitational collapse of 10^5 super-particle clouds at 1 AU, with the stated goal of predicting the mutual orbits, spin states, and shapes of newly formed planetesimals. Two simulation suites are run: one varying the initial cloud rotation scale factor f at fixed contact parameters, and one varying the coefficients of restitution and friction at two rotation rates. The authors report that SSDEM reproduces the binary mass ratios and accretion efficiencies of earlier perfect-merger models while additionally producing resolved spins, shapes, and tighter binary orbits. The simulated binary orbit properties are compared with observed trans-Neptunian binaries, and the spin and shape distributions are compared with asteroid and Kuiper Belt observations. The paper concludes that gravitational collapse can directly produce the observed ranges of binary orbits and that the simulated spin and shape statistics are the first direct predictions of these properties from this formation channel.","tokens_in":37084,"tokens_out":4666,"duration_ms":51010,"significance":"If the results hold, this is a substantial methodological advance: it is the first SSDEM treatment of pebble-cloud gravitational collapse that resolves planetesimal spins, shapes, and tight binary orbits, and it makes falsifiable predictions that can be compared with relict Kuiper Belt and asteroid populations. The work has concrete strengths: it explicitly validates against prior perfect-merger simulations, uses a well-established SSDEM code lineage, and compares rather than fits to observations, so there is no circularity in the central comparison. The paper also correctly identifies and discusses several limitations, including the lack of very tight observed binaries in the simulations and the absence of collisional evolution in the spin predictions. However, the external validity of the spin and shape predictions depends on an untested environmental assumption (1 AU vs 30 AU), and some of the headline reporting is internally inconsistent; these issues need to be addressed before the observational comparisons can be taken as secure.","major_comments":[{"comment":"The claim that \"the choice of the cloud's heliocentric orbit does not have a significant effect on the accreted planetesimal properties\" is load-bearing for the paper's comparisons to Kuiper Belt observations in Figures 8, 9, 11, and 16, yet it is supported only by the scaling argument that collision velocities scale with surface escape velocity. No simulation at 30 AU or with a controlled variation of heliocentric distance is presented, and the model omits gas drag and the changing ratio of the heliocentric tidal field to the cloud's self-gravity. Because the novel spin and shape predictions are the central contribution, this assumption needs a dedicated numerical test, even a single 30 AU run or a controlled variation of a_sun, or the observational comparisons must be explicitly qualified as conditional on the environment being dynamically equivalent.","section":"Section 2"},{"comment":"The abstract states that \"Newly-formed planetesimals exhibit 10-hr rotation periods on average,\" but Section 3.5 reports mean rotation periods of 12.9 hr (sigma = 8.8 hr) and 13.3 hr (sigma = 8.5 hr) for the two simulation suites. The 10-hr value appears to describe only the largest planetesimals, or the comparison populations, rather than the full simulated sample. The abstract should be corrected so that the headline number matches the reported statistics.","section":"Abstract and Section 3.5"},{"comment":"The six shape categories and the frequencies in Table 4 are assigned \"by eye\" (stated twice in the text), with no quantitative decision rule and no reproducibility check. Because the shape frequencies are a novel headline result, the classification should be made reproducible—for example, by defining thresholds on the measured oblateness gamma and prolateness beta plus an additional metric for ridge or egg-like features—or by releasing the labeled dataset and classification code alongside the processed data.","section":"Section 3.6"},{"comment":"The normal spring constant is chosen so that maximum particle overlaps are about 20%, deliberately relaxing the <1% criterion used in prior SSDEM work (Schwartz et al. 2012). Since the paper's new spin and shape results emerge directly from contact-force resolution, a convergence test with a stiffer spring, at least for one representative rotation rate and contact-physics case, is needed to rule out an artifact of the soft contact model.","section":"Section 2.1"}],"minor_comments":[{"comment":"The text refers to \"286958 Arrokoth\" but the correct Minor Planet Center designation is 486958 Arrokoth.","section":"Figure 16"},{"comment":"There is a small typo in the introduction: \"so that the they can rest upon each other\" should read \"so that they can rest upon each other.\"","section":"Section 1"},{"comment":"The statement that the largest simulated planetesimals have a \"similar mean spin period\" to observed populations is based on a subset of the simulated population; please clarify in the text that the 10-hr figure refers to the largest bodies rather than the full sample.","section":"Section 3.5"}],"recommendation":"major_revision","confidential_remarks":"The main unresolved issue is the environmental extrapolation from 1 AU to the Kuiper Belt. The paper's own scaling argument is plausible but is not a substitute for a numerical test, and the abstract's rotation-period number should be fixed. If the authors add a 30 AU control run or substantially temper the observational claims, and make the shape classification reproducible, the paper would be a strong contribution to the planetesimal formation literature."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this one is worth reading. The authors apply SSDEM in PKDGRAV to run 54 gravitational collapse simulations of 100-km-class clouds, and for the first time they resolve spin states, shapes, and tight binaries directly out of the collapse. The binary orbits and mass ratios agree well with the older perfect-merger results from Nesvorný and Robinson, which validates the new method. That is a real achievement, and the tight-binary population (a/RHill around 0.02–0.15) is a genuinely new result that perfect-merger models cannot produce. The spin and shape distributions, though preliminary, are direct predictions from a plausible formation channel and give observers something to test.\n\nThe soft spots are real but mostly addressable. The abstract says \"10-hr rotation periods on average\" but the reported means are 12.9 and 13.3 hours; that mismatch needs fixing. The shape taxonomy is by eye, which is fine for a first pass but should be backed by quantitative criteria. Data are only promised for Zenodo, not deposited, so the paper is not independently reproducible yet. The biggest substantive concern is the collapse at 1 AU used as a proxy for the Kuiper Belt. The paper's scaling argument is that collision outcomes depend on relative velocities, which scale with surface escape speed, not on orbital radius. That is plausible for the internal dynamics, but gas drag is far stronger at 1 AU and the tidal field relative to self-gravity changes with radius; neither is tested with a 30-AU run or a controlled comparison. Without that test, the comparisons of spins and shapes to Kuiper Belt objects rest on an unverified assumption. It is a moderate flaw, not a fatal one, because the contact-physics results likely transfer, but the paper should either add a 30-AU run or damp the observation-matching language.\n\nI would send this to a serious referee. The method is sound, the new observables are important, and the flaws are fixable in revision. A good referee can push for the missing scaling test and better quantitative shape metrics. I'd cite it once the data are out.","headline":"First resolved spin/shape/tight-binary outcomes from SSDEM collapse are worth reading; the 1-AU-to-Kuiper-Belt scaling is the main caveat.","tokens_in":37566,"tokens_out":2659,"would_cite":true,"duration_ms":28017,"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":"Simulating planetesimal collapse with contact physics yields first predictions of spins and shapes.","keywords":["Planetesimal formation","Gravitational collapse","Soft-sphere discrete element method","Binary planetesimals","Spin states","Shapes","Kuiper Belt objects","Numerical simulation"],"falsifier":"A direct test would be to run the same gravitational-collapse simulations at a heliocentric distance of 30 AU (or with a different gas drag environment) and compare the resulting binary orbits, spins, and shapes to the 1 AU simulations; if the outcomes differ significantly, the assumption that orbit choice does not affect accreted properties is invalidated.","tokens_in":1618,"feed_emoji":"🪐","tokens_out":2587,"duration_ms":41952,"temperature":0.7,"pith_summary":"This paper argues that the gravitational collapse of pebble clouds, a leading route to planetesimal formation, can now be modeled with realistic contact physics using the soft-sphere discrete element method (SSDEM) in PKDGRAV. The authors show that the method reproduces the binary-formation results of earlier perfect-merger models while additionally resolving planetesimal spins, shapes, and tight binary orbits. If the approach is correct, the simulated populations—with mean rotation periods near 10 hours and a variety of shapes from spherical to prolate—provide the first direct predictions of these properties from the gravitational-collapse formation channel, testable against observed asteroids and trans-Neptunian objects.","feed_headline":"Collapse simulations now resolve planetesimal spins and shapes","feed_subtitle":"Soft-sphere contacts produce 10-hour spins and six shape classes, matching relict asteroids and Kuiper Belt objects.","key_machinery":"The soft-sphere discrete element method (SSDEM) within the PKDGRAV N-body integrator is the central tool. Instead of treating collisions as perfect mergers of inflated particles, SSDEM simulates contact forces with spring-dashpot models, friction, and restitution, allowing super-particles to rest, roll, bounce, and aggregate into realistic, volume-filling planetesimals. This enables the resolution of shapes, spins, and tight binary orbits, and the modeling of both accretion and decretion during collapse. The simulations use a small timestep of 4.7 s at 1 AU, with the cloud collapse timescale set by the heliocentric distance; they show that the choice of orbit does not affect outcomes because collision velocities scale with surface escape velocity.","core_discovery":"The central claim is that SSDEM can model the collapse of a cloud of super-particles into planetesimals while tracking mutual orbits, spin states, and shapes, and that the resulting systems match and extend earlier perfect-merger results. Simulations produce many binary systems per cloud, with the most massive ones on tight, low-inclination orbits (a/RHill ~ 0.02–0.15, i ≲ 15°, e ≲ 0.40), while less massive systems span a wider range. Newly formed planetesimals spin with a mean period near 10 hours, with the largest objects rotating faster when formed from slowly rotating clouds. Six shape classes emerge—spherical, oblate, top-shaped, flattened, egg-shaped, and prolate—with the most massive planetesimals typically spherical or oblate. The shapes correspond to low internal friction angles (ϕ ≲ 10°), indicating relaxation toward a low-strength state during collapse.","pith_inferences":["The paper's use of a single cloud mass and orbit implies that spin and shape predictions depend primarily on cloud angular momentum, not location or mass; extending to a range of cloud masses and orbits would test whether the ~10-hour rotation scale is universal.","The absence of simulated rotators with periods under about 6 hours, while observed small asteroids spin down to ~2.4 hours, could indicate a formation limit or a missing process; targeted simulations with different initial angular momentum distributions could distinguish these.","The shape taxonomy relies on by-eye classification; a quantitative, unsupervised classification of the full 793-object sample would reveal whether the shape classes are discrete or continuous, sharpening comparisons to observed asteroids.","The comparison to trans-Neptunian binaries implicitly assumes alignment of the mutual orbital plane with the cloud's angular momentum; real clouds may have a spread in orientations, which would broaden the inclination distributions."],"forward_implications":["If the method is accepted, the simulated spin and shape distributions become the first direct predictions of planetesimal properties from gravitational collapse, testable against asteroid and Kuiper Belt observations.","The match between simulated large-planetesimal spins (7–70 hours) and observed large asteroid and trans-Neptunian spins supports the hypothesis that these spins are primordial, set at formation.","The finding that small simulated planetesimals spin slower than comparable asteroids by about 5 hours implies that collisional evolution in the Main Belt has accelerated rotation over 4.5 Gyr.","The prevalence of spherical and oblate shapes with low internal friction angles suggests that collapse assembles bodies near relaxed, low-strength states, informing theories of rubble-pile structure.","The ability to form many tight binaries per cloud, inaccessible to perfect-merger models, strengthens the gravitational-collapse explanation for the binary fraction in the cold classical Kuiper Belt."],"supporting_citations":[{"why":"The baseline perfect-merger collapse model whose binary formation rates and orbit properties the SSDEM simulations are designed to reproduce and extend.","marker":"Nesvorný et al. (2010)"},{"why":"Independent perfect-merger collapse simulations whose binary masses, orbits, and accretion efficiencies are directly compared to the SSDEM results, serving as the primary numerical benchmark.","marker":"Robinson et al. (2020)"},{"why":"The source of the SSDEM formulation and contact-physics parameter conventions used in PKDGRAV, supplying the spring-dashpot and friction model.","marker":"Schwartz et al. (2012)"},{"why":"Provides the Mohr–Coulomb yield criterion and tri-axial ellipsoid stress equations used to infer minimum internal friction angles for the simulated shapes.","marker":"Holsapple (2001)"},{"why":"The compilation of 41 observed trans-Neptunian binary mutual orbits used as the observational comparison set for orbit size, eccentricity, and inclination distributions.","marker":"Grundy (2023)"},{"why":"The asteroid spin-period dataset and non-Maxwellian analysis for large objects, used to argue that the largest simulated planetesimal spins are primordial.","marker":"Pravec et al. (2002)"},{"why":"Provides trans-Neptunian object rotation periods and the binary vs. non-binary spin distinction used to compare simulated spin distributions.","marker":"Thirouin et al. (2014)"},{"why":"The observed Kuiper Belt binary radius-ratio distribution used to validate the simulated binary size ratios.","marker":"Noll et al. (2008)"}],"fun_headline_variants":["Soft-sphere collapse yields planetesimal spins and shapes","Contact physics reveals 10-hour spins and six shape classes","Binary planetesimals form with realistic spins and shapes","Gravitational collapse with contacts matches asteroid shapes","SSDEM resolves spins and shapes in planetesimal collapse"],"cache_read_input_tokens":39808,"weakest_assumption_plain":"The load-bearing premise is that the exact heliocentric orbit of the collapsing cloud does not matter because collision velocities scale with surface escape velocity, so collapse simulations run at 1 AU correctly represent collapse in the Kuiper Belt at 30 AU.","fun_headline_variants_meta":{"raw":{"variants":["Soft-sphere collapse yields planetesimal spins and shapes","Contact physics reveals 10-hour spins and six shape classes","Binary planetesimals form with realistic spins and shapes","Gravitational collapse with contacts matches asteroid shapes","SSDEM resolves spins and shapes in planetesimal collapse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00018,"raw_usage":{"total_tokens":1342,"prompt_tokens":1022,"completion_tokens":320,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":638,"completion_tokens_details":{"reasoning_tokens":244}},"tokens_in":638,"tokens_out":320,"duration_ms":3952,"temperature":1.0,"reasoning_tokens":244,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:03:14.991162+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to run the same gravitational-collapse simulations at a heliocentric distance of 30 AU (or with a different gas drag environment) and compare the resulting binary orbits, spins, and shapes to the 1 AU simulations; if the outcomes differ significantly, the assumption that orbit choice does not affect accreted properties is invalidated.","supporting_citations":[{"cited_title":"E., Fraser , W","cited_arxiv_id":null,"evidence_quote":"Independent perfect-merger collapse simulations whose binary masses, orbits, and accretion efficiencies are directly compared to the SSDEM results, serving as the primary numerical benchmark."}],"review_version":1}