{"id":"1a2f1a65-7af8-439b-a220-6f851e83275c","arxiv_id":"2411.14241","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Giant planets in the low-density CO gas of debris disks can open gaps that fall below ALMA's detection threshold, making planetary gaps observable in the gas and offering a new indirect exoplanet detection channel.","lead":"Simulations show that giant planets can carve gaps in the thin carbon-monoxide gas around old debris disks, gaps deep enough to show up in ALMA images, something that rarely works in younger gas-rich disks. These gaps could become a new way to find cold giant planets that direct imaging misses.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 0.5 M_J threshold rests entirely on the assumed alpha=1e-3; at the alpha~0.1 allowed for beta Pic, the same scaling requires roughly 5-10 M_J, so the headline and HD138813 prediction are not robust.","rationale":"The reader's CONDITIONAL verdict is supported by this stress test. The central idea, that CO gas in low-column-density debris disks can be empty enough that the gap falls below ALMA's 3-sigma floor while the surrounding disk is detected, is physically plausible and is demonstrated by the pipeline for the chosen parameters. The independent support includes the use of standard public codes (FARGO3D, RADMC-3D, CASA) and a clean, falsifiable HD138813 prediction. The vulnerability is not an internal inconsistency; it is that all quantitative thresholds are evaluated at one point in (alpha, h) space, with alpha fixed at a medium value while the one existing constraint (beta Pic, Sect. 5.2.1) allows alpha up to 0.1. Since gap depth depends as alpha^-1 and q^2, the allowed alpha range shifts the minimum planet mass by roughly an order of magnitude. The paper's Eq. 9 validation (Fig. 8) is partly circular for this question because the images and the criterion use the same hydrodynamical surface densities; it cannot constrain the alpha dependence. A re-run at alpha=0.1 is cheap relative to the claim and would settle whether the 0.5 M_J threshold is physical or an artifact of the assumed viscosity. Until then, the correct verdict remains CONDITIONAL, with the quantitative thresholds stated as valid only for low-alpha disks.","tokens_in":25011,"tokens_out":7124,"duration_ms":69549,"concrete_test":"Re-run the fiducial low-mass-disk case (mp=0.5 M_J, rp=50 AU, mdisk=10^-3 M_Earth, 40 pc) through the FARGO3D->RADMC-3D->CASA pipeline with alpha=0.1 instead of 10^-3, keeping h and all other settings fixed. Compute Fgap/(3Fsens) from the new Sigma_gap and produce moment-zero images. If the ratio exceeds 1 or the gap is not visually detected, the 0.5 M_J threshold must be replaced by a mass-vs-alpha curve; if it remains observable, the alpha concern is retired for this case. An analytic version of the same check is to insert alpha=0.1, h=0.027, q=0.0005 into Eq. 10 and then Eq. 9, but the simulation is the decisive test because Fig. 11 shows analytic criteria can be off by more than 10x in this regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The most load-bearing assumption is Sect. 2.2's fixed alpha=10^-3 and h=0.01(r/1AU)^0.25, because the quantitative claims ('planets as small as 0.5 M_J', the HD138813 Jupiter-mass gap) are computed from FARGO3D surface densities at that single point. Sect. 5.2.1 explicitly admits that beta Pic's turbulence can be as high as alpha~0.1 and that higher alpha requires larger planetary masses, but no simulation or recalculated threshold is given for alpha>10^-3. The observability criterion (Eq. 9) uses Sigma_gap from those simulations; increasing alpha by a factor of 100 changes the gap-depth parameter K=q^2 h^-5 alpha^-1 (Eq. 10) by the same factor. For the 0.5 M_J planet at 50 AU, this takes a very deep gap (K~10^4) to a shallow partial gap (K~10^2, Sigma_gap/Sigma_0~0.1-0.3). Recovering the original depth requires q^2 to rise by ~100, i.e. mp~5 M_J. Thus the 0.5 M_J headline and the 1 M_J HD138813 prediction are not robust to the observationally allowed viscosity. The paper is internally consistent, but the central quantitative claim is a one-parameter extrapolation; Sect. 5.3's Fig. 11 also shows that published gap-depth formulas disagree with the simulations by over an order of magnitude in low-mass cases, so using Eq. 10 as a cross-check is uncertain and strengthens the need for an alpha sweep.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript combines FARGO3D hydrodynamic simulations, RADMC-3D radiative transfer, and CASA synthetic ALMA observations to ask when a giant planet embedded in a low-mass CO debris disk produces an observable gas gap. It explores planet masses 0.5, 1, and 5 M_J at 10, 50, and 100 AU; disk masses 10^-5, 10^-3, and 10^-1 M_Earth; and distances 40, 100, and 130 pc. A one-beam LTE flux estimate (Eq. 9) is compared with the 3-sigma ALMA sensitivity to define gap observability, and the criterion is checked against the synthetic images (Fig. 8). The same machinery is applied to HD 138813, predicting that a Jupiter-mass planet near 70 AU would open a resolvable gap in roughly 7 hours of ALMA C-7 observations. The central claim is that, unlike protoplanetary disks, debris-disk gas is so tenuous that residual gas in planetary gaps falls below ALMA's detection threshold, making planets as small as 0.5 M_J indirectly detectable.","tokens_in":25240,"tokens_out":5793,"duration_ms":56630,"significance":"If the quantitative claims are robust, the paper opens a practical new channel for detecting cold giant planets around old stars, and the gap/kink dichotomy is a clean observable distinction. The modeling pipeline is standard and the synthetic-observation setup is realistic; the HD 138813 prediction is a genuine, falsifiable target that can be checked with new ALMA data. The transparent analytical criterion in Eq. 9 is a useful tool even if its absolute calibration needs refinement. However, the headline mass threshold and the target prediction are tied to two fixed disk parameters, alpha=10^-3 and h=0.01(r/1AU)^0.25, and the validation of the criterion in Fig. 8 uses the same hydrodynamical surface densities that produce the synthetic images, so the evidence is not yet as strong as the abstract's phrasing suggests.","major_comments":[{"comment":"The fixed viscosity alpha=10^-3 in Sect. 2.2 is load-bearing: the abstract's '0.5 M_J' claim and the Sect. 5.1 HD 138813 Jupiter-mass prediction are computed at this single value. The paper itself notes that beta Pic allows alpha up to 0.1 (Sect. 5.2.1). Because the gap-depth parameter scales as K=q^2 h^-5 alpha^-1 (Eq. 10), an increase of alpha by two orders of magnitude moves a 0.5 M_J planet from the deep-gap regime to a shallow partial gap; recovering the same depth requires roughly an order of magnitude more planet mass. The manuscript acknowledges this in words but provides no simulation or re-derived threshold for higher alpha, leaving the headline claim as a one-parameter extrapolation. Please add an alpha sweep (at least alpha=10^-2 and 10^-1 at the fiducial planet locations) or, if this is too costly, an explicit analytic bound with the caveat prominently reflected in the abstract and conclusions.","section":"Sect. 5.2.1, Eq. (10)"},{"comment":"The validation of the observability criterion is partly by construction. The F_gap/(3F_sens) map in Fig. 8 is computed (Eq. 9) from the same FARGO3D Sigma_gap fields that generate the synthetic images used for the red/green dots, so the agreement measures internal consistency rather than predictive power. The application to HD 138813 is independent and therefore valuable, but the claim in Sect. 4 that the criterion is 'an efficient way to determine whether gaseous gaps can be observed' should be supported by an out-of-sample test, for example a leave-one-configuration-out check or a comparison with an independently calibrated gap-depth relation. In light of Fig. 11, the latter option requires acknowledging that published relations can differ by more than an order of magnitude.","section":"Fig. 8, Eq. (9)"},{"comment":"Fig. 11 shows that the Kanagawa, Fung, and Pichierri gap-depth formulas disagree with the simulations by more than an order of magnitude in low-mass cases. This has two consequences for the paper's logic. First, Eq. 10 cannot be used as a reliable quantitative cross-check for the alpha sensitivity discussed in my first comment, although the direction of the effect is not in doubt. Second, because the simulations are the only source for Sigma_gap, the authors should state more explicitly that the numerical results, not the analytic scaling, set the quantitative thresholds; this strengthens the need for a documented resolution/convergence check and for additional hydrodynamical parameters (viscosity, aspect ratio) before the mass threshold is quoted as a general result.","section":"Sect. 5.3, Fig. 11"}],"minor_comments":[{"comment":"The article contains many typographical errors from spacing ('di fferent', 'wether', 'critera', 'us thus') and a copyediting pass is needed before publication.","section":"Throughout"},{"comment":"The text defines the gap width by w_gap=4r_H and then cautions that Eq. 6 should use the width at the bottom of the gap; since observers will measure an intensity profile, please state explicitly how the measured full width and the inclination projection enter Eq. 6.","section":"Eq. (6), Appendix B"},{"comment":"The dot colors in Fig. 8 are described in the caption and text, but the orange and blue cases are only explained later; adding a legend to the figure would improve readability.","section":"Fig. 8"}],"recommendation":"major_revision","confidential_remarks":"None beyond the report. I would advise the editor that the paper is within A&A's scope and the target application is a good selling point, but the mass threshold and the HD 138813 prediction should not be advertised before an alpha-sensitivity run is available."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Core idea is original and worth engaging: in the low-column-density CO gas of debris disks, a giant planet's gap can be so empty that the residual gas falls below ALMA's detection threshold, making gas gaps observable in a regime where protoplanetary disks hide them. That inverts the usual logic and leads to a clean, falsifiable prediction for HD138813. The paper does this with a standard, reproducible-in-principle pipeline (FARGO3D, RADMC-3D, CASA), a broad parameter study, and an honest discussion of what is uncertain. The kink-versus-gap behavior as a function of disk mass is a nice result, and their Eq. 9 criterion is physically sensible.\n\nThe soft spots are real but not fatal. First, the validation of the criterion in Fig. 8 is partly circular: the images are synthesized from the same hydrodynamical surface densities that feed Eq. 9, so the agreement is a consistency check of the imaging pipeline, not an independent test of the criterion. That doesn't undermine the observational channel; it just means the criterion's predictive power outside this simulation suite is not yet demonstrated. Second, the quantitative thresholds (0.5 MJ, the HD138813 Jupiter-mass gap) rest on a single viscosity choice, alpha=1e-3, with a thin disk. The paper explicitly notes beta Pic's turbulence can be as high as alpha~0.1 but does not propagate that into the claims. Eq. 10 (and their own Fig. 11, which shows analytic gap-depth scalings can be off by an order of magnitude) implies that a 100x higher alpha would require roughly 10x more massive planet to open the same gap, so the 0.5 MJ headline is a one-parameter extrapolation. An alpha sweep would have settled this. Third, no code or data are shipped; the CASA setup is 'available upon request,' which is weak but minor.\n\nI read the stress-test note as landing on the paper; the central claim survives, but the headline number does not. The paper's own admission in Sect. 5.2.1 confirms the concern, so I don't think the reader's conditional verdict is too harsh.\n\nBottom line: this deserves a serious referee. The idea is new, the prediction is testable, and a good referee can push on the viscosity sensitivity without sending the authors back to step zero. I'd bring it to reading group and would cite it if I worked in this area.","headline":"A genuinely new observable channel for planets in debris-disk gas, with a testable HD138813 prediction, but the headline planet masses rest on a single viscosity assumption.","tokens_in":25926,"tokens_out":3007,"would_cite":true,"duration_ms":26926,"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":"In the thin CO gas of old debris disks, planets as small as half a Jupiter-mass should open gaps that ALMA can see.","keywords":["debris disks","planet-disk interactions","CO gas","ALMA observations","planetary gaps","kinks","exoplanet indirect detection","hydrodynamical simulations"],"falsifier":"Observe HD138813's CO $J=2$-$1$ emission with ALMA in the C-7 configuration (beam $\\sim$0.09 arcseconds) for about 7 hours. If the moment-zero map shows no surface-brightness depression centered near 70 AU with width $\\sim$0.14 arcseconds, the paper's headline prediction for a Jupiter-mass planet fails. A second check: measure the turbulent $\\alpha$ of any debris-disk gas, as attempted for $\\beta$ Pictoris; a value near 0.1 would imply that the half-Jupiter threshold is too optimistic by orders of magnitude at the assumed scale height.","tokens_in":24666,"feed_emoji":"🪐","tokens_out":10775,"duration_ms":86871,"temperature":0.7,"pith_summary":"Debris disks around old stars carry far less gas than the planet-forming disks of their youth, and this paper argues that the difference makes planet-carved gaps visible for the first time in gas emission. Combining hydrodynamical simulations, radiative transfer, and simulated ALMA observations, the authors show that in disks with a few $10^{-3}\\,M_\\oplus$ of CO, a giant planet empties its gap so thoroughly that the leftover gas radiates below ALMA's $3\\sigma$ sensitivity while the surrounding disk stays bright. Planets as small as 0.5 Jupiter masses can do this in the lightest disks. The paper also derives a simple flux criterion (Eq. 9) that predicts when a gap is observable, and applies it to HD138813, where a Jupiter-mass planet at about 70 AU should produce a resolvable gap in a roughly 7-hour ALMA observation. If correct, the result turns debris-disk gas into a new indirect exoplanet detector, reaching cold giant planets that direct imaging cannot see.","feed_headline":"Half-Jupiter planets may open visible gaps in old debris disks","feed_subtitle":"Sparse CO gas makes gap interiors fall below ALMA's detection floor, offering a new way to find cold giant planets.","key_machinery":"The load-bearing object is an observability criterion built from a geometric resolution condition and a flux condition. The gap width is tied to the planet's Hill sphere, $w_{\\rm gap}=4r_{\\rm H}=4r_p(q/3)^{1/3}$, which the simulations find matches the width at the bottom of the gap. The flux inside the gap is $F_{\\rm gap}=h\\nu_{u,l}A_{u,l}x_u S_{\\rm beam}\\Sigma_{\\rm gap}/(4\\pi d^2 m_{\\rm mol})$ (Eq. 9), with $\\Sigma_{\\rm gap}$ taken directly from the hydrodynamical runs rather than from an analytic gap-depth formula. A gap is observable when it is resolved and $F_{\\rm gap}<3F_{\\rm sens}$. The mechanism works because the absolute surface density of debris-disk gas is low: the same fractional gap depth that hides gaps in protoplanetary disks pushes debris-disk gaps below the sensitivity floor. The half-Jupiter threshold and the HD138813 prediction both rest on this criterion.","core_discovery":"The central claim is that planet-gas interactions in the late, H2-poor stage of disk evolution are observable in a way that the protoplanetary stage is not. In a protoplanetary disk the gas surface density is high enough that even a deep gap remains luminous, so gaps are hard to see in CO; in a debris disk the same relative gap depth leaves an absolute column inside the gap so small that its emission falls below the ALMA detection threshold. The paper demonstrates this for planets of 0.5, 1, and 5 Jupiter masses at 10, 50, and 100 AU, in disks of $10^{-5}$, $10^{-3}$, and $10^{-1}$ Earth masses, and finds that gap observability is governed by the absolute amount of gas remaining in the gap, not by the relative depth. It condenses this into a criterion: the gap must be resolved, with its width set by four Hill radii of the planet, and the flux from one beam inside the gap, computed from the residual surface density, must lie below $3F_{\\rm sens}$. Applied to HD138813, the criterion predicts that a Jupiter-mass planet at roughly 70 AU would be detectable.","pith_inferences":["Because the gap depth scales as $\\Sigma_{\\rm gap}/\\Sigma_0 \\simeq 1/(1+0.04K)$ with $K=q^2h^{-5}\\alpha^{-1}$, the paper's half-Jupiter threshold is tied to $\\alpha=10^{-3}$ and a thin flaring disk; if independent measurements show $\\alpha\\sim0.1$, the required mass rises by orders of magnitude, so a search should prioritize disks with low inferred turbulence.","The criterion could be inverted: a measured gap flux and gap width in a debris disk would jointly constrain planet mass and local viscosity, complementing dust-gap analyses that cannot easily locate the planet.","An efficient observational test would be a high-resolution ALMA survey of disks with CO masses between $10^{-4}$ and $10^{-2}\\,M_\\oplus$, looking for surface-brightness depressions before attempting deep direct imaging.","Because kinks and gaps are visible in complementary disk-mass regimes, observing both in one disk could distinguish a low-mass planet in a low-mass disk from a massive planet in a massive disk, though the paper does not work out that degeneracy."],"forward_implications":["In a $10^{-3}\\,M_\\oplus$ CO disk at 40 pc, simulated ALMA images show gaps for planets as small as $0.5\\,M_{\\rm J}$ when the planet sits at 50-100 AU; at 10 AU the gap is not resolved.","For $10^{-1}\\,M_\\oplus$ disks, gaps become too luminous to see, but kinks in the channel maps become observable for planets more massive than about $1\\,M_{\\rm J}$.","The flux criterion is disk-agnostic: it applies to any line and any disk, including low-mass protoplanetary disks, provided the gas outside the gap is detected and the gap is resolved.","HD138813, HD121191, and HD156623 are identified as ideal targets; for HD138813 a Jupiter-mass planet at 70 AU would produce an easily resolvable $\\sim$0.14 arcsecond gap in a 7-hour ALMA C-7 observation.","Observed gaps would permit a planet-mass estimate through the Hill-sphere relation $m_p=3M_*(w_{\\rm gap}/4r_p)^3$, probing giant planets too cold for direct imaging."],"supporting_citations":[{"why":"Supplies the FARGO3D code used for all 2D hydrodynamical gap-opening simulations.","marker":"Benítez-Llambay & Masset 2016"},{"why":"Supplies the RADMC-3D radiative transfer code that converts density and velocity fields into CO line emission.","marker":"Dullemond et al. 2012"},{"why":"Supplies the observation simulator used to add realistic ALMA noise and beams and produce the synthetic images.","marker":"CASA Team et al. 2022"},{"why":"Provides the gap-opening framework and the Hill-sphere width scaling used to define resolvable gaps.","marker":"Crida et al. 2006"},{"why":"Provides the analytic gap-depth scaling $K=q^2h^{-5}\\alpha^{-1}$ used to discuss how turbulence and scale height set the planet mass needed.","marker":"Kanagawa et al. 2015"},{"why":"Provides the second-generation gas model, temperature profile, and disk-mass expectations that set the simulated disk structure.","marker":"Kral et al. 2019"},{"why":"Provides the LTE excitation formalism (Eqs. 7-8) used to convert gas mass inside the gap into CO $J=2$-$1$ flux.","marker":"Matrà et al. 2015"},{"why":"Supplies the HD138813 CO flux, mass, and inner dust edge at about 70 AU used for the specific target prediction.","marker":"Hales et al. 2019"}],"fun_headline_variants":["Sparse gas in old disks reveals planets through gaps","Planet-carved gaps in debris gas offer new exoplanet sign","ALMA can see gaps carved by cold Jupiters in debris disks","Half-Jupiter planets expose themselves in debris disk gas","Old disks' thin gas makes planet gaps detectable"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The disk is modeled as a thin, weakly turbulent single-fluid disk with a fixed turbulence parameter $\\alpha=10^{-3}$ and scale height $h=0.01\\,(r/1\\,\\mathrm{AU})^{0.25}$; if real debris-disk gas is substantially more turbulent or thicker, the same planets would carve shallower gaps and the predicted observable masses would rise.","fun_headline_variants_meta":{"raw":{"variants":["Sparse gas in old disks reveals planets through gaps","Planet-carved gaps in debris gas offer new exoplanet sign","ALMA can see gaps carved by cold Jupiters in debris disks","Half-Jupiter planets expose themselves in debris disk gas","Old disks' thin gas makes planet gaps detectable"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000273,"raw_usage":{"total_tokens":1734,"prompt_tokens":1139,"completion_tokens":595,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":755,"completion_tokens_details":{"reasoning_tokens":512}},"tokens_in":755,"tokens_out":595,"duration_ms":5838,"temperature":1.0,"reasoning_tokens":512,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:24:17.595568+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe HD138813's CO $J=2$-$1$ emission with ALMA in the C-7 configuration (beam $\\sim$0.09 arcseconds) for about 7 hours. If the moment-zero map shows no surface-brightness depression centered near 70 AU with width $\\sim$0.14 arcseconds, the paper's headline prediction for a Jupiter-mass planet fails. A second check: measure the turbulent $\\alpha$ of any debris-disk gas, as attempted for $\\beta$ Pictoris; a value near 0.1 would imply that the half-Jupiter threshold is too optimistic by orders of magnitude at the assumed scale height.","supporting_citations":[{"cited_title":"2006, , 181, 587","cited_arxiv_id":null,"evidence_quote":"Provides the gap-opening framework and the Hill-sphere width scaling used to define resolvable gaps."},{"cited_title":"D., Tanaka , H., Muto , T., Tanigawa , T., & Takeuchi , T","cited_arxiv_id":null,"evidence_quote":"Provides the analytic gap-depth scaling $K=q^2h^{-5}\\alpha^{-1}$ used to discuss how turbulence and scale height set the planet mass needed."},{"cited_title":"C., Kama , M., & Matr \\`a , L","cited_arxiv_id":null,"evidence_quote":"Provides the second-generation gas model, temperature profile, and disk-mass expectations that set the simulated disk structure."},{"cited_title":"S., Gorti , U., Carpenter , J","cited_arxiv_id":null,"evidence_quote":"Supplies the HD138813 CO flux, mass, and inner dust edge at about 70 AU used for the specific target prediction."}],"review_version":1}