{"id":"9512964f-659f-4feb-b4c4-afb7c2e8d3de","arxiv_id":"2507.14645","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Cohesion increasingly controls ricochet, roll-out, and full-stop behavior in low-speed granular impacts as gravity decreases, and additive Froude-Bond scaling cannot predict the simulated trajectories.","lead":"Simulations of low-speed impacts into granular beds show that cohesion between grains matters far more under asteroid gravity than under Earth or Moon gravity. The paper shows that current scaling laws based on Froude and Bond numbers fail to capture this effect, pointing to a need for a new dimensionless parameter.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Figure 7's negative DRFT result uses cohesion-independent α and λ, so the claim that additive Bond scaling is incongruent—and the call for a new dimensionless number—is not yet demonstrated.","rationale":"Read in good faith, the paper's strongest supported result is the behavior-map trend: cohesion reshapes impact outcomes far more at Bennu gravity than at Earth or Moon. This is visible in Fig. 5 and does not depend on the DRFT comparison. The reader's conditional verdict and weakest assumption coincide with the key issue: the negative DRFT result is not a test of the full additive model. The manuscript even includes the needed calibration data (Table 3) but does not close the loop by applying those α and λ values to the trajectory calculation. If closing that loop restores Bond-number collapse, the abstract's conclusion that a new dimensionless parameter is required would be an overreach; if it fails, the conclusion would be substantially strengthened. Because the paper presents no code or out-of-sample validation and the fit forms in Eqs. 6–9 are flexible, the appropriate verdict remains conditional pending this specific check. The authors' explicit acknowledgment of the uncalibrated coefficients is the basis for the requested test, and no ad hominem is intended.","tokens_in":12241,"tokens_out":5912,"duration_ms":75766,"concrete_test":"Use the measured α_x, α_z, λ from Table 3 for cohesive Bennu (γs=0.2) and Earth (γs=1.0), with interpolated values for Moon, and integrate the DRFT equations (Eqs. 13–15) for the same Froude number and angle as Fig. 7, keeping the cohesion coefficient c scaled by Bond number. Compare the predicted behavior (full-stop/roll-out/ricochet) with the DEM outcomes. If the calibrated additive model reproduces the DEM divergence, the additive Bond-number scaling is viable and the new-dimensionless-number conclusion is unsupported; if it still does not, the paper's negative claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central negative claim is that DEM trajectories are incongruent with an additive DRFT model using Froude/Bond scaling, motivating a new dimensionless parameter. That claim depends on the comparison in Fig. 7, but the derivation in Eqs. 13–15 assumes α and λ are unaffected by cohesion, and the authors state in §3.2 that 'neither of the experimental value α(β, γ) or the shape factor λ is calibrated for this example.' The later plate-drag measurements (Table 3) directly show that α is strongly cohesion- and gravity-dependent: for the cohesive Bennu bed α_x rises 441%. Thus Fig. 7 tests only a constant-coefficient additive model, not the additive model itself. The authors notice this and suggest coefficients must be modified, but they never integrate the measured α, λ back into Eqs. 13–15. Until that is done, the 'incongruence' can be read as a calibration artifact, and the proposed new dimensionless number lacks its main evidence. The qualitative low-gravity enhancement of cohesion is unaffected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript studies low-speed oblique impacts of a disc into a 2D polydisperse granular bed using DEM simulations under Earth, Moon, and Bennu gravity, with four levels of cohesion. Impact outcomes are classified as full-stop, roll-out, or ricochet and presented as behavior maps (Figure 5). The authors introduce a modified Froude number π5 with a cohesion term (Eq. 4) and fitted separation lines (Eqs. 6-9), and they test an additive cohesive extension of Dynamic Resistive Force Theory (Eqs. 10-15). The central positive finding is that cohesion has a larger effect on impact behavior at lower gravity, and the central negative finding is that DEM trajectories do not collapse under Bond-number scaling of cohesion, leading the authors to call for a new dimensionless parameter.","tokens_in":12457,"tokens_out":4284,"duration_ms":49861,"significance":"If the qualitative finding holds, it is relevant for spacecraft-surface interactions on small bodies such as Bennu and Ryugu, where cohesion is significant relative to gravity. The study has notable strengths: a large simulation campaign (1716 cases), calibration of bulk tensile strength via tensile tests (Eq. 2), and direct measurement of the DRFT coefficients α and λ for cohesive and cohesionless beds (Table 3). The negative result against a constant-coefficient additive Bond scaling is potentially useful, but as discussed below it does not yet establish that additive scaling fails once the coefficients are recalibrated for cohesion. The quantitative fitting apparatus in Eqs. 4-9 is currently descriptive rather than predictive, since the constants are fit to the same data they are used to describe.","major_comments":[{"comment":"The central negative claim—that DEM results are incongruent with additive DRFT and hence a new dimensionless parameter is needed—is not yet demonstrated, because the comparison fixes α and λ at their cohesionless values. The authors state in §3.2 that neither α(β,γ) nor λ is calibrated for this example, and Table 3 shows that α_x increases by 441% for the cohesive Bennu bed while α_z changes by 611%. Figure 7 therefore tests only a constant-coefficient additive model, not the additive model itself. The authors' own plate-drag measurements provide the missing coefficients; inserting them into Eqs. 13-15 would either restore or robustly falsify the additive Bond scaling. As written, the incongruence can be read as a calibration artifact, and the appeal for a new dimensionless number rests on this unresolved point.","section":"§3.2, Eqs. (13)-(15), Fig. 7"},{"comment":"The modified Froude number π5 is introduced with a fitted dimensionless constant c_g=100,000, and the separation-line coefficients A_rc, B_rc, A_ro, and B_ro are fitted to the same behavior maps they are used to describe. No out-of-sample test or uncertainty estimate is provided. The qualitative observation that cohesion shifts behavior more strongly at Bennu gravity does not depend on these fits, but the claim that π5 'successfully matched' the simulations is weakened by the fitting-to-the-data circularity. I ask the authors to either validate the scaling on held-out conditions or present Eqs. 4-9 as descriptive fits with appropriate caveats.","section":"§3.1, Eq. (4) and Eqs. (6)-(9)"},{"comment":"The behavior maps are computed as averages over only three impactor locations, and no measure of run-to-run variability or statistical uncertainty is reported. The Earth-row anomaly at 40° is explicitly attributed to one specific impact point, indicating sensitivity to local packing. Without error bars or a variance map, the strength of the qualitative low-gravity cohesion trend and the claimed critical-angle shifts (e.g., 50° to 60° on Bennu) are hard to evaluate. Please report per-location outcomes or a variance measure for at least the key comparisons.","section":"§2, Fig. 5"}],"minor_comments":[{"comment":"The phrase 'granular granular surfaces' should be corrected to 'granular surfaces'.","section":"Abstract and §1"},{"comment":"'as a ration of ρga2/γs' should read 'as a ratio of ρga2/γs'.","section":"§1"},{"comment":"The text reports packing fractions 0.83, 0.82, and 0.82 in one place, while Figure 5 reports 0.812 and 0.811; please reconcile these values or clarify which quantity is being reported.","section":"§2 and Fig. 5"},{"comment":"Writing 'Bond number of ∞' for the cohesionless case is notationally confusing; consider writing 'no cohesion (Bo → ∞)' instead.","section":"Fig. 7 caption"},{"comment":"No data or code availability statement is provided; adding one would improve reproducibility of the DEM workflow and fitting procedures.","section":"Availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is within scope for a granular-physics journal and the simulation campaign is substantial. My main concern is that the headline negative result is underdetermined by the current comparison: Figure 7 tests only a constant-coefficient additive model, despite Table 3 showing strong cohesion- and gravity-dependence of the coefficients. I would be comfortable with acceptance after the authors either integrate the measured α and λ into the DRFT comparison or substantially soften the claim about a new dimensionless parameter. The fitted π5 scaling should also be framed as descriptive unless out-of-sample tests are added."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know about this paper. First, the qualitative result—that cohesion has a stronger effect on impact outcomes at low gravity—is convincing and well-supported by the behavior maps. That is the paper's real contribution. Second, the negative claim against additive Bond-number scaling is not yet demonstrated; the DRFT comparison in Fig. 7 uses cohesion-independent α and λ, and the authors' own Table 3 shows α changes by hundreds of percent under cohesion at Bennu gravity. So the 'incongruence' is plausibly a calibration artifact, and the call for a new dimensionless number is premature.\n\nWhat the paper does well: it is a systematic simulation campaign—1716 DEM runs across three gravity levels and four cohesion values—with a sensible tensile-test calibration to convert grain-scale surface energy into a bulk cohesive strength. The cohesionless baseline matches the experimental separation lines of Wright et al. The behavior maps (Fig. 5) show a clear, monotonic expansion of the ricochet/roll-out regions with cohesion at Bennu gravity, and almost no effect at Earth gravity. That trend is robust and physically plausible: cohesion becomes relevant when it competes with gravitational pressure.\n\nThe soft spots are real but not fatal to the qualitative take. The quantitative scaling is post-hoc: the modified Froude number π5 uses a fitted constant c_g=100,000, and the separation-line equations have multiple fitted coefficients with no out-of-sample validation. Error bars are missing despite only three repeated impactor locations. No code or data is provided, which limits reproducibility but is not unusual for a DEM study.\n\nThe bigger issue is the DRFT comparison. Equations 13–15 assume α and λ are unchanged by cohesion, and the text admits they are not calibrated for this example. Table 3 then shows that α_x increases by 441% under Bennu gravity with cohesion. That undercuts the conclusion that additive scaling 'fails'—it only shows that a constant-coefficient additive model fails. The authors could reasonably argue that this motivates a new parameter, but they have not shown that recalibrating α and λ would not restore the collapse. Until that is tested, the negative result is suggestive, not conclusive.\n\nWho is this for? The planetary science and impact-engineering community will find the behavior maps and the qualitative gravity-cohesion trend useful. The paper deserves a serious referee, but it needs a major revision: either recalibrate α and λ for cohesive beds and retest the scaling, or clearly frame the negative result as a limitation of the current DRFT parametrization rather than evidence for a new dimensionless number. I would send it to review and let the referees push on that point.","headline":"A solid DEM study showing cohesion matters more at low gravity, but the negative result against Bond scaling is undercut by uncalibrated DRFT coefficients.","tokens_in":13015,"tokens_out":2469,"would_cite":true,"duration_ms":27513,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["45.70.-n"],"model":"deepseek-v4-flash","headline":"Grain cohesion dominates low-speed impact outcomes on low-gravity granular surfaces, and the standard additive Froude–Bond scaling fails to predict it.","keywords":["granular impact","cohesion","discrete element method","dynamic resistive force theory","Bond number","Froude number","low gravity","asteroid regolith"],"falsifier":"Run the same DEM protocol but measure $\\alpha(\\beta,\\gamma)$ and $\\lambda$ by plate-drag tests at each cohesion level and gravity condition, insert those recalibrated values into the additive dynamic resistive force theory, and check whether trajectories at equal Bond number collapse across Earth, Moon, and Bennu gravity. If they collapse, the additive Bond-number scaling is not falsified; if they still diverge, the claimed cohesion–friction coupling is supported.","tokens_in":11988,"feed_emoji":"☄️","tokens_out":8836,"duration_ms":97117,"temperature":0.7,"pith_summary":"Grain-by-grain simulations of a disc striking a loose granular bed across Earth, Moon, and Bennu gravity show that inter-particle cohesion changes the outcome of a low-speed oblique impact—full stop, roll-out, or ricochet—far more strongly when gravity is weak. The paper establishes that cohesion acts through its ratio to gravity, $\\gamma_s/g$, rather than as an independent additive stress, and that behavior maps can be fitted by a modified Froude number. It then tests an analytic extension of dynamic resistive force theory in which a cohesive traction is simply added to frictional and inertial forces and scaled by Bond number; the simulations do not follow that prediction. The authors conclude that cohesion and friction are coupled, so a new dimensionless parameter is needed. If the paper is right, impact trajectories on asteroid surfaces carry a direct signal of surface cohesion, and existing scaling laws for low-speed impacts are incomplete.","feed_headline":"At asteroid gravity, grain cohesion controls the projectile's fate","feed_subtitle":"Simulations show projectile fate on Bennu-like surfaces is set by grain stickiness, and additive scaling laws break down.","key_machinery":"Two complementary tools carry the argument. The first is discrete element method simulation of roughly 6900 polydisperse soft discs in a 2D bed, with a Hertzian normal repulsion, a JKR/DMT cohesive traction proportional to surface energy density $\\gamma_s$, and a damping term; impacts are classified as full-stop, roll-out, or ricochet and summarized as color-coded behavior maps over Froude number and impact angle. The paper converts grain-scale $\\gamma_s$ into a bulk tensile strength $\\sigma$ by simulated tensile tests, giving the four cohesion levels used in the maps. The second tool is dynamic resistive force theory, in which the force on the intruding disc is written as a quasi-static traction $\\alpha(\\beta,\\gamma)H(-\\bar{z})|\\bar{z}|$ plus an inertial traction $\\lambda \\rho |V|^2 \\hat{n}$; the paper's proposed extension adds a cohesive traction $c\\hat{n}$ and predicts trajectory collapse across gravities when $c$ is scaled by $g/g^*$, i.e., by Bond number. The mismatch between that analytic prediction and the DEM trajectories is the evidence for the paper's call for new dimensionless parameters.","core_discovery":"The paper's central positive claim is that cohesion becomes a dominant control on low-speed oblique impact outcomes as gravity decreases. Under Earth and Moon gravity, the boundaries between ricochet, roll-out, and full stop barely move as bulk cohesive strength rises from 0 to 26.9 Pa, but under Bennu gravity those boundaries expand sharply: the ricochet and roll-out regions grow, and at the highest cohesion the full-stop region disappears from the studied range. The behavior is organized by a modified Froude number, $\\pi_5=\\rho v^2 a/(\\rho g a^2 + c_g \\gamma_s)$, which reduces to the ordinary Froude number when $\\gamma_s=0$. The paper's central negative claim is that dynamic resistive force theory, extended by adding a cohesive traction and scaling it across gravities with Bond number, fails to reproduce the simulated trajectories: with Bond-number-scaled cohesion, Earth and Moon runs ricochet while the Bennu run rolls out. The authors attribute this failure to cohesion acting as an inter-granular pressure that raises friction, so that cohesion cannot be added independently to frictional and inertial resistance.","pith_inferences":["The paper's negative result would not rule out an additive model whose coefficients are re-fitted per cohesion level; the new dimensionless parameter is one possible repair, not the only one.","A direct testable consequence is that for any regolith there should be a crossover gravity where the full-stop/roll-out boundary starts shifting with cohesion; locating that crossover experimentally would confirm the maps.","Because the modified Froude number uses a single empirical constant $c_g$, the same functional form could be transferred to 3D beds and natural polydisperse regolith if $c_g$ is rescaled by particle and impactor radii.","The local surface void seen in the high-cohesion Bennu runs suggests cohesive beds can fail internally rather than at the impact point, which would alter how crater dimensions are interpreted on asteroid surfaces."],"forward_implications":["At Bennu-level gravity, a given cohesive strength changes the projectile from full stop to roll-out to ricochet at low Froude numbers, so the observed trajectory can serve as a measure of the surface's cohesive strength.","Future impact scaling laws must weight cohesion by gravity, e.g. through $\\gamma_s/g$, rather than treating cohesion as an environment-independent additive stress.","Dynamic resistive force theory requires cohesion-dependent $\\alpha(\\beta,\\gamma)$ and $\\lambda$; a constant additive cohesive traction is insufficient.","Design of landing and sampling operations on rubble-pile asteroids should treat even weak surface cohesion as a dominant factor at those gravities.","The fitted separation lines based on the modified Froude number $\\pi_5$ give a practical engineering criterion for classifying impact outcomes on cohesive regolith."],"supporting_citations":[{"why":"Supplies the experimental low-velocity grazing-impact behavior categories (full-stop, roll-out, ricochet) and the separation-curve form the DEM maps are compared against.","marker":"[4]"},{"why":"Establishes the Froude-number scaling for impact trajectories in cohesionless granular media across gravitational conditions that this paper extends to cohesive beds.","marker":"[25]"},{"why":"Provides the dynamic resistive force theory (DRFT) with quasi-static and inertial tractions that the paper augments with a cohesive traction.","marker":"[40]"},{"why":"Prior extension of granular RFT to cohesive powder-scale media; the source of the expectation that cohesion raises $\\alpha(\\beta,\\gamma)$.","marker":"[41]"},{"why":"Relates tensile strength to grain size for regolith, used to convert grain-level surface energy into bulk cohesive stress.","marker":"[30]"},{"why":"JKR contact model used for the surface-energy cohesive force in the DEM normal contact law.","marker":"[34]"},{"why":"DMT adhesion contact model contributing the cohesive force term in the DEM force law.","marker":"[24]"},{"why":"LAMMPS is the simulation engine in which the DEM beds and impacts are run.","marker":"[16]"}],"fun_headline_variants":["Cohesion dominates impact outcomes at asteroid gravity","Low gravity makes grain stickiness the key impact factor","Additive scaling fails for cohesive impacts on asteroids","Grain cohesion alters projectile fate only at low gravity","Stickiness controls ricochet and roll-out on small worlds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The negative result rests on leaving the two coefficients of the resistive-force model, $\\alpha(\\beta,\\gamma)$ and $\\lambda$, at their cohesionless values instead of recalibrating them for cohesive beds; if recalibration restores the trajectories at equal Bond number, the paper's conclusion that a fundamentally new dimensionless parameter is needed would not follow.","fun_headline_variants_meta":{"raw":{"variants":["Cohesion dominates impact outcomes at asteroid gravity","Low gravity makes grain stickiness the key impact factor","Additive scaling fails for cohesive impacts on asteroids","Grain cohesion alters projectile fate only at low gravity","Stickiness controls ricochet and roll-out on small worlds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000235,"raw_usage":{"total_tokens":1535,"prompt_tokens":1014,"completion_tokens":521,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":630,"completion_tokens_details":{"reasoning_tokens":445}},"tokens_in":630,"tokens_out":521,"duration_ms":6687,"temperature":1.0,"reasoning_tokens":445,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:51:27.295012+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same DEM protocol but measure $\\alpha(\\beta,\\gamma)$ and $\\lambda$ by plate-drag tests at each cohesion level and gravity condition, insert those recalibrated values into the additive dynamic resistive force theory, and check whether trajectories at equal Bond number collapse across Earth, Moon, and Bennu gravity. If they collapse, the additive Bond-number scaling is not falsified; if they still diverge, the claimed cohesion–friction coupling is supported.","supporting_citations":[{"cited_title":"Quillen, Juliana South, Randal C","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental low-velocity grazing-impact behavior categories (full-stop, roll-out, ricochet) and the separation-curve form the DEM maps are compared against."},{"cited_title":"Froude number scaling unifies impact trajectories into granular media across gravitational conditions","cited_arxiv_id":"2307.10998","evidence_quote":"Establishes the Froude-number scaling for impact trajectories in cohesionless granular media across gravitational conditions that this paper extends to cohesive beds."},{"cited_title":"Surprising simplicity in the modeling of dynamic granular intrusion.Science Advances, 7(17):eabe0631, 2021","cited_arxiv_id":null,"evidence_quote":"Provides the dynamic resistive force theory (DRFT) with quasi-static and inertial tractions that the paper augments with a cohesive traction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Prior extension of granular RFT to cohesive powder-scale media; the source of the expectation that cohesion raises $\\alpha(\\beta,\\gamma)$."},{"cited_title":"S´ anchez and D","cited_arxiv_id":null,"evidence_quote":"Relates tensile strength to grain size for regolith, used to convert grain-level surface energy into bulk cohesive stress."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"JKR contact model used for the surface-energy cohesive force in the DEM normal contact law."},{"cited_title":"Effect of contact deformations on the adhesion of particles.Journal of Colloid and Interface Science, 53(2):314–326, 1975","cited_arxiv_id":null,"evidence_quote":"DMT adhesion contact model contributing the cohesive force term in the DEM force law."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"LAMMPS is the simulation engine in which the DEM beds and impacts are run."}],"review_version":1}