{"id":"954a033e-7b20-4ec4-bbf9-571874c44c9e","arxiv_id":"2505.23902","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"For τ_s=0.1 dust, stratified and unstratified streaming-instability simulations agree on midplane gas and dust dynamics across four pressure gradients, validating the unstratified approximation before strong clumping.","lead":"The authors ran 2D simulations of the streaming instability in disks both with and without vertical gravity, across four pressure-gradient strengths, for dust grains with stopping time τ_s=0.1. They found that the cheaper gravity-free simulations reproduce the midplane dust and gas behavior of the full stratified disks until strong clumping begins, which supports their continued use in planet-formation studies.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Fiducial L_x = 4ΠH is below the paper's own convergence requirement; COM-frame midplane comparison may be affected by unconverged vertical undulations.","rationale":"I read the paper as claiming a quantitative bridge for τ_s = 0.1: saturated midplane dynamics in stratified disks are statistically indistinguishable from unstratified AB runs before strong clumping. The evidence is substantial: dispersion scaling (Figure 9), CDFs (Figure 10), diffusion coefficients (Figure 11), and power spectra (Figure 16) all match or track B24. I also credit the explicit caveats and the Appendix A convergence tests. The weakest point is not the physics but the domain size: the authors themselves state L_x > 4ΠH is required. Their rebuttal is plausible because the quantities they rely on for the comparison—dispersions and density distributions—are converged across L_x. But the mean vertical flows and COM velocities that carry the interpretive weight in Section 3.4 are not checked at L_x = 8ΠH, and the L_x = 4ΠH runs show measurable asymmetries (nonzero gas u_z at the midplane, shifted z̄). Since the central assertion is precisely that stratification's main effect is vertical momentum redistribution and that unstratified models capture everything else, leaving the redistribution itself at an unconverged domain size is a load-bearing gap. It does not overturn the paper; it makes the conditional verdict appropriate. I therefore keep the reader's conditional disposition, so the verdict is unchanged.","tokens_in":31378,"tokens_out":7941,"duration_ms":74944,"concrete_test":"Re-run the Π = 0.05 stratified case with L_x = 8ΠH (and, if computationally feasible, Π = 0.01 and Π = 0.1) using the same L_z = 8ΠH and N_ΠH = 256, and recompute over t/T = 10–50 the COM velocities U_x,CM(z) and U_y,CM(z) (Equations 10–11), the midplane dispersions in Table 2, and the time-averaged mean dust position z̄. If the midplane COM radial velocity or any dispersion shifts by more than the quoted 1σ temporal variability, or if agreement with B24's AB models worsens, then L_x = 4ΠH is inadequate for the central claim; if these quantities remain within 1σ, the concern is retired.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing weakness is the fiducial radial domain size. The paper's own convergence test concludes that L_x > 4ΠH is required to adequately resolve the vertical undulations of the dust layer (Appendix A, Figures A.3–A.4), yet all main runs use L_x = 4ΠH. In these runs the gas vertical velocity crosses zero away from the midplane and the time-averaged mean dust position is shifted by ≈0.02ΠH, while the L_x = 8ΠH test has zero crossing at z = 0 and z̄ ≈ 8×10⁻⁴ΠH. The authors argue the midplane comparison is unaffected because density and velocity dispersions and dust density distributions converge (Figure A.1). However, those converged diagnostics are not sensitive to the mean/vertical-flow asymmetry that enters the COM velocities and vertical momentum redistribution, which is exactly the mechanism invoked to explain why unstratified and stratified simulations agree (Section 3.4). No L_x = 8ΠH COM-frame analysis is shown, so the central claim that unstratified simulations reliably predict midplane dynamics rests on a numerical ingredient the paper itself flags as not converged. This is a genuine but resolvable gap; the paper is transparent about it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents 2D axisymmetric, vertically stratified shearing-box simulations of the streaming instability for τ_s=0.1, with four radial pressure gradients Π=0.01–0.1 and dust surface density ratios Z chosen to keep Z/Π=0.24 and a midplane dust-to-gas density ratio ϵ≈1. The authors compare the saturated state (t/T=10–50) with the unstratified AB models of B24 and report close agreement in midplane velocity dispersions, density distributions, dust diffusion coefficients, and COM-frame velocities. They argue that vertical gravity redistributes momentum vertically but does not change midplane SI dynamics before strong clumping.","tokens_in":31644,"tokens_out":5382,"duration_ms":54009,"significance":"If correct, the result is useful: it gives quantitative evidence that vertically unstratified SI simulations can approximate midplane dynamics in stratified disks for at least one stopping time, and it clarifies the role of vertical gravity and the COM frame. The paper is transparent and thorough in its diagnostics, includes explicit caveats and convergence tests in Appendix A, and uses a well-documented comparison to B24. The scope is narrow (one τ_s, 2D, fixed Z/Π), and the domain-size caveat below tempers the strength of the central claim.","major_comments":[{"comment":"The fiducial runs use L_x=4ΠH, yet Appendix A concludes that L_x>4ΠH is required to adequately resolve the vertical undulations of the dust layer. For L_x=4ΠH, the gas vertical velocity crosses zero away from the midplane and the mean particle position is shifted by ≈0.02ΠH, whereas the L_x=8ΠH run has zero crossing at z=0 and z̄≈8×10^-4ΠH. The authors state this does not affect the midplane comparison because density and velocity dispersions and density distributions converge (Figure A.1), but those diagnostics are not the ones used for the central COM-frame argument in Figures 5–6. Since the mechanism invoked to explain the agreement is vertical momentum redistribution, the COM-frame quantities (Equations 10–11) should be verified at L_x=8ΠH, or an explicit argument connecting the converged diagnostics to the COM-frame result must be given.","section":"Appendix A, Figures A.3–A.4; Section 3.4"},{"comment":"The chosen saturated-state window t/T=10–50 is used for all four runs, but the Π=0.01 run begins a sharp density increase at t/T≈50 and reaches a second saturated state with a maximum density of ≈300ρ_g0. The time averages in Figures 8–9 and Table 2 for Π=0.01 may therefore be contaminated by the approach to strong clumping. Please show that the reported agreement is robust to using a shorter window, e.g., t/T=10–30, or otherwise justify the window for that run.","section":"Section 3.2, Figure 3; Section 4.2"}],"minor_comments":[{"comment":"The phrase \"unstratified simulations represents well\" should be \"unstratified simulations represent well,\" and similar grammatical fixes are needed in a few other places.","section":"Abstract; Summary"},{"comment":"There is a typo in \"c_s is isotermal sound speed\" — should be \"isothermal.\"","section":"Section 2"},{"comment":"The text contains \"the the axisymmetric KHI\" and \"possiblity\"; these should be corrected.","section":"Section 4.1"},{"comment":"The sentence beginning \"The equation is the analytical solution...\" is ambiguous about which equation it refers to; please clarify that it refers to Equation (16).","section":"Section 3.7"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest and well-structured, and the main gap is fixable: an L_x=8ΠH COM-frame analysis (or an explicit diagnostic linking the converged dispersions to the COM-frame result) should be requested in revision. No concerns about novelty or attribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The punchline: this is a solid, quantitative numerical comparison that mostly delivers on its central claim. For τ_s=0.1, the saturated-state midplane dispersions, COM-frame velocities, density distributions, and diffusion coefficients in stratified runs line up with B24's unstratified runs across Π=0.01–0.1. That is a genuinely useful result because unstratified boxes are far cheaper, and the paper is the first to make this comparison systematically and across pressure gradients.\n\nWhat is new: four new stratified runs with Z/Π held fixed, COM-frame analysis that isolates the role of vertical gravity in redistributing radial momentum, Richardson-number checks for shear instabilities during settling and saturation, and power-law scaling comparisons with B24. The analysis is careful: multiple diagnostics, temporal variability, convergence tests for grid resolution and domain size, explicit caveats about 2D axisymmetry, single τ_s, and the late-time strong clumping in the Π=0.01 run. The circularity concern is minor: matching Z to reach ϵ≈1 is a sensible design choice, and the quantities compared (dispersions, distributions, diffusion coefficients) are not fixed by that choice and could have disagreed.\n\nThe important soft spot is the radial domain size. All main runs use L_x=4ΠH, but the paper's own Appendix A concludes that L_x>4ΠH is required to adequately resolve the vertical undulations of the dust layer. At 4ΠH the gas vertical velocity crosses zero away from the midplane and the time-averaged mean dust position is shifted by ~0.02ΠH, whereas at 8ΠH both are essentially zero. The authors argue this does not affect the midplane comparison because density and velocity dispersions and dust density distributions converge across L_x. That argument is reasonable for those diagnostics, but the COM-frame velocities — exactly the quantities used to show that vertical gravity is the only extra ingredient — are sensitive to the mean vertical flow asymmetry. No L_x=8ΠH COM-frame analysis is shown. This is a genuine but resolvable gap, and the paper is transparent about it. A referee should ask for either the 8ΠH COM-frame panel or a more explicit defense of why the mean-flow shift cannot bias the central comparison.\n\nMinor things: some factor-of-two differences in vertical diffusion coefficients are acknowledged and discussed; the abstract has a grammar slip ('represents well'); the comparison to B24 uses different domain sizes and resolutions, but the convergence tests make that acceptable. None of this undercuts the main conclusion.\n\nWho it is for: anyone using unstratified SI simulations as a proxy for midplane dynamics, and anyone planning large parameter surveys. It deserves a serious referee. My recommendation: send it to review, with a request to address the L_x=8ΠH COM-frame check.","headline":"A careful, quantitative numerical bridge for τ_s=0.1 SI; the main caveat is that the fiducial radial box is smaller than the paper's own convergence criterion, and the missing L_x=8ΠH COM-frame check should be addressed before this becomes the standard reference.","tokens_in":32195,"tokens_out":2423,"would_cite":true,"duration_ms":23231,"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":"For $\\tau_s=0.1$, stratified and unstratified streaming-instability simulations agree at the midplane before strong clumping begins.","keywords":["streaming instability","planetesimal formation","protoplanetary disks","dust-gas dynamics","stratified simulations","unstratified simulations","pressure gradient","dust settling"],"falsifier":"Re-run the $\\Pi=0.05$ comparison at $L_x=8\\Pi H$ (the size used in Appendix A) and recompute the center-of-mass midplane radial and azimuthal velocities; if they move by more than the quoted uncertainties relative to the $L_x=4\\Pi H$ run, the claim that the fiducial box captures the midplane dynamics would fail.","tokens_in":31195,"feed_emoji":"🪐","tokens_out":9580,"duration_ms":83490,"temperature":0.7,"pith_summary":"This paper tackles a standing question in planetesimal-formation simulations: can the computationally cheap setup that ignores vertical gravity and stratification reproduce the behavior of a real stratified disk? For grains with stopping time $\\tau_s=0.1$, sampled at four radial pressure-gradient strengths, the authors say yes near the midplane before strong clumping begins. They find that vertical gravity mainly redistributes momentum vertically, and that once gas and dust velocities are measured in the center-of-mass frame, midplane velocities, density and velocity dispersions, density distributions, and dust diffusion coefficients agree with the unstratified simulations. The practical payoff is that unstratified runs can act as a reliable predictor of streaming-instability-driven dynamics at the midplane, while stratified runs add settling and shear-driven instabilities on top of the same local turbulence.","feed_headline":"For τ_s=0.1, unstratified runs match midplane dust-gas flow","feed_subtitle":"Vertical gravity only shifts momentum; local streaming-instability turbulence is the same.","key_machinery":"The central object is the density-weighted center-of-mass velocity profile $U_{\\mathrm{CM}}(z)=(\\rho_g \\mathbf{u}+\\rho_p\\mathbf{v})/(\\rho_g+\\rho_p)$ and its vertical structure. By comparing gas and dust velocities in the COM frame rather than the lab frame, the analysis removes the vertical-gravity-driven mean flow and isolates the fluctuation part of the streaming turbulence, which is what matches the unstratified models at the midplane. The supporting diagnostics are the radial and azimuthal Richardson numbers, used to check the Kelvin-Helmholtz and symmetric instabilities; the kinetic-energy power spectrum, which locates the SI at wavenumber $k_x\\Pi H/(2\\pi)\\sim1$; and measured radial and vertical dust diffusion coefficients, which follow $\\propto\\Pi^2$.","core_discovery":"The central claim is that for $\\tau_s=0.1$ and a fixed ratio $Z/\\Pi=0.24$ of dust surface density to pressure-gradient parameter, the saturated streaming instability at the midplane of a vertically stratified disk is the same physical phenomenon as in a vertically unstratified box across $\\Pi = 0.01,0.02,0.05,0.1$. The paper supports this by comparing stratified runs to the unstratified AB models: dust filaments form during settling with similar morphologies; vertical gravity creates a center-of-mass radial velocity and vertical gradients that peak near $\\pm H_p$; subtracting that mean flow brings midplane radial velocities into agreement, while azimuthal velocities differ by only $2{-}3\\%$ of $\\Pi c_s$; density and velocity dispersions match inside one dust scale height; and the dust density distributions and radial diffusion coefficients are nearly identical, with vertical diffusion somewhat weaker because gravity constrains the vertical random walk. The authors conclude that unstratified simulations represent midplane dust-gas dynamics well before strong clumping, and that the streaming turbulence in stratified disks is fundamentally the same as in unstratified disks for these parameters.","pith_inferences":["The authors' own Appendix A indicates that $L_x>4\\Pi H$ is required to capture the dust layer's vertical undulations; a natural next check, which they did not report, is the COM-frame midplane comparison at $L_x=8\\Pi H$.","If the agreement holds for other stopping times and in three dimensions, unstratified simulations could serve as calibrated subgrid sources for pebble diffusion and momentum feedback in larger disk-evolution calculations.","The $\\Pi=0.01$ late clumping suggests that, despite equal $Z/\\Pi$, weaker pressure gradients may be especially conducive to planetesimal formation; longer runs with larger domains would test whether $Z_{\\mathrm{crit}}/\\Pi$ loses predictive power."],"forward_implications":["For $\\tau_s=0.1$, unstratified simulations give a quantitatively reliable estimate of midplane velocity dispersions, density distributions, and diffusion coefficients in stratified disks, so cheaper unstratified runs can be used for these diagnostics.","The main stratification effect is vertical momentum redistribution; after removing the COM motion, the residual streaming turbulence is statistically indistinguishable from the unstratified case.","During the saturated state both Richardson numbers stay above their thresholds, indicating that the KHI and SymI are not active and the SI dominates turbulence in the explored parameter range.","The late strong clumping in the $\\Pi=0.01$ run, despite the same $Z/\\Pi$ as the other runs, indicates that $Z/\\Pi$ alone is not a sufficient predictor of strong clumping."],"supporting_citations":[{"why":"Provides the linear streaming instability and its characteristic scale ~$\\Pi H$ that the paper takes as the baseline.","marker":"A. N. Youdin & J. Goodman 2005"},{"why":"Provides the unstratified AB models that are the direct comparison set for all midplane diagnostics.","marker":"B24"},{"why":"Gives the NSH equilibrium used for the initial horizontal velocities and the theoretical COM azimuthal profile.","marker":"Y. Nakagawa et al. 1986"},{"why":"Motivates fixing $Z/\\Pi$ constant to keep the midplane density ratio near unity across $\\Pi$.","marker":"M. Sekiya & I. K. Onishi 2018"},{"why":"Supplies the Richardson-number criteria and the KHI/SymI interpretation used to explain stratification-specific turbulence.","marker":"D. Sengupta & O. M. Umurhan 2023"},{"why":"Supplies the strong-clumping threshold and prior stratified simulations that the late clumping in the $\\Pi=0.01$ run is measured against.","marker":"R. Li & A. N. Youdin 2021"},{"why":"Provides earlier unstratified saturated-state dynamics and diffusion measurements that this work extends.","marker":"A. Johansen & A. Youdin 2007"},{"why":"Provides the particle-gas coupling method and the particle-resolution guidance used in the numerical setup.","marker":"X.-N. Bai & J. M. Stone (2010b)"}],"fun_headline_variants":["Unstratified and stratified SI agree for tau_s=0.1 grains","Streaming instability: stratified runs match unstratified midplane","For tau_s=0.1, midplane flow is same with or without gravity","SI bridging: stratified disks match unstratified boxes at tau_s=0.1","tau_s=0.1 SI: vertical gravity doesn't change midplane dynamics"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the fiducial radial box width $L_x=4\\Pi H$ is wide enough not to bias the midplane comparison, even though the paper's own convergence test found that $L_x>4\\Pi H$ is required to capture the dust layer's vertical undulations; if that width changes the mean vertical dust position or the COM velocities, the central claim would be affected.","fun_headline_variants_meta":{"raw":{"variants":["Unstratified and stratified SI agree for tau_s=0.1 grains","Streaming instability: stratified runs match unstratified midplane","For tau_s=0.1, midplane flow is same with or without gravity","SI bridging: stratified disks match unstratified boxes at tau_s=0.1","tau_s=0.1 SI: vertical gravity doesn't change midplane dynamics"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000864,"raw_usage":{"total_tokens":3816,"prompt_tokens":1084,"completion_tokens":2732,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":700,"completion_tokens_details":{"reasoning_tokens":2629}},"tokens_in":700,"tokens_out":2732,"duration_ms":19375,"temperature":1.0,"reasoning_tokens":2629,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T12:38:56.713125+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-run the $\\Pi=0.05$ comparison at $L_x=8\\Pi H$ (the size used in Appendix A) and recompute the center-of-mass midplane radial and azimuthal velocities; if they move by more than the quoted uncertainties relative to the $L_x=4\\Pi H$ run, the claim that the fiducial box captures the midplane dynamics would fail.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Motivates fixing $Z/\\Pi$ constant to keep the midplane density ratio near unity across $\\Pi$."}],"review_version":1}