{"id":"3d5a4220-7b2a-4119-8de1-5703925248b0","arxiv_id":"2506.13592","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Using global 3D two-fluid simulations, the authors report the first numerical evidence that planetary waves from a low-mass planet can trigger resonant drag instabilities and produce filamentary dust structures in a protoplanetary disc.","lead":"Simulations of a planet embedded in a dusty gas disc show that planetary waves, the Rossby waves launched near the planet's horseshoe region, can couple with the streaming motion of dust and gas and drive a resonant drag instability. The resulting filamentary dust structures appear within a few orbits and may offer a new route to dust clumping and planetesimal formation around growing planets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The RDI diagnosis is never tested: no dispersion relation or resonance condition k·w_s = ω_gas is measured, so the observed stripes could be a forced or advected response rather than a resonant instability.","rationale":"The reader's weakest assumption identifies exactly the same gap: the paper does not compute the wave dispersion relation or the resonance condition, so the RDI interpretation is inferred rather than demonstrated. I agree with the conditional verdict because the feedback/no-feedback comparison, the finite-dust-layer experiment, and the softened-potential run in Appendix B provide meaningful controls that support the involvement of dust back-reaction. However, those controls cannot distinguish a genuinely resonant instability from a forced, feedback-amplified response, which is the central claim of the paper. A convergence test would address numerical concerns, but it would not settle the mechanism; the missing resonance check is the more load-bearing issue. The proposed Fourier-based test is feasible with the existing simulations and would either confirm the RDI identification or require the conclusions to be weakened. Since the concern is addressable but unresolved, the reader's CONDITIONAL verdict is appropriate, and my assessment does not move it.","tokens_in":16064,"tokens_out":4271,"duration_ms":51416,"concrete_test":"In the gas-only RDISt0Mp3 run, extract the time series of the vertical vorticity (or radial velocity) along a ring just outside the horseshoe separatrices and perform a spatiotemporal Fourier transform to obtain the planetary-wave dispersion branch ω_gas(k_phi), with frequency resolution set by the 20-orbit run length. In the matching dusty run RDISt0_2Mp3, measure the local dust-gas relative velocity w_s in the same region, averaged over azimuth and over a few orbits. Then take the dominant stripe wavenumber k from a two-dimensional Fourier transform of the dust-density map and test whether |k·w_s − ω_gas(k)| is smaller than the numerical linewidth. If no branch satisfies the resonance condition, the stripes are not a resonant drag instability but a forced or sheared response.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the stripe-like disturbances along the horseshoe separatrices are intrinsic planetary-wave modes whose frequency satisfies the RDI resonance condition k·w_s = ω_gas(k). The paper never measures ω_gas(k) from the simulations; the identification rests on visual similarity with Masset & Benítez-Llambay (2016) and on the morphological resemblance of the dust filaments to SI. Section 3.6 states that the pattern arises from a similarity between the propagation velocity of the gas vorticity waves and the dust-gas relative velocity components, but no quantitative comparison is presented. The Appendix A experiment does not close this gap: Eq. (A1) injects a Gaussian vortex-like velocity perturbation, and its shearing into two stripes is expected in Keplerian shear even without any resonant wave. Furthermore, an instability requires exponential growth of a mode, yet no growth-rate measurement or mode-amplitude time series is reported; disturbances spreading spatially at finite amplitude can be a linear response to the planet's potential in the presence of dust feedback. The no-feedback and finite-dust-layer controls usefully show that dust back-reaction is involved, but they do not establish that the resonance condition is met. Since the abstract claims the first numerical evidence of RDIs driven by planetary waves, the missing dispersion-relation check is the load-bearing gap.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports 3D high-resolution two-fluid simulations of a dusty protoplanetary disc with an embedded low-mass planet (0.3–3 Earth masses) and dust Stokes numbers in 0.01–0.5. The central claim is that planetary waves (also called Rossby waves) propagating along the horseshoe-region separatrix resonate with the dust-gas relative streaming motion and thereby activate a resonant drag instability (RDI), producing global filamentary dust structures. The authors also claim that dust feedback enables buoyancy resonances in an otherwise isothermal disc. The evidence consists of: gas vorticity maps showing planetary-wave-like disturbances and stripes, dust-density maps showing filamentary structures when dust feedback is included, a suite of control runs (no feedback, finite-thickness dust layer, softened potential, mass tapering), and a 2D experiment with an injected Gaussian perturbation. The paper does not report a quantitative test of the RDI resonance condition, nor measurements of wave dispersion or growth rates.","tokens_in":16290,"tokens_out":6353,"duration_ms":63061,"significance":"If validated, the claim that low-mass planets can trigger RDIs via planetary waves would be an important new pathway for dust clumping and planetesimal formation, extending RDI theory to non-axisymmetric, global disc settings. The simulations are ambitious and carefully parameterized, and the inclusion of both-feedback and no-feedback controls, a finite-thickness dust layer, potential softening, and mass tapering is a strength. The use of the public FARGO3D code and the high resolution are additional positive features. However, the absence of any direct check of the RDI resonance condition and the lack of growth-rate measurements mean that the principal conclusion is not yet quantitatively established.","major_comments":[{"comment":"The RDI resonance condition k·w_s = ω_gas(k) stated in the Introduction is never evaluated using the simulation data. The paper does not measure the frequencies or wavenumbers of the planetary waves, nor does it compute the local dust-gas relative velocity w_s. Section 3.6 merely states that the filaments arise from a 'similarity' between the gas-wave propagation velocity and the dust-gas relative velocity, without any quantitative comparison. As a result, the identification of the observed structures as an RDI is inferred from morphology and from the presence of dust feedback, which are necessary but not sufficient conditions. The disturbance could be a non-resonant forced response of the dust to the planet-induced gas flow modified by feedback. This is the central gap in the evidence for the paper's main claim.","section":"Section 1 and Section 3.6"},{"comment":"No growth-rate measurement or mode-amplitude time series is presented. An instability requires exponential growth of a mode amplitude, but the paper shows only spatial spreading of finite-amplitude disturbances over 5–20 orbits. The 'rapid onset' of the structures could equally be interpreted as a fast-propagating linear response. A quantitative growth-rate analysis (e.g., tracking Fourier-mode amplitudes in time) is essential to distinguish a genuine instability from an advected or forced pattern.","section":"Sections 3.3, 5 (Figs. 5, 6, 9)"},{"comment":"The buoyancy-resonance interpretation of the colored rays is circular as presented. The overlay in Fig. 4 uses z=2H_g in Eq. (17), but footnote 4 states that the perturbations do not appear for values of z at which the dust has settled, and no independent determination of z is given. Choosing the height to match the observed pattern means Eq. (17) is not being used as a predictive test. This substantially weakens the secondary claim that dust feedback enables buoyancy resonances in a globally isothermal disc.","section":"Section 3.4, Eq. (17), footnote 4"},{"comment":"The two-dimensional experiment in Appendix A is not a valid test of RDI activation. Injecting a localized Gaussian velocity perturbation into a Keplerian disc will inevitably shear into two inclined stripes, regardless of any resonance. This experiment demonstrates only that a finite-amplitude perturbation can produce filamentary dust density variations, not that a resonant drag instability is at work. It therefore does not close the gap left by the absence of a resonance-condition check in the main simulations.","section":"Appendix A, Eq. (A1)"}],"minor_comments":[{"comment":"There are several typographical errors: 'deveolps' (Introduction), 'mow' (Section 4.3), 'Once can see' (Section 3.3), 'bouyancy' (Sections 3.4, 3.5), 'adibatic' (Section 5), and the rendering of 'sin' as 's i n' in Eq. (A1). These should be corrected.","section":"Throughout"},{"comment":"The claim that the instability 'does not depend on the dust-to-gas mass ratio' is too strong given that only two values are tested (ε=0.01 in the main simulations and ε=0.001 in the finite-layer run, which also has a different Stokes number and dust scale-height setup). A systematic variation of ε would be needed to support this conclusion.","section":"Section 4.2"},{"comment":"The identification of the waves as 'planetary waves (also known as Rossby waves)' is not justified in detail; the disturbances shown in Fig. 1 are localized to the horseshoe region and may be a different mode of vortensity wave. The authors should clarify the relationship to the classical Rossby-wave literature or temper the terminology.","section":"Section 3.2 and Fig. 1"},{"comment":"The description of the computational mesh reports (N_r, N_θ, N_φ) = (3200, 100, 12560), which gives over 4×10^9 cells. While this is stated as high resolution, the authors do not discuss how this extreme resolution is handled computationally or whether any grid-convergence tests were performed beyond the softening/taper experiments in Appendix B.","section":"Section 2.3"}],"recommendation":"major_revision","confidential_remarks":"The core idea is interesting and the simulations are extensive, but the manuscript currently makes a strong physical claim (first numerical evidence of RDI driven by planetary waves) without providing the quantitative analysis that would substantiate it. I would advise the editor to ask the authors to either (i) directly measure the wave dispersion relation and the dust-gas drift in the simulations and check the resonance condition, including a growth-rate measurement, or (ii) substantially soften the claims to describe the observed filamentary dust structures as a dust-feedback-driven response to planet-induced gas flows, leaving the RDI interpretation as a hypothesis. The former would make the paper much stronger; the latter would avoid overstating the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You should know this paper has a good core idea: planetary waves launched by a low-mass planet along the horseshoe separatrices might satisfy the RDI resonance condition and drive dust clumping. The authors support it with a set of control simulations that are genuinely thoughtful. They show filamentary dust structures appear only when dust back-reaction is on, that they survive changes in softening and mass tapering, that a finite-thickness dust layer with St=0.01 and low dust-to-gas ratio still produces the same pattern, and that a planet-free run with forced white noise does not. Those controls rule out several mundane alternatives, and that is the paper's real strength.\n\nThe new thing here is the non-axisymmetric trigger: previous RDI work in discs is mostly axisymmetric, and this is, as far as I know, the first numerical evidence that a planet's waves can excite an RDI-like instability. The runs span a useful range of Stokes numbers and planetary masses, and the dust-to-gas ratio stays well below unity where the instability appears, which matters for planetesimal formation.\n\nNow the soft spots, and they are not minor. The paper never verifies the RDI resonance condition. It never measures the dispersion relation omega_gas(k) of the planetary waves, never computes k·w_s, and never reports a mode growth rate or amplitude versus time. All we see is spatial spreading of stripes, which can happen as a linear response to the planet's potential even without a resonant instability. The buoyancy-ray overlay in Section 3.4 uses z=2Hg, and the text says the positions 'match,' but no quantitative metric is given; that parameter looks picked to fit the simulation. Appendix A does not close the gap: injecting a Gaussian vortex into a Keplerian shear flow will shear into two stripes whether or not any resonant wave exists. The grid is also extraordinary—3200x100x12560 cells—and there is no convergence test or code version reported. Input files are only on request, which is a reproducibility weakness, though the code itself is public.\n\nNone of this makes me think the paper is wrong. The central claim could well be true, and the controls are better than most papers at this level. But the claim 'resonate with the streaming motion and trigger the RDI' is asserted, not demonstrated. A referee can fix this: ask for a dispersion-relation measurement, a resonance-condition check, and a growth-rate curve. If those come out positive, this becomes a solid result.\n\nWho gets value from it: people working on planet-disc interactions, dust dynamics, and planetesimal formation. It deserves serious peer review, not a desk reject, and it should go back for major revision with the quantitative RDI diagnosis as the required addition.","headline":"Plausible but not proven: the paper isolates dust-feedback-driven filament formation with well-designed controls, but it never measures the RDI resonance condition or a growth rate, so the central mechanism remains an inference.","tokens_in":16845,"tokens_out":1705,"would_cite":false,"duration_ms":18819,"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":"Planetary (Rossby) waves launched by a low-mass planet can resonate with the dust–gas streaming motion and trigger a resonant drag instability that produces global filaments in a protoplanetary disc.","keywords":["resonant drag instability","planetary waves","Rossby waves","protoplanetary discs","dust-gas two-fluid dynamics","streaming instability","planet-disc interaction","buoyancy resonance"],"falsifier":"Take the gas-only simulations, Fourier-transform the vorticity perturbations along the downstream separatrices to measure the wave frequency $\\omega_{\\mathrm{gas}}(\\mathbf{k})$, and compare it with the dust–gas drift $\\mathbf{w}_s$ measured in the dusty runs; if no wavenumber in the filament-forming region satisfies $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$, the RDI resonance is not what drives the instability.","tokens_in":15817,"feed_emoji":"🪐","tokens_out":6523,"duration_ms":67190,"temperature":0.7,"pith_summary":"The paper sets out to show that planetary (Rossby) waves, launched by a low-mass planet in a protoplanetary disc, can act as the gas wave mode that the resonant drag instability (RDI) needs. Using global 3D two-fluid simulations of an inviscid, globally isothermal disc with 0.3–3 Earth-mass planets and dust with Stokes number 0.01–0.5, the authors find filamentary dust structures that grow on dynamical timescales and spread beyond the planet's horseshoe region. They argue that because these structures appear only when dust feeds back on gas, and because they resemble streaming-instability filaments, the waves must be resonating with the dust–gas streaming motion. If true, planet-induced RDIs give a pathway to dust concentration and planetesimal formation that works at low dust-to-gas ratios and does not depend strongly on grain size or vertical settling.","feed_headline":"Planet-driven waves ignite resonant drag instability in discs","feed_subtitle":"3D simulations grow dust filaments within a few orbits, even at low dust-to-gas ratios—a fast route to planetesimals.","key_machinery":"The load-bearing machinery is the resonant drag instability condition together with the wave source that satisfies it. The RDI occurs when a dust–gas streaming velocity $\\mathbf{w}_s$ obeys $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$, meaning that dust drifting through gas resonantly drives an intrinsic gas wave at wavenumber $\\mathbf{k}$. In this disc the intrinsic waves are planetary (Rossby) waves propagating along the downstream separatrices of the horseshoe region, excited by the vortensity structure the planet sets up. Dust feedback is what couples the streaming motion to those waves, and the inclined 'pw-stripes' seen in the vorticity and dust density are the visible signature of the resonance being driven. Vertical dust settling is present but is not the trigger; the planet's waves are.","core_discovery":"On the paper's own terms, the discovery is that the vortensity disturbances previously seen in the horseshoe region of a planet in a 3D isothermal disc—stripe-like vertical-oscillation features and the planetary waves propagating along downstream separatrices—become the seed of a resonant drag instability once dust aerodynamic feedback is included. The relevant waves are planetary (Rossby) waves propagating along the separatrices; the dust–gas drift provides the streaming motion; and the resonance condition $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$ turns the two-fluid mixture unstable. The authors report global, streaming-instability-like filamentary dust concentrations that develop within about five to twenty orbits, spread beyond the horseshoe region, and appear across the full range of Stokes numbers and planet masses tried, even when the midplane dust-to-gas ratio stays far below unity. They also report buoyancy-resonance ray patterns in the vertical velocity that emerge in a globally isothermal disc only because dust feedback is present. Their conclusion is that this is the first numerical evidence of RDI activation driven by planetary waves.","pith_inferences":["The paper leaves open a direct spectral test: measuring the wave dispersion relation in the simulations and checking that some wavenumber satisfies $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$ in the region where filaments first appear would close the gap between morphology and mechanism.","The mechanism need not be limited to planets: any localised vortensity source that launches Rossby waves, such as a gap edge or an eccentric vortex, could seed the same RDI, widening the relevance to observed rings and asymmetries.","If real discs support these waves, the instability could erase the asymmetric dust signatures expected near low-mass planets within tens of orbits; searching for low-contrast extended filaments beyond the horseshoe region at low dust-to-gas ratios would be a concrete observational test."],"forward_implications":["Low-mass planets of 0.3–3 $M_\\oplus$ can trigger dust clumping on dynamical timescales, so planetesimal formation may begin near an embryo well before classical streaming instability would act.","The instability operates at midplane dust-to-gas ratios $\\epsilon\\ll1$, so dust-rich regions are not required for it to start; a wave source is the key ingredient.","Because growth is fast and weakly dependent on Stokes number, it can organise a broad grain-size spectrum into coherent filaments and suppresses the planet-localised dust voids and asymmetric structures that dust-only dynamics would produce.","In a globally isothermal disc, dust feedback mimics buoyancy resonances normally associated with adiabatic discs, so wave patterns in isothermal dusty discs cannot be read using a purely gas adiabatic model."],"supporting_citations":[{"why":"Establishes the RDI resonance condition $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$ that the paper invokes as the activation mechanism.","marker":"Squire & Hopkins 2018a"},{"why":"Places the streaming instability within the RDI family and supplies the taxonomy used to identify the filamentary patterns as SI-like.","marker":"Squire & Hopkins 2020"},{"why":"Defines the streaming instability whose filamentary dust concentrations the paper compares its results to.","marker":"Youdin & Goodman 2005"},{"why":"Found the stripe-like vortensity patterns and planetary-wave propagation in a 3D isothermal disc that this paper extends by adding dust feedback.","marker":"Masset & Benítez-Llambay 2016"},{"why":"Provides the analogy between isothermal dusty discs and adiabatic gas discs used to interpret the buoyancy-resonance rays.","marker":"Lin & Youdin 2017"},{"why":"Supplies the SI stability criterion used in the finite-thickness dust-layer model to rule out classical SI and DSI.","marker":"Chen & Lin 2020"}],"fun_headline_variants":["Planetary waves ignite resonant drag instability in dust-gas discs","Rossby waves trigger resonant drag instabilities that build planetesimals","Dust feedback lets planetary waves ignite drag instabilities","Planetary waves trigger fast dust filaments via RDI"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the wave-like disturbances travelling along the horseshoe separatrices are genuine planetary (Rossby) waves whose frequency can satisfy the RDI resonance condition $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$ against the local dust–gas drift; the paper infers this from morphology and from the dependence on dust feedback rather than measuring the resonance directly.","fun_headline_variants_meta":{"raw":{"variants":["Planetary waves ignite resonant drag instability in dust-gas discs","Rossby waves trigger resonant drag instabilities that build planetesimals","Dust feedback lets planetary waves ignite drag instabilities","Planetary waves trigger fast dust filaments via RDI"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000692,"raw_usage":{"total_tokens":3200,"prompt_tokens":1081,"completion_tokens":2119,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":697,"completion_tokens_details":{"reasoning_tokens":2052}},"tokens_in":697,"tokens_out":2119,"duration_ms":16898,"temperature":1.0,"reasoning_tokens":2052,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:58:47.115904+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the gas-only simulations, Fourier-transform the vorticity perturbations along the downstream separatrices to measure the wave frequency $\\omega_{\\mathrm{gas}}(\\mathbf{k})$, and compare it with the dust–gas drift $\\mathbf{w}_s$ measured in the dusty runs; if no wavenumber in the filament-forming region satisfies $\\mathbf{k}\\cdot\\mathbf{w}_s=\\omega_{\\mathrm{gas}}(\\mathbf{k})$, the RDI resonance is not what drives the instability.","supporting_citations":[],"review_version":2}