{"id":"d2af41aa-76c1-4453-b68a-7f909203c70f","arxiv_id":"2411.11742","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Collisions between dust particles inside planet-induced spiral waves in 2D simulations reach velocities far above fragmentation thresholds, even for similar-sized particles and sub-thermal-mass planets.","lead":"Particles drifting through a protoplanetary disk can collide at very high speeds inside the spiral waves stirred up by an embedded planet, often fast enough to shatter them. The result suggests planet-induced spirals may act as additional fragmentation sites, shrinking dust grains near the planet.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"3D weakening of spiral wave velocities could invalidate the quantitative 'far exceed' claim; a 3D test is needed.","rationale":"The paper is a clean proof-of-concept that spiral-wave crossings can produce high relative velocities between dust particles of different Stokes numbers in a 2D shearing-sheet model. The method is internally consistent: the 'constant stream' assumption makes the intersection-velocity calculation a valid estimate of the collisional velocity between coexisting particle streams, and the authors are explicit that collision frequencies are not computed. The strongest claim in the abstract, however, is quantitative ('far exceed the typical fragmentation threshold'), and it is not certified against the dominant known physical simplification: the razor-thin 2D geometry. The 3D reduction in density contrast is acknowledged, but the velocity perturbation—the actual driver of particle collision speeds—is not quantified for the relevant parameters. Because the marginal small-ΔSt collision velocities are only slightly above the fragmentation threshold in the outer disk, a factor-of-two reduction in velocity would make the claim fail at large radii. This is a real, load-bearing uncertainty rather than a manufactured concern. Nevertheless, the paper's own framing is appropriately cautious in Section 4, and the central mechanism is plausible and novel. I therefore agree with the reader's weakest assumption and see no reason to change the verdict; the paper should remain accepted as an idealized demonstration, with the quantitative caveat clearly flagged.","tokens_in":7984,"tokens_out":18823,"duration_ms":178512,"concrete_test":"Run a 3D, vertically stratified local shearing-box simulation with M_p=0.5 M_th and α=0.01 (matching the 2D run), extract the equilibrium gas velocity field, and post-process particle trajectories for St=0.05 and 0.01 to compute the distribution of collision velocities at intersections inside the spiral. If the median vcoll/cs is reduced by more than a factor of 2 relative to the 2D run, or if the median converted collision velocity falls below the 10 m/s fragmentation threshold at r>50 au, then the abstract's 'far exceed' claim is not robust to 3D effects and should be qualified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim rests on the amplitude of the 2D, razor-thin spiral wave. The authors concede in Section 4 that density contrasts can be roughly halved in 3D (Tanaka et al. 2002) and that it is 'unclear' whether the velocity perturbation is reduced by as much. The collision velocities that drive the conclusion are computed from particle responses to this 2D velocity field. For the smallest ΔSt=0.025 pairs, reported collision velocities are 5–10% of the sound speed; when converted to SI (Fig. 3), these values are comparable to the top of the nominal fragmentation threshold only in the outer disk. If 3D velocity perturbations are also reduced by roughly a factor of two, the collision velocities in the outer disk could fall below the fragmentation threshold, directly undermining the abstract's 'far exceed' claim. The assertion that the 3D velocity perturbation 'can remain high' (Zhu et al. 2015; Rabago & Zhu 2021) is not quantified for the specific M_p/St values used here, leaving the quantitative applicability of the claim unverified.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses local 2D shearing-sheet hydrodynamical simulations of a gap-opening planet to study how planet-induced spiral waves affect the collisional velocities of dust particles. The gas field is evolved with the PENCIL CODE, and particle trajectories for Stokes numbers 0.01-0.1 are integrated in the post-processed equilibrium state. Collisional velocities are computed at spatial intersections of trajectories of particles with different Stokes numbers, ignoring the time dimension and assuming continuous streams of particles. The authors find that collisions occurring inside the spiral wave have much higher relative velocities than collisions elsewhere, with values increasing with planet mass, local gas density, and Stokes-number difference, reaching up to 35% of the sound speed near the gap edge and 5-10% of the sound speed even for small Stokes-number differences. They interpret these velocities as likely to exceed typical fragmentation thresholds and discuss implications for gap filtering, pebble accretion, and planetesimal formation. The paper is concise, the numerical checks (gap-depth comparison with Kanagawa et al. 2015, comparison of post-processed trajectories with PENCIL runs, validation of radial drift against analytical estimates) are appropriate, and the limitations of the 2D model are openly acknowledged in Section 4.","tokens_in":8192,"tokens_out":5091,"duration_ms":53873,"significance":"If the reported velocity enhancement carries over to more realistic 3D disks, the paper identifies a previously underappreciated fragmentation channel inside planet-induced spiral waves, with consequences for the size distribution of solids near growing planets and for the leakiness of planetary gaps. The study's strengths include the use of an independent hydrodynamic simulation, external calibration against published gap-depth fits, explicit comparison with analytical drift velocities, and the demonstration that the gas velocity field, rather than density variations, is the dominant driver of the enhanced collision velocities. The central qualitative trend is robust within the 2D model. The main weakness is that the quantitative 'far exceed' claim in the abstract is tied to the amplitude of 2D, razor-thin spiral waves, and the paper's own Section 4 acknowledges that 3D density contrasts can be reduced by roughly a factor of two without establishing whether the velocity perturbations (and hence the derived collision velocities) are similarly reduced for the specific parameter range explored.","major_comments":[{"comment":"The abstract's central quantitative claim that collision velocities 'far exceed' the fragmentation threshold is not robust to the 3D uncertainty that the authors themselves acknowledge. The text states that the 3D density contrast can be reduced by about a half (Tanaka et al. 2002) and that it is 'unclear' whether the velocity perturbation is reduced by as much. For the smallest Stokes-number difference considered, ΔSt=0.025, the reported collision velocities are 5-10% of the sound speed; when converted to SI units in Fig. 3, these values are comparable only to the upper end of the 1-10 m/s fragmentation threshold in the outer disk. A factor-of-two reduction of the velocity perturbation in 3D would place a substantial fraction of these collisions below the threshold, directly undermining the unsupported 'far exceed' wording. Please provide a quantitative estimate from existing 3D simulations (Zhu et al. 2015; Rabago & Zhu 2021) specifically for the M_p/M_th and St values used here, or perform a targeted 3D test, or alternatively qualify the claim throughout the abstract and conclusions as applying to the 2D razor-thin model and soften 'far exceed' accordingly.","section":"Section 4 and Fig. 3"},{"comment":"The entire quantitative analysis rests on a single deterministic trajectory per Stokes number per planetary mass. The statement 'we limit our study to one integrated particle trajectory per combination of St and M_p' means that all intersection-based collisional velocities are drawn from one pair of streamlines. The trajectory shape depends on the initial radial and azimuthal position (particles are introduced just interior of the radial damping zone with random azimuthal positions), and a different starting azimuth would lead to a different crossing geometry through the spiral wave and potentially different collision velocities. Because the paper describes velocities as 'commonly obtained' rather than as properties of a single path, please demonstrate robustness by integrating a small ensemble of trajectories with different initial azimuthal positions for at least one representative case, or otherwise justify that the chosen trajectory is representative of the particle population.","section":"Sections 2.2 and 3"},{"comment":"No numerical convergence test is presented for the velocity field that drives the particle collisions. The paper fixes the resolution at 32 cells per H_g and does not compare against a coarser or finer run. The highest collision velocities (exceeding 35% of the sound speed near the gap edge) occur in regions where the spiral perturbation is strong and potentially shock-like, and the amplitude of such perturbations can be resolution-dependent even with high-order schemes. A resolution study with, e.g., 16 and 64 cells per H_g for one planetary mass would establish whether the quantitative collision velocities reported in Figs 2-4 are converged, and whether the 'far exceed' conclusion is stable.","section":"Section 2.1 and Figs 2-4"}],"minor_comments":[{"comment":"The caption states that the colored areas indicate the region between the 10th and 90th percentiles, but the figure legend simply labels 'collisional velocity'; please add a sentence in the caption clarifying that the shaded band is a percentile range, not a measurement uncertainty.","section":"Section 3, Fig. 3 caption"},{"comment":"The sentence 'Observations of protoplanetary disks generally suggest turbulent parameters on the order of α_turb ~ 10^-5 - 10^-3' would benefit from a brief mention that the corresponding turbulent collision velocities are for the specific particle sizes/Stokes numbers considered here, since Eq. (7) depends on St.","section":"Section 4, text after Eq. (7)"},{"comment":"The keyword list includes 'planets and satellites: general'; given the paper's focus on dust evolution, consider adding 'protoplanetary discs' as a keyword.","section":"Introduction, first paragraph"},{"comment":"The gravitational potential term is written as -∇Φ without a subscript; for clarity, please denote the planetary potential as Φ_p or state explicitly that Φ is the smoothed planetary potential from Eq. (1).","section":"Section 2.2, Eq. (5)"},{"comment":"The name 'Kruijer' in the reference list uses a ligature that may not render correctly in all bibliographic styles; 'Kruijer' is the standard ASCII rendering.","section":"Reference list"}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of MNRAS and is a solid idealized proof-of-concept. My main concern is that the abstract overstates the certainty of a 2D result: the authors themselves list the 3D uncertainty and the absence of collision probabilities as limiting factors, yet the abstract states that collisions 'far exceed' the fragmentation threshold without qualification. A major revision that either supplies a quantitative 3D argument or tempers the claim would resolve this. I do not see any circularity or fitted-constant issues; the comparison with external fragmentation thresholds and the internal checks (Kanagawa et al. 2015 gap depths, analytical drift) are appropriate. The paper would also benefit from a brief resolution study, as the quantitative amplitudes are the core of the claimed effect."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a solid proof-of-concept that spiral waves launched by a planet can act as a fragmentation site for pebbles. The new result is that collisions between particles with different Stokes numbers crossing the spiral have relative velocities well above typical fragmentation thresholds, with the trend increasing with planet mass and ΔSt. That is distinct from the gap-edge fragmentation story of Drążkowska et al. (2019) and worth having on the record.\n\nThe paper does several things right. The hydro setup follows Yang & Zhu (2020), the gap depths match Kanagawa et al. (2015) within 10%, and the particle drift far from the planet agrees with the analytical drift rate. The authors are explicit that they are not computing collision frequencies, and they flag the 2D limitation and the mass-transfer ambiguity. That is the right level of honesty for a short paper.\n\nSoft spots: the central velocity estimate rests on one deterministic trajectory per Stokes number. That is fine for a proof of concept, but it means the spatial distribution of 'collisions' is limited to the intersections of a handful of trajectories, and the percentile ranges in Fig. 3 should not be read as a statistical sample. There are no resolution or convergence tests, which is a minor gap but not fatal given the checks they do. The 3D question is the one that matters for the quantitative claim: spiral waves are weaker in 3D, and if the velocity perturbation is also roughly halved, the outer-disk collision velocities could drop below the fragmentation threshold. The authors acknowledge this and cite Zhu et al. (2015) and Rabago & Zhu (2021) for the velocity staying high, but they do not quantify it for their parameter range. That is a real uncertainty, not a misstep: the paper's contribution is the mechanism, not a precise fragmentation map.\n\nWho this is for: anyone working on dust evolution near planets, gap filtering, or pebble accretion. It deserves a serious referee; it is not a desk reject. The main thing I would ask for in revision is a short convergence check and a more explicit caveat in the abstract that the 'far exceed' claim is 2D and conditional on collision frequency.","headline":"A clean, honest proof-of-concept that planet-induced spiral waves can fragment dust; the 3D caveat is real but already flagged.","tokens_in":8712,"tokens_out":1974,"would_cite":true,"duration_ms":18726,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Planet-launched spiral waves drive dust collisions fast enough to fragment pebbles.","keywords":["protoplanetary discs","planet-disc interactions","dust fragmentation","spiral density waves","pebble dynamics","shearing sheet simulations","planetesimal formation"],"falsifier":"A 3D local or global simulation of the same planet masses ($M_{\\rm p}/M_{\\rm th}=0.25$–$1$) measuring the gas velocity jump across the spiral and the resulting particle collisional velocities would settle the claim: if peak collisional velocities inside the spiral drop below about 1–10 m/s for a Stokes-number difference of 0.025, the central claim would fail. A laboratory measurement showing the fragmentation threshold is substantially higher than 10 m/s for the relevant pebbles would also remove the effect.","tokens_in":7787,"feed_emoji":"🌀","tokens_out":8515,"duration_ms":75082,"temperature":0.7,"pith_summary":"The paper argues that the spiral density waves a growing planet excites in its protoplanetary disk are sites where dust pebbles collide fast enough to fragment, not just stick. Using local 2D shearing-sheet hydrodynamical simulations with a gap-opening planet, the authors track particles of different sizes and find that their trajectories bend sharply where they cross the spiral. Collisional velocities at those crossings far exceed the typical fragmentation threshold of $1$–$10\\,\\mathrm{m\\,s^{-1}}$, even for collisions between nearly equal-sized particles and for planet masses below the pebble isolation mass. If such collisions occur often enough, they would grind particles into progressively smaller sizes closer to the planet, with consequences for dust crossing gaps, pebble accretion, and planetesimal formation.","feed_headline":"Planet spiral waves shatter pebbles that should stick","feed_subtitle":"Even similar-sized particles collide fast enough inside the spiral to fragment, and the grinding intensifies near the planet.","key_machinery":"The central object is the planet-induced spiral wave in a razor-thin local shearing-sheet model. The wave is a bending of gas streamlines driven by the planet's gravitational potential; particles with different Stokes numbers (a dimensionless measure of drag coupling to the gas) respond differently to that bend, so their trajectories cross. At each crossing the paper computes a collisional velocity from the difference in particle velocities and compares it with the fragmentation threshold $v_{\\rm frag}$. The Stokes-number difference is the lever: larger $\\Delta\\mathrm{St}$ means more divergent trajectories and faster collisions.","core_discovery":"The central claim is that planet-induced spiral waves are a fragmentation site: the velocity perturbation across the spiral bends gas streamlines and particle trajectories, and intersections between trajectories of particles with different Stokes numbers occur at high relative velocity inside the spiral. In the simulations, collisional velocities there reach up to 20–35% of the sound speed near the gap edge, and even for a Stokes-number difference of only $\\Delta\\mathrm{St}=0.025$ they exceed 5–10% of the sound speed, above typical fragmentation thresholds in most of the disk. The effect increases with planetary mass and local gas density, and it is driven by the gas velocity field rather than by density enhancements: holding the Stokes number constant changes collision velocities by only about 20%. The authors conclude that with sufficiently frequent collisions, the spiral produces progressively smaller particles with decreasing distance from the planet.","pith_inferences":["Editorial extension: a direct 3D simulation test is the natural next step; the paper notes 3D spirals have roughly halved density contrasts, but if the velocity perturbation stays high, the fragmentation claim could survive in 3D disks.","Editorial extension: coupling these trajectory crossings to a collision clock (orbital-phase information) would convert the velocity map into a fragmentation rate, the missing ingredient for coagulation models.","Editorial extension: dust back-reaction is neglected here; in dust-rich regions, feedback on the gas could alter the spiral's velocity field and either moderate or amplify the collision speeds."],"forward_implications":["Inside the spiral, collision speeds computed from turbulence alone would be at most a few percent of the sound speed for typical $\\alpha_{\\rm turb}$, so the spiral adds a fragmentation channel that standard dust-growth models miss.","Collisional velocities increase with planetary mass and local gas density and decrease with distance from the planet, so the spiral acts as a radial grinder that makes particle sizes decrease inward.","Smaller particles cross planetary gaps more easily, so the result supports leaky gap models and weakens the case that Jupiter's early core isolated the inner and outer Solar System reservoirs.","Smaller particles drift more slowly and are less efficiently accreted by pebble accretion, and they are poorer streaming-instability planetesimal precursors."],"supporting_citations":[{"why":"Provides the numerical framework for the shearing-sheet gas-plus-particle model, including the planet potential and damping zones, which this work adapts.","marker":"Yang & Zhu (2020)"},{"why":"Gives the empirical gap-depth scaling used to validate the simulated gap depths.","marker":"Kanagawa et al. (2015)"},{"why":"Earlier finding that fragmentation at the gap edge enhances small particles; the baseline for the new spiral-interior mechanism.","marker":"Drążkowska et al. (2019)"},{"why":"Basis for the radial drift picture and for the gas-density-dependent stopping time used in the particle integrator.","marker":"Weidenschilling (1977)"},{"why":"Supplies the turbulent relative-velocity formula used as the baseline for comparison with spiral-induced collision speeds.","marker":"Ormel & Cuzzi (2007)"},{"why":"Supplies the typical fragmentation velocity range of 1–10 m/s used to judge whether collisions fragment.","marker":"Blum & Wurm (2008)"},{"why":"Temperature profile used to convert simulated velocities into physical units (m/s).","marker":"Chiang & Goldreich (1997)"},{"why":"Source of the caveat that 3D spirals have weaker density contrasts, defining the main limitation of the 2D result.","marker":"Tanaka et al. (2002)"},{"why":"Cited as evidence that 3D velocity perturbations can remain high in sound-speed units, countering the 3D caveat.","marker":"Zhu et al. (2015)"}],"fun_headline_variants":["Spiral waves shatter pebbles near planets","Planet spirals grind pebbles into dust","Spiral wave collisions shatter pebble growth","High-speed spiral collisions fragment pebbles"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the spiral-wave velocity perturbations seen in the 2D shearing sheet are representative of real 3D disks; if 3D effects weaken the velocity jumps as much as they weaken density contrasts, collision speeds could fall below the fragmentation threshold.","fun_headline_variants_meta":{"raw":{"variants":["Spiral waves shatter pebbles near planets","Planet spirals grind pebbles into dust","Spiral wave collisions shatter pebble growth","High-speed spiral collisions fragment pebbles"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000512,"raw_usage":{"total_tokens":2457,"prompt_tokens":879,"completion_tokens":1578,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":495,"completion_tokens_details":{"reasoning_tokens":1518}},"tokens_in":495,"tokens_out":1578,"duration_ms":11380,"temperature":1.0,"reasoning_tokens":1518,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T18:12:08.790524+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A 3D local or global simulation of the same planet masses ($M_{\\rm p}/M_{\\rm th}=0.25$–$1$) measuring the gas velocity jump across the spiral and the resulting particle collisional velocities would settle the claim: if peak collisional velocities inside the spiral drop below about 1–10 m/s for a Stokes-number difference of 0.025, the central claim would fail. A laboratory measurement showing the fragmentation threshold is substantially higher than 10 m/s for the relevant pebbles would also remove the effect.","supporting_citations":[],"review_version":1}