{"id":"b734a0c9-25ed-472a-bfc0-2573e5ad79f9","arxiv_id":"2412.11507","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Long-lived disk vortices excite density waves that shock, open gas gaps, and concentrate dust into rings, with elongated vortices producing rings at larger separations.","lead":"Simulations show that vortex-shaped gas flows in young planetary disks can create the same ring-and-gap patterns that astronomers usually attribute to hidden planets. If real vortices behave this way, many observed disk rings may need a new explanation, and the paper shows vortices can be far more efficient than small planets at making distant rings.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'two orders of magnitude' AMF claim rests on a periodic shearing box in which the vortex spans the y-domain and its two waves overlap (footnote 1); no box-length or isolation test is provided, so the comparison to planet waves is not yet secure.","rationale":"Good-faith reading: the paper is a serious simulation study. It establishes the causal chain (vortex launches waves, waves shock, AMF decays, vortensity jumps at the same radius as the gas gap, pressure bumps form dust rings) with internally consistent diagnostics and a resolution-convergence test in Appendix A. The synthetic image and the discussion of missing physics are honest. I am not objecting to the basic mechanism. The single load-bearing concern is the isolation of the vortex in the AMF measurement. The shearing-box setup deliberately truncates the azimuthal domain to 1/m* and keeps one vortex; because the vortex is as long as the domain, its two waves form a coherent, periodically repeated pattern rather than the two arms of a compact perturber. The FJ integral in Eq. 13 and the shock location both come from that pattern. The comparison planet is a point-like Lindblad perturber, whose wave flux is given by standard isolated-planet formulas, so this is not a matched comparison. If the periodic coherence boosts FJ or moves the shock, the two-orders result is not transferable. The paper itself flags the overlap in footnote 1 but does not quantify its effect. I considered whether the empirical x_ring-χ_v relation is instead the weakest point, because τ20 varies by a factor of roughly six between models and the gap widens with time; that is a genuine secondary issue and should be checked, but it is separate from the central AMF comparison. The AMF ratio is the most striking quantitative claim and is the one most directly threatened by the acknowledged geometry. The proposed test (an isolated-vortex run in a larger box, or an FJ integral restricted to the wave's physical support) would settle whether the effect is a box artifact. Until that is done, the paper should remain conditional rather than being accepted on its headline numbers. This matches the reader's verdict, so no change is needed.","tokens_in":16155,"tokens_out":18977,"duration_ms":187443,"concrete_test":"Run Model A4 again with the identical initial density bump, resolution, and m*=3 perturbation, but double the azimuthal box length (L_y = 2 x 20π/3 h) while keeping a single vortex; if the vortex does not remain isolated in the larger box, instead initialize an isolated elliptical vortex with the same aspect ratio χ_v=8.28 and size as in Model A4 at τ20. Measure FJ(x) via Eq. 13 and the shock location via the AMF-decay/vortensity jump. If FJ per vortex or x_shock changes by more than ~30% relative to the current full-domain measurement, the periodic overlap is biasing the AMF comparison and the 'two orders' claim must be re-evaluated; if they are unchanged, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's flashiest quantitative result is that vortex-induced density waves carry over two orders of magnitude more angular momentum flux than planet-induced waves that shock at the same location (abstract, Fig. 8). That comparison uses the shearing-box FJ integral in Eq. 13, measured in the truncated azimuthal domain (1/m*). Footnote 1 states that 'since the vortex is as long as the simulation domain, the two waves overlap.' With periodic y-boundaries, the box therefore contains not one isolated vortex but a periodic train of touching vortices, and the measured FJ and the shock location (from AMF decay, Fig. 2c) are properties of the overlapping superposition, not of an isolated vortex. The planet-wave flux, by contrast, is evaluated with standard isolated-planet formulas (Eqs. 16-17). If the y-coherence inflates FJ or shortens the shock length, the claimed two-orders-of-magnitude advantage could shrink or disappear in a global disk where a vortex has finite azimuthal extent. No sensitivity test to L_y or to the 1/m* truncation is reported, and L_y changes between models (Table 1), so the same geometry could also feed the x_ring-χ_v correlation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper investigates whether Rossby Wave Instability (RWI) vortices can open gaps and produce dust rings in protoplanetary disks through the shock dissipation of vortex-excited density waves. The authors run 2D inviscid shearing-box simulations with gas and a single dust fluid (St = 0.1) using Athena++, initializing four models with different Gaussian density-bump amplitudes to produce vortices of different elongations. In the lowest-amplitude model (A4), the vortex excites density waves whose angular momentum flux (AMF) stays roughly constant until x = 13.6h, where the AMF begins to decay, a vortensity jump appears, and a gas gap opens; dust drifts to pressure bumps and forms two rings. Over the four models, more elongated vortices produce rings at larger separations, fit by the empirical relation x_ring/h = 2.07 chi_v - 0.2 (Eq. 15). Using analytic scalings for planet-disk interactions (Goodman & Rafikov 2001; Dong et al. 2011a), the authors compare vortex-induced and planet-induced waves and claim the former carry over two orders of magnitude more AMF when their shock locations coincide. A synthetic ALMA image from a global simulation is presented as resembling HD 135344B.","tokens_in":16342,"tokens_out":4816,"duration_ms":46700,"significance":"If the central claims hold, the paper establishes a genuinely non-planetary pathway to the ring-and-gap substructures observed in protoplanetary disks, and it makes a strong, falsifiable statement that vortex-induced density waves are far more efficient angular-momentum carriers than planet-induced waves with the same shock radius. The paper has clear strengths: the causal chain from wave excitation to shock, AMF decay, gas-gap formation, and dust-ring placement is checked with multiple independent diagnostics (AMF profiles, vortensity jumps, and surface-density morphology); the convergence test in Appendix A explicitly motivates the 256 cells/h resolution; and the model parameters and fitting are reported in tabular form. However, the headline quantitative results rest on a small number of models and on an idealized shearing-box geometry, so the significance is conditional on additional robustness tests.","major_comments":[{"comment":"The most serious concern is that the shearing box is periodic in y and the vortex spans the full azimuthal domain, so the two density waves on opposite sides overlap (footnote 1, p. 4). Consequently, the AMF measured with Eq. (13) and the shock location inferred from the AMF decay in Fig. 2(c) are properties of a superposition of overlapping waves, not of an isolated vortex. The comparison with planet-induced waves in §3.4 uses analytic scalings derived for an isolated planet (Eqs. 16-17), so the claimed 'two orders of magnitude' advantage may be an artifact of the periodic geometry. No sensitivity test to the azimuthal domain size L_y or to the 1/m* truncation is reported, and L_y differs among the models in Table 1 (10π/4 vs. 10π/3), which could also feed the x_ring-chi_v correlation. The authors should add simulations with larger L_y (or otherwise isolated vortices) and demonstrate that FJ, the shock location, and the ring location do not depend on the box width.","section":"§2.3, footnote 1, and Fig. 2/3"},{"comment":"The comparison in Fig. 8 mixes directly simulated vortex data with analytic planet-wave predictions. It should be stated explicitly that the blue points are not from simulations but from Eq. (17) applied to the Mp,shock values from Eq. (16). As written, the reader cannot tell whether the factor-of-100 difference is a robust physical result or an artifact of concatenating two approximate scalings. Moreover, Eq. (17) is calibrated for weak, weakly nonlinear planet waves; the authors should justify that it remains valid for the planet masses inferred here, and ideally verify the comparison with a direct planet-disk simulation at the same resolution and box size.","section":"§3.4, Fig. 8, and Eqs. 16-17"},{"comment":"The linear correlation in Eq. (15) and Fig. 6 is based on four models, and the ring position x_ring is evaluated at a different time for each model (tau_20 ranges from 149 to 860 orbits). Section 3.1 states that 'the gas gap deepens and widens with time, pushing the two dust rings further apart,' so x_ring is time-dependent. The fit therefore conflates vortex elongation with evolution time. The authors should either evaluate all models at a common epoch, demonstrate that the ordering of x_ring is stable over the simulated interval, or include time as an explicit variable in the scaling.","section":"§3.3, Table 1, and §3.1"},{"comment":"The synthetic ALMA image that supports the observational claim 'detectable by ALMA' is produced from a global simulation whose maximum resolution near the vortex is approximately 33 cells per scale height (Appendix C.1), far below the 256 cells/h convergence requirement established in Appendix A for capturing weak-wave propagation and shock locations. The authors should either increase the global resolution to a level consistent with their own convergence criterion, or soften the claim to a qualitative morphological match rather than a resolved prediction of detectability.","section":"§4.1 and Appendix C"}],"minor_comments":[{"comment":"The text 'HLLE Reimann solver' and 'the default Reimann solver' should read 'HLLE Riemann solver' and 'Riemann solver'; also 'trancation errors' should be 'truncation errors'.","section":"§2.1"},{"comment":"The sentence 'where 4 are the gas continuity and momentum equation, respectively' should be 'where Equations (3) and (4) are the gas continuity and momentum equations, respectively.'","section":"§2.1"},{"comment":"The caption does not explain the dotted and solid arrows that are referenced in the text; please describe them explicitly in the caption.","section":"Figure 5 caption"},{"comment":"The fitted formula is rendered with garbled characters in the figure label ('𝑥!\"#$ℎ=2.07𝜒v−0.2'); please fix the typesetting.","section":"Figure 6"},{"comment":"The statement that the 'planet induced ones have a full-width-half-magnitude ~70% smaller' would be clearer if the comparison were made at matched wave amplitude or matched shock distance, since Fig. 3(a) and 3(b) use different normalizations and different shock locations.","section":"§3.2 and Figure 3"},{"comment":"The table header contains the typo 'T able 1'; also, the unit for x_ring is given as '[h]' but the text in Eq. (15) reports x_ring/h, so please make the notation consistent.","section":"Table 1"}],"recommendation":"major_revision","confidential_remarks":"The core mechanism is credible and the diagnostic chain is well validated, but the headline 'two orders of magnitude' claim and the x_ring-chi_v correlation need substantial hardening. In my view, the paper can become acceptable if the authors add box-size/isolation tests, address the time-dependence of ring locations, and clarify the analytic basis of the planet comparison. I do not see a reason for rejection if these load-bearing points are resolved."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, what's new: they actually demonstrate a vortex opening a gas gap and producing dust rings in a shearing box, with the causal chain checked – the vortensity jump, AMF decay, and ring location all line up at the same radius. That is a concrete result I hadn't seen before. The linear correlation between vortex aspect ratio and ring location is neat, though it's an empirical fit to four models.\n\nThe flashier claim – that vortex waves carry over two orders of magnitude more AMF than planet waves shocking at the same place – is built on indirect ground. The measured FJ comes from a box in which the vortex spans the full azimuthal extent, so the two waves overlap (their footnote 1). That means they're measuring a periodic train of touching vortices, not an isolated vortex. No test of box length or truncation is reported, so I'm not ready to trust the two-orders-of-magnitude number as a global statement. The planet comparison comes from analytic formulas (Goodman & Rafikov; Dong et al.), which are standard, but it's not a matched simulation. None of this breaks the central gap-opening mechanism, but it does mean the quantitative headline needs more support.\n\nOther soft spots are minor: the 4-point fit for Eq. 15, a single dust species at St=0.1, no dust feedback, inviscid 2D. They're honest about these in §4.2. The synthetic ALMA image for HD 135344B is a nice touch, though it's one example.\n\nThe paper deserves a serious referee. The right referee will ask for matched planet simulations, a domain-length sensitivity test, and ideally a global disk run or at least a justification that the overlapping-wave geometry doesn't inflate the AMF. If those hold up, this is a viable non-planetary route to rings. If not, the qualitative conclusion still stands but the 'two orders of magnitude' needs to be softened.\n\nI'd send it to review. It's a genuinely useful simulation result with honest limitations, and the authors seem to know what they don't include.","headline":"A genuinely new gap-opening mechanism that deserves refereeing, but the 'two orders of magnitude' AMF advantage over planets is not yet secure because the simulation vortex is a periodic train, not an isolated one.","tokens_in":16927,"tokens_out":2603,"would_cite":true,"duration_ms":24610,"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":"Long-lived vortices in protoplanetary disks can open the dust rings and gaps ALMA sees, and their density waves carry over two orders of magnitude more angular momentum than planet waves shocking at the same radius.","keywords":["protoplanetary disks","dust rings and gaps","Rossby wave instability","vortices","density waves","angular momentum flux","planet-disk interaction","hydrodynamic simulations"],"falsifier":"A global, full-azimuth simulation with a finite vortex (no $1/m_*$ truncation) that measures the angular momentum flux of the vortex wave at its shock radius would settle the central quantitative claim: if the flux is not more than two orders of magnitude larger than a planet wave shocking at the same radius, or if the outer dust ring fails to track $x_{\\mathrm{ring}}/h = 2.07\\,\\chi_v - 0.2$, the result fails.","tokens_in":15881,"feed_emoji":"🌀","tokens_out":9446,"duration_ms":78701,"temperature":0.7,"pith_summary":"The paper sets out to establish that vortices, the same anticyclonic eddies thought to produce crescent-shaped dust asymmetries, can also create the concentric dust rings and gaps that ALMA observes in protoplanetary disks, with no planet required. It reports two-dimensional hydrodynamic simulations, in a small co-rotating patch of an inviscid disk with dust, of a single long-lived vortex formed by the Rossby Wave Instability. The vortex excites spiral density waves; those waves steepen into shocks, deposit their angular momentum, and open a gas gap whose pressure maxima gather dust into rings. More elongated vortices put the outer dust ring farther away, following the fitted relation $x_{\\mathrm{ring}}/h = 2.07\\,\\chi_v - 0.2$ across the four simulated models. The paper further claims that vortex-driven density waves carry over two orders of magnitude more angular momentum flux than planet-driven waves that shock at the same location, making vortices potentially far more effective at opening distant gaps and rings.","feed_headline":"Disk vortices carve ALMA-visible rings and gaps without planets","feed_subtitle":"Simulations show vortex-driven shock waves can move over 100 times the angular momentum of planet waves.","key_machinery":"The engine of the result is the vortex itself, produced by the Rossby Wave Instability, an instability that amplifies a radial pressure bump into a persistent anticyclonic eddy. The vortex excites a pair of spiral density waves, and the paper tracks the wave's angular momentum flux $F_J(x) = \\int \\Sigma v_x\\,\\delta v_y\\,dy$, which remains constant until the wave shocks and then decays as angular momentum is deposited into the gas. Shock locations are identified independently by jumps in vortensity, $\\zeta = (\\nabla\\times\\mathbf{v})_z/\\Sigma$, which is conserved in inviscid barotropic flow except at shocks. The quantitative payload is carried by the empirical relation between vortex aspect ratio $\\chi_v$ and outer ring position, and by comparing $F_J$ for vortex waves against planet waves at equal shock distance.","core_discovery":"The central discovery is that a single long-lived vortex in an inviscid protoplanetary disk acts as a gap-opening agent through its own density waves, reproducing the ring-and-gap morphology usually attributed to planets. The vortex excites a pair of spiral density waves whose angular momentum flux stays roughly constant until the wave shocks; after the shock the flux decays and the lost angular momentum opens a gas gap, with dust rings forming at the gap edges. Quantitatively, the paper finds that the vortex aspect ratio $\\chi_v$ and the location of the outer dust ring obey $x_{\\mathrm{ring}}/h = 2.07\\,\\chi_v - 0.2$, and that vortex waves carry more than two orders of magnitude more angular momentum flux than planet waves that shock at the same distance, implying a substantially larger gap-opening capacity for the vortex in the far disk.","pith_inferences":["If the angular momentum flux excess survives in global simulations with finite vortices, vortex-driven wave transport could compete with planet-driven transport in setting the angular momentum budget of low-viscosity disks.","A testable extension is to compare disks that have crescents with disks that have rings but no detected planets; the vortex mechanism predicts the outer ring radius should correlate with the crescent's elongation, unlike planet-driven ring spacing.","Because the simulations use a single dust Stokes number and neglect dust feedback, the predicted ring contrast and location may depend on grain size, and multi-species dust runs would show whether the rings remain ALMA-detectable across realistic grain populations.","The shearing-box geometry makes the vortex span the whole simulated azimuthal strip, so the two emitted waves overlap; the reported two-orders-of-magnitude excess and the ring-location correlation could shift in a global disk where the vortex is finite."],"forward_implications":["A disk showing a crescent-shaped vortex should also tend to show a dust ring outside it, with the ring's separation controlled by how elongated the vortex is.","Rings and gaps in very young disks, where planets may not have had time to form, could be produced by vortices instead of requiring rapid planet formation.","Because vortex waves carry over two orders of magnitude more angular momentum than equal-shock planet waves, vortices may dominate gap opening in disks that host them, especially at large radii.","The fitted relation $x_{\\mathrm{ring}}/h = 2.07\\,\\chi_v - 0.2$ offers observers a quantitative check: measure the vortex aspect ratio from a crescent and predict where the outer ring should sit."],"supporting_citations":[{"why":"Supplies the wave-shocking theory and shock-length relation used to locate gap opening and to convert vortex shock radii into equivalent planet masses.","marker":"Goodman & Rafikov 2001"},{"why":"Provides the angular momentum flux relation for planet-induced waves used in the comparison with vortex waves.","marker":"Dong et al. 2011a"},{"why":"Gives the Gaussian-bump initial condition and linear Rossby Wave Instability mode selection used to form the simulated vortices.","marker":"Ono et al. 2018"},{"why":"Supplies the pressureless dust module and the wave-damping boundary treatment used in the simulations.","marker":"Huang & Bai 2022"},{"why":"Gives the angular momentum flux integral used to measure wave transport and identify shock locations.","marker":"Miranda & Rafikov 2020"},{"why":"Establishes the gap-opening mechanism by which angular momentum deposited by dissipating waves opens gaps in the disk.","marker":"Lin & Papaloizou 1986"},{"why":"Supplies the ALMA observations of HD 135344B used for the synthetic-image comparison.","marker":"Cazzoletti et al. 2018"},{"why":"Documents the co-existing crescent and inner ring in HD 135344B that motivates the vortex-ring connection.","marker":"van der Marel et al. 2016"}],"fun_headline_variants":["Vortices, not planets, carve disk gaps and rings","Disk vortices create ALMA-visible rings and gaps","Vortex waves outdo planets at carving disk rings","Single vortex can explain protoplanetary disk rings"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that a vortex stretched across the whole simulated azimuthal strip, whose two density waves overlap coherently, represents a real finite vortex in a disk, because the measured angular momentum flux, shock radius, and ring-position correlation all come from that overlapping geometry.","fun_headline_variants_meta":{"raw":{"variants":["Vortices, not planets, carve disk gaps and rings","Disk vortices create ALMA-visible rings and gaps","Vortex waves outdo planets at carving disk rings","Single vortex can explain protoplanetary disk rings"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000302,"raw_usage":{"total_tokens":1700,"prompt_tokens":868,"completion_tokens":832,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":484,"completion_tokens_details":{"reasoning_tokens":776}},"tokens_in":484,"tokens_out":832,"duration_ms":7865,"temperature":1.0,"reasoning_tokens":776,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T14:51:40.231769+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A global, full-azimuth simulation with a finite vortex (no $1/m_*$ truncation) that measures the angular momentum flux of the vortex wave at its shock radius would settle the central quantitative claim: if the flux is not more than two orders of magnitude larger than a planet wave shocking at the same radius, or if the outer dust ring fails to track $x_{\\mathrm{ring}}/h = 2.07\\,\\chi_v - 0.2$, the result fails.","supporting_citations":[{"cited_title":"2018, The Astrophysical Journal, 864, 70, doi: 10.3847/1538-4357/aad54d","cited_arxiv_id":null,"evidence_quote":"Gives the Gaussian-bump initial condition and linear Rossby Wave Instability mode selection used to form the simulated vortices."}],"review_version":1}