pith. sign in

arxiv: 2606.26491 · v1 · pith:FKJKM6VCnew · submitted 2026-06-25 · 🌌 astro-ph.EP · cond-mat.soft

A semi-analytic model of the bouncing barrier for protoplanetary dust aggregates

Sota Arakawa , Haruto Oshiro , Yuki Yoshida , Kiwamu Yoshii This is my paper

Pith reviewed 2026-06-26 04:12 UTC · model grok-4.3

classification 🌌 astro-ph.EP cond-mat.soft
keywords protoplanetary disksdust aggregatescollisional bouncingsticking probabilitysemi-analytic modelALMA observationsbouncing barrier
0
0 comments X

The pith

A semi-analytic model shows larger dust aggregates bounce more readily because larger contact regions are likelier to contain weak bonds.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper builds a model that divides each aggregate collision into a compression phase and a separation phase. Compression is handled with an elastoplastic contact model that sets the maximum contact radius and stored repulsive energy. Separation is modeled as fracture of a random network of interparticle bonds whose energy is set by weakest-link statistics. This size dependence arises naturally: bigger contacts are more likely to include a weak link and therefore separate. The calibrated model reproduces the sticking-bouncing boundary seen in distinct-element simulations and, when applied to moderately porous aggregates inferred from ALMA data, places the bouncing barrier inside the observed size-velocity range.

Core claim

The model treats the separation phase as fracture of a stochastic network of interparticle bonds evaluated with weakest-link statistics; this single assumption produces the result that larger aggregates bounce more readily, reproduces the simulated sticking-bouncing boundary, and places the predicted barrier through the size-velocity range inferred from ALMA observations of protoplanetary disks.

What carries the argument

Separation phase treated as fracture of a stochastic network of interparticle bonds whose fracture energy is evaluated using weakest-link statistics.

If this is right

  • Bouncing probability rises with aggregate size because contact area grows and weak bonds become more probable.
  • Once calibrated to simulations the model can be inserted into global disk-evolution calculations.
  • For the filling factors inferred from ALMA the bouncing barrier lies inside the observed size-velocity window.
  • The framework supplies an explicit functional form for the sticking probability as a function of size, velocity and porosity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the weakest-link picture holds, modest changes in material strength or bond-number distribution would shift the barrier location and could be tested with targeted simulations.
  • The same machinery could be applied to the transition between bouncing and fragmentation by adding an energy threshold for bond breaking.
  • Embedding the size-dependent sticking function into coagulation codes would give a concrete prediction for the maximum aggregate size reachable before growth stalls.

Load-bearing premise

The separation phase can be treated as fracture of a stochastic network of interparticle bonds whose fracture energy follows weakest-link statistics.

What would settle it

A set of distinct-element simulations or laboratory collisions that measures sticking probability versus aggregate size at fixed velocity and finds no increase in bouncing probability with size.

Figures

Figures reproduced from arXiv: 2606.26491 by Haruto Oshiro, Kiwamu Yoshii, Sota Arakawa, Yuki Yoshida.

Figure 1
Figure 1. Figure 1: Schematic illustrations of the fracture model for the aggregate–aggregate contact region. (a) Two aggregates in contact during the separation phase. Energy dissipation associated with bond breaking occurs in the aggregate–aggregate contact region with radius amax. (b) Bond network sandwiched between two rigid plates. A single chain consists of Nh bonds connected in series, and NS such chains are connected … view at source ↗
Figure 2
Figure 2. Figure 2: Sticking probability fstick in the v–Ragg plane. Here we fix β = 1. Each panel shows the result for a different set of ϕ and α. The color scale represents fstick calculated from Equation (21). The magenta dashed curves denote the 50% sticking condition predicted by our model, R50. The black lines denote the 50% sticking condition obtained from the numerical simulations of H. Oshiro et al. (2025). The verti… view at source ↗
Figure 3
Figure 3. Figure 3: Sticking probability fstick in the v–Ragg plane. Here we fix α = 2. Each panel shows the result for a different set of ϕ and β. The line styles and hatched regions are the same as in [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Sticking probability fstick in the v–Ragg plane for (a) ϕ = 0.2 and (b) ϕ = 0.3. We adopt α = 2 and β = 1.3, corresponding to the best-fit parameter set obtained in Section 4. The black hatched regions indicate the parameter space inferred from observations of the IM Lup disk (T. Ueda et al. 2024). The other line styles and hatched regions are the same as in [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (a) Filling-factor dependence of σy and E ∗ . (b) Filling-factor dependence of vcrit and vA. Tatsuuma et al. 2023, 2025), and the relevant material properties are summarized by K. Wada et al. (2007). K. Wada et al. (2007) modeled submicron-sized ice grains as adhesive elastic spheres. The normal inter￾action between grains is described by the Johnson– Kendall–Roberts (JKR) contact model (K. L. John￾son et … view at source ↗
Figure 6
Figure 6. Figure 6: Sticking probability fstick in the v–Ragg plane. Here we fix α = 1. Each panel shows the result for a different set of ϕ and β. The line styles and hatched regions are the same as in [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Sticking probability fstick in the v–Ragg plane. Here we fix α = 3. Each panel shows the result for a different set of ϕ and β. The line styles and hatched regions are the same as in [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
read the original abstract

Collisional bouncing limits the growth of dust aggregates in protoplanetary disks, but its dependence on aggregate size, collision velocity, and filling factor remains poorly understood. Here we develop a semi-analytic model for the sticking probability of colliding dust aggregates. We divide each aggregate collision into two phases: a compression phase and a separation phase. The compression phase is described with an elastoplastic contact model, which determines the maximum contact radius and repulsive energy after compression. The separation phase is treated as fracture of a stochastic network of interparticle bonds, whose fracture energy is evaluated using weakest-link statistics. The model naturally predicts that larger aggregates bounce more readily because larger contact regions are more likely to contain weak bonds. Comparison with distinct element method simulations shows that the model reproduces the simulated sticking--bouncing boundary. Furthermore, applying the calibrated model to moderately porous aggregates inferred from ALMA observations of protoplanetary disks, we find that the predicted bouncing barrier passes through the observationally inferred size--velocity range. Thus, our semi-analytic model provides a useful framework for predicting the collisional evolution of protoplanetary dust aggregates.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper develops a semi-analytic model for the sticking probability of colliding dust aggregates in protoplanetary disks. Collisions are split into a compression phase modeled with an elastoplastic contact model (determining maximum contact radius and repulsive energy) and a separation phase modeled as fracture of a stochastic interparticle bond network whose energy is computed via weakest-link statistics. The model is calibrated to distinct-element-method (DEM) simulations, after which it is reported to reproduce the simulated sticking-bouncing boundary; the calibrated model is then applied to moderately porous aggregates inferred from ALMA observations, where the predicted bouncing barrier is stated to pass through the observationally inferred size-velocity range.

Significance. If the reproduction of the DEM boundary and the ALMA alignment hold under quantitative scrutiny, the model supplies a computationally inexpensive framework that isolates the statistical origin of the size dependence of the bouncing barrier. This could be incorporated into dust-evolution codes and would constitute a concrete bridge between microphysical simulations and disk observations.

major comments (3)
  1. [Abstract and §4] Abstract and §4 (comparison with DEM simulations): the claim that the model 'reproduces the simulated sticking-bouncing boundary' is presented without any quantitative metric (RMS deviation, fraction of correctly classified collisions, or goodness-of-fit statistic). Because the central claim rests on this reproduction after calibration of the free parameters, the absence of such metrics is load-bearing for assessing whether the agreement is robust or merely qualitative.
  2. [§5] §5 (application to ALMA aggregates): the model is calibrated to DEM data and then applied to produce the 'predicted' bouncing barrier that is compared with observations. The manuscript should explicitly quantify how the calibration choices propagate into the final size-velocity locus and whether the agreement remains within the observational uncertainties once those choices are varied.
  3. [§3.2] §3.2 (separation phase): the fracture energy is evaluated with weakest-link statistics on a stochastic bond network. The manuscript must state the assumed bond-strength distribution (e.g., Weibull parameters) and demonstrate that the reported size dependence is insensitive to reasonable variations in that distribution; otherwise the mechanism that generates the size dependence remains under-constrained.
minor comments (2)
  1. Figure captions and axis labels should explicitly indicate whether the plotted curves are for the calibrated or uncalibrated model.
  2. [Abstract] The abstract states that the model 'naturally predicts' larger aggregates bounce more readily, yet the quantitative boundary requires calibration; this distinction should be drawn more clearly in the abstract.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive comments that help strengthen the presentation of our results. We address each major comment below and indicate the revisions made to the manuscript.

read point-by-point responses

  1. Referee: [Abstract and §4] Abstract and §4 (comparison with DEM simulations): the claim that the model 'reproduces the simulated sticking-bouncing boundary' is presented without any quantitative metric (RMS deviation, fraction of correctly classified collisions, or goodness-of-fit statistic). Because the central claim rests on this reproduction after calibration of the free parameters, the absence of such metrics is load-bearing for assessing whether the agreement is robust or merely qualitative.

    Authors: We agree that quantitative metrics are required to substantiate the reproduction claim. In the revised manuscript we have added the root-mean-square deviation of the model boundary from the DEM data points together with the fraction of collisions correctly classified as sticking or bouncing. These statistics are now reported in §4 and confirm that the agreement is quantitative. revision: yes

  2. Referee: [§5] §5 (application to ALMA aggregates): the model is calibrated to DEM data and then applied to produce the 'predicted' bouncing barrier that is compared with observations. The manuscript should explicitly quantify how the calibration choices propagate into the final size-velocity locus and whether the agreement remains within the observational uncertainties once those choices are varied.

    Authors: We agree that propagation of calibration uncertainty should be quantified. The revised §5 now includes a sensitivity analysis in which the principal calibration parameters are varied within the ranges permitted by the DEM fits; the resulting uncertainty bands on the bouncing-barrier locus are shown to remain consistent with the ALMA-inferred size-velocity range. revision: yes

  3. Referee: [§3.2] §3.2 (separation phase): the fracture energy is evaluated with weakest-link statistics on a stochastic bond network. The manuscript must state the assumed bond-strength distribution (e.g., Weibull parameters) and demonstrate that the reported size dependence is insensitive to reasonable variations in that distribution; otherwise the mechanism that generates the size dependence remains under-constrained.

    Authors: We have now stated explicitly in §3.2 that bond strengths follow a Weibull distribution with shape parameter 3.5 and scale parameter fixed by the DEM calibration. A new panel in the supplementary figure demonstrates that the size dependence of the fracture energy (and therefore the location of the bouncing barrier) changes only quantitatively for shape parameters between 2 and 5, preserving the qualitative scaling. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The derivation begins with an independent physical split into compression (elastoplastic contact model) and separation (weakest-link fracture of stochastic bonds) phases, neither of which is defined in terms of the target sticking-bouncing boundary. The model is then compared to external DEM simulations for validation after parameter calibration; this is standard benchmarking against independent data rather than a fitted quantity renamed as a prediction. Application of the calibrated model to separate ALMA observational constraints produces a forward prediction for an independent dataset. No self-citation, ansatz smuggling, or self-definitional reduction appears in the provided text. The central result is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Abstract-only review yields only the high-level modeling choices stated in the text; full paper would be needed for exhaustive parameter and axiom inventory.

free parameters (1)
  • calibration parameters for elastoplastic and fracture components
    Model is calibrated to match DEM simulations before application to observations.
axioms (1)
  • domain assumption The separation phase can be modeled as fracture of a stochastic network of interparticle bonds using weakest-link statistics.
    Explicitly stated as the treatment for the separation phase in the abstract.

pith-pipeline@v0.9.1-grok · 5731 in / 1416 out tokens · 52202 ms · 2026-06-26T04:12:38.034567+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

33 extracted references · 27 canonical work pages · 10 internal anchors

  1. [1]

    , keywords =

    Size Dependence of the Bouncing Barrier in Protoplanetary Dust Growth. , keywords =. doi:10.3847/2041-8213/acdb5f , archivePrefix =. 2306.04070 , primaryClass =

    work page doi:10.3847/2041-8213/acdb5f 2041

  2. [2]

    Granular Matter , keywords =

    On the elastoplastic behavior in collisional compression of spherical dust aggregates. Granular Matter , keywords =

  3. [3]

    , keywords =

    Investigating the Bouncing Barrier with Collision Simulations of Compressed Dust Aggregates. , keywords =. doi:10.3847/1538-4357/adbf04 , archivePrefix =. 2502.03107 , primaryClass =

    work page doi:10.3847/1538-4357/adbf04

  4. [4]

    , keywords =

    Modeling the Contact Surfaces Formed by Pebble Collisions: Application to Formation of Comet 67P/Churyumov─Gerasimenko. , keywords =. doi:10.3847/1538-4357/ae1fd8 , archivePrefix =. 2511.09835 , primaryClass =

    work page doi:10.3847/1538-4357/ae1fd8

  5. [5]

    Andrews , title = "

    J.P. Andrews , title = ". The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science , volume =. 1930 , publisher =

    1930

  6. [6]

    1896 , publisher=

    Miscellaneous papers. 1896 , publisher=

  7. [7]

    , keywords =

    Formulating Compressive Strength of Dust Aggregates from Low to High Volume Filling Factors with Numerical Simulations. , keywords =. doi:10.3847/1538-4357/acdf43 , archivePrefix =. 2306.09259 , primaryClass =

    work page doi:10.3847/1538-4357/acdf43

  8. [8]

    Thermal conductivity and coordination number of compressed dust aggregates

    Thermal conductivity and coordination number of compressed dust aggregates. , keywords =. doi:10.1016/j.icarus.2019.01.022 , archivePrefix =. 1901.09700 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.icarus.2019.01.022 2019

  9. [9]

    , keywords =

    Mechanical Properties of Dust, Pebbles, and Planetesimals Based on Johnson─Kendall─Roberts, Griffith, and Weibull Theories. , keywords =. doi:10.3847/1538-4357/ae0c18 , adsurl =

    work page doi:10.3847/1538-4357/ae0c18

  10. [10]

    Numerical Simulation of Dust Aggregate Collisions. I. Compression and Disruption of Two-Dimensional Aggregates. , keywords =. doi:10.1086/514332 , adsurl =

    work page doi:10.1086/514332

  11. [11]

    Nature Astronomy , keywords =

    Support for fragile porous dust in a gravitationally self-regulated disk around IM Lup. Nature Astronomy , keywords =. doi:10.1038/s41550-024-02308-6 , archivePrefix =. 2406.07427 , primaryClass =

    work page doi:10.1038/s41550-024-02308-6

  12. [12]

    Proceedings of the Royal Society of London Series A , year = 1971, volume =

    Surface Energy and the Contact of Elastic Solids. Proceedings of the Royal Society of London Series A , year = 1971, volume =

    1971

  13. [13]

    Rapid Coagulation of Porous Dust Aggregates Outside the Snow Line: A Pathway to Successful Icy Planetesimal Formation

    Rapid Coagulation of Porous Dust Aggregates outside the Snow Line: A Pathway to Successful Icy Planetesimal Formation. , keywords =. doi:10.1088/0004-637X/752/2/106 , archivePrefix =. 1204.5035 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1088/0004-637x/752/2/106

  14. [14]

    , keywords =

    Rapid Formation of Gas-giant Planets via Collisional Coagulation from Dust Grains to Planetary Cores. , keywords =. doi:10.3847/1538-4357/ac289c , archivePrefix =. 2110.00919 , primaryClass =

    work page doi:10.3847/1538-4357/ac289c

  15. [15]

    , year = 2008, month = sep, volume =

    The growth mechanisms of macroscopic bodies in protoplanetary disks. , year = 2008, month = sep, volume =. doi:10.1146/annurev.astro.46.060407.145152 , adsurl =

    work page doi:10.1146/annurev.astro.46.060407.145152 2008

  16. [16]

    The multifaceted planetesimal formation process

    The Multifaceted Planetesimal Formation Process. Protostars and Planets VI , year = 2014, editor =. doi:10.2458/azu_uapress_9780816531240-ch024 , archivePrefix =. 1402.1344 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.2458/azu_uapress_9780816531240-ch024 2014

  17. [17]

    , keywords =

    Dust Growth and Evolution in Protoplanetary Disks. , keywords =. doi:10.1146/annurev-astro-071221-052705 , archivePrefix =. 2312.13287 , primaryClass =

    work page doi:10.1146/annurev-astro-071221-052705

  18. [18]

    The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier

    The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? II. Introducing the bouncing barrier. , keywords =. doi:10.1051/0004-6361/200912976 , archivePrefix =. 1001.0488 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/200912976

  19. [19]

    , keywords =

    The bouncing barrier revisited: Impact on key planet formation processes and observational signatures. , keywords =. doi:10.1051/0004-6361/202347716 , archivePrefix =. 2312.06000 , primaryClass =

    work page doi:10.1051/0004-6361/202347716

  20. [20]

    Bouncing Behavior of Microscopic Dust Aggregates

    Bouncing behavior of microscopic dust aggregates. , keywords =. doi:10.1051/0004-6361/201220946 , archivePrefix =. 1301.3629 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201220946

  21. [21]

    , keywords =

    The Rebound Condition of Dust Aggregates Revealed by Numerical Simulation of Their Collisions. , keywords =. doi:10.1088/0004-637X/737/1/36 , adsurl =

    work page doi:10.1088/0004-637x/737/1/36

  22. [22]

    The Physics of Protoplanetesimal Dust Agglomerates. II. Low-Velocity Collision Properties. , keywords =. doi:10.1086/525841 , archivePrefix =. 0711.2148 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1086/525841

  23. [23]

    The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals?. I. Mapping the zoo of laboratory collision experiments. , keywords =. doi:10.1051/0004-6361/200912852 , archivePrefix =. 0910.4251 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/200912852

  24. [24]

    A Statistical Theory of the Strength of Materials. Ingeni

  25. [25]

    , keywords =

    The tensile strength of dust aggregates consisting of small elastic grains: constraints on the size of condensates in protoplanetary discs. , keywords =. doi:10.1093/mnras/staa1641 , archivePrefix =. 2006.05107 , primaryClass =

    work page doi:10.1093/mnras/staa1641 2006

  26. [26]

    , keywords =

    Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization. , keywords =. doi:10.3847/1538-4357/ab45f0 , archivePrefix =. 1907.00189 , primaryClass =

    work page doi:10.3847/1538-4357/ab45f0 1907

  27. [27]

    Free collisions in a microgravity many-particle experiment. III. The collision behavior of sub-millimeter-sized dust aggregates. , keywords =. doi:10.1016/j.icarus.2013.02.034 , archivePrefix =. 1302.5532 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1016/j.icarus.2013.02.034 2013

  28. [28]

    Collisions between sintered icy aggregates

    Collisions between Sintered Icy Aggregates. , keywords =. doi:10.3847/1538-4357/aa6fad , archivePrefix =. 1705.04778 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.3847/1538-4357/aa6fad

  29. [29]

    Erosion and the limits to planetesimal growth

    Erosion and the limits to planetesimal growth. , keywords =. doi:10.1051/0004-6361/201425222 , archivePrefix =. 1412.3593 , primaryClass =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.1051/0004-6361/201425222

  30. [30]

    1986 , publisher=

    Theory of Elasticity (Third Edition). 1986 , publisher=

    1986

  31. [31]

    , keywords =

    Collisional properties of cm-sized high-porosity ice and dust aggregates and their applications to early planet formation. , keywords =. doi:10.1093/mnras/stab3348 , archivePrefix =. 2111.09141 , primaryClass =

    work page doi:10.1093/mnras/stab3348

  32. [32]

    , keywords =

    Multi-frequency analysis of the ALMA and VLA high resolution continuum observations of the substructured disc around CI Tau: Preference for submillimetre-sized low-porosity amorphous carbon grains. , keywords =. doi:10.1051/0004-6361/202452986 , archivePrefix =. 2507.08797 , primaryClass =

    work page doi:10.1051/0004-6361/202452986

  33. [33]

    , keywords =

    Dust Scattering Albedo at Millimeter Wavelengths in the TW Hya Disk. , keywords =. doi:10.3847/1538-4357/ad9f31 , archivePrefix =. 2412.10731 , primaryClass =

    work page doi:10.3847/1538-4357/ad9f31