{"id":"a73889fa-f538-4eae-9914-38f0c99644e8","arxiv_id":"2508.20070","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Fully resolved Keplerian turbulence concentrates pebbles into clumps and reduces their radial drift, potentially bypassing the drift and fragmentation barriers to planetesimal formation at canonical dust-to-gas ratios.","lead":"Simulations that resolve, rather than average over, the turbulence in a protoplanetary disk show that pebble-sized dust gathers into tight clumps inside swirling eddies and that the dust's fall toward the star slows or stops. If the effect survives in three dimensions, it would offer a single-mechanism route to planet formation starting from the normally observed low dust-to-gas ratio.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Clustering amplitude is never translated into a volume density: the paper does not show that the point clusters exceed the Roche density, so the 'unique path' to gravitational collapse is unsupported.","rationale":"The reader's weakest assumption is that the 2D forced-turbulence model captures real disk dynamics and that clustering is sufficient for collapse. I focus on the latter half of that assumption because it is the most directly load-bearing for the paper's central claim: even if the 2D clustering is quantitatively reproduced in 3D, the paper still does not show that the resulting overdensities reach the Roche density needed for gravitational collapse. The abstract's 'unique path' is an overstatement without this demonstration. The paper does have strengths: it uses resolved, low-dissipation turbulence rather than a diffusion coefficient, explores a wide Stokes-number range, and explicitly acknowledges the 2D restriction and the reliance on unpublished 3D work. These factors support a conditional verdict, not outright rejection. The proposed concrete test—computing the volume density of the point clusters and comparing to the Roche density—would directly settle whether the mechanism can produce planetesimals. Because this gap is already reflected in the reader's conditional verdict, I leave the verdict unchanged.","tokens_in":4477,"tokens_out":3608,"duration_ms":45330,"concrete_test":"Post-process the point-cluster snapshots from the 2D simulations (Fig. 1, 4th panel). Assign each particle a physical radius from its stopping time using the Epstein drag law at 5 AU (gas surface density 100 g/cm^2, solid density 3 g/cm^3). Compute the local dust volume density in a sphere of radius equal to the cluster's effective radius (e.g., the typical particle spacing or the mean particle radius) and compare to the Roche density ρ_R = 3 M_sun / (4π a^3) at 5 AU. If the cluster density does not exceed ρ_R (or a Toomre-style criterion for the dust layer), then the mechanism cannot trigger gravitational collapse, and the 'unique path' to planetesimal formation is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that Keplerian turbulence concentrates solids into point clusters that can overcome the barriers to planetesimal formation at canonical dust-to-gas ratios. The paper identifies the point-cluster regime (§3.1) and shows drift reduction (§3.2), but it never quantifies the resulting dust density enhancement or compares it to the threshold for gravitational collapse. In the simulations, particles are Lagrangian points in 2D; 'all the particles of a family reach the same position' is a mathematical concentration, not a physical density. The cluster's physical volume is set by the particle radius and stopping time, and the mass within a Hill sphere may be far below the Roche density. The paper mentions that concentrations of 'up to four orders of magnitude' are needed (Introduction) but provides no measurement or estimate that the clusters achieve this. It also explicitly neglects back-reaction and self-gravity, so the step from a passive overdensity to a collapsing cloud is absent. Even if 3D simulations (Gerosa et al., in prep) confirm columnar vortices, the current paper does not demonstrate the final collapse, making the advertised 'unique path' an extrapolation rather than a result.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This short proceedings paper reports two-dimensional incompressible shearing-box simulations of forced Keplerian turbulence with Lagrangian dust particles. For low turbulent intensity (α ~ 10^-3–10^-4) and pebble-sized stopping times, it identifies a 'point-cluster' regime in which particles of a given stopping time converge to the same position and velocity inside anticyclonic eddies, and a companion set of simulations in which a constant azimuthal force is applied to the particles to model pressure-gradient drift. The authors show that the radial drift can be reduced or halted in this regime. They conclude that Keplerian turbulence alone can overcome the drift and fragmentation barriers and provide a unique path to planetesimal formation at canonical dust-to-gas ratio. The manuscript is explicitly a synthesis and re-interpretation of the authors' earlier work (Gerosa et al. 2023, 2024), with a 3D extension cited as in preparation.","tokens_in":4763,"tokens_out":8021,"duration_ms":88005,"significance":"The proposed mechanism, if quantitatively confirmed, would be a valuable addition to planetesimal-formation scenarios because it would show that a weakly turbulent disk can concentrate solids without invoking pressure bumps or the streaming instability. The paper has real strengths: it uses direct numerical simulations of the flow rather than diffusion models; it classifies particle behavior with a Lyapunov-dimension diagnostic; it recovers the known high-turbulence behavior (diffusion, turbulent concentration, enhanced drift), which gives confidence in the numerical setup; and it reports a parameter sweep over Stokes and Rossby numbers. The central limitation is that the step from Lagrangian clustering to a physical, gravitationally collapsible density enhancement is not made. Because of that missing quantitative link, the significance currently rests on an extrapolation rather than on a demonstrated result.","major_comments":[{"comment":"The point-cluster regime is not converted into the quantity that actually matters for the paper's central claim, namely the local dust volume density or dust-to-gas ratio. The text states that 'all the particles of a family reach the same position with the same velocity' (Sec. 3.1), but in a Lagrangian simulation this is a mathematical concentration in phase space. The Introduction quotes a required enhancement of 'up to four orders of magnitude,' yet no measurement of the density contrast in the clusters, no estimate of the mass within a Hill radius, and no comparison to the Roche density are given. Without this, 'a unique path to planetesimal formation' is not established by the data presented.","section":"§3.1, Fig.1"},{"comment":"The drift-reduction result has no associated error bars, ensemble statistics, or convergence study. The verbal summary 'reduced by a factor between 0.5 and 1.0' is ambiguous and does not specify whether the ratio plotted is v_turb/v_laminar or 1 - v_turb/v_laminar, nor how it depends on Stokes number and α. The statement that drift is 'fully arrested' is only supported by the point-cluster coincidence with zero relative velocity. Please provide the measured values, their uncertainties, and the exact definition of the plotted quantity, since this is one of the two pillars of the claim that the drift barrier is overcome.","section":"§3.2, Fig.2"},{"comment":"The transition from this idealized setup to protoplanetary disks is not yet documented. The simulations are 2D, driven by an external forcing at a single wavenumber, neglect particle back-reaction, and do not include stratification or self-gravity. The justification for 2D relies on a paper in preparation (Gerosa et al., in prep). Moreover, if the clusters actually achieve the factor of 10^4 concentration quoted in the Introduction, the local dust-to-gas ratio would exceed unity and the neglect of back-reaction would be inconsistent at the cluster scale. Please provide a published 3D check or a quantitative argument for why these approximations do not change the conclusion, and estimate the local dust-to-gas ratio implied by the clustering.","section":"§2, last paragraph; §4, last paragraph"},{"comment":"The title/abstract claim of a 'unique' solution or 'unique path' is stronger than what is shown. The manuscript demonstrates one possible mechanism in a simplified model; it does not compare the efficiency or robustness of this path with pressure-bump plus streaming-instability scenarios, nor does it prove that point clusters form from global disk initial conditions. I recommend either adding a comparative discussion or softening the uniqueness claim (e.g., 'a direct path' or 'a single mechanism').","section":"Abstract and §4"}],"minor_comments":[{"comment":"The abstract says turbulence is 'fully captured rather than modeled with a turbulent diffusion or turbulent viscosity parameter,' but §2 states that the onset of turbulence is not captured and that a turbulent state is sustained by an imposed external forcing. Please revise to avoid overstatement.","section":"§2"},{"comment":"The text refers to '2nd panel,' '3rd panel,' and '4th panel,' but the figure panels are not labeled; add panel letters and color bars with a vorticity scale.","section":"§3.1, Fig.1"},{"comment":"The conversion from turbulent Stokes number to Keplerian Stokes number is given by reference to Sengupta et al. (2024) but no formula is provided; please state the relation used.","section":"§3.1"},{"comment":"Please state the Reynolds number and the resolution relative to the forcing scale and the Kolmogorov scale, so the reader can assess whether the flow is fully resolved.","section":"§2"},{"comment":"The invocation of Taylor (1922) to justify 2D columnar vortices in 3D disks is an oversimplification; the Taylor-Proudman theorem is modified by stratification, vertical boundaries, and shear. Please qualify.","section":"§4"},{"comment":"The entry 'Simon, J. B., Birnstiel, J. B., & Nesvorný, D. 2024' appears to have a wrong author name ('J. B. Birnstiel'?) and the arXiv number is missing; please check.","section":"References"},{"comment":"The manuscript would benefit from a sentence stating which simulations are new to this paper and which are re-analyzed from Gerosa et al. (2023, 2024), given that the text says 'We present numerical simulations.'","section":"§1 and §4"}],"recommendation":"major_revision","confidential_remarks":"This is a proceedings-style paper with a very broad claim. The quality of the numerics appears adequate for the qualitative phase diagram, but the advertised 'unique path to planetesimal formation' is not backed by a density-collapse calculation. I recommend major revision with a request for a quantitative estimate of the clustering amplitude and a more careful framing. It may be appropriate for a conference proceedings after those additions; I would not accept the uniqueness claim in its current form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a conference proceedings paper that repackages the authors' own published simulations (Gerosa et al. 2023, 2024) under a new astrophysical framing. The new element is the identification of a 'Keplerian turbulence' regime at low Rossby number, where point clusters form in anticyclonic eddies and the radial drift of pebbles is reduced by a factor of 0.5-1.0, or even fully arrested for the clustered particles. That is a genuinely interesting result, and the paper is honest that it is a 'novel perspective' rather than new simulations. The 2D assumption is justified with reference to Taylor-Proudman columnar vortices and to unpublished 3D simulations (Gerosa et al., in prep.), which is a reasonable though unverifiable position for a proceedings paper.\n\nThe main soft spot is exactly what the stress-test note flags: the paper never converts the observed clustering into a volume density or mass. It says that gravitational collapse requires 'up to four orders of magnitude' concentration, but it never measures whether the point clusters achieve that. The clusters are mathematical points in the Lagrangian sense; their physical density depends on particle radius and stopping time, which are not modeled. The paper also neglects back-reaction and self-gravity, so the step from a passive overdensity to a collapsing cloud is absent. Combined with the lack of error bars on the drift reduction and the 2D setup, the advertised 'unique path' is an extrapolation, not a demonstrated result.\n\nThat said, this is a short SF2A contribution, and the authors clearly state their simplifications and point to their prior work for the full method. The physics is plausible and the parameter-space map (Lyapunov dimension classification) is a useful organizing principle. The paper does not commit the sin of hiding its limitations; it openly says that streaming instability and back-reaction are neglected, and that 3D confirmation is forthcoming.\n\nMy verdict: worth a serious referee if expanded into a full journal paper with a density calculation, but as written it is a conference abstract with an overconfident title. I'd bring it to a reading group, and I'd cite the earlier Gerosa papers rather than this one. If it lands on my editorial desk, I would send it out—on the strength of the underlying simulations and the potential importance of the regime—but I would expect the referee to require the quantitative collapse threshold and either the 3D results or a clear statement that they are still preliminary.","headline":"A plausible, clearly presented synthesis of the authors' own simulations, identifying a low-turbulence regime where pebbles cluster in anticyclonic eddies and radial drift is reduced or halted; but the leap to a 'unique path' to planetesimal formation is not backed by a density or collapse calculation.","tokens_in":5245,"tokens_out":1213,"would_cite":false,"duration_ms":16022,"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":"Weakly turbulent disks alone may gather pebbles into dense clusters and halt their drift, offering a single-process path to planetesimal formation at ordinary dust-to-gas ratio.","keywords":["protoplanetary disks","planetesimal formation","Keplerian turbulence","radial drift","anticyclonic eddies","dust clustering","shearing box","pebbles"],"falsifier":"A 3D stratified shearing-box simulation with self-gravity and particle back-reaction, using the same forcing and parameter range, that shows the peak dust-to-gas ratio inside anticyclonic eddies never reaches order unity would disprove the claimed route; equivalently, a 2D run with back-reaction included that destroys the point clusters would do the same.","tokens_in":4387,"feed_emoji":"🪐","tokens_out":5802,"duration_ms":59600,"temperature":0.7,"pith_summary":"This paper proposes that the turbulence present in protoplanetary disks, when captured in full rather than approximated by a diffusion coefficient, can itself solve the two classic barriers to planetesimal formation: inward radial drift and insufficient solid concentration. In two-dimensional shearing-box simulations of Keplerian turbulence at low turbulent intensity, centimeter-to-meter solids clump into point clusters inside anticyclonic eddies, and their drift toward the star is slowed or completely stopped. The authors argue this 'Keplerian turbulence clustering' is a distinct mechanism from ordinary turbulent concentration, and that it could allow gravitational collapse to form planetesimals without invoking pressure bumps or the streaming instability. A sympathetic reader would care because it suggests a unique, self-contained route to planet formation in the standard low-dust-to-gas-ratio disk.","feed_headline":"Weak disk turbulence traps pebbles and stops their drift","feed_subtitle":"Low-intensity Keplerian turbulence alone can gather solids into clusters, opening a new planetesimal route.","key_machinery":"The central mechanism is 'Keplerian turbulence clustering': long-lived anticyclonic vortices in weakly turbulent disks concentrate solid grains toward their centers via the Coriolis force, forming point clusters. The same elongated eddy geometry that traps particles also explains the drift reduction, because the eddies' small radial extent limits the radial acceleration of particles and their long azimuthal extent blocks inward motion. A complementary tool, the Lyapunov dimension of the particle set in phase space, is used to classify the three dynamical behaviors: diffusion, filamentary structures, and point clusters.","core_discovery":"The paper reports that when protoplanetary disk turbulence is treated as fully developed Keplerian turbulence with Rossby number below one, three dynamical regimes appear depending on turbulent intensity and particle stopping time. At low turbulent intensity (alpha around 10^-3 to 10^-4), pebble-sized solids concentrate into point clusters inside anticyclonic eddies, where the Coriolis force overcomes centrifugal expulsion; all particles of a given size reach the same position and velocity and evolve as a single particle. In the same regime, the radial drift of centimeter-to-meter solids is reduced by a factor of 0.5 to 1.0 relative to the laminar case, and can be fully arrested when the sol","pith_inferences":["A testable implication the paper leaves unexplored is the actual density contrast inside the point clusters: if the local dust-to-gas ratio in these clusters reaches the order-unity threshold, gravitational collapse would follow directly; the paper does not report that ratio explicitly.","If the mechanism survives vertical stratification, as the authors' unpublished 3D simulations suggest, then low-viscosity 'dead zones' of disks would become preferred sites for planetesimal formation rather than quiet regions where formation stalls.","The toy model of elongated eddies implies a quantitative prediction: the drift reduction should scale with eddy aspect ratio, so simulations or observations that constrain eddy shapes in low-Rossby-number turbulence could test the mechanism independently of cluster formation.","Including particle back-reaction could go either way: the streaming instability might add to the concentration, or the feedback might disrupt the eddies; the robustness of the mechanism to back-reaction remains the main untested margin."],"forward_implications":["If this mechanism operates, planetesimals could form in weakly turbulent disk regions without any external pressure bump, removing a common requirement in current formation scenarios.","The radial drift barrier would be relaxed for the pebble sizes most vulnerable to loss, since their inward velocity is reduced by up to a factor of two and can vanish entirely in the cluster regime.","The collision/fragmentation barrier may be mitigated because particles inside a point cluster share the same velocity and position, reducing relative velocities within the cluster.","The mechanism reproduces standard turbulent concentration at higher turbulent intensities, so it extends rather than replaces existing results, and its clustering is distinct from that of Hogan & Cuzzi (2001).","The formation of point clusters at arbitrarily low dust-to-gas ratio, in the absence of back-reaction, implies that the concentration step does not need a preconditioning of the dust-to-gas ratio."],"supporting_citations":[{"why":"Establishes the concentrating effect of anticyclonic vortices that the paper invokes for its point-cluster mechanism.","marker":"Barge & Sommeria 1995"},{"why":"Shows how Coriolis force drives particle concentration in vortices, a basis for the cluster formation.","marker":"Tanga et al. 1996"},{"why":"Provides the standard turbulent concentration model that the paper contrasts with its new Keplerian clustering regime.","marker":"Hogan & Cuzzi 2001"},{"why":"Supplies the laminar radial drift solution against which the measured drift reduction is compared.","marker":"Nakagawa et al. 1986"},{"why":"Defines the streaming instability, the clustering mechanism the paper deliberately excludes and does not need.","marker":"Youdin & Goodman 2005"},{"why":"Gives the definitions of turbulent and Keplerian Stokes numbers and Rossby number used to map particles to sizes.","marker":"Sengupta et al. 2024"},{"why":"Provides the shearing-box framework used for the simulations.","marker":"Hawley et al. 1995"},{"why":"Justifies the incompressible description of subsonic disk turbulence at scales smaller than the scale height.","marker":"Lesur & Longaretti 2005"},{"why":"Contains the detailed methodology and prior results that this paper builds on and interprets astrophysically.","marker":"Gerosa et al. 2024"},{"why":"Quantifies the collision and fragmentation barrier that the clustering mechanism is claimed to help overcome.","marker":"Ormel & Cuzzi 2007"}],"fun_headline_variants":["Weak turbulence traps pebbles, halts drift","New planetesimal route: low-intensity turbulence","Keplerian turbulence alone clusters pebbles","Pebble drift stopped by weak disk turbulence","Turbulence enables planetesimals at low dust"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The path to planetesimal formation assumes that the point clusters seen in 2D, externally forced, incompressible turbulence without back-reaction or self-gravity correspond to real dense clumps in a 3D stratified disk that can gravitationally collapse; if 3D effects dilute the clusters or the concentration never reaches the collapse threshold, the claim fails.","fun_headline_variants_meta":{"raw":{"variants":["Weak turbulence traps pebbles, halts drift","New planetesimal route: low-intensity turbulence","Keplerian turbulence alone clusters pebbles","Pebble drift stopped by weak disk turbulence","Turbulence enables planetesimals at low dust"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000186,"raw_usage":{"total_tokens":1118,"prompt_tokens":660,"completion_tokens":458,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":404,"completion_tokens_details":{"reasoning_tokens":386}},"tokens_in":404,"tokens_out":458,"duration_ms":5335,"temperature":1.0,"reasoning_tokens":386,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T15:12:59.499140+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A 3D stratified shearing-box simulation with self-gravity and particle back-reaction, using the same forcing and parameter range, that shows the peak dust-to-gas ratio inside anticyclonic eddies never reaches order unity would disprove the claimed route; equivalently, a 2D run with back-reaction included that destroys the point clusters would do the same.","supporting_citations":[{"cited_title":"1996, Icarus, 121, 158","cited_arxiv_id":null,"evidence_quote":"Shows how Coriolis force drives particle concentration in vortices, a basis for the cluster formation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the standard turbulent concentration model that the paper contrasts with its new Keplerian clustering regime."},{"cited_title":"1986, Icarus, 67, 375","cited_arxiv_id":null,"evidence_quote":"Supplies the laminar radial drift solution against which the measured drift reduction is compared."},{"cited_title":"Length and Velocity Scales in Protoplanetary Disk Turbulence","cited_arxiv_id":"2402.15475","evidence_quote":"Gives the definitions of turbulent and Keplerian Stokes numbers and Rossby number used to map particles to sizes."},{"cited_title":"& Longaretti, P.-Y","cited_arxiv_id":null,"evidence_quote":"Justifies the incompressible description of subsonic disk turbulence at scales smaller than the scale height."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Quantifies the collision and fragmentation barrier that the clustering mechanism is claimed to help overcome."}],"review_version":1}