{"id":"908a0409-d4d9-458f-8327-fdd7ea57e93e","arxiv_id":"2504.20036","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Rotation-curve measurements of 15 protoplanetary disks show that the majority of dust rings and gaps sit at gas pressure maxima and minima, supporting pressure-driven dust trapping.","lead":"Using ALMA observations of 15 planet-forming disks, astronomers measured the rotation speed of the gas and found small wiggles that line up with the dusty rings and gaps in each disk. The pattern suggests that gas pressure bumps and dips, possibly carved by young planets, are the main organizing force behind most dust substructures.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 75%/80% co-location statistic lacks a well-defined denominator and a random-co-location baseline; because δυϕ substructures are ubiquitous, the headline percentage is not yet interpretable as evidence for pressure-driven ring/gap formation.","rationale":"The reader's verdict is CONDITIONAL with moderate confidence, and I agree that the paper is careful and that the central claim is plausible but not yet quantitatively secured. However, I do not think the single most load-bearing weakness is the velocity-decomposition assumption emphasized in the reader's weakest_assumption field. Appendix D directly tests that assumption for a planet-driven spiral and shows that although amplitudes can change by ~25%, the sign pattern of δυϕ is preserved in that simulation; the paper also flags beam-smearing and warp/radial-flow caveats in Sec. 6.3. The statistical interpretation of the co-location fraction is more directly load-bearing because it concerns the headline number itself. The abstract presents 75%/80% as the key evidence for pressure variations being the dominant mechanism, but the fraction is computed without a stated null expectation and with selection criteria that may encode the expected sign. The paper does include Table 2 listing non-aligning and unresolved features, which is good reporting, and the VADP release would make a re-analysis feasible. The absence of a random-co-location baseline and an exact denominator means the quantitative strength of the conclusion is currently overstated; a conditional verdict with a request for a baseline and sensitivity analysis is appropriate. I do not see an internal inconsistency that would justify rejection, and the projection-effect appendix provides some independent support for the sign-based method.","tokens_in":36557,"tokens_out":10109,"duration_ms":105559,"concrete_test":"Using the public VADP radial profiles and the Curone et al. continuum catalog, recompute the co-location fraction under two variants. First, build a null model by randomly shifting each ring/gap radius (preserving the number of features, their widths, and the observed δυϕ profiles) and recomputing the sign-of-gradient match 10^4 times; report the shuffled baseline and its 95% interval. Second, perform a sensitivity run that drops selection criterion (4) and, when possible, criterion (5), counting all strong axisymmetric rings/gaps with at least one independent beam of extracted δυϕ data, with flat or reversed gradients scored as non-aligned. If the shuffled baseline approaches ~50% or the relaxed co-location fraction drops well below the claimed 75%/80%, the headline percentage must be re-stated with an explicit denominator and a significance level.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim rests on the co-location percentage: \"More than 75% of the rings and 80% of the gaps in the dust continuum emission resolved in δυϕ are co-located with gas pressure maxima and minima.\" The denominator of this statistic is not fully specified, and the selection procedure in Sec. 3.4 may embed the expected sign into the sample. Criterion (4) requires that a continuum substructure \"must exhibit a radially increasing (at dust gaps) or decreasing (at dust rings) δυϕ-profile\" to be reported in Table 1. Section 5.2 then says that non-aligning features are listed in Table 2 and that features that could not be assessed are \"not included in the total count of non-aligning δυϕ-substructures considered.\" If criterion (4) is effectively used to define the measured sample, the 75%/80% match is partly definitional rather than a test of the co-location hypothesis. If it is not, the paper should state explicitly which continuum features enter the denominator and which are excluded as unresolved, non-axisymmetric, or low-contrast. In addition, no random-co-location baseline is provided. The paper itself reports that δυϕ substructures are ubiquitous (Sec. 5.1), so even if pressure bumps are unrelated to dust rings, a substantial fraction of rings will sit on a local δυϕ extremum by chance. The conclusion that pressure variations are \"likely the dominant mechanism\" requires a null expectation for the co-location fraction, not just a high raw percentage. The Appendix D projection study and the public VADP release are real supporting assets, but they address velocity decomposition errors, not the statistical interpretation of the headline fraction.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents azimuthally averaged rotation curves and Keplerian-subtracted velocity deviations (δυϕ) for 15 protoplanetary disks from the exoALMA Large Program, measured in 12CO and 13CO J=3–2 emission. It reports vertical thermal stratification inferred from the differing rotation of the two CO isotopologues, ubiquitous small- and large-scale δυϕ substructures, and a co-location statistic between dust continuum rings/gaps and gas pressure maxima/minima identified from the sign of the radial derivative of δυϕ. The paper also derives midplane logarithmic pressure gradients for a subset of disks, discusses systematic uncertainties including beam smearing, projection effects, and stellar-mass offsets, provides rough planet-mass estimates for some gaps, and releases the extracted radial profiles as a Value-Added Data Product. The central claim is that more than 75% of dust rings and 80% of dust gaps are co-located with pressure maxima and minima, implying that gas pressure variations are the likely dominant mechanism for continuum substructure formation.","tokens_in":36896,"tokens_out":5705,"duration_ms":58236,"significance":"If the co-location claim holds, the paper would supply the strongest observational support to date for pressure-bump dust trapping as the leading explanation of ring and gap formation in protoplanetary disks, extending earlier smaller-sample results (e.g., Izquierdo et al. 2023) to a homogeneous ALMA Large Program sample. The paper's strengths include a careful and transparent velocity extraction with documented systematics, a sign-based diagnostic that is independent of the absolute stellar-mass offset, and a public release of the kinematic data products. The midplane pressure derivative derivations, though limited to a subset of sources, are a promising new observable. However, the headline co-location statistic is not yet on a sound statistical footing: its denominator is not fully defined, and no random-co-location baseline is established, so the significance of the percentages relative to a null hypothesis is currently unclear.","major_comments":[{"comment":"The co-location fraction reported in the abstract and Section 5.2 is computed on a sample that is partly defined by the outcome. Criterion (4) in Section 3.4 restricts entries in Table 1 to continuum substructures with a δυϕ gradient of the theoretically expected sign, and Section 5.2 explicitly states that substructures that could not be assessed are not included in the total count of non-aligning features. The reader cannot reconstruct the denominator of '16 out of 21 rings and 10 out of 12 gaps' from the text and tables alone. Please specify a fixed, pre-registered sample—all axisymmetric continuum rings and gaps with contrast ID/IB < 0.8 that are beam-resolved, including those with flat, reversed, or unassessable δυϕ gradients—and tabulate the co-location fraction over that full sample for each tracer separately.","section":"Section 3.4 and Section 5.2"},{"comment":"The headline statistic lacks a null expectation. The paper itself reports that δυϕ substructures are ubiquitous (Section 5.1), so a substantial fraction of continuum features would coincide with a local sign change of ∂δυϕ/∂R even if pressure bumps were unrelated to dust rings and gaps. To support the conclusion that gas pressure variations are 'likely the dominant mechanism' for ring and gap formation, provide a random-co-location baseline or permutation test, for example comparing the observed sign at continuum features with the distribution of signs at randomly drawn radii in the same disks.","section":"Section 5.1 and Section 5.2"},{"comment":"The sign-based diagnostic is robust to the stellar-mass offset because a constant offset shifts δυϕ vertically without changing the sign of its radial derivative, but the paper's assertion that 'the shape of the δυϕ-profile is hardly affected' by projection effects is only demonstrated for one planet-driven spiral simulation. Section 6.3 itself identifies beam smearing at steep intensity gradients, warps, and low-SNR outer-disk biases as sources of radial-dependent velocity errors that could, in principle, alter the sign of ∂δυϕ/∂R at some features. Because the co-location claim rests entirely on signs, please add tests that bound the fraction of sign flips under these systematics, for example using existing hydrodynamic simulations with radial or vertical flows, or by masking the innermost and outermost beams where biases are strongest, or temper the 'dominant mechanism' conclusion to the level of a consistency check.","section":"Appendix D and Section 6.3"}],"minor_comments":[{"comment":"The phrase 'resolved in δυϕ' should be defined explicitly; the current wording leaves open whether a continuum substructure is counted only when the δυϕ sign is measurable, which is precisely the selection effect at issue in the major comments.","section":"Abstract and Section 5.2"},{"comment":"The note that unassessable substructures are 'not included in the total count of non-aligning δυϕ-substructures considered' should state how the total count of considered features is built and where that list is published, so that the denominator of the co-location fraction is transparent.","section":"Table 2 note"},{"comment":"The sentence 'Out of 17 continuum gaps with a positive CO δυϕ radial gradient' is hard to reconcile with the earlier counts (10 gaps for 12CO and 8 for 13CO); please clarify whether this is the sum over both tracers and over which gap sample, and align the numbering with Table 1 and Table 2.","section":"Section 6.2"},{"comment":"The caption describes both the High Resolution Images and the High Surface Brightness Sensitivity Images as 'red'; one of the colors in the caption is presumably a typo.","section":"Figure E.1 caption"}],"recommendation":"major_revision","confidential_remarks":"The central statistical concern about the co-location statistic is fixable without new observations, through a rigorous definition of the denominator and the addition of a null baseline. If the authors address this, I would expect the paper to be publishable. I also note that the acknowledgments thank an anonymous referee, which is standard and not an issue for this report."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear colleague,\n\nQuick take: this is a solid, carefully written exoALMA paper that delivers genuinely new measurements — full rotation curves and δυϕ profiles for 15 disks in 12CO and 13CO, plus the first observational estimate of the midplane pressure derivative. The vertical stratification result (13CO rotates faster than 12CO beyond the isothermal expectation) is convincing and well documented. The public VADP release is a real asset for the community.\n\nThe central claim — that >75% of rings and >80% of gaps are co-located with pressure maxima/minima — is plausible but the quantitative framing is weaker than it should be. The stress-test concern about the denominator is half right. The paper does include non-aligning features: Table 2 explicitly lists reversed and flat δυϕ gradients, so the 16/21 and 10/12 counts are not purely definitional. But the exact denominator is not clearly enumerated. The text says features that could not be assessed are excluded, and the reader has to reconstruct the full sample from Curone et al. (2024) plus the selection criteria. That should be stated explicitly. More importantly, there is no random-co-location baseline. Given that δυϕ substructures are ubiquitous (Sec. 5.1), a raw 75–80% alignment is not enough to establish that pressure variations dominate ring/gap formation. The paper would be stronger if it showed what fraction would coincide by chance, even approximately.\n\nOther soft spots are minor: the midplane pressure derivative error budget omits covariances from the multi-level fits (acknowledged), and the planet mass estimates in Appendix B are rough. The projection-effect study in Appendix D is a nice piece of due diligence.\n\nWho is this for? Planet formation people, especially those working on disk kinematics and dust trapping. It is a valuable catalog paper that extends earlier MAPS work to the exoALMA sample. I would send it to a serious referee — the data and analysis deserve scrutiny, and the co-location statistic should be tightened before publication.\n\nRecommendation: engage with it. The authors should be asked to make the denominator explicit and add a null expectation, but this is a good paper.","headline":"A careful, transparent survey of rotation-curve perturbations in 15 disks; the headline co-location statistic is plausible but needs a clearer denominator and a null baseline before it can carry the 'dominant mechanism' conclusion.","tokens_in":37693,"tokens_out":3405,"would_cite":true,"duration_ms":32141,"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":"Rotation curves of 15 planet-forming disks show that over 75% of dust rings and 80% of dust gaps coincide with gas pressure maxima and minima, making pressure variations the likely dominant cause of the observed continuum substructures.","keywords":["protoplanetary disks","planet formation","rotation curves","gas pressure gradients","dust ring and gap formation","molecular line kinematics","ALMA observations","disk vertical stratification"],"falsifier":"A decisive test is to fit the full three-dimensional velocity field—including radial and vertical motions—from the same CO cubes in a few disks and recompute the pressure extrema from the sign of the corrected radial velocity gradient; any dust ring or gap where the sign reverses under this correction would falsify the co-location claim. A more targeted version is to examine disks with known warps or cavity flows, such as J1604 and HD 143006, and check whether their apparent pressure bumps survive when the axisymmetric-rotation assumption is relaxed.","tokens_in":36348,"feed_emoji":"🪐","tokens_out":8636,"duration_ms":82580,"temperature":0.7,"pith_summary":"This paper uses ALMA observations of 15 planet-forming disks to test whether the rings and gaps seen in the dust are carved by variations in gas pressure. It measures how the gas rotation speed departs from the velocity expected from the star's gravity alone, and shows that in over 75% of dust rings the gas rotates more slowly just inside the ring and faster just outside it—the signature of a pressure maximum—while over 80% of dust gaps show the reverse pattern, the signature of a pressure minimum. The same rotation curves reveal that the hotter upper gas layer ($^{12}$CO) rotates more slowly than the cooler lower layer ($^{13}$CO), evidence that the disks are vertically stratified. From these curves the paper derives, for the first time from observations, the radial pressure gradient at the disk midplane, and finds it consistent with the dust substructure positions. If the interpretation is right, gas pressure bumps become the leading explanation for the rings and gaps seen in planet-forming disks.","feed_headline":"Over 75% of dust rings trace gas pressure maxima","feed_subtitle":"Rotational wobbles in 15 exoALMA disks tie most rings and gaps to gas pressure peaks and troughs.","key_machinery":"The load-bearing object is the azimuthally averaged deviation from Keplerian rotation, $\\delta\\upsilon_{\\phi} = \\upsilon_{\\phi} - \\upsilon_k$, extracted from molecular-line centroid velocity maps. In centrifugal balance the deviation is tied to the pressure gradient by $\\partial\\ln P/\\partial\\ln R \\approx 2\\upsilon_k^2/c_s^2\\, \\delta\\upsilon_{\\phi}/\\upsilon_k$, so the paper reads off pressure maxima and minima from the sign of $\\partial\\delta\\upsilon_{\\phi}/\\partial R$ at each dust feature rather than from absolute velocities. A toy model with an imposed Gaussian gap in the pressure calibrates this sign mapping. Supporting pieces are a vertical-structure integral that converts the pressure gradient at the CO emitting height into a midplane pressure gradient, and a projection study of a planet-driven spiral wave showing that radial-flow contamination can change the amplitude of $\\delta\\upsilon_{\\phi}$ by up to 25% while preserving the sign pattern of its radial profile.","core_discovery":"The central claim is that most dust continuum rings and gaps in the sample are co-located with gas pressure maxima and minima, so that gas pressure variations are likely the dominant mechanism forming these substructures. The evidence is the sign of the radial derivative of the deviation from Keplerian rotation, $\\delta\\upsilon_{\\phi}$: a pressure maximum makes the rotation speed decrease with radius, and a pressure minimum makes it increase. In $^{12}$CO, 16 of 21 continuum rings and 10 of 12 gaps show the expected sign, and in $^{13}$CO the counts are 14 of 17 rings and 8 of 10 gaps. Because the diagnostic uses only the sign of the velocity gradient, the inferred pressure extrema do not depend on the exact stellar mass. The paper also reports vertical thermal stratification across most of the sample, sub-Keplerian rotation in the outer disks, and midplane pressure derivative profiles for a subset of sources.","pith_inferences":["The sign-based diagnostic is robust to the 25% amplitude projection error the paper quantifies, but not to sign flips from unmodeled radial flows; if some of the substructures classed as 'not accessible' were counted as failures, the 75% and 80% rates would be lower.","The pressure substructures seen beyond the dust continuum predict that deep continuum imaging at those radii should find little or no trapped dust; a targeted search is a straightforward test of whether those bumps are leaky, short-lived, or dust-free.","Applying the same rotation-curve technique to multiple isotopologues at different heights can separate density-driven from temperature-driven pressure variations, potentially turning the co-location statistic into a diagnostic of the physical origin of each bump.","If the co-location rates hold across a larger sample, continuum ring catalogs alone would become statistical tracers of gas pressure maxima, letting population studies probe the frequency of pressure bumps without needing kinematic data."],"forward_implications":["Most dust rings and gaps are not independent dust phenomena: they mark gas pressure maxima and minima, so models of ring and gap formation should reproduce pressure bumps at those radii.","The midplane pressure derivative is now measurable from rotation curves and temperature maps, allowing direct confrontation with predicted dust drift speeds and dust-trapping efficiencies.","Vertical stratification is a general property of these disks, so single-molecule rotation curves cannot be translated into stellar masses without correcting for pressure support at the emitting height.","Sub-Keplerian rotation in the outer disk speeds up the inward drift of pebbles, potentially replenishing the inner disk with dust from large radii.","Gas pressure substructures beyond the dust continuum imply that outer pressure bumps either trap dust inefficiently or are short-lived, sharpening the question of what sets the outer dust radius."],"supporting_citations":[{"why":"Establishes the azimuthal-averaging velocity extraction and previously found nine of eleven continuum rings co-located with pressure maxima; the method this paper extends to the full sample.","marker":"(Izquierdo et al. 2023)"},{"why":"Introduces the sign-of-the-radial-derivative diagnostic that identifies pressure minima and maxima without requiring an exact stellar mass.","marker":"(Rosotti et al. 2020)"},{"why":"Provides the dust continuum ring, gap, and crescent catalog whose locations are the reference for the co-location statistics.","marker":"(Curone et al. 2024)"},{"why":"Supplies the 2D temperature structure and emission-height profiles used to convert surface-layer pressure gradients to the midplane.","marker":"(Galloway-Sprietsma et al. 2024)"},{"why":"Provides self-gravity-corrected rotation curve fits and stellar mass estimates used in the midplane pressure derivative calculation.","marker":"(Longarini et al. 2024)"},{"why":"Describes the imaging products and spatial/spectral resolution sets that make the small-scale velocity measurements possible.","marker":"(Teague et al. 2024)"},{"why":"Establishes that dust radial drift is set by the midplane pressure gradient, the theoretical link between pressure bumps and dust substructures.","marker":"(Whipple 1972; Weidenschilling 1977)"},{"why":"Introduces the line-centroid and surface-fitting machinery used to produce the velocity maps.","marker":"(Izquierdo et al. 2021)"},{"why":"Provides the simulation-based relation that converts measured gap widths in the velocity perturbation into planet mass estimates used in the appendix.","marker":"(Yun et al. 2019)"}],"fun_headline_variants":["Pressure peaks align with most disk dust rings","Gas pressure variations shape planetary disk rings","Rotation deviations expose pressure-driven disk rings","Dust rings and gaps trace gas pressure extrema","Most protoplanetary disk rings trace pressure maxima"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument assumes the measured gas motion is almost entirely circular rotation, so that any remaining small-scale velocity variation can be attributed to a pressure gradient; if unresolved radial flows or warps flip the sign of the velocity gradient at a dust feature, the pressure extremum inferred there is an artifact.","fun_headline_variants_meta":{"raw":{"variants":["Pressure peaks align with most disk dust rings","Gas pressure variations shape planetary disk rings","Rotation deviations expose pressure-driven disk rings","Dust rings and gaps trace gas pressure extrema","Most protoplanetary disk rings trace pressure maxima"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000324,"raw_usage":{"total_tokens":1867,"prompt_tokens":1044,"completion_tokens":823,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":756}},"tokens_in":660,"tokens_out":823,"duration_ms":8076,"temperature":1.0,"reasoning_tokens":756,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T05:37:29.519148+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A decisive test is to fit the full three-dimensional velocity field—including radial and vertical motions—from the same CO cubes in a few disks and recompute the pressure extrema from the sign of the corrected radial velocity gradient; any dust ring or gap where the sign reverses under this correction would falsify the co-location claim. A more targeted version is to examine disks with known warps or cavity flows, such as J1604 and HD 143006, and check whether their apparent pressure bumps survive when the axisymmetric-rotation assumption is relaxed.","supporting_citations":[{"cited_title":"2019, , 884, 142, 10.3847/1538-4357/ab3fab","cited_arxiv_id":null,"evidence_quote":"Provides the simulation-based relation that converts measured gap widths in the velocity perturbation into planet mass estimates used in the appendix."}],"review_version":1}