{"id":"c661f5a1-18eb-46a4-a812-8ceb64d4a6c3","arxiv_id":"2508.05858","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"In weakly turbulent particle-laden midplane layers of planet-forming disks, small-scale dust clumping is produced by turbulent concentration, with particles gathered in high-strain regions and absent from coherent vortices.","lead":"The authors show that small-scale dust clumping in their simulations of planet-forming disk midplanes is produced by turbulent concentration, the vortex-strain mechanism that clusters droplets in turbulent flow, rather than mainly by the streaming instability. If right, previously reported density spikes in streaming-instability simulations are turbulent concentration operating inside large-scale filaments.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Claimed St_omega collapse to 0.4-0.7 is contradicted by Table 4: StK=0.01 peak is 0.35, outside the stated range.","rationale":"The reader's stated weakest assumption was the representativeness of the 2D axisymmetric turbulent state, which is certainly important and explicitly acknowledged by the authors. However, the most immediately load-bearing concern is the internal contradiction between the abstract/Fig. 6 caption and Table 4 regarding the St_omega 'collapse.' This contradiction affects the paper's central quantitative claim and is directly checkable without new simulations. The reader did note the tension ('the Fig. 6 caption contradicts Table 4') but did not elevate it to the weakest-assumption level. I agree with the reader's overall CONDITIONAL verdict: the qualitative TC evidence (strain-vorticity correlation, RDF trends) is convincing and prevents a REJECT, but the headline St_omega range and the implied universality need correction or additional support. The 3D caveat remains a secondary limitation rather than an internal flaw; the paper already frames its conclusions as axisymmetric. Correcting the St_omega reporting would align the paper's claims with its data, and would be a precondition for any stronger acceptance. The proposed test—recomputing the peak values from the existing data—is cheap and decisive. If the peak is genuinely 0.35 for StK=0.01, the 'collapse' narrative is misleading and the comparison with H17's critical range (0.4-0.7) loses its anchor.","tokens_in":25185,"tokens_out":5634,"duration_ms":56720,"concrete_test":"Independently recompute the St_omega PDF peaks from the simulation data (or read them from the histograms in Fig. 6) and compare with Table 4. Specifically, verify the peak for StK=0.01 in the zlim=h and zlim=0.02H panels. If the peak is indeed 0.35, the abstract's '0.4-0.7' range is unsupported; then either correct the abstract and Fig. 6 caption, or provide multi-snapshot/time-averaged PDFs that shift the peak into the claimed range. Also recompute the peaks for all StK from the same snapshots to confirm whether the spread is genuine.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The abstract's headline quantitative result—'effective Stokes number within particle voids, St_omega, collapses to values in the range ~0.4-0.7'—is not supported by the paper's own Table 4. For the two zlim choices that Fig. 6's caption claims show convergence (zlim=h and zlim=0.02H), Table 4 lists St_omega(peak) = 0.35 for StK=0.01, 0.44/0.43 for StK=0.02, and 0.72/0.75 for StK=0.04. Thus the peaks span 0.35-0.75, not 0.4-0.7, and the lowest-Stokes case sits well outside the claimed range. The Fig. 6 caption asserts '0.4<St_omega(peak)<0.75' for the left/middle panels, directly contradicting the table. This is not a minor typo: the 'collapse' to a universal range is a central piece of evidence for the TC interpretation and for the comparison with H17's critical Stokes numbers. If the StK=0.01 peak is 0.35, the data show a monotonic increase with StK rather than a collapse, weakening the claim of universality. The paper offers no time-averaged statistics or error bars to reconcile this, and each run contributes only a single snapshot. This internal inconsistency must be resolved before the central claim can be accepted.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper analyzes three axisymmetric shearing-box simulations of settled particle layers in a globally laminar protoplanetary disk, with particle Stokes numbers St_K = 0.01, 0.02, and 0.04 and midplane metallicities Z = 0.01–0.02. It reports small-scale particle clumping attributed to turbulent concentration (TC): particle-rich filaments align with high gas strain rate and enclose gas-only voids with coherent vorticity. The authors introduce a void-based effective Stokes number St_omega = t_s omega and claim its peak collapses to ~0.4–0.7 across the St_K range. They further compare measured radial distribution functions (RDFs) with the Hartlep et al. (2017) cascade model using a generalized spectral slope n ≈ -2, compare SI linear growth rates with turbulent overturn frequencies, and suggest Symmetric Instability rather than SI drives the turbulence. They conclude that TC is a persistent, self-generated mechanism and that small-scale Roche-exceeding fluctuations in SI filament simulations may be TC operating on elevated filament densities.","tokens_in":25421,"tokens_out":3777,"duration_ms":41124,"significance":"If the central claim holds, the paper provides the first direct demonstration of turbulent concentration under nebula conditions and challenges the default assumption that small-scale clumping in settled, low-metallicity midplane layers is dominated by the streaming instability. The work is strengthened by direct, non-inferential diagnostics: the spatial anti-correlation between particles and vorticity and correlation with strain rate (Fig. 7), a generalized spectral rescaling of the RDF (Eq. 25b, Fig. 8), and an explicit first-principles treatment of particle–gas velocity coupling in Appendices A–B. The authors are also candid about the axisymmetric nature of the simulations. However, the quantitative centerpiece of the paper — the claimed St_omega collapse — is internally inconsistent with the paper's own Table 4, and the statistical support from single snapshots is thin. These issues must be repaired before the headline claim can be accepted.","major_comments":[{"comment":"The abstract states that St_omega 'collapses to values in the range ~0.4–0.7', and the Fig. 6 caption states that for z_lim = h and 0.02H the peak satisfies 0.4 < St_omega(peak) < 0.75. Table 4, however, lists St_omega(peak) = 0.35 for St_K = 0.01 at z_lim = 0.01H, h, and 0.02H. The values are 0.35, 0.44, and 0.72 for St_K = 0.01, 0.02, and 0.04 respectively at z_lim = h. The lowest-Stokes case therefore lies outside the stated range, and the sequence is monotonic with St_K rather than a collapse. Because this claimed universality is a central piece of evidence for TC and for the connection to the critical St ≈ 0.3 scale, this inconsistency must be resolved. Either the stated range must be corrected, or additional time-averaged statistics must be provided to demonstrate that the peak values are robust and not single-snapshot fluctuations.","section":"§4.2 / Table 4"},{"comment":"The simulations are 3D-axisymmetric: azimuthal gradients are absent, and the turbulence spectra, strain/vorticity statistics, void structure, and RDFs are all computed in the (R,z) plane. The paper itself concedes in §6 that 'Whether this behavior persists with similar clarity in fully 3D simulations remains to be seen.' This is a load-bearing limitation because the title and abstract generalize to planet-forming disks, and because the claimed universality of St_omega near 0.4–0.7 and the attribution of small-scale clumping to TC depend on the 2D turbulent state. Without a 3D test—or at minimum a careful statement that the conclusions are provisional to axisymmetric dynamics—the central claim outruns the presented evidence.","section":"§6"},{"comment":"Each simulation contributes a single snapshot: Table 2 footnote c states that the analyzed timestamp is a single time for each run. All statistics—spectra, PDFs, St_omega peaks, RDFs, and particle–gas velocity comparisons—are therefore derived from one instant per simulation. The paper repeatedly describes TC as a 'persistent' feature and stresses robustness, but no time averaging, multiple snapshots, or uncertainty estimates are given for the St_omega peaks or RDF inflection points. A single snapshot could produce a favorable configuration by chance. At minimum, the analysis should be repeated over several snapshots in the turbulent phase to support the persistence claim.","section":"Table 2 / §4"}],"minor_comments":[{"comment":"The text says the particle density field is normalized by the mean vertical particle density profile '(equation 16)', but the relevant definition appears to be Eq. (12) or the Gaussian fit Eq. (14). The equation reference should be corrected.","section":"§3.6"},{"comment":"The right panel of Fig. 8 is described as collapsing the RDFs, but for St_K = 0.01 and 0.02 the simulations do not resolve scales at or below ell_* and the S-curve inflection is not captured. The qualitative agreement with the cascade model is limited to the onset location. The text should more precisely state that the collapse is only partial and that the lowest-Stokes cases provide weaker evidence.","section":"§4.3"},{"comment":"The St_K = 0.01 and 0.02 runs use a 0.4H box with 4096^2 resolution, while the St_K = 0.04 run uses a 0.2H box with 2048^2 resolution. Box size and resolution differ across the three simulations; this complicates direct comparisons of clustering statistics and void scales. A brief comment on how this affects the claimed universality would be helpful.","section":"Table 2"},{"comment":"The empirical expression for sigma_SI in the epsilon < 1 regime is stated without derivation or error estimate. Given that the SI-vs-turbulence timescale comparison is central to ruling out SI as the driver, a reference to the source of this fit and its range of validity would strengthen the argument.","section":"Eq. (27b)"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses an important question and contains several genuinely direct diagnostics, but the headline St_omega 'collapse' is contradicted by the paper's own Table 4, and the single-snapshot axisymmetric setup leaves the generalized disk claim under-supported. I would be willing to reconsider after the stated range is corrected, time-averaged statistics are added, and the conclusions are more carefully scoped to axisymmetric simulations. The paper may also benefit from a slightly less categorical title if the 3D caveat remains."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The key thing to know: this is the first paper I know of that shows turbulent concentration (TC) actually operating in the self-generated turbulence of a settled midplane layer at St_K 0.01-0.04, and the core evidence is direct and persuasive. But the paper's headline number--the St_omega collapse to 0.4-0.7--is contradicted by its own Table 4, and that needs to be fixed before the quantitative claim is taken at face value.\n\nWhat's genuinely new and good: the strain-vorticity anti-correlation in Fig. 7 is about as direct a TC fingerprint as you can get in a disk context. The RDF analysis with the generalized cascade scaling (Eq. 25b) is a sensible way to handle the non-Kolmogorov spectrum, and the comparison with H17's S-curves is informative even though the low-St runs never resolve the St_ell=0.3 scale. I also credit the appendix that actually tests the traversal-time concern quantitatively rather than hand-waving it away. The paper is honestly hedged about the 3D question and about SymI.\n\nThe soft spots are real but not fatal to the qualitative conclusion. The stress-test note is right: Table 4 lists St_omega(peak)=0.35 for St_K=0.01, which sits outside the abstract's 0.4-0.7 range, and the Fig. 6 caption claims 0.4-0.75 while the table says 0.35. That's an internal inconsistency, not a typo, because the \"collapse\" is central to the universal clustering claim. Also, each run contributes one snapshot with no error bars, and the z_lim choice shifts the peaks noticeably. The SI rule-out depends on laminar linear growth rates applied to a clumpy turbulent flow; it's suggestive, not decisive. And the abstract's Roche-density reinterpretation outruns the more careful \"seems plausible\" wording in Sec. 6.\n\nMy take: the qualitative story is solid and worth taking seriously. The quantitative universality claim is oversold as stated. The authors need to reconcile Table 4 with the abstract/caption, provide multi-snapshot statistics or at least say plainly that these are single-snapshot values, and scale back the abstract to match what the paper actually demonstrates.\n\nFor a specialist in planetesimal formation or particle-turbulence interactions, this is a useful and genuinely new result. It deserves a serious referee--the issues are fixable with revision, not fatal. I'd send it to review.","headline":"Genuine and direct evidence for turbulent concentration in settled midplane layers, but the St_omega 'collapse' claim is contradicted by the paper's own Table 4 and needs revision before the quantitative story holds up.","tokens_in":26101,"tokens_out":2680,"would_cite":true,"duration_ms":26654,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper shows that small dust grains in settled midplane layers of planet-forming disks clump by turbulent concentration rather than by the streaming instability, with a universal effective Stokes number in the voids near 0.4–0.7.","keywords":["turbulent concentration","streaming instability","symmetric instability","protoplanetary disks","planetesimal formation","Stokes number","radial distribution function","particle clumping"],"falsifier":"Run the same three parameter sets in a fully 3D shearing-box simulation and measure the peak of the $\\mathrm{St}_\\omega$ distribution inside particle voids, the spectral slope of the gas kinetic energy, and the linear growth rate of the fastest axisymmetric SI mode. The central claim fails if the $\\mathrm{St}_\\omega$ peak does not fall near 0.4–0.7, if voids no longer contain coherent vorticity bounded by high-strain filaments, or if the fastest-growing SI mode grows faster than the eddy turnover frequency at $k_{\\mu\\epsilon}$.","tokens_in":24922,"feed_emoji":"🪐","tokens_out":9206,"duration_ms":86719,"temperature":0.7,"pith_summary":"This paper argues that in settled, weakly turbulent midplane layers of planet-forming disks, small dust grains spontaneously clump by turbulent concentration (TC): particles are flung out of coherent vortices and pile up in high-strain filaments. Working from axisymmetric shearing-box simulations with particle Stokes numbers $\\mathrm{St}_K = 0.01$–$0.04$, it finds that the effective Stokes number measured inside particle voids, $\\mathrm{St}_\\omega$, collapses to roughly 0.4–0.7, close to the value known to maximize TC clustering. A timescale comparison shows that streaming instability growth rates are one to two orders of magnitude slower than eddy turnover frequencies, so SI is not driving the turbulence there; the paper proposes the symmetric instability (SymI) as the driver. If true, TC is a persistent, self-generated clumping mechanism in exactly the regimes where planetesimals must form, and many small-scale Roche-density-exceeding clumps attributed to SI are actually TC operating on top of elevated filament densities.","feed_headline":"Turbulent concentration clumps grains, not streaming instability","feed_subtitle":"Small dust grains clump in settled disk layers even where the streaming instability cannot act.","key_machinery":"The analysis is carried by the cascade-model radial distribution function (RDF) of Hartlep et al. 2017, generalized to the non-Kolmogorov spectra of these axisymmetric runs. The scale-local Stokes number is written $\\mathrm{St}_\\ell = \\mathrm{St}_L (\\ell/L)^{(-3-n)/2}$ with $n \\approx -2$ instead of the Kolmogorov $-5/3$, and the clustering onset scale is $\\ell_* = \\mathrm{St}_L^{2/(3+n)} L$; rescaling separations by $\\ell_*$ collapses the measured RDFs onto the predicted S-curve. Complementary diagnostics are the void/filament grid classification ($X^{\\{vn\\}}$, $X^{\\{fn\\}}$) used to compute vorticity inside particle voids and the corresponding $\\mathrm{St}_\\omega = t_s \\omega$, plus timesca","core_discovery":"In axisymmetric (radial–vertical plane) shearing-box simulations of a settled particle layer, the authors observe that the particle field develops filament–void complexes: particle-rich filaments bound gas-only voids filled with coherent vorticity. Particles are expelled from vortex cores and collect in regions of high gas strain rate, the hallmark of TC. Across runs with $\\mathrm{St}_K = 0.01$, $0.02$, and $0.04$, the effective Stokes number $\\mathrm{St}_\\omega = t_s \\omega$ evaluated on particle-free void grid points peaks in the range 0.4–0.7, independent of $\\mathrm{St}_K$. The radial distribution functions, when separation is rescaled by $\\ell_* = \\mathrm{St}_L^{2/(3+n)} L$ with the mea","pith_inferences":["If the $\\mathrm{St}_\\omega$ collapse persists in fully 3D turbulence, TC could be identified in simulations by the strain-vorticity alignment of the gas field rather than by particle density alone, offering a cleaner observable diagnostic.","A testable extension is to run the same $\\mathrm{St}_K$ and metallicity parameters in a fully 3D shearing box: if the kinetic energy spectral slope reverts to Kolmogorov $-5/3$, the RDF collapse should use the standard scaling and the predicted clustering onset would shift accordingly.","The SymI attribution remains suggestive; a direct measurement of SymI eigenmode growth from the simulated mean azimuthal shear profile, or identification of its energy injection band, would settle whether SymI or a SymI–SI hybrid drives the turbulence.","If TC alone can push local densities past the Roche threshold in higher-resolution runs of these small Stokes numbers, planetesimal formation may not require SI for the smallest particles, changing the expected sizes and locations of first planetesimals."],"forward_implications":["In settled midplane layers with $\\mathrm{St}_K \\lesssim 0.04$ and midplane dust-to-gas ratio below unity, self-generated turbulence alone can produce persistent small-scale particle clumping via TC, independent of the streaming instability.","The effective void Stokes number $\\mathrm{St}_\\omega \\approx 0.4$–$0.7$ offers a scale-free diagnostic: whenever particles inside coherent vortex voids meet this condition, TC expulsion should be strongest regardless of the nominal $\\mathrm{St}_K$.","The streaming instability cannot be assumed to be the turbulence driver or the small-scale clumping mechanism in this parameter regime, because its linear growth is overwhelmed by eddy overturn at the relevant scales.","Numerical simulations seeking TC-driven Roche-density exceedance must resolve the $\\mathrm{St}_\\ell = O(1)$ scale; simulations that do not will under-report clustering, explaining resolution-dependent clumping results in the literature.","Small-scale density fluctuations exceeding the Roche density inside large-scale SI filaments are reinterpreted as TC acting on top of the locally elevated particle density, rather than as direct SI nonlinear structure."],"supporting_citations":[{"why":"Provides the axisymmetric simulations analyzed here and the identification of KHI, SymI, and SI as the midplane layer's turbulence mechanisms.","marker":"SU23"},{"why":"Supplies the cascade model for TC clustering whose S-curve RDF and $\\mathrm{St}_\\ell \\sim 0.3$ criterion are the comparison baseline.","marker":"H17"},{"why":"Supplies the methods for measuring particle scale height, $\\alpha$, $R'$, and $\\mathrm{St}_L$, including the correction that $\\Omega_L$ typically exceeds $\\Omega_K$.","marker":"S+24"},{"why":"Gives the linear streaming-instability growth rates used in the timescale comparison bounding $\\sigma_{\\rm SI}$.","marker":"Youdin & Goodman 2005"},{"why":"Provides the resonance criterion used to estimate $k_{\\mu\\epsilon}$, the fastest-growing SI radial wavenumber.","marker":"Squire & Hopkins 2018"},{"why":"Supplies the analytic particle-gas velocity difference model used in Appendix A to show particles remain coupled to the void vortices.","marker":"Cuzzi & Hogan 2003"},{"why":"Establishes the centrifuge mechanism by which inertial particles are expelled from vortex cores, the physical basis of TC.","marker":"Maxey 1987"},{"why":"Shows that maximum TC intermittency requires resolving specific scales, supporting the paper's resolution argument.","marker":"Hartlep & Cuzzi 2020"},{"why":"Reports no strong clumping for $\\mathrm{St}_K = 0.01$ at low metallicity, providing the resolution-dependent contrast the paper reinterprets.","marker":"Li & Youdin 2021"},{"why":"Reports Roche-exceeding small-scale clumps in higher-resolution SI-filament simulations, which the paper attributes to TC on top of elevated filament densities.","marker":"Lim et al. 2025b"}],"fun_headline_variants":["Turbulent concentration clumps dust, not streaming instability","Dust clumps arise from turbulence, not streaming instability","Vortices drive dust clumping in planet-forming disks","Dust clumping: turbulent concentration over streaming instability","Streaming instability not the driver: turbulent concentration is"],"cache_read_input_tokens":2816,"weakest_assumption_plain":"The conclusions assume that the axisymmetric, two-dimensional turbulent state seen in the simulations faithfully represents the real three-dimensional midplane layer; if 3D dynamics change the strain-vorticity statistics, shorten vortex lifetimes, or let SI filaments grow, the TC attribution and the $\\mathrm{St}_\\omega$ collapse lose their footing.","fun_headline_variants_meta":{"raw":{"variants":["Turbulent concentration clumps dust, not streaming instability","Dust clumps arise from turbulence, not streaming instability","Vortices drive dust clumping in planet-forming disks","Dust clumping: turbulent concentration over streaming instability","Streaming instability not the driver: turbulent concentration is"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000991,"raw_usage":{"total_tokens":4100,"prompt_tokens":872,"completion_tokens":3228,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":616,"completion_tokens_details":{"reasoning_tokens":3149}},"tokens_in":616,"tokens_out":3228,"duration_ms":23083,"temperature":1.0,"reasoning_tokens":3149,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T23:07:58.200386+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same three parameter sets in a fully 3D shearing-box simulation and measure the peak of the $\\mathrm{St}_\\omega$ distribution inside particle voids, the spectral slope of the gas kinetic energy, and the linear growth rate of the fastest axisymmetric SI mode. The central claim fails if the $\\mathrm{St}_\\omega$ peak does not fall near 0.4–0.7, if voids no longer contain coherent vorticity bounded by high-strain filaments, or if the fastest-growing SI mode grows faster than the eddy turnover frequency at $k_{\\mu\\epsilon}$.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the centrifuge mechanism by which inertial particles are expelled from vortex cores, the physical basis of TC."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that maximum TC intermittency requires resolving specific scales, supporting the paper's resolution argument."}],"review_version":1}