{"id":"5d124fc0-5b53-4b90-92cc-b3f171b6aefb","arxiv_id":"2505.14856","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":8.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For small polytropic galaxies with a central point mass, the gravitational force from linear perturbations decays as (1+t)^{-b}, with decay order set by initial-data and steady-state regularity.","lead":"This paper proves the first quantitative decay rates for the linearized gravitational Vlasov-Poisson system around nontrivial, compactly supported galaxy equilibria with a central point mass. It shows the gravitational force decays in time like a power law determined by the regularity of the initial data and the steady state.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Decay claim hinges on imported no-eigenvalue theorem and unproved global invertibility of the (y,z)-foliation; both need verification.","rationale":"The reader's weakest assumption identifies exactly the load-bearing dependency: the proof requires the spectral and geometric input that the linearized operator has no embedded eigenvalues and that the (y,z) foliation is a valid global change of variables. My stress-test pass found no internal contradiction in the main resolvent and decay argument, but it confirmed that these inputs are imported rather than proved in the manuscript. Proposition 2.5 is stated as a consequence of [22,44], and the proof of Lemma 3.6 establishes only nonvanishing of the Jacobian, not global injectivity; the latter is plausible from the monotonicity of y in L along constant-z curves, but the argument is not written. Since the central quantitative decay claim collapses if either input fails, and since the paper does not itself verify the citation's exact hypotheses, the correct verdict is conditional acceptance: the mathematical framework is coherent and the novel contributions appear sound, but the theorem's full range is contingent on external results. The reader's conditional verdict therefore remains appropriate, and I recommend no change to it. My concrete test would settle the concern by checking the cited theorems apply verbatim and by supplying the missing global injectivity proof for (3.36).","tokens_in":62376,"tokens_out":23070,"duration_ms":206559,"concrete_test":"Verify that [44, Theorems 6.4.1 and 6.5.5] apply to the exact operator L=T L on the weighted space H with weight |phi'| for all mu>2, nu>1 and 0<=eta<eta_0, re-deriving Proposition 2.5 including absence of embedded eigenvalues. Also check global injectivity of (3.36) by computing d_L y = P_R omega > 0 on each z-level set of J_R; if the cited spectral proof has an unhandled parameter range, or the global diffeomorphism requires an extra condition, the proof of Theorem 1.1 is incomplete.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Theorem 1.1 rests on the unproven-in-this-paper spectral input of Proposition 2.5: sigma(L)=sigma_ess(L)=sigma(T) with no embedded eigenvalues, together with uniform monotonicity T'(I)>=1/c0 and |d_L omega|=O(eta) in (2.12)-(2.13). These feed the Stone formula (2.38), the limiting absorption principle, and the invertibility of the Plemelj operator F^{±,epsilon}_lambda in Proposition 6.12. If the cited [44, Theorems 6.4.1 and 6.5.5] do not apply verbatim to the weighted space H in (1.12) with polytropic exponents mu>2, nu>1, or if an embedded eigenvalue exists, the resolvent bounds (6.48) and (6.57) can acquire point-spectrum contributions and the decay (1.20) fails. The paper does not reproduce or verify this input, so it is a black box. Additionally, Lemma 3.6 only shows the Jacobian determinant of (3.36) is nonzero; it does not explicitly prove the global diffeomorphism property used to define J_R and to justify the change of variables in all Plemelj integrals. For each fixed R this follows from monotonicity of y=omega(z+Psi_L(R),L) in L since d_L y = P_R omega > 0, but the proof is omitted.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies quantitative decay for the linearized gravitational Vlasov-Poisson system around compactly supported polytropic steady states with a central point mass. Theorem 1.1 claims that, for small steady-state amplitude eta, initial data satisfying the orthogonality condition (1.19), and regularity parameter K = min{mu-1, nu, k}, the gravitational force decays as (1+t)^{-b} for 0 <= b <= K-1, with the decay rate improving with the regularity of the data. The proof combines action-angle variables, a Birkhoff-Poincare normal form, a new two-dimensional foliation of the action support, a quantitative pure-transport decay theorem, Stone's formula, and uniform resolvent bounds with Plemelj-type singular integral estimates. The paper is self-contained from Section 3 onward in the sense of proving the new estimates, but it relies on prior results for the spectral structure of the linearized operator and for regularity of action-angle variables.","tokens_in":62673,"tokens_out":4491,"duration_ms":48817,"significance":"If correct, the result is a significant advance: it gives the first quantitative linear decay theorem around nontrivial inhomogeneous gravitational equilibria with trapping, and it introduces a genuinely new mechanism (the (y,z)-foliation and the normal-form cancellations) that is likely to be reusable in related problems such as massive perturbations of the Schwarzschild solution. The central claim is a genuine theorem, not an identity: no parameter is fitted, and the decay rate is derived rather than imposed. The authors are also honest about the main external inputs: the absence of embedded eigenvalues and the spectral description are quoted from [22,44], and several bootstrap steps are compressed. My assessment is that the main architecture is coherent and the new contributions are substantial, but two load-bearing points need to be made fully explicit before the theorem can be considered established: the exact scope of the imported spectral theorem, and the global invertibility of the foliation change of variables.","major_comments":[{"comment":"Proposition 2.5 is load-bearing: the claims sigma(L)=sigma_ess(L)=sigma(T) with no embedded eigenvalues are used in Stone's formula (2.38), in the definition of sigma_L in (2.39), and in every resolvent bound, in particular Theorems 6.13 and 6.14. The proof of Proposition 2.5 consists of a citation to [22,21,44], but the manuscript does not verify that the theorems of [44] apply verbatim to the weighted space H in (1.12), to the polytropic exponents mu>2, nu>1, and to the small-eta regime used here. This is not merely a cosmetic gap: if an embedded eigenvalue existed, the resolvent bounds (6.48) and (6.57) would acquire point-spectrum contributions and the decay (1.20) would fail. Please either state precisely the theorems being cited and check their hypotheses, or give a self-contained argument for the spectral statement in this setting.","section":"Section 2.3, Proposition 2.5 and Eq. (2.29)"},{"comment":"Lemma 3.6 proves only that the Jacobian determinant in (3.38) is nonzero at each point, which gives local invertibility. The global statement that the map (3.36) is a diffeomorphism onto J_R is asserted in Remark 3.10 ('the map ... is a bijection') and is used throughout the paper to justify the change of variables in the Plemelj integrals, e.g. in (5.5), (6.7), and (6.11). The missing step is global injectivity together with a description of the boundary of J_R. For fixed R, injectivity should follow from monotonicity of y=omega(E,L) in L, since d_L y = P_R omega > 0, but the proof is not written. Because all subsequent estimates are formulated in (R,y,z) coordinates, this global fact must be proved explicitly rather than left as an assertion.","section":"Section 3.3, Lemma 3.6 and Remark 3.10"},{"comment":"The bound (7.14) is introduced with the phrase 'by the same argument as in Theorem 6.14' and is then used to justify integration by parts in lambda and the passage to the limit in (7.18). This is a load-bearing estimate: it supplies the |lambda|^{-2} decay that produces the (1+t)^{-b} factor for b>=1. It should be proved explicitly, because the objects being estimated, R^{+,lambda,epsilon}_{a,b}[U^+_epsilon] - R^{-,lambda,epsilon}_{a,b}[U^-_epsilon], mix differences of resolvents and differences of Plemelj kernels. In the same lemma, the inductive proof of (7.11) says 'arguing as in (7.8)' and omits the details of the epsilon -> 0 limit for the higher lambda-derivatives; since (7.11) is used to kill boundary terms in the integration by parts, the argument needs to be written out.","section":"Section 7.1, Eq. (7.14), and Lemma 7.2"}],"minor_comments":[{"comment":"The notation in (6.27)-(6.28), written as omega(I)+lambda+i epsilon/(2 pi m), is ambiguous: it should be omega(I) + (lambda + i epsilon)/(2 pi m), with the parentheses displayed explicitly. The same ambiguity appears in several later formulas.","section":"Throughout, esp. Eq. (6.27)"},{"comment":"The proof of Theorem 5.1 passes an infinite sum over m through the integration by parts in y without an explicit summability statement. The estimate (5.9) suggests the sum is finite, but a sentence justifying the interchange would improve readability.","section":"Theorem 5.1, Eq. (5.9)"},{"comment":"In the definition of Sigma^{pm}, the indicator chi_{[omega(R,z), omega(R,z)]}(y) appears to denote a degenerate interval and is not used in the displayed trace formula. This should be corrected or deleted.","section":"Remark 7.4"},{"comment":"The symbol C is reused with different meanings in Proposition 4.6, e.g. in 'C_{k,j} C^{|m|} |m|^{k+j}'. The constants are independent of m, but the exposition would be clearer if the various constants were labeled distinctly or if the dependence on k,j was stated once.","section":"Section 4.3, Proposition 4.6"},{"comment":"The abstract says 'first linear asymptotic stability result around such equilibria', while Theorem 1.1 proves decay of the gravitational force and density, and Corollary 7.3 provides H-scattering only under the extra assumption K>=3. The wording in the abstract is stronger than the proved statements; I suggest qualifying it as decay of macroscopic quantities rather than full asymptotic stability.","section":"Abstract and Introduction"}],"recommendation":"major_revision","confidential_remarks":"The paper is substantial and the main strategy is credible. My main concern is external verification: the spectral input from [44] is a doctoral thesis and [22] is an accepted paper, but the current manuscript does not spell out how those results apply to exactly the weighted space and polytropic exponents used here. Given the length of the proof and the fact that several bootstraps are compressed, I would advise the editor to request that the authors either verify the cited hypotheses in detail or state the relevant theorems verbatim. The global invertibility issue in Section 3.3 is also correctable and should be addressed. This is a strong candidate for publication after a careful revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth your time. This is the first quantitative decay theorem for the full linearized gravitational Vlasov-Poisson flow around nontrivial, compactly supported polytropic equilibria with a central point mass. The pure transport decay alone beats the previous t^{-2} bound, and the full result requires genuinely new machinery: the Birkhoff-Poincare normal form near the elliptic trapping set, a new (y,z) foliation of the action support, and uniform Plemelj/resolvent estimates. The regularity theory at trapping, especially the product rule with Bell polynomials and the cancellation in Proposition 4.9, is the technical heart and it holds up. The theorem is honest: no fitted parameters, and the decay rate follows exactly from the regularity exponents.\n\nThe soft spots are two, and both are manageable. First, Proposition 2.5, which asserts purely essential spectrum with no embedded eigenvalues, is imported from the authors' earlier work [44]. If those theorems do not apply verbatim to the weighted space H with mu>2, nu>1, the Stone formula and the limiting absorption principle do not close. I have no concrete reason to doubt the transfer, but a referee should verify it carefully. Second, the global invertibility of the (y,z) foliation is a bit terse in Lemma 3.6. The Jacobian computation plus monotonicity in L does give it, but the paper does not spell that out. That is a minor omission, not a defect. A few bootstrap estimates, notably (7.14), are only sketched as 'by the same argument'; I would want those details in a final version, but I found no contradiction.\n\nThis is a paper for kinetic theory and spectral analysis specialists. It deserves a serious referee and a strong journal. I would send it out, asking the referee to focus on the imported spectral theorem and the compressed bootstrap steps.","headline":"A novel and technically demanding proof of quantitative decay for the linearized gravitational Vlasov-Poisson system; the main soft spot is the imported spectral input that needs referee verification.","tokens_in":63182,"tokens_out":2939,"would_cite":true,"duration_ms":37592,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35Q83","35B40","37J40","47A10"],"pacs":[],"model":"deepseek-v4-flash","headline":"Small smooth perturbations of compact galactic equilibria in the gravitational Vlasov-Poisson system decay at quantified rates: the gravitational force decays as a power of time set by the regularity of the data.","keywords":["Vlasov-Poisson system","gravitational relaxation","linear Landau damping","action-angle variables","Birkhoff-Poincaré normal form","limiting absorption principle","stable trapping","resolvent estimates"],"falsifier":"An eigenvalue inside the essential spectrum of the linearised operator for any steady state in the stated family would be a direct obstruction: the Stone-formula representation would then contain a non-decaying oscillatory component, and the gravitational force could not satisfy the claimed $(1+t)^{-b}$ bound uniformly over the full range of $b\\le K-1$. A concrete way to look for it is to scan polytropic parameters near the thresholds $\\mu=2$, $\\nu=1$ or at small but nonzero $\\eta$ for resonant modes.","tokens_in":62196,"feed_emoji":"🌌","tokens_out":13900,"duration_ms":115183,"temperature":0.7,"pith_summary":"This paper proves quantitative decay for the linearised gravitational Vlasov-Poisson system around a family of compactly supported polytropic equilibria, with a fixed point mass at the centre: smooth perturbations whose angular average vanishes produce a gravitational force that decays like $(1+t)^{-b}$ for every $0\\le b\\le K-1$, where $K=\\min\\{\\mu-1,\\nu,k\\}$ is set by the regularity of the steady state and of the initial datum. The result is the first quantitative linear asymptotic stability theorem around such inhomogeneous galactic equilibria, and is framed as the gravitational analogue of linear Landau damping in the radial setting. The proof must overcome two obstructions that do not appear in spatially homogeneous plasma damping: stable trapping, where particles oscillate around elliptic points and resist phase mixing, and frequency mixing, where the gravitational response at one Fourier mode depends on all other modes. A reader should care because the decay is not a formal heuristic but is obtained from a spectral representation plus resolvent bounds, and the same machinery is meant to be the analytic core for nonlinear asymptotic stability.","feed_headline":"Gravitational force decays at quantified rates near galaxy equilibria","feed_subtitle":"Linear stability with quantified decay around compact galaxies, despite trapping and frequency mixing.","key_machinery":"Two coordinate changes and one decomposition carry the argument. Action-angle variables $(\\theta,I)$ convert the transport operator into $\\omega(I)\\partial_\\theta$ with $\\omega(I)=1/T(I)$, so pure transport is solved by oscillating phases $e^{-2\\pi i m t\\omega(I)}$. The new foliation $(R,y,z)=(R,\\omega(I),E-\\Psi_L(R))$ separates the singular behaviour at the trapping set $E=E_L^{\\min}$ from the vacuum boundary, and is invertible when $\\partial_L\\omega=O(\\eta)$. Near the elliptic points, the Birkhoff-Poincaré normal form (a canonical Hamiltonian normal form around a nondegenerate elliptic point) gives explicit action-angle expansions, and the Green's function decomposes as $S_m(\\theta(R,I))=E^{-|m|/2}H_m(R,y,z)$ with $H_m$ analytic; this factorisation cancels the singular factor in the Fourier coefficients of the data and yields the quantitative derivative bound (1.32). The full linearised decay is obtained from Stone's formula by proving uniform resolvent estimates, based on a limiting absorption principle and on Plemelj-type estimates split into near-resonant and non-resonant frequencies.","core_discovery":"The paper's central claim is Theorem 1.1. For polytropic steady states of the form (1.5) with $\\mu>2$, $\\nu>1$ and amplitude $\\eta$ sufficiently small, any initial datum $f_0\\in C^k(\\Omega)$ satisfying the orthogonality condition (1.19) produces a unique solution of the linearised problem (1.13)–(1.14) whose gravitational force obeys $\\|\\partial_R U_{|\\varphi'|f}(t,\\cdot)\\|_{C^{K-1-b}}\\le C\\|f_0\\|_{C^K(\\Omega)}(1+t)^{-b}$ for $0\\le b\\le K-1$, with $K=\\min\\{\\mu-1,\\nu,k\\}$. This bound applies to the full linearised flow, not only to the pure transport part: the gravitational back-reaction is shown to contribute a strictly lower-order term, so the transport mechanism dictates the decay. The theorem also implies, for $k\\ge 2$, decay of the spatial density in $C^{K-2-b}$, and, for $K\\ge 3$, existence of a scattering profile in the natural weighted Hilbert space $H$.","pith_inferences":["An extension the paper leaves implicit is that the decay mechanism should be stable under small perturbations of the steady state that preserve the single-gap spectrum: the proof's central spectral input is absence of embedded eigenvalues, so nearby equilibria in the same regularity class should satisfy the same bounds with modified constants.","The factored Green's function identity suggests a sharper heuristic: the effective transfer function is nearly diagonal in the new $(y,z)$ foliation, so a numerical scheme that discretises $(y,z)$ rather than $(E,L)$ should resolve the decay with fewer modes.","A testable extrapolation is that violating the smallness of $\\eta$ or the monotonicity of $\\omega$ will destroy or slow the decay; constructing an equilibrium where $\\partial_L\\omega$ is not small and measuring the gravitational force would probe the sharpness of the foliation condition.","The method points toward inhomogeneous plasma analogues, where BGK-type equilibria with a trapping region could be treated by the same combination of normal form and resolvent estimates, provided an absence-of-eigenvalues statement is available."],"forward_implications":["The theorem gives explicit, regularity-determined decay rates: smoother steady states and smoother initial data yield faster algebraic decay, with the exponent limited by $K=\\min\\{\\mu-1,\\nu,k\\}$.","For data with $k\\ge 2$, the spatial density inherits the decay in $C^{K-2-b}$, so macroscopic observables relax at quantified rates rather than merely oscillating.","For $K\\ge 3$, the linearised solution scatters to a well-defined profile in the weighted space $H$, which is the linear precursor of a relaxation picture for nearby galaxies.","The result is new even in the pure transport case $\\eta=0$; the mechanism that converts regularity into decay is tied to the point-mass potential and to the three-dimensional geometry, so the same rates should not be expected in one-dimensional confining models.","Quantitative linear decay is an essential stepping stone to nonlinear asymptotic stability, since the same cancellations near trapping control the quadratic terms in a bootstrap argument."],"supporting_citations":[{"why":"Establishes existence, uniform-in-eta regularity, and the single-gap spectral structure of the polytropic steady states, used in Proposition 2.5 and Lemma 2.1.","marker":"[44]"},{"why":"Proves that the linearised operator has purely essential spectrum and no embedded eigenvalues, the spectral input on which the Stone-formula and resolvent analysis rests.","marker":"[22]"},{"why":"Supplies the quantitative pure-transport decay and the Green's function identity (2.42) that the paper extends to the full linearised flow.","marker":"[23]"},{"why":"Provides the Birkhoff-Poincaré canonical normal form with analytic twist used to derive the action-angle expansions and the factored Green's function (Theorem 3.1).","marker":"[35]"},{"why":"Gives the decomposition of the state space as the range of the transport operator and its orthogonal complement, making the angular-average condition (1.19) the correct invariant subspace.","marker":"[21]"}],"fun_headline_variants":["Gravitational decay rates proven for linear Vlasov-Poisson","First linear stability result for compact galaxies with trapping","Gravitational analog of Landau damping in radial setting","Quantitative decay for gravitational potential near equilibria","Overcoming stable trapping: gravitational decay theorem"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the spectral input that the linearised operator has no embedded eigenvalues, together with the geometric input that the frequency map $\\omega(E,L)$ is strictly monotone in $E$ with $\\partial_L\\omega=O(\\eta)$, so that the $(y,z)$-foliation is a true coordinate change; if an embedded eigenvalue existed or the foliation degenerated, the Stone-formula and Plemelj machinery would not deliver the decay estimate.","fun_headline_variants_meta":{"raw":{"variants":["Gravitational decay rates proven for linear Vlasov-Poisson","First linear stability result for compact galaxies with trapping","Gravitational analog of Landau damping in radial setting","Quantitative decay for gravitational potential near equilibria","Overcoming stable trapping: gravitational decay theorem"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000581,"raw_usage":{"total_tokens":2713,"prompt_tokens":903,"completion_tokens":1810,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":519,"completion_tokens_details":{"reasoning_tokens":1733}},"tokens_in":519,"tokens_out":1810,"duration_ms":12772,"temperature":1.0,"reasoning_tokens":1733,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T15:28:21.926619+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An eigenvalue inside the essential spectrum of the linearised operator for any steady state in the stated family would be a direct obstruction: the Stone-formula representation would then contain a non-decaying oscillatory component, and the gravitational force could not satisfy the claimed $(1+t)^{-b}$ bound uniformly over the full range of $b\\le K-1$. A concrete way to look for it is to scan polytropic parameters near the thresholds $\\mu=2$, $\\nu=1$ or at small but nonzero $\\eta$ for resonant modes.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes existence, uniform-in-eta regularity, and the single-gap spectral structure of the polytropic steady states, used in Proposition 2.5 and Lemma 2.1."},{"cited_title":"Damping versus oscillations for a gravitational Vlasov-Poisson system","cited_arxiv_id":"2301.07662","evidence_quote":"Proves that the linearised operator has purely essential spectrum and no embedded eigenvalues, the spectral input on which the Stone-formula and resolvent analysis rests."},{"cited_title":"Quantitative phase mixing for Hamiltonians with trapping","cited_arxiv_id":"2405.17153","evidence_quote":"Supplies the quantitative pure-transport decay and the Green's function identity (2.42) that the paper extends to the full linearised flow."},{"cited_title":"Differential Equations 179, 625–646 (2002)","cited_arxiv_id":null,"evidence_quote":"Provides the Birkhoff-Poincaré canonical normal form with analytic twist used to derive the action-angle expansions and the factored Green's function (Theorem 3.1)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the decomposition of the state space as the range of the transport operator and its orthogonal complement, making the angular-average condition (1.19) the correct invariant subspace."}],"review_version":1}