{"id":"0f380441-2e03-42cf-9313-02d313064d21","arxiv_id":"2411.16071","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For small-energy Hénon-Heiles orbits, the slow variables u and w rotate with explicit formulas (17) and (58) that stay accurate for n ε^(5/2) << 1, far beyond the ε^(-2) validity of naive perturbation series.","lead":"The paper derives explicit asymptotic formulas for the slow drift of the Hénon-Heiles Hamiltonian at small energies, with rigorously controlled errors over times much longer than standard perturbation theory allows. A smart generalist would read it as a case study in how two-scale dynamics can turn a famous non-integrable system into a tractable slow-oscillation problem.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Theorem 3.1(iii) rests on a bootstrap whose only source of uniform O(ε³) remainder control is Lemma 4.2, and Lemma 4.2 is a two-sentence proof sketch; the uniformity in n up to N ε₀³ ≤ K₀ is the unverified crux. Keep conditional.","rationale":"I read the proof in detail rather than scanning. The theorem's meaningful content is the rotation law (13)–(16) for the slow variables over n ~ ε⁻³ loops, and the chain is: complex-contour Poincaré map (§4.1) → ε-expansion with explicit v[1], w[1], v[2], w[2] and remainder equations (31) → one-loop contraction (Lemma 4.1) → per-loop remainder bounds (Lemma 4.2) → recurrence estimates (Lemma 4.3) → bootstrap (Lemma 4.4) → Theorem 3.1(iii). I verified as much of the algebra as is checkable by hand: (78) follows from (19)–(21); (82) differentiates to (78); the one-loop increments in §4.2.2 follow from the arcsin-2π monodromy and the trivial monodromy of (v₀−y²)^{3/2} after winding both branch points; the Tₙ recurrence (62), the bound (65), and the phase computation (68)–(72) are algebraically consistent up to absorbed constant factors. Lemma 4.1's contraction proof is detailed and sound. The bootstrap Lemma 4.4 is logically sound conditional on Lemma 4.2. What remains unproved is Lemma 4.2 — the uniform bound on the ε³ remainders along every loop, with constants independent of the cycle index — and the reader's weakest_assumption identifies exactly this same point. I do not see a contradiction or a concrete failure mode; the numerics (Figure 4, at nε³ = 1.024, ε = 0.01) provide independent support that the actual remainders are far smaller than the allowed O(1) bound. So the concern does not warrant rejection or an unverdictable status; it warrants keeping the CONDITIONAL verdict: Theorem 3.1(iii) should be accepted conditional on a completed proof of Lemma 4.2, or on a verified uniform remainder bound of the kind tested above. One additional presentation issue that does not affect (iii): Theorem 3.1(i) claims uₙ = O(nε⁴) when ẏ₀ = 0, but then u₀ = h ≠ 0, so the statement as printed is inconsistent; presumably vₙ = O(nε⁴) or uₙ − h = O(nε⁴) was intended, matching Appendix D.","tokens_in":21669,"tokens_out":47676,"duration_ms":382901,"concrete_test":"Extract per-loop remainders numerically. Fix h = 0.1, section y(0) = 0, and solve (4) with the §5 protocol (WorkingPrecision 32, MaxSteps 10⁸) for ε = 10⁻² and 5×10⁻³, at T₀ = u₀²+w₀² = 0.0085 (the Figure 4 case), 0.005, and 0.0099 (near the boundary h² = 0.01). At each return to y = 0 compute Δvₙ = v_{n+1} − vₙ and Rₙ = [Δvₙ − (14π/3)ε²wₙ√(2hvₙ−vₙ²−wₙ²)]/ε³, and Sₙ similarly, for n up to nε³ ≈ 1. Accept Lemma 4.2 only if maxₙ|Rₙ| and maxₙ|Sₙ| are bounded by one constant across both ε values and all three T₀ values. Growth with n, or divergence like ε⁻¹ as ε → 0, would show the uniform remainder bound in (49) — the crux of Lemma 4.4's bootstrap — does not hold as stated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim Theorem 3.1(iii) is proved in §4.7.3 by combining Lemmas 4.3 and 4.4. Lemma 4.4 is a bootstrap whose sole source of the uniform bound |Rₙ|, |Sₙ| ≤ M in the recurrence (49) is Lemma 4.2, applied at every cycle (uₙ, wₙ) in the annulus c₁ ≤ Tₙ ≤ c₀². Lemma 4.2's proof, however, is a sketch: 'the structure of the integral operators (31) is similar to that of (40). One difference is the appearance of arcsin(y/√v₀) ... but this is regular on C_{v0}. The same arguments go through straightforwardly.' The sketch never establishes the specific quantitative content the bootstrap needs: (i) sup-norm bounds on the O(1) sources f₁, g₁ of (31), which contain v[2], w[2] from (84)–(85) with arcsin(y/√v₀) and (v₀−y²)^{±3/2} terms that are ramified along the two-circle contour; (ii) a contraction estimate for (I₁, I₂) on an O(1)-ball — the remainders R, S are O(1), not small in ε — requiring the Taylor remainder (30) to remain analytic throughout that ball, a genuinely different operator family from J in Lemma 4.1; (iii) uniformity in the cycle index n up to N ε₀³ ≤ K₀ and in (uₙ, wₙ) over the annulus. Because Lemma 4.4 applies Lemma 4.2 about ε₀⁻³ times, any hidden dependence of the Lemma 4.2 constants on the cycle would break the bootstrap and with it Theorem 3.1(iii). Independent checks support the rest of the proof: formula (82) for v[1] differentiates to (78), the monodromy bookkeeping in §4.2.2 (arcsin gains 2π, (v₀−y²)^{3/2} has trivial monodromy around both branch points) is consistent, and the algebra of Lemma 4.3 (62)–(72) checks out up to constant-tracking sloppiness (a missing 1/√T₀ factor in the error term of (71), absorbed into 'depends only on c₀, c₁'). The gap is localized and plausible, but it is load-bearing.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies the Hénon–Heiles system in the small-energy regime after rescaling, and derives long-time asymptotic formulas for the slow variables u = h - (y^2 + \\dot y^2) and w = \\dot x\\dot y + xy. The main result, Theorem 3.1(iii), states that for initial data with u_0^2 + w_0^2 < h^2, the iterated Poincaré map of the slow variables satisfies u_n + i w_n = \\sqrt{u_0^2+w_0^2} e^{i\\varphi_n}(1 + n\\varepsilon^3\\eta_n) with an explicit phase \\varphi_n and uniformly bounded error, for n up to \\varepsilon_0^{-3}. For n\\varepsilon^{5/2}\\ll 1 the phase reduces to a uniform rotation. The proof uses a complex-plane Poincaré map, an \\varepsilon-expansion with remainders, and a bootstrap argument; the paper also includes explicit second-order expansions and a numerical comparison.","tokens_in":22135,"tokens_out":22999,"duration_ms":198205,"significance":"If the proof gap discussed below is filled, the result is significant: it gives rigorous asymptotic control of the slow dynamics on a timescale n\\sim\\varepsilon^{-3}, which goes well beyond the standard perturbation horizon n\\varepsilon^2=O(1). The derivation is self-contained and has no fitted parameters; the explicit formulas in Appendices B and C and the independent numerical check in Figures 4 and 6 are genuine strengths. The main bottleneck is the proof of Lemma 4.2, which is sketched rather than proved and is load-bearing for Theorem 3.1(iii).","major_comments":[{"comment":"The proof of Lemma 4.2 is only a two-sentence sketch ('the structure of the integral operators ... is similar', 'the same arguments go through'), but this lemma is the sole source of the uniform remainder bounds |R_n|, |S_n|\\le M used in the bootstrap Lemma 4.4. A complete proof needs to establish several quantitative facts that are not immediate from Lemma 4.1: (i) uniform sup-norm bounds on the source terms f_1,g_1 in (31), which contain v[2], w[2] from (84)-(85) with arcsin(y/\\sqrt{v_0}) and (v_0-y^2)^{\\pm 3/2} along the two-circle contour C_{v_0}; (ii) a contraction estimate for (I_1,I_2) in an O(1)-ball of the Banach space B, noting that R,S are O(1), not small in \\varepsilon, so the Taylor remainder Q_F,Q_G in (30) must be controlled throughout that ball; (iii) uniformity of all constants in the initial data (v_0,w_0) ranging over the annulus c_1\\le (h-v_0)^2+w_0^2\\le c_0^2, since Lemma 4.4 re-applies Lemma 4.2 about \\varepsilon_0^{-3} times. Without this uniformity the bootstrap in Lemma 4.4, and with it Theorem 3.1(iii), does not follow as written.","section":"4.5, Lemma 4.2"},{"comment":"The statement 'if \\dot y_0=0, then w_n=O(n\\varepsilon^3) and u_n=O(n\\varepsilon^4)' is inconsistent with the definition u=h-v. Since \\dot y_0=0 gives v_0=0 and hence u_0=h, the claim u_n=O(n\\varepsilon^4) fails already at n=0 for every h>0. Appendix D likewise gives v(t)=O(\\varepsilon^2)+O(t\\varepsilon^4), so the correct statement is of the form u_n=h+O(\\varepsilon^2)+O(n\\varepsilon^4), not u_n=O(n\\varepsilon^4). Please correct the theorem statement and the corresponding proof summary in \\S4.7.1.","section":"3.2, Theorem 3.1(i)"}],"minor_comments":[{"comment":"The paragraphs 'Solutions along one loop' and 'Solutions Along One Loop' are duplicated nearly verbatim; one copy should be removed.","section":"4.4-4.5"},{"comment":"Both subsections are titled 'Proof of (iii)'; the second should refer to the simplified phase formula (16) or to the n\\varepsilon^{5/2} regime, not repeat the same title.","section":"4.7.3-4.7.4"},{"comment":"The displayed recursion for B_n is typeset in a way that is hard to parse: the product index, the \\sqrt{T_0} factor, and the relation between \\tilde C_\\ell and C_\\ell should be written out explicitly, because the telescoping identity \\prod \\tilde C_\\ell=(\\sqrt{T_0}/\\sqrt{T_n})\\prod C_\\ell is used implicitly to pass to (71).","section":"4.6, around Eq. (69)"},{"comment":"The caption reports n=1024000 and \\varepsilon=0.01, giving n\\varepsilon^3=1.024; since the theorem guarantees the estimates for N\\varepsilon_0^3\\le K_0, please state the value of K_0 for these parameters or note explicitly that the plotted data lie in the asymptotic (not necessarily proven) regime.","section":"5, Figure 4"},{"comment":"In the sentence immediately after (16), 'Formula (14) simplifies to (16)' is correct, but the line 'For n slightly smaller, such that n\\varepsilon^{5/2}\\ll 1' should be reconciled with the qualitative claim in Remark 3.2 that this range can still be much larger than the secular range n\\varepsilon^2=O(1); a sentence clarifying the ordering would help.","section":"3.2, Eq. (14)"}],"recommendation":"major_revision","confidential_remarks":"The main gap is the missing proof of Lemma 4.2; given that the bootstrap in Lemma 4.4 applies this lemma repeatedly, the paper needs a complete proof of Lemma 4.2, or at least a detailed appendix establishing the uniform bounds and contraction estimates, before the central claim can be accepted. The rest of the proof structure is coherent and the final formulas appear correct modulo the theorem-statement typo in item (i), so I view this as repairable rather than rejectable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The real news here is Theorem 3.1(iii): for small ε, the slow variables (u,w) rotate at a slowly varying rate with amplitude fixed to O(nε^3) over n ~ ε^{-3} periods. Formulas (13)-(17) and the simplified (58) are explicit, parameter-free, and I don't see them anywhere in the previous literature. The numerical agreement in Figures 4 and 6 is real evidence, not a fit.\n\nThe proof strategy is honest: pass to the Poincaré map in y, expand to ε^2 with remainders, then bootstrap the remainders along iterations. The algebra in Appendices B and C is heavy; I spot-checked some terms and they are consistent. The monodromy bookkeeping (arcsin picking up 2π, (v0-y^2)^{3/2} unramified) is careful.\n\nThe soft spot is Lemma 4.2. Its proof is two sentences: 'similar to Lemma 4.1', plus a remark that arcsin is regular on Cv0. This lemma is the only source of the uniform O(1) bounds on the remainders R,S that Lemma 4.4 applies ε0^{-3} times. The stress-test is right that the sketch never establishes the needed quantitative content: sup-norm bounds on f1,g1 (which have (v0-y^2)^{-3/2} terms), contraction on an O(1)-ball, and uniformity in the cycle index n across the annulus. If any constant in Lemma 4.2 depends on the base point (vn,wn), the bootstrap collapses. I want to be clear: this looks fixable — the operators really are similar — but as written the central theorem leans on a placeholder. A referee should ask for a full proof.\n\nMinor issues: Theorem 3.1(i) likely has a typo (un = O(nε^4) looks suspect), and equation (71) drops a 1/√T0 factor into the constant. Neither affects the main claim.\n\nWho gets value: people in Hamiltonian perturbation theory and celestial mechanics, and anyone teaching slow-fast systems. The paper deserves a serious referee, not a desk rejection. I'd send it out with a request to expand Lemma 4.2 and check the typo. After that, I'd cite it.","headline":"A serious, likely correct paper on long-time Hénon-Heiles asymptotics; the one load-bearing lemma is a sketch and needs a real proof before the main theorem can be fully trusted.","tokens_in":22777,"tokens_out":3884,"would_cite":false,"duration_ms":36180,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["34E10","37J40","70H05"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves an explicit Poincaré-map formula showing that the Hénon–Heiles slow variables rotate with nearly constant amplitude for times of order ε^{-3}.","keywords":["Hénon–Heiles system","slow variables","Poincaré map","long-time asymptotics","adiabatic invariants","small-energy dynamics","secular perturbation theory","two-scale dynamics"],"falsifier":"Take $\\varepsilon=0.01$, $h=0.1$, and the initial conditions of Figure 4; integrate with a high-accuracy solver and record $u_n,w_n$ at every return to $y=0$. The theorem predicts $\\sqrt{T_n}=\\sqrt{T_0}(1+O(n\\varepsilon^3))$, so the relative deviation after $n=10^5$ loops must be bounded by a fixed constant times $0.1$; if instead it grows linearly in $n\\varepsilon^2$, reaching order $10$ at that $n$, the claim is falsified.","tokens_in":21453,"feed_emoji":"🔄","tokens_out":13893,"duration_ms":104846,"temperature":0.7,"pith_summary":"The paper proves that at small energies the Hénon–Heiles system has a pair of slow observables whose motion is an almost uniform rotation over a time of order $\\varepsilon^{-3}$. Specifically, the iterated Poincaré map on the section $y=0$ satisfies $u_n+iw_n=\\sqrt{u_0^2+w_0^2}\\,e^{i\\varphi_n}(1+n\\varepsilon^3\\eta_n)$ with bounded $\\eta_n$, for all $n$ up to order $\\varepsilon^{-3}$. This extends far beyond the time horizon $t=O(\\varepsilon^{-2})$ at which the ordinary perturbation series develops secular terms. The result gives rigorous asymptotic constants of motion that explain the numerically observed slow drift and filling of the phase-space region. If correct, it provides an explicit formula for long-time evolution in a benchmark non-integrable Hamiltonian system.","feed_headline":"Hénon–Heiles slow variables rotate for ε^{-3} loops","feed_subtitle":"An explicit map fixes the amplitude of u+iw to within O(n ε^3), far past the usual perturbation limit.","key_machinery":"The carrying object is the iterated Poincaré map of the slow variables with respect to the section $y=0$, computed by integrating the complexified equations around the contour $C_{v_0}$ that encircles the two branch points $\\pm\\sqrt{v_0}$ of $\\sqrt{v-y^2}$. The one-loop map produces the exact second-order coefficients displayed in (33), and the reduced recurrence $u_{n+1}=u_n-\\varepsilon^2 w_n\\sqrt{h^2-u_n^2-w_n^2}-\\varepsilon^3 R_n$, $w_{n+1}=w_n+\\varepsilon^2 u_n\\sqrt{h^2-u_n^2-w_n^2}+\\varepsilon^3 S_n$ carries the slow dynamics. The recursion is a rotation by angle $\\varepsilon^2\\sqrt{h^2-T_n}$ plus a small remainder, and the bootstrap in Lemma 4.4 keeps the remainders $R_n,S_n$ bounded uniformly in $n$ as long as $n\\varepsilon^3\\le K_0$; this uniform remainder control is what extends validity beyond the perturbation-series horizon.","core_discovery":"The central discovery is that, for initial conditions with $y(0)=0$, $\\dot y_0\\ne0$, and $u_0^2+w_0^2<h^2$, the slow variables $u=h-(y^2+\\dot y^2)$ and $w=\\dot x\\dot y+xy$ do not drift chaotically on short timescales. Instead, after each return to the section $y=0$ they rotate: $u_n+iw_n=\\sqrt{T_0}\\,e^{i\\varphi_n}(1+n\\varepsilon^3\\eta_n)$, where $T_0=u_0^2+w_0^2$, $|\\eta_n|\\le M$, and $\\varphi_n$ is given by $\\varphi_0+\\frac{14\\pi}{3}\\varepsilon^2\\sum_{k=0}^{n-1}\\sqrt{h^2-T_k}$. For $n\\varepsilon^{5/2}\\ll1$ the phase simplifies to $\\varphi_0+n\\varepsilon^2\\sqrt{h^2-T_0}+O(n^2\\varepsilon^5)$, so the radius $\\sqrt{T_n}$ is conserved up to $O(n\\varepsilon^3)$ while the angle advances at a slowly varying rate. The proof runs through a Poincaré map in complex time, integrating once around the branch points $y=\\pm\\sqrt{v_0}$; the one-loop change of the slow variables is exactly of order $\\varepsilon^2$ with coefficients $\\frac{14\\pi}{3}w_0\\sqrt{2hv_0-v_0^2-w_0^2}$ and $\\frac{14\\pi}{3}(h-v_0)\\sqrt{2hv_0-v_0^2-w_0^2}$, and the $\\varepsilon^3$ remainder is controlled uniformly for $N\\varepsilon^3\\le K_0$ by a bootstrap.","pith_inferences":["A natural next step, not taken in the paper, is to carry the expansion one order higher: if the remainders can be controlled at $\\varepsilon^4$, the same contour method should yield validity up to $n\\varepsilon^4=O(1)$ with a corrected phase.","Because the proof only uses the rational structure of the integrands in $x$, $\\sqrt{v-y^2}$, and $S$, the same two-scale Poincaré-map construction may apply to other two-degree-of-freedom Hamiltonians with cubic nonlinearities and a resonant harmonic limit.","The conservation of $T_n$ suggests that the level sets $T=\\text{constant}$ act as approximate invariant tori for small $\\varepsilon$; a direct numerical test would be to measure how long a trajectory stays within $O(\\varepsilon)$ of such a level set.","The phase formula (14) depends only on $T_k$, so it can be iterated as a cheap map even when the simplified formula (16) is no longer accurate; comparing the two against numerics would locate the exact $n\\varepsilon^{5/2}$ threshold."],"forward_implications":["The standard $\\varepsilon$-expansion of $x(t),y(t)$ develops secular terms when $t\\varepsilon^2=O(1)$; the theorem implies the slow-variable rotation remains accurate for $n$ of order $\\varepsilon^{-3}$, extending the reliable time horizon by a factor of $\\varepsilon^{-1}$.","For $n\\varepsilon^{5/2}\\ll1$, formulas (16) and (58) give explicit trigonometric expressions for $u_n$ and $w_n$, so long-time evolution can be predicted without solving the differential equations over that entire interval.","The squared radius $T_n=u_n^2+w_n^2$ is an approximate adiabatic invariant: it changes by at most $O(n\\varepsilon^3)$, and in the simplified regime the amplitude of the slow oscillation is constant to that accuracy.","In the degenerate cases $\\dot y_0=0$ or $u_0=w_0=0$, the slow variables stay $O(n\\varepsilon^3)$, so the system exhibits no slow rotation beyond the trivial drift.","The numerical comparisons in Section 5 show the formulas tracking $v_n$ and $w_n$ over more than a million loops at $\\varepsilon=0.01$, in agreement with the predicted error size."],"supporting_citations":[{"why":"introduces the Hénon–Heiles model and the long-time drift and mixing observed numerically, which the present theorem makes rigorous.","marker":"[14]"},{"why":"introduces the approximate-adiabatic-invariant method with complex-time Poincaré maps that the paper adapts to this system.","marker":"[10]"},{"why":"companion source of the same method, cited as the origin of the approach used here.","marker":"[11]"},{"why":"standard theory of the Poincaré map used to define the iterated returns of the slow variables.","marker":"[8]"}],"fun_headline_variants":["Explicit map steers Hénon–Heiles slow variables past ε^{-3}","Hénon–Heiles rotation law: invariants hold for ε^{-3} loops","Slow variables in Hénon–Heiles don't drift; they rotate","New asymptotic invariants for small-energy Hénon–Heiles","Explicit Poincaré map extends Hénon–Heiles perturbation validity"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the $\\varepsilon^3$ remainder terms $R_n,S_n$ in the slow-variable recurrence stay bounded by a constant depending only on the initial distance from the boundary $h^2$, uniformly over all loops with $N\\varepsilon^3\\le K_0$; the whole theorem rests on that bootstrap succeeding.","fun_headline_variants_meta":{"raw":{"variants":["Explicit map steers Hénon–Heiles slow variables past ε^{-3}","Hénon–Heiles rotation law: invariants hold for ε^{-3} loops","Slow variables in Hénon–Heiles don't drift; they rotate","New asymptotic invariants for small-energy Hénon–Heiles","Explicit Poincaré map extends Hénon–Heiles perturbation validity"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000795,"raw_usage":{"total_tokens":3562,"prompt_tokens":1069,"completion_tokens":2493,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":685,"completion_tokens_details":{"reasoning_tokens":2388}},"tokens_in":685,"tokens_out":2493,"duration_ms":17264,"temperature":1.0,"reasoning_tokens":2388,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T13:35:12.540014+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take $\\varepsilon=0.01$, $h=0.1$, and the initial conditions of Figure 4; integrate with a high-accuracy solver and record $u_n,w_n$ at every return to $y=0$. The theorem predicts $\\sqrt{T_n}=\\sqrt{T_0}(1+O(n\\varepsilon^3))$, so the relative deviation after $n=10^5$ loops must be bounded by a fixed constant times $0.1$; if instead it grows linearly in $n\\varepsilon^2$, reaching order $10$ at that $n$, the claim is falsified.","supporting_citations":[{"cited_title":"and Heiles, C., 1964","cited_arxiv_id":null,"evidence_quote":"introduces the Hénon–Heiles model and the long-time drift and mixing observed numerically, which the present theorem makes rigorous."},{"cited_title":"and Huang, M., 2016","cited_arxiv_id":null,"evidence_quote":"introduces the approximate-adiabatic-invariant method with complex-time Poincaré maps that the paper adapts to this system."},{"cited_title":"and Huang, M., 2015","cited_arxiv_id":null,"evidence_quote":"companion source of the same method, cited as the origin of the approach used here."},{"cited_title":"and Teichmann, T., 1956","cited_arxiv_id":null,"evidence_quote":"standard theory of the Poincaré map used to define the iterated returns of the slow variables."}],"review_version":1}