{"id":"a0a69115-1d0b-4c5c-b552-6eecb4aa6fda","arxiv_id":"2607.08750","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":4,"one_line_summary":"Secular perturbations from an inner planet create sharp, marching peaks and caustic discontinuities in a debris disk's azimuthally-averaged surface density, with radial locations following a predictable power-law pattern.","lead":"This paper shows that a planet orbiting inside a debris disk creates a predictable pattern of peaks and gaps in the disk's radial density profile as gravity slowly reshapes particle orbits over millions of years. A smart generalist might read it to understand how unseen planets can be detected and characterized just from the wavy structure they imprint on a disk of dust and rock.","discovery_kind":"unclear","skeptic_critique":{"model":"glm-5.2","headline":"Post-hoc Gaussian convolution of ASD is not equivalent to computing ASD with an eccentricity distribution; the detectability claims in Section 8 rest on an unvalidated approximation that DebrisPy could test exactly but does not.","rationale":"The paper's core analytical contribution — the identification of ASD peaks at eccentricity nulls, jumps at caustics, and the radial spacing pattern a_k ∝ (t/k)^{2/7} — is well-supported. The mathematics in Appendix C is detailed and self-consistent, the Monte Carlo verification shows excellent agreement (2.88% RMS deviation), and the key observational prediction (equation 29) is falsifiable and robust to random eccentricity since the e-nulls are set by the deterministic secular evolution. The DebrisPy software is publicly released, adding reproducibility. The concern I identify is narrower than the reader's: it is not that random eccentricity exists (the paper acknowledges this), but that the specific approximation used to assess its impact — post-hoc Gaussian convolution of the output ASD — is not equivalent to the exact calculation, and the discrepancy is unquantified in the regime where it matters most. This affects the detectability claims but not the existence or pattern of the features themselves. The paper is appropriately cautious in its language ('can still be detectable in high-resolution observations') and acknowledges the limitation. The verdict of ACCEPT is appropriate; the concern is a gap in validation rather than a flaw in the central argument. A future study using DebrisPy's distribution capability would close this gap cleanly.","tokens_in":39975,"tokens_out":2795,"duration_ms":94417,"concrete_test":"Use DebrisPy's eccentricity-distribution capability (equation A1 with ψ_e as a Rayleigh distribution centered on e(a,t) with dispersion σ_e,0/a) to recompute the ASD profiles shown in Figure 8 for e_p = 0.2, t = 15t_sec and t = 40t_sec, with σ_e,0 = 0.1a_p and 0.25a_p. Compare the exact result to the post-hoc Gaussian convolution shown in the figure. If the peak contrasts in the exact calculation drop by more than ~30% relative to the convolved profiles, the detectability thresholds are materially optimistic and the observational claims of Section 9.1 need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The reader correctly identifies the unique-e(a) assumption as the weakest link. I would sharpen the concern as follows. In Section 8, the paper models the effect of random free eccentricity by convolving the exact ASD (computed with deterministic e(a)) with a 1D Gaussian of width σ_e,0·r. This is not the same as computing ASD from first principles using an eccentricity distribution ψ_e(e,a) — which equation (A1) and DebrisPy can in principle handle. The mapping from e(a) to Σ̄(r) in equation (14) is nonlinear (it involves a Heaviside constraint and a square root), so convolving the output ASD is not generally equivalent to integrating over a distribution of inputs. The discrepancy matters most precisely when σ_e is comparable to e_f — the regime where features are most threatened. If the exact calculation shows features are washed out more aggressively than the convolution predicts, the detectability thresholds in Section 7–8 (e.g., the claim that 2–3 peaks survive for e_p ≥ 0.2 at σ = 0.25a_p) would be optimistic. The paper itself notes (Section 9.3) that DebrisPy can handle eccentricity distributions, so the tool to test this exists but is not exercised. This does not undermine the core analytical results (the peak-location pattern of equation 19 is robust since random eccentricity doesn't shift the deterministic e-nulls), but it does affect the observational detectability claims that form part of the paper's stated contribution.","agreement_with_reader":"agree"},"referee_report":{"model":"glm-5.2","summary":"This manuscript studies the azimuthally-averaged surface density (ASD) of a debris disc secularly perturbed by an inner planet. Building on the exact analytical framework of Rafikov (2023), the authors show that secular evolution of the particle eccentricity profile e(a) produces sharp ASD features: weakly singular (logarithmically divergent) peaks at eccentricity nulls and finite discontinuous jumps at caustic (tangent) points. These features march outward through the disc as it ages, with their radial locations following the well-defined pattern a_k proportional to (t/k)^{2/7}. The authors develop detailed asymptotic analysis (Appendix C) for the conditions under which peaks form (parameter zeta_k > 1) versus mere bumps, and for the amplitudes and widths of the caustic jumps. They verify all analytical predictions against a new numerical tool, DebrisPy, which is made publicly available. The paper explores parameter dependence (planetary eccentricity, free eccentricity fraction, disc mass profile), discusses observational detectability under PSF convolution and random eccentricity, and illustrates how detection or non-detection of secular ASD features can constrain planetary mass, semi-major axis, and eccentricity.","tokens_in":40857,"tokens_out":1360,"duration_ms":237745,"significance":"The paper makes a solid contribution to the debris-disc dynamics literature. The analytical results are derived from first principles: the peak-location pattern (Eq. 19), the peak-to-bump transition criterion (Eq. 22), the jump amplitudes (Eq. C21), and the pedestal overlap threshold (Eq. C33) are all falsifiable predictions verified against DebrisPy calculations. The public release of DebrisPy as a reproducible computational tool is a genuine strength. The proposed observational test (Eq. 29) for identifying secular features via their radial spacing pattern is concrete and immediately applicable to ALMA/JWST data. The framework for constraining planetary parameters from both detection and non-detection (Section 9.1-9.2, Figure 9-10) is well-constructed and practical.","major_comments":[{"comment":"Section 8: The effect of random free eccentricity on ASD is modeled by convolving the exact ASD (computed with deterministic e(a)) with a 1D Gaussian of width sigma_e,0 * r. This is not equivalent to computing ASD from first principles using an eccentricity distribution psi_e(e,a), which equation (A1) and DebrisPy can in principle handle. The mapping from e(a) to Sigma_bar(r) in equation (14) is nonlinear (it involves a Heaviside constraint and a square root), so convolving the output ASD is not generally equivalent to integrating over a distribution of inputs. The discrepancy matters most when sigma_e is comparable to e_f, which is precisely the regime where features are most threatened. The paper itself notes (Section 9.3) that DebrisPy can handle eccentricity distributions, so the tool to test this exists but is not exercised. This does not undermine the core analytical results (the e","section":null}],"minor_comments":[{"comment":"Section 2.2, equation (9): The general eccentricity solution is presented here with f_free defined as e_free/e_f, but the specific case f_free = 1 is not adopted until Section 4. It would help the reader to state upfront in Section 2.2 that f_free = 1 will be the default for most of the paper.","section":null},{"comment":"Section 4.1, Figure 1: The Monte Carlo comparison uses N_p = 10^8 particles and reports 2.88% RMS fractional deviation. It would be useful to state the computational time for this comparison to give readers a sense of DebrisPy's efficiency advantage.","section":null},{"comment":"Section 7: The 1D Gaussian convolution used to model PSF effects is an approximation, since the actual observational procedure involves 2D convolution of the sky image followed by azimuthal averaging. The authors should briefly note this approximation and its expected validity for moderately inclined discs.","section":null},{"comment":"Section 9.2, Figure 10: The threshold psi_e = 2 AU is described as 'chosen rather arbitrarily and only for illustration.' Given that this parameter determines the boundary of the excluded grey region, a brief justification of why 2 AU is a reasonable order-of-magnitude choice (or a note that it should be calibrated per system) would help readers gauge sensitivity.","section":null},{"comment":"Equation (32): The notation [(i+k_0)^2 r_i^{7/2}]^{1/2} is slightly ambiguous. It should be clarified whether the exponent 1/2 applies to the product (i+k_0)^2 * r_i^{7/2} or just to r_i^{7/2}. Based on the derivation, it appears to be (i+k_0) * r_i^{7/2}.","section":null},{"comment":"Appendix C1.1, equation (C11): The constant C(a_1, a_2, zeta_k) depends on integration endpoints a_1, a_2 that are defined only implicitly. A brief note on how these are determined (or a reference to where this is specified) would be helpful.","section":null},{"comment":"Section 6.3, Figure 7: For the narrow Gaussian ring models (center and right columns), the radial range is different from the left column, which can cause confusion. Adding a panel showing the underlying e(a) profile for the narrow ring case, or at least noting the different radial scale, would improve clarity.","section":null},{"comment":"The paper uses both 'ASD' and 'Sigma_bar(r)' interchangeably. While defined in Section 3, a brief reminder at first use in later sections would aid readability.","section":null}],"recommendation":"minor_revision","confidential_remarks":"The stress-test concern about the post-hoc Gaussian convolution in Section 8 is valid and is the most substantive technical issue in the paper. However, it affects the observational detectability claims rather than the core analytical results, which are mathematically rigorous and well-verified. The peak-location pattern (Eq. 19) is robust to random eccentricity since it depends only on the deterministic secular precession rate. I recommend minor revision with the request that the authors either validate the convolution approximation with at least one exact calculation using DebrisPy's eccentricity-distribution capability, or clearly scope the detectability claims as approximate pending such validation. The paper is otherwise well-executed and appropriate for publication."},"author_rebuttal":null,"desk_editor":{"model":"glm-5.2","letter":"The main thing to know: this paper takes the exact ASD framework from Rafikov (2023) and works out, in considerable analytical detail, what happens when a debris disc is secularly perturbed by an inner planet. The result is a set of sharp ASD features — logarithmically divergent peaks at eccentricity nulls, finite jumps at caustic points — whose radial locations follow a_k ∝ (t/k)^{2/7}. That spacing pattern (equation 29) is the paper's most valuable output: it's a clean, falsifiable diagnostic that could distinguish secular features from gaps carved by embedded planets. The asymptotic analysis in Appendix C is genuinely careful, with predictions for peak amplitudes (C11), jump amplitudes (C21), and pedestal widths (C31) all verified against DebrisPy calculations. The code is public, the Monte Carlo cross-check is convincing, and the non-detection framework (Section 9.2) is a nice bonus that turns a null result into useful constraints on M_p a_p^2 and e_p a_p. The application to HD 107146 and the ARKS sample, even though it yielded no detections, shows the method is ready for real data. The stress-test concern about Section 8 lands. The paper models random free eccentricity by convolving the deterministic ASD with a 1D Gaussian, which is not equivalent to computing ASD from an eccentricity distribution ψ_e(e,a). The ASD integral is nonlinear — it has a Heaviside constraint and a square root — so output convolution ≠ input distribution integration. This matters most when σ_e ~ e_f, which is precisely the regime where features are most threatened. The irony is that DebrisPy can handle eccentricity distributions natively (equation A1), so the exact test is available but not performed. That said, this gap affects the detectability thresholds in Sections 7–8, not the core analytical results. The peak-location pattern of equation 19 is robust because random eccentricity doesn't shift the deterministic nulls. The observational claims would be more credible with one figure showing the exact distribution calculation alongside the convolution approximation. This is a well-executed paper with a real new result, a public tool, and one soft spot that's fixable with a calculation the authors already have the machinery to do. It deserves a serious referee who should ask them to run the exact ψ_e calculation in Section 8.","headline":"Solid analytical work on debris disc substructure with one testable gap in the observational claims","tokens_in":40768,"tokens_out":1481,"would_cite":true,"duration_ms":70120,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"glm-5.2","headline":"Planet-induced density waves in debris discs trace unseen planets","keywords":["debris disc","secular perturbation","planetary dynamics","surface density","exoplanet detection","celestial mechanics","debris disc substructure","planet-disc interaction"],"falsifier":"If the random eccentricity dispersion of debris particles is comparable to or larger than the forced eccentricity from the planet, the sharp ASD peaks and jumps predicted by the theory would be washed out below the detection threshold, making the predicted (t/k)^(2/7) spacing pattern unobservable and undermining the planet-constraining claims.","tokens_in":40260,"feed_emoji":"🪐","tokens_out":1232,"duration_ms":181201,"temperature":0.7,"pith_summary":"When a planet on an eccentric orbit sits inside a debris disc, its gravity slowly winds up the eccentricities of debris particles over millions of years. This paper shows that this process imprints a characteristic pattern of sharp peaks and steps into the disc's azimuthally-averaged surface density (ASD) — the radial density profile that observers routinely extract from images. The peaks appear at radii where particle eccentricities momentarily drop to zero, and the steps appear where particle orbits pile up at their closest or farthest points from the star. Crucially, the radial locations of these peaks follow a precise geometric progression scaling as (t/k)^(2/7), where t is the system age and k is an integer counting peaks from the outside in. This pattern is a fingerprint of secular planetary perturbation and nothing else, so detecting it in a debris disc would point to an unseen planet and immediately constrain the combination of planetary mass and semi-major axis. The authors build an analytical theory of how these features form and evolve, verify it with a new public numerical tool, and show that even after the smoothing effects of finite telescope resolution and random particle eccentricities, the features should remain detectable in high-resolution observations. They further argue that even a non-detection is informative: it rules out planets in a well-defined region of parameter space.","feed_headline":"Density peaks in debris discs trace unseen planets","feed_subtitle":"A planet's slow gravitational winding of debris orbits creates a fingerprint pattern of density peaks whose spacing reveals the planet'smass","key_machinery":"The central mechanism is the secular precession of debris particle eccentricities driven by an interior planet. Each particle's eccentricity oscillates in time at a rate proportional to a^(-7/2), passing through zero at semi-major axes a_k that march outward as the system ages. The ASD is computed via an exact integral (Equation 14) that sums contributions from all particle orbits crossing a given radius, with the integrand controlled by the interplay between the eccentricity profile e(a) and a geometric function κ(r,a) = |1 − r/a|. Peaks arise where e(a) vanishes (eccentricity nulls) and the parameter ζ_k — the ratio of eccentricity gradient to κ-gradient at the null — exceeds unity. Jumps,","core_discovery":"Secular perturbation by an inner planet on an eccentric orbit produces a train of density peaks in a debris disc's azimuthally-averaged surface density at radii a_k = a_in × (t / t_sec × k)^(2/7), and the detection or non-detection of this pattern constrains the perturbing planet's mass, semi-major axis, and eccentricity.","pith_inferences":["If multiple planets orbit interior to a debris disc, their combined secular perturbations could produce eccentricity profiles with more complex null patterns, potentially breaking the clean (t/k)^(2/7) spacing rule — the absence of this rule in observed discs might itself signal multi-planet architecture.","The pedestal overlap phenomenon at high null order suggests that very old or very massive planets could produce discs with featureless inner regions but structured outer regions, creating a radial gradient of secular signature visibility that evolves with system age.","If future high-resolution surveys of debris discs systematically fail to find the predicted peak spacing pattern, this would either imply that most debris discs are not secularly perturbed by interior planets on eccentric orbits, or that collisional and stirring processes randomize eccentricities faster than secular oscillations can organize them."],"forward_implications":["Observing a sequence of density peaks in a debris disc following the (t/k)^(2/7) spacing rule would provide a mass-times-semi-major-axis-squared measurement of an unseen planet, independent of direct imaging.","Combining the secular ASD constraint on M_p × a_p^2 with a stellar radial-velocity or astrometric acceleration measurement (which constrains M_p / a_p^2) would uniquely determine the planet's mass and orbital distance.","A non-detection of secular ASD features in a well-resolved disc excludes planets in a calculable region of the (M_p × a_p^2, e_p × a_p) parameter space, guiding target selection for direct imaging campaigns.","Narrow debris rings observed with ALMA or JWST that show skewed or multi-peaked radial profiles may be explained by a single interior planet rather than requiring multiple distinct planetesimal belts.","The publicly released DebrisPy tool enables forward-modeling of ASD profiles for any assumed eccentricity distribution, allowing observers to test planetary hypotheses against measured radial density profiles."],"fun_headline_variants":["Debris disc density peaks mark hidden planets","Secular planet perturbations carve density peaks in debris discs","Debris disc peak spacing encodes planet mass and orbit","Marching density peaks in debris discs trace perturbing planets","Planet's slow warp of debris orbits leaves detectable peak pattern"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The analytical theory assumes that every debris particle at a given orbital distance has a single, deterministic eccentricity value set by secular evolution alone, with no random scatter. Real debris discs have eccentricity dispersions from self-stirring, collisions, and gravitational stirring by embedded bodies, which the paper handles only by post-hoc smoothing of the idealized result rather than by incorporating the dispersion into the core calculation.","fun_headline_variants_meta":{"raw":{"variants":["Debris disc density peaks mark hidden planets","Secular planet perturbations carve density peaks in debris discs","Debris disc peak spacing encodes planet mass and orbit","Marching density peaks in debris discs trace perturbing planets","Planet's slow warp of debris orbits leaves detectable peak pattern","Density peak positions in debris discs constrain unseen planets","Debris discs record planet perturbations as density peak trains","Peak spacing in debris discs constrains undetected planets","Density peaks march outward as planets perturb debris discs","Secular debris disc features fingerprint perturbing planets"]},"model":"glm-5.2","effort":"low","cost_usd":0.0,"raw_usage":{"total_tokens":1335,"prompt_tokens":558,"completion_tokens":777,"prompt_tokens_details":null},"tokens_in":558,"tokens_out":777,"duration_ms":45105,"temperature":1.0,"reasoning_tokens":679,"cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-10T01:49:18.924446+00:00","model_set":{"reader":"glm-5.2"},"falsifier":"If the random eccentricity dispersion of debris particles is comparable to or larger than the forced eccentricity from the planet, the sharp ASD peaks and jumps predicted by the theory would be washed out below the detection threshold, making the predicted (t/k)^(2/7) spacing pattern unobservable and undermining the planet-constraining claims.","supporting_citations":[],"review_version":1}