{"id":"0b868739-70c5-4e15-b57d-3b6f40e31521","arxiv_id":"2509.02884","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Fomalhaut's debris ring has a forced eccentricity that falls with semi-major axis as a^-1.75 +/- 0.16, the first reported eccentricity gradient in a debris disk.","lead":"Astronomers fitted ALMA images of Fomalhaut's dust ring and found the ring's oval shape becomes rounder with distance from the star, a gradient never measured before for such a disk. The finding gives a new way to infer hidden planets and the early history of planetary systems.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Free-eccentricity subspace is not explored: the steep npow=-1.75 is inferred with ep forced to zero, while the full ep+npow model is explicitly deferred; a hidden degeneracy could mimic the gradient.","rationale":"The reader's conditional verdict is appropriate. The strongest evidence is the large Delta-chi^2 against constant-eccentricity models and the visibility-domain check, so the sign is likely real. However, the specific npow value is not secure: the model with ep=0 may absorb width/brightness asymmetries that a free-eccentricity population or a warp would produce, and the paper's own Proper 2 and Verify numbers bracket a shallower gradient. The manuscript is also internally inconsistent in Section 6, reversing the stability result of Section 5.1.2 ('only initially circular disks are stable' vs. the opposite), which weakens the origin discussion but not the detection. My recommended action remains CONDITIONAL, i.e., no change to the reader's verdict: the paper should fit the fully general free+forced eccentricity model and re-derive the model comparison with a valid non-nested statistic before the steepness is adopted.","tokens_in":25751,"tokens_out":7179,"duration_ms":88387,"concrete_test":"Refit the high-resolution image with the full model: ef0, ep, omega_p, npow, and the usual geometric parameters all free, using the same mask and likelihood. Then examine the marginalized posterior of npow and the joint (npow, ep) profile likelihood. If npow remains about -1.75 with ep consistent with 0, the concern is resolved; if the posterior shifts toward -1 or broadens past about -1.3, the steep gradient is an artifact of forcing ep=0. Report a non-nested comparison (AIC/BIC or k-fold cross-validation) because the General and Proper models are not nested.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The signature of the gradient is read out of a parametric model that fixes the free eccentricity to zero (ep=0). Section 3.4.2 fits 'Proper' models that include ep only at fixed npow=0 and -1, and explicitly defers the general case with both ep and npow free. This is the load-bearing gap because ep and wr are degenerate (the paper itself notes the wide wr posteriors), and a small free-eccentricity population can produce the same qualitative signatures as a negative forced-eccentricity gradient: broader pericenter and brighter apocenter. The two anchor points leave the magnitude unsettled: Proper 2 (npow=-1, ep~2%) is only Delta-chi^2=21.5 worse than General (about 2.6 sigma by their Eq. 7), and the low-resolution Verify fit gives npow=-1.16 +/- 0.16, about 2.6 sigma shallower than -1.75. The negative sign of the gradient is probably robust, but the headline steepness, and hence the quantitative 'first gradient' claim, is conditional on an untested part of parameter space.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper fits high-resolution ALMA 1.32 mm continuum images of Fomalhaut with parametric RADMC-3D models in which the forced eccentricity varies as e_f ∝ a^{n_pow}, and reports n_pow = -1.75 ± 0.16 in the 'General' model. The model is preferred over constant-eccentricity models with free eccentricity (Proper 1, Δχ²=70) and over n_pow=-1 with free eccentricity (Proper 2, Δχ²=21.5), with a visibility-domain check. The paper interprets the preferred model as evidence for 'Eccentric Velocity Divergence,' proposes single-planet 'Gap' and 'Resonant' scenarios, and uses REBOUND simulations to argue that an initially eccentric disk survives 440 Myr.","tokens_in":26101,"tokens_out":5023,"duration_ms":58457,"significance":"If the steep negative gradient is real, this is the first reported eccentricity gradient in a debris disk and a potentially valuable diagnostic of planet-disk interactions. The parametric model is physically motivated, the EVD kinematic identity in §2 is clearly derived, the public code release supports reproducibility, and the MCMC/visibility comparisons are extensive. However, the headline steepness is not yet established at the claimed confidence: the free-eccentricity subspace is not explored, the lower-resolution verification gives a shallower slope, and the Δχ² significance conversion is statistically questionable. The sign of the gradient is likely robust; the magnitude and the quantitative 'first gradient' claim need further work.","major_comments":[{"comment":"The General model fixes e_p=0; Proper 1 and Proper 2 fix n_pow=0 or -1. The general case with both e_p and n_pow free is explicitly deferred to future work. Since e_p and w_r are degenerate (the Proper fits show wide w_r posteriors), and a free-eccentricity population can produce qualitatively similar pericenter broadening and apocenter brightening, the fitted n_pow=-1.75 may be biased. This is the central load-bearing gap. Fit the full model with e_p free and n_pow free, or provide injection-recovery tests showing that n_pow is unbiased under e_p>0, before claiming a quantitative gradient.","section":"§3.4.2, Eq. (4), Table 1"},{"comment":"The quoted 6.7σ and 2.6σ preferences use a chi-square CDF with Npar=9 degrees of freedom. General, Proper 1, and Proper 2 are not nested models and have the same number of free parameters, so this p-value conversion is not valid. Report an appropriate model-comparison criterion (e.g., BIC/AIC, cross-validation, or a genuinely nested sequence with a full model) instead. Without this, the statistical preference for the General model is overstated.","section":"§3.4.3, Eq. (7)"},{"comment":"The Verify fit to the lower-resolution data gives n_pow=-1.16±0.16, about 2.6σ shallower than the high-resolution -1.75. Section 3.5 acknowledges the high-resolution value is 4.6σ from -1 while the low-resolution value is consistent with -1. Thus the sign of the gradient is robust, but the magnitude is not established. The abstract and §6 present -1.75 as the headline result; they should be conditioned on resolution and on the unmodeled e_p subspace.","section":"§3.4.1 vs Appendix B and Table 1"},{"comment":"The fitting mask selects primary beam response >0.66 and projected radius <200 au, i.e., the bright ansae. The residual maps show the model is fainter than the data along the minor axes, which are excluded from the fit. Appendix A tests only one alternative mask and reports a worse apocentre fit without providing Δχ² or the resulting parameter changes. Demonstrate that n_pow and the model preference are robust to mask choice, or quantify how the mask affects the inference.","section":"§3.3, Appendix A, Eq. (6)"}],"minor_comments":[{"comment":"The abstract cites 'Lynch & Lovell 2022' while the body and reference list use Lynch & Lovell (2021). Unify the citation.","section":"Abstract and References"},{"comment":"Define the incomplete gamma function and use standard regularized notation. The exponent Npar/2 with Npar=9 appears ad hoc for comparing non-nested models.","section":"Eq. (7)"},{"comment":"The second width in the footnote is written 'wr, in = 15.7±0.7 au'; this should presumably be the outer width, wr, out.","section":"Footnote 5"},{"comment":"The sentence 'only scenarios with an initially circular disk are stable/finalise as eccentric disks' contradicts §5.1.2, which finds that circular disks are disrupted. It should read 'initially eccentric'.","section":"§6"},{"comment":"The extracted figure annotation appears to show 'npow = 1.1651' with no minus sign, while the axis runs from -1.5 to -0.75; check the sign convention in the figure. Also, Table 1 lists Proper 2 e_p as '<3.6' without the nominal 1.9±0.9% given in the text; make consistent.","section":"Fig. 6 and Table 1"}],"recommendation":"major_revision","confidential_remarks":"The paper is well-written and the data analysis is substantial; the central claim is plausible but the quantitative gradient is conditional on an untested e_p+n_pow model and a statistically questionable model comparison. I would not reject it, but the authors should be asked to close the free-eccentricity gap and re-derive the significance before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the negative eccentricity gradient in Fomalhaut is real in sign, but the headline value npow = -1.75 is not yet closed. The paper deserves a careful referee.\n\nWhat's actually new: the first fit of a freely varying power-law eccentricity index to Fomalhaut's ALMA ring, and a serious model-comparison effort. The General model beats constant-eccentricity plus free-eccentricity models by large Delta-chi^2 in the image domain, and a visibility-domain check gives the same ordering. The authors also ship the code and do a lower-resolution verification fit, both of which strengthen reproducibility.\n\nSoft spots, in order of importance. First, the General model sets free eccentricity ep=0. The full model with both ep and npow free is explicitly deferred. Since ep and wr are degenerate (the paper notes the wide wr posteriors), a small free-eccentricity population could mimic part of the signature; Proper 2 (npow=-1, ep~2%) is only Delta-chi^2=21.5 worse, about 2.6 sigma by their Eq. 7. The sign is robust - Proper 1 (constant e, ep free) is far worse - but the magnitude is conditional. Second, the verification fit to lower-resolution data gives npow=-1.16+-0.16, about 2.6 sigma shallower than -1.75. The authors are transparent about this, but it means the steep value is not yet pinned down. Third, the fitting mask drops the minor-axis regions where the model is fainter than the data; the residual there is acknowledged and looks noise-dominated, but it does mean the fit is driven by the two ansae. Fourth, there's a sentence in Section 6 that inverts the N-body result: it says only initially circular disks are stable and finalise as eccentric, when the simulations show the opposite. That's a typo, but it will confuse readers. Finally, the N-body simulations start with an initially eccentric disk using the same power law they fit; that is a stability test, not an independent formation test, and the paper mostly says so.\n\nWho benefits: debris disk observers and anyone working on planet-disk secular interactions. This is a well-executed parametric modeling paper with honest caveats. My recommendation: send it to peer review. The core claim about a negative gradient is likely correct; the paper should be asked to either fit the joint ep+npow model or state clearly why that is infeasible, and to soften the 'first gradient' claim until the magnitude is more secure.","headline":"Fomalhaut's ring really needs a negative eccentricity gradient, but the steepness (-1.75) is not yet closed - the sign looks robust, the magnitude is conditional.","tokens_in":26604,"tokens_out":3760,"would_cite":true,"duration_ms":43772,"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":"This paper claims that Fomalhaut's debris disk is best described by a forced eccentricity that falls steeply with semi-major axis, ef ∝ a^-1.75±0.16, the first eccentricity gradient reported in a debris disk.","keywords":["circumstellar disks","debris disks","eccentricity","eccentricity gradient","Fomalhaut","planet-disk interactions","Eccentric Velocity Divergence","millimeter imaging"],"falsifier":"Measure the local radial width and surface brightness separately at the two ansae in a deeper millimeter image that also resolves the minor axis. The EVD model with npow = -1.75 predicts a pericenter width about 1.2 times the apocenter width and a fainter, broader pericenter; observing the opposite width asymmetry, or minor-axis emission that no apse-aligned power-law model can reproduce, would rule out the claimed gradient.","tokens_in":25643,"feed_emoji":"🪐","tokens_out":9792,"duration_ms":104923,"temperature":0.7,"pith_summary":"The paper argues that the eccentricity of Fomalhaut's main debris belt is not constant but drops steeply with distance from the star, following a power law ef ∝ a^-1.75 ± 0.16. This is presented as the first reported eccentricity gradient in any debris disk. The gradient produces, through mass continuity, exactly the asymmetries seen in the millimeter images: a broader, fainter pericenter and a narrower, brighter apocenter. If correct, the belt's shape encodes the presence and orbit of an unseen planet, and the eccentricity likely originated during the protoplanetary disk phase.","feed_headline":"Fomalhaut's disk eccentricity drops steeply with radius","feed_subtitle":"Millimeter-wave images favor a -1.75 power law, the first eccentricity gradient seen in a debris disk.","key_machinery":"The load-bearing object is a single-Gaussian-ring parametric disk model whose forced eccentricity follows a power law in semi-major axis, ef(a) = ef,0(a/a0)^npow, with aligned pericenters (∂ωf/∂a = 0). The essential identity is the Jacobian j = [1 - e(e + a∂e/∂a)]/√(1 - e²) · (1 - q cos E), with q = (a∂e/∂a)/[1 - e(e + a∂e/∂a)], which converts the eccentricity gradient into surface-density and width asymmetries. This 'Eccentric Velocity Divergence' is what lets the single parameter npow explain both the broad-and-faint pericenter and the narrow-and-bright apocenter seen by the millimeter observations.","core_discovery":"On the paper's own terms, the central discovery is that Fomalhaut's main belt is modeled better by a forced eccentricity profile ef(a) = ef,0(a/a0)^npow with npow = -1.75 ± 0.16 than by any constant-eccentricity model, with or without a free-eccentricity component. The gradient is not a mere curve-fit extra: it simultaneously reproduces two observed asymmetries that a constant eccentricity cannot, namely a pericenter that is broader and fainter and an apocenter that is narrower and brighter. The mechanism behind this, called Eccentric Velocity Divergence, is mass continuity in an apse-aligned eccentric disk: a negative eccentricity gradient converts the orbit packing into exactly those width","pith_inferences":["The EVD relation predicts a quantitative link between the sign of a disk's width asymmetry and the sign of its brightness asymmetry; a survey of resolved eccentric debris disks could measure a population of eccentricity-gradient slopes and test whether values near -1 (single planet) or steeper (self-gravity) dominate.","Because the fit masks the minor axis and assumes zero warp, an alternative reading of the data is a disk with a mild warp or a wider free-eccentricity distribution; a visibility-domain fit that includes minor-axis emission and allows a warp could separate these.","If the 'born eccentric' conclusion holds, young debris disks around A-type stars should show steeper eccentricity gradients than old disks, since long-term planet-disk interactions would tend to flatten or re-process the primordial profile; this is a testable age trend with deeper imaging.","The small apocenter over-subtraction noted in the residuals suggests the true profile may be even steeper or non-Gaussian; a two-sided radial-width model on deeper data could determine whether the -1.75 slope is absorbing an asymmetric radial profile."],"forward_implications":["Future models of Fomalhaut's main belt must treat eccentricity as a radially varying parameter; a single forced eccentricity underfits the data by a large margin.","The measured slope is steeper than the classical a^-1 expected for a massless disk forced by an internal planet, so either the disk is massive enough for self-gravity to steepen the profile, or the perturber sits closer than the simple gap-carving geometry suggests.","The 440-Myr N-body integrations show that initially circular disks are disrupted, while initially eccentric disks survive, implying the ring's lopsidedness was present since the protoplanetary disk phase.","A single unseen planet, either a gap-carver at roughly 109-120 au or a 2:1 resonant-clearing planet at 70-75 au, can sculpt the belt's inner edge and gaps while remaining below current detection limits, making the eccentricity gradient a probe of such planets.","The lower-resolution archival data also require a negative gradient (npow = -1.16 ± 0.16), so the result does not depend on the highest-resolution image alone."],"supporting_citations":[{"why":"Supplies the parametric eccentricity-gradient disk models (npow = 0, -1, +1) that this work generalizes to arbitrary npow.","marker":"Lynch & Lovell (2021)"},{"why":"Provides the calibrated, high-resolution millimeter image that is the dataset fitted here.","marker":"Chittidi et al. (2025)"},{"why":"Prior fit establishing the pericenter-broader width asymmetry; the model here reproduces and supersedes its residuals.","marker":"Kennedy (2020)"},{"why":"Earlier millimeter observations and modeling used as the low-resolution verification dataset.","marker":"MacGregor et al. (2017)"},{"why":"Provides the free-eccentricity modeling method used to test whether a proper-eccentricity population can mimic the gradient.","marker":"Lovell & Lynch (2023)"},{"why":"Classical secular theory giving the ef ∝ a^-1 slope for an internal planet, the baseline against which -1.75 is compared.","marker":"Murray & Dermott (1999)"},{"why":"Shows disk self-gravity can steepen forced eccentricity profiles toward a^-4.5, the physics invoked to explain a gradient steeper than -1.","marker":"Sefilian (2024)"},{"why":"Establishes that eccentric protoplanetary disks require eccentricity gradients to avoid differential precession, supporting the primordial-origin argument.","marker":"Teyssandier & Ogilvie (2016)"},{"why":"Infrared observations of the intermediate belt whose edges and eccentricities constrain the two single-planet scenarios.","marker":"Gáspár et al. (2023)"},{"why":"The N-body integrator used for the 440-Myr stability simulations.","marker":"Rein & Liu (2012)"}],"fun_headline_variants":["Disk eccentricity gradient discovered at Fomalhaut","Fomalhaut's lopsided ring explained by eccentricity gradient","Eccentricity gradient shapes Fomalhaut's debris disk","Fomalhaut's ring eccentricity drops steeply outward"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The model assumes the belt is a single apse-aligned Gaussian ring with no warp and negligible free eccentricity, while masking out the minor-axis regions where the model is faintest; if a warp or free-eccentricity population is present, the inferred npow ≈ -1.75 could be biased.","fun_headline_variants_meta":{"raw":{"variants":["Disk eccentricity gradient discovered at Fomalhaut","Fomalhaut's lopsided ring explained by eccentricity gradient","Eccentricity gradient shapes Fomalhaut's debris disk","Fomalhaut's ring eccentricity drops steeply outward"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000659,"raw_usage":{"total_tokens":2905,"prompt_tokens":855,"completion_tokens":2050,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":599,"completion_tokens_details":{"reasoning_tokens":1978}},"tokens_in":599,"tokens_out":2050,"duration_ms":20930,"temperature":1.0,"reasoning_tokens":1978,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T11:19:34.794252+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the local radial width and surface brightness separately at the two ansae in a deeper millimeter image that also resolves the minor axis. The EVD model with npow = -1.75 predicts a pericenter width about 1.2 times the apocenter width and a fainter, broader pericenter; observing the opposite width asymmetry, or minor-axis emission that no apse-aligned power-law model can reproduce, would rule out the claimed gradient.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the calibrated, high-resolution millimeter image that is the dataset fitted here."},{"cited_title":"D., & Dermott , S","cited_arxiv_id":null,"evidence_quote":"Classical secular theory giving the ef ∝ a^-1 slope for an internal planet, the baseline against which -1.75 is compared."}],"review_version":1}