{"id":"ebb84f13-02bf-490d-96b7-5bb76cf0ac64","arxiv_id":"2505.03325","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Sub-catastrophic impacts into rotating porous 100 km asteroids can directly place deep subsurface debris onto stable high-pericenter orbits, forming satellites, when the impact leaves the primary elongated and fast-spinning.","lead":"Using impact simulations that track the spin and shape of the leftover asteroid, this paper finds a direct way for large asteroids to form small moons: an off-center impact stretches and spins up the body, and debris lifted from the elongated side can land in stable orbits. The result matches the observed pattern that satellites orbit mostly fast-spinning, elongated, dark asteroids, and it explains why giant asteroid families do not correlate with having moons.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Post-hoc spin boost at handoff may not reproduce the impact-induced distortion that the formation pathway depends on; the 1 km appendix test does not cover the 100 km regime.","rationale":"The reader identified the same load-bearing weakness: the post hoc rotation approximation. I agree and sharpen it. The mechanism being proposed is not simply that an already-elongated body sheds material; it is that the impact itself, acting on a rotating target, distorts the body and launches deep subsurface material onto stable orbits. The numerical shortcut replaces the coupled impact-plus-rotation dynamics with a rigid velocity boost applied after the shock has fully propagated and the early ejecta field has been set. This is precisely the phase in which the 'immediate distortion' in Fig. 2 originates, so the central claim is not actually simulated in the paper. The appendix test at 1 km / 10 m impactor is reassuring for integrated quantities (shape, spin, remnant mass) but not for the provenance-depth statistics or the high-pericenter mass that the headline result depends on. The paper does not release code or data, and each parameter setting is a single realization, so there is no internal check on the sensitivity of the result to this approximation. Other concerns (neglected debris collisions, no error bars, no released data) are real but secondary; they would affect quantitative yields, not the existence of the pathway. The lacking test is feasible with the authors' own SPH code and would directly settle whether the pathway survives the removal of the approximation. Until then, the paper should be accepted only conditionally on that validation.","tokens_in":18668,"tokens_out":5705,"duration_ms":60902,"concrete_test":"Run the key case (13 km impactor, 60 deg, 5 km/s, 100 km target) with SPH initialized with a uniformly rotating target at 10 hr and, separately, at 4 hr, using the same porosity, resolution, and handoff time. Then hand off to pkdgrav with identical settings and compare with the paper's post-hoc-boost outcomes for those spin periods: (i) largest remnant axis ratios and spin at 38 hr; (ii) the provenance depth distribution of temporary satellites with q > 1.5 a_pri; (iii) total mass on high-pericenter orbits. The approximation is validated only if all three agree within the spread seen between the two 6 hr treatments in Appendix A.1. If the high-pericenter mass or provenance depth changes by more than ~50%, the claimed direct pathway is not established by the current simulations.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim requires that the immediate distortion and elongation of a rotating target launches deep subsurface material onto high-pericenter orbits. But every SPH impact is run with a non-rotating target, and rotation is imposed only at the SPH-to-pkdgrav handoff (400 s) by adding to each particle the velocity difference between rotating and non-rotating targets at t=0 (Appendix A.1). This assumes that pre-impact rotation is dynamically decoupled from shock propagation, compaction, and early crater growth. In the real coupled problem, rotation changes the strain field during shock passage and the direction and speed of early ejecta, so the elongated shape and the 10-20 km provenance may be artifacts of the prescription rather than a direct consequence of the impact. The appendix validation (Figs. 9-10) compares only a 1 km target struck by a 10 m impactor, and only checks remnant c/a, spin period, and accumulated mass; it does not test the provenance-depth distribution, the high-pericenter debris mass, or the 100 km target / 3-18 km impactor regime where deformation is strong. Because Fig. 2 shows the elongation developing at 0.1 h, just at the handoff time, the very feature central to the mechanism is the one being approximated rather than simulated.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a direct formation pathway for satellites around large (D>100 km) asteroids: sub-catastrophic impacts into pre-rotating, porous, low-density targets. Using SPH impact simulations handed off to pkdgrav granular N-body reaccumulation and then to REBOUND stability mapping, the authors argue that impact-induced elongation of a fast-spinning primary launches deep subsurface material (10–20 km below the surface) onto eccentric orbits with pericenters large enough to avoid prompt re-impact. They connect this mechanism to the observed preference for fast-spinning, elongated, non-S-type primaries and to the lack of a strong correlation between satellite-hosting asteroids and large collisional families.","tokens_in":18922,"tokens_out":6935,"duration_ms":70741,"significance":"If correct, this is a valuable single-mechanism explanation for several observed demographic features of large-asteroid satellites, and it is one of the first collisional models to track both spin and shape continuously through impact and reaccumulation. The paper's strengths are its end-to-end numerical pipeline, the systematic impact-parameter table, the explicit medium-term stability maps, and the direct comparison with observed rotation-period and lightcurve-amplitude distributions. The central quantitative conclusion, however, rests on the post-hoc rotation-boost approximation, whose validation is limited and does not cover the most important outputs. The mechanism is plausible and the observational correlations are suggestive, but the load-bearing numerical shortcut needs stronger support before the central claim can be accepted.","major_comments":[{"comment":"The load-bearing step is the addition of pre-impact rotation only at the SPH-to-pkdgrav handoff. This procedure assumes that pre-impact rotation is dynamically decoupled from shock propagation, compaction, and early crater growth. Figure 2 shows that the elongation that drives the proposed mechanism is already developing at 0.1 h (360 s), essentially coincident with the nominal 400 s handoff time, so the very distortion central to the mechanism is not computed for a genuinely rotating target. The appendix validation (Figs. 9–10) tests only a 1 km target struck by a 10 m impactor and compares the remnant c/a, spin period, and accumulated mass; it does not compare the provenance-depth distribution, the mass of high-pericenter debris, or the 100 km / 3–18 km impactor regime in which the deformation is strong. Please provide SPH simulations with initially rotating 100-km targets (or a validated equivalence argument) and compare the q>1 a_pri debris mass and depth provenance; without this, the deep-subsurface high-pericenter debris could be an artifact of imposing rotation on an already shocked and displaced particle field.","section":"§2.1 and Appendix A.1 (Figs. 2, 9–10)"},{"comment":"The abstract states that 'the specific energy and resultant total mass loss in satellite-forming collisions are not constraining,' but §3 reports that the normalized satellite mass is correlated with the rotation-dependent normalized energy Q/Q*_RD, and Fig. 4d displays this trend. If the intended claim is that un-normalized Q or M_rem/M_target is not constraining while Q/Q*_RD is, the wording must be corrected; as written, the abstract contradicts the paper's own principal quantitative correlation and weakens the interpretation of the observed lack of family correlation.","section":"Abstract vs. §3 and Fig. 4d"},{"comment":"The claim of a 'sharp transition' in satellite incidence at post-impact spin periods of 10 h or less is stronger than the tabulated outcomes warrant. Table 2 contains fast-rotating cases with zero or very few high-pericenter satellites, e.g., Simulation 13 (P_rem = 4.22 h, 0 with q>1 a_pri) and Simulation 19 (P_rem = 4.45 h, 1). The trend in Fig. 4b is real but the narrative of a sharp observational match should be qualified to mass-weighted or per-impact measures, not raw incidence.","section":"§2.3 and Table 2 (Fig. 4b)"}],"minor_comments":[{"comment":"The column header repeats 'b/a' for the second axis ratio; the final axis ratio column should be 'c/a'.","section":"Table 2"},{"comment":"The main text gives the handoff time as 400 s after impact, while the Fig. 9 caption says the spin adjustment is applied 100 s after impact; these times should be reconciled.","section":"Appendix A.1 and §2.1"},{"comment":"The phrase 'from (38)' appears to be an unresolved citation; a reference should be supplied.","section":"§2.2"},{"comment":"The definition of Q*_RD is not given in the text; please provide the explicit functional form or a precise reference so the x-axis of Fig. 4d is reproducible.","section":"Fig. 4d and §3"},{"comment":"The text says impactors are 10–36 km in diameter, but Table 2 lists impactor radii of 3, 5, 7, 13, and 18 km, i.e., diameters of 6–36 km; the range should be corrected.","section":"§2.2 and Table 2"},{"comment":"The color description 'pericenter less than zero' for hyperbolic orbits is physically inaccurate; unbound orbits have positive pericenter distances, so the caption should instead say 'unbound (hyperbolic) trajectories'.","section":"Figs. 2, 11, 12"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for the journal and I see no citation or novelty concerns. My reservation is specifically the validation of the rotation-handoff approximation, which is load-bearing for the central mechanism; the requested rotating-target SPH tests are feasible and would settle whether the deep-provenance, high-pericenter debris is robust."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Your instinct to send this my way was right. The paper is a credible simulation study that for the first time tracks the spin and shape of the largest remnant through the reaccumulation phase, and it uses that to explain three observed features of the large-asteroid satellite population: the preference for fast-spinning, elongated primaries, the absence of S-type satellites, and the lack of correlation with large families. The mechanism—immediate distortion of a rotating target launching deep subsurface material onto high-pericenter orbits—is concrete and visible in the simulations. That is a real step beyond earlier impact models.\n\nWhat the paper does well: it is honest about what is not tested. The authors state explicitly that they only explore prograde rotation, spherical targets, a small set of impact geometries, and no debris-debris collisions. The REBOUND stability maps are a sensible way to connect the short-term reaccumulation to long-term survival. The comparison with observations is qualitative, but the observed correlations are clear enough that a qualitative match carries meaning.\n\nThe main soft spot is the post-hoc rotation prescription. The SPH impacts are all non-rotating; rotation is added by hand at the 400 s handoff. The appendix validation uses a 1 km target and a 10 m impactor, which is not the same regime as a 100 km target with multi-km impactors. Because the mechanism hinges on the distorted, elongated shape forming at around the handoff time, there is a real question whether a genuinely rotating target would behave the same way during shock propagation and early crater growth. This is not a disqualifying issue—the authors are transparent about the approximation, and the validation gives some confidence—but a referee should ask for a demonstration in the target regime, or at least an estimate of the error. Secondary issues: single realizations per geometry, no released code or data, and resolution limits around 1 km particles are also worth noting but are not disqualifying.\n\nBottom line: the paper is a solid contribution that deserves a serious referee. The central mechanism is plausible and the simulation suite supports it, but the handoff approximation leaves enough uncertainty that the paper should be accepted with a request for additional validation rather than taken as the last word. I would cite it if I worked on asteroid satellite formation.","headline":"First impact model to track post-impact spin and shape of the largest remnant, connecting impact geometry to observed satellite demographics; the post-hoc rotation prescription is the main caveat.","tokens_in":19460,"tokens_out":4570,"would_cite":true,"duration_ms":43729,"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":"This paper claims that sub-catastrophic impacts can form asteroid satellites directly, by stretching a fast-spinning target into an elongated shape that flings subsurface debris onto wide, stable orbits.","keywords":["asteroid satellites","satellite formation","sub-catastrophic impacts","asteroid elongation","fast rotation","porous asteroids","reaccumulation","main-belt asteroids"],"falsifier":"Run a shock-physics simulation of a 13 km impactor into a 100 km, 50-percent-porous target that is already rotating with a 4 or 6 hour period from time zero, without the handoff-time momentum boost, and compare the resultant largest remnant axis ratio, final spin, and depth of origin of debris with pericenter above the primary's long semi-axis; if the rotating-target run does not reproduce the elongated shape and deep-sourced high-pericenter debris, the direct pathway is an artifact of the post-hoc spin treatment.","tokens_in":18440,"feed_emoji":"☄️","tokens_out":10561,"duration_ms":94193,"temperature":0.7,"pith_summary":"This paper tries to show that the small moons around the largest asteroids can be made in a single impact, without needing later collisions or tides to arrange their orbits. The key is that a sub-catastrophic strike on an already fast-spinning, porous body stretches the body into an elongated shape, and that distortion flings subsurface debris onto orbits wide enough to miss the asteroid on the next pass. If true, the mechanism connects satellite formation to the observed spin and shape of the primary, explains why the impact energy that forms satellites is not linked to family size, and offers a reason no large S-type asteroid is known to host a moon. The paper reaches this conclusion by chaining shock-physics impact simulations, gravitational reaccumulation, and long-term orbital stability calculations.","feed_headline":"Impact-stretched asteroids can launch their own moons","feed_subtitle":"New simulations show subsurface debris tossed onto wide orbits when a strike elongates a rotating primary.","key_machinery":"The load-bearing mechanism is the immediate distortion of the rotating target into an elongated figure, which acts as a launch pad: debris shed from the stretched long axis is placed onto eccentric orbits with pericenter above the primary's longest semi-axis, bypassing the re-impact that ordinary ballistic ejecta suffer. The paper tracks this through a three-stage modeling chain: shock propagation through a porous, low-density target; handoff of the largest remnant's shape into a granular N-body code that lets it reaccumulate, reshape, and spin over 38 hours; and finally a stability map using the remnant's J2 and C22 gravity terms to identify debris that survives hundreds of orbits. The dynamically equivalent equal-volume ellipsoid of the remnant is the shape quantity that ties the simulations to observations of lightcurve amplitude and spin period.","core_discovery":"The paper's central claim is that the satellites around large main-belt asteroids can be made directly by sub-catastrophic impacts into rotating, porous parent bodies. When the impact is energetic enough to distort the target but not to disrupt it, the target is momentarily stretched into an elongated shape; material torn from the stretched long axis, mostly from 10 to 20 km below the surface, travels on eccentric orbits whose pericenters are already above the primary's long axis, so it does not crash back down. This converts what would otherwise be re-impacting ballistic ejecta into a stable orbiting population. Because the satellite-bearing outcome depends on the post-impact shape and spin rather than on total mass loss, the mechanism naturally explains why the observed primaries are fast rotators and elongated, why satellite presence does not track membership in large collision families, and possibly why no large S-type primaries have moons.","pith_inferences":["If the launch mechanism is generic, it should also operate when a fast-spinning small rubble pile is struck by a small projectile, offering a satellite-formation route for 10 to 100 km primaries that does not depend on YORP spin-up; the paper only notes the mass-shedding analogy for near-Earth asteroids.","The depth provenance result suggests a compositional test the paper does not run: material launched from 10 to 20 km depth could be less space-weathered or less porous than surface regolith, so the moons' spectra might differ subtly from their primaries' surfaces.","The abundance of temporary satellites in every run implies a selection step the paper leaves open; we infer that collisional relaxation into a disk, followed by re-accretion, is the most likely way to go from hundreds of clumps to the one-to-three moons typically observed.","The spin-boost shortcut could be checked without new observational data: rerunning a few 100 km impacts with rotation present from time zero in the shock code would confirm whether the elongation and 10 to 20 km provenance survive."],"forward_implications":["Asteroids with satellites should be preferentially fast-rotating and elongated, because both properties are produced by the same sub-catastrophic impact that launches the satellites.","The specific impact energy and total mass lost in a collision should not predict whether a satellite forms, which matches the observed absence of a strong correlation with membership in large asteroid families.","The absence of known satellites around large S-type asteroids can be explained by their higher density shrinking the stable orbital region around an elongated fast spinner, even if their shape and spin distributions resemble C-types.","Extremely elongated, fast-spinning bodies such as Kleopatra could acquire young satellites from a later small impact, because their existing shape and spin raise debris pericenters without needing a large family-forming event.","Temporary satellites formed this way start on orbits with pericenter above the primary's longest axis, so they should be stable for hundreds of orbits; observationally, known large-asteroid satellites sit at 3 to 14 primary radii with low eccentricity."],"supporting_citations":[{"why":"Provides the porous low-density asteroid shock model with compaction used for all impact simulations.","marker":"(Jutzi et al. 2019)"},{"why":"Supplies the shock-physics impact code used in the handoff and the rotation-dependent disruption threshold metric.","marker":"(Jutzi 2015)"},{"why":"Supplies the alpha-shape handoff that carries the largest remnant's irregular shape into the N-body phase.","marker":"(Ballouz et al. 2019)"},{"why":"Provides the gravitational N-body integrator used for reaccumulation.","marker":"(Richardson et al. 2000)"},{"why":"Adds soft-sphere granular interactions that let the remnant reshape and spin over 38 hours.","marker":"(Schwartz et al. 2012)"},{"why":"Provides the hierarchical search that identifies bound debris clumps as temporary satellites.","marker":"(Leinhardt & Richardson 2005)"},{"why":"Supplies the universal catastrophic-disruption law used to define Q_star_RD.","marker":"(Leinhardt & Stewart 2012)"},{"why":"Defines the rotation-dependent disruption threshold that puts the simulations on a common energy scale.","marker":"(Ballouz et al. 2014)"},{"why":"Supplies direct shape models of the largest asteroids that establish the elongation correlation among primaries.","marker":"(Vernazza et al. 2021)"},{"why":"Supplies the lightcurve database of rotation periods and amplitudes used for the spin and shape statistics.","marker":"(Warner et al. 2021)"}],"fun_headline_variants":["Impact-stretched asteroids fling their own moons","Asteroid moons born from spin and stretch, not crash size","Stretched rotating asteroids launch subsurface moons","How big impacts skip the family and still make moons","Fast-spinning elongated asteroids can spawn satellites"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing shortcut is to add the target's rotation only after the impact has already been simulated on a non-rotating body, and the test of that shortcut used a target one-hundredth the size, so the elongation and deep-sourced debris that make the pathway work could be artifacts of adding spin after the fact.","fun_headline_variants_meta":{"raw":{"variants":["Impact-stretched asteroids fling their own moons","Asteroid moons born from spin and stretch, not crash size","Stretched rotating asteroids launch subsurface moons","How big impacts skip the family and still make moons","Fast-spinning elongated asteroids can spawn satellites"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000145,"raw_usage":{"total_tokens":1164,"prompt_tokens":913,"completion_tokens":251,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":529,"completion_tokens_details":{"reasoning_tokens":178}},"tokens_in":529,"tokens_out":251,"duration_ms":3111,"temperature":1.0,"reasoning_tokens":178,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:54:23.527265+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run a shock-physics simulation of a 13 km impactor into a 100 km, 50-percent-porous target that is already rotating with a 4 or 6 hour period from time zero, without the handoff-time momentum boost, and compare the resultant largest remnant axis ratio, final spin, and depth of origin of debris with pericenter above the primary's long semi-axis; if the rotating-target run does not reproduce the elongated shape and deep-sourced high-pericenter debris, the direct pathway is an artifact of the post-hoc spin treatment.","supporting_citations":[{"cited_title":"J., Richardson, D","cited_arxiv_id":null,"evidence_quote":"Supplies the alpha-shape handoff that carries the largest remnant's irregular shape into the N-body phase."},{"cited_title":"M., & Richardson, D","cited_arxiv_id":null,"evidence_quote":"Provides the hierarchical search that identifies bound debris clumps as temporary satellites."},{"cited_title":"D., Harris, A","cited_arxiv_id":null,"evidence_quote":"Supplies the lightcurve database of rotation periods and amplitudes used for the spin and shape statistics."}],"review_version":1}