pith. sign in

arxiv: 2605.01864 · v1 · submitted 2026-05-03 · 🧮 math.NA · cs.NA

Numerical Construction of Elliptic Lower-Dimensional Quasi-Periodic Solutions with a Priori Bound

Pith reviewed 2026-05-09 16:34 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords quasi-periodic solutionsnumerical constructionalternating schemeMelnikov conditionsGevrey decaylower-dimensional toriresolvent identitynearly integrable systems
0
0 comments X

The pith

Alternating numerical scheme extends to elliptic lower-dimensional quasi-periodic solutions while controlling inverses via resolvent identity.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper extends an alternating numerical scheme, previously used for full-dimensional quasi-periodic solutions, to construct elliptic lower-dimensional ones in nearly integrable systems. It incorporates the first and second Melnikov conditions to remove small divisors arising from normal frequencies and resonances. The perturbation is handled as a general real-valued function with Gevrey decay, without assuming a Hankel structure. By applying the resolvent identity to simplify multi-scale analysis, the global inverse is expressed as a linear combination of local inverses, with weak short-range and strong long-range interactions that keep the decay uniformly controlled. Numerical tests on the Hénon-Heiles and Fermi-Pasta-Ulam models show the scheme produces solutions with a priori bounds.

Core claim

The authors establish that the alternating scheme extends to elliptic lower-dimensional quasi-periodic solutions by using the resolvent identity to represent the global inverse of the perturbation operator as a linear gluing of local inverses. This yields a regime-dependent interaction structure—weak at short range, strong at long range—that preserves Gevrey decay uniformly, even for real-valued perturbations lacking Hankel structure. With the Melnikov conditions eliminating small divisors, the method delivers numerical constructions with a priori bounds, as verified on the Hénon-Heiles and FPU models.

What carries the argument

The resolvent identity applied in multi-scale analysis, which linearly combines local inverses into a global one whose Gevrey decay stays controlled through short-range weak and long-range strong interactions.

If this is right

  • The scheme applies to general real-valued perturbations, not only polynomials.
  • Numerical solutions are obtained for the Hénon-Heiles model and the Fermi-Pasta-Ulam model with demonstrated effectiveness.
  • Gevrey decay of the inverse remains uniformly controlled despite the lack of Hankel structure.
  • Small divisors are eliminated once the first and second Melnikov conditions are met.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The linear gluing representation may simplify similar inversion problems in other KAM constructions for non-polynomial perturbations.
  • The approach could be tested on additional nearly integrable systems with known lower-dimensional tori to check uniformity of the decay control.
  • Practical use in problems like the restricted three-body problem becomes feasible if the Melnikov conditions can be verified numerically.

Load-bearing premise

The frequencies must satisfy the first and second Melnikov conditions and the perturbation operator must exhibit Gevrey decay for the uniform control of the inverse to hold.

What would settle it

A long-time numerical integration of the constructed Hénon-Heiles solution that shows secular drift or violation of the a priori bound would falsify the claim that the scheme produces valid elliptic lower-dimensional quasi-periodic solutions.

Figures

Figures reproduced from arXiv: 2605.01864 by Bin Shi, Mingwei Fu.

Figure 1
Figure 1. Figure 1: Geometric illustration of resonance regions. In (b), the solid lines denote the boundaries view at source ↗
Figure 2
Figure 2. Figure 2: Phase space trajectories of the H´enon-Heiles model with markers indicating specific time view at source ↗
Figure 3
Figure 3. Figure 3: Convergence profiles across iterations for the H´enon-Heiles model, where the top and view at source ↗
Figure 4
Figure 4. Figure 4: Phase space trajectories of the FPU model with markers indicating specific time points view at source ↗
Figure 5
Figure 5. Figure 5: Convergence profiles across iterations for the FPU model. The panels represent different view at source ↗
Figure 6
Figure 6. Figure 6: Phase space trajectories of the FPU model with markers indicating specific time points view at source ↗
read the original abstract

A numerical framework for constructing full-dimensional quasi-periodic solutions in nearly integrable systems was recently developed by Fu and Shi[2026]. Based on an alternating scheme, this approach effectively overcomes the secular drift in angle variables, a fundamental limitation of symplectic integrators. However, in many applications, such as the restricted three-body problem, lower-dimensional quasi-periodic solutions hold greater significance. The construction of these solutions is considerably more challenging due to the presence of normal frequencies, leading to intricate resonance phenomena. Beyond the subspace resonance, one must also account for the first and second Melnikov conditions to eliminate small divisors. In this study, we extend the proposed alternating numerical scheme to compute the elliptic lower-dimensional quasi-periodic solutions. Numerical experiments are presented for the H\'{e}non-Heiles model and the Fermi--Pasta--Ulam (FPU) model, demonstrating the effectiveness of the proposed method. Furthermore, we emphasize that the perturbation is not merely a polynomial with real coefficients but is a real-valued function. As a result, the associated perturbation operator exhibits Gevrey decay without possessing a Hankel structure. Meanwhile, we further simplify the multi-scale analysis by exploiting the resolvent identity, showing that the global inverse can be expressed linearly in terms of local inverses via the gluing procedure. This representation reveals a regime-dependent interaction structure: weak interactions dominate at short range, while strong interactions emerge at long range. This balance ensures that the Gevrey decay of the inverse remains uniformly controlled. Moreover, within this linear representation, the inversion conditions provide a clearer characterization of the localization properties.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper extends the alternating numerical scheme introduced by Fu and Shi (2026) for full-dimensional quasi-periodic solutions to the elliptic lower-dimensional case in nearly integrable Hamiltonian systems. It incorporates the first and second Melnikov conditions to handle normal frequencies and small divisors, presents numerical experiments on the Hénon-Heiles and Fermi-Pasta-Ulam models, and analyzes the perturbation operator for real-valued (non-polynomial) functions exhibiting Gevrey decay without Hankel structure. The multi-scale analysis is simplified via the resolvent identity and a gluing procedure that expresses the global inverse linearly in terms of local inverses, revealing regime-dependent weak-to-strong interactions that preserve uniform Gevrey decay.

Significance. If the numerical effectiveness and a priori bound hold, the work provides a practical extension of KAM-type numerical methods to lower-dimensional tori, which are central in applications such as the restricted three-body problem. The emphasis on real-valued perturbations and the resolvent-based gluing for controlling inverses without Hankel structure offers a technical advance over prior polynomial-focused schemes; reproducible code or explicit bounds would strengthen its utility for the field.

major comments (2)
  1. The abstract and title claim an 'a priori bound' and 'effectiveness' for the extended scheme, yet no explicit derivation, theorem statement, or quantitative error measures (e.g., residual norms, convergence rates, or comparison to the bound) are referenced; this is load-bearing for the central numerical claim and requires a dedicated section or theorem to verify uniform control of the Gevrey decay under the gluing procedure.
  2. The weakest assumption—that the perturbation is a general real-valued function with Gevrey decay (no Hankel structure) and that Melnikov conditions hold—is stated but not verified numerically or analytically for the chosen models; § on numerical experiments should include a check that the first and second Melnikov conditions are satisfied for the computed solutions to support the small-divisor elimination.
minor comments (2)
  1. The abstract contains a minor formatting inconsistency in 'Hénon-Heiles' (appears with escaped LaTeX in one place); ensure consistent rendering throughout.
  2. The description of 'regime-dependent interaction structure' (weak at short range, strong at long range) is conceptually useful but would benefit from a brief equation or diagram illustrating the linear representation of the inverse via gluing.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript extending the alternating numerical scheme to elliptic lower-dimensional quasi-periodic solutions. We address the major comments point by point below and will revise the manuscript to strengthen the presentation of the a priori bounds and numerical validations.

read point-by-point responses
  1. Referee: The abstract and title claim an 'a priori bound' and 'effectiveness' for the extended scheme, yet no explicit derivation, theorem statement, or quantitative error measures (e.g., residual norms, convergence rates, or comparison to the bound) are referenced; this is load-bearing for the central numerical claim and requires a dedicated section or theorem to verify uniform control of the Gevrey decay under the gluing procedure.

    Authors: We agree that an explicit theorem statement would improve clarity. The multi-scale analysis in the theoretical sections derives the a priori bound on the Gevrey decay of the inverse operator via the resolvent identity and gluing procedure, establishing uniform control through the regime-dependent weak-to-strong interactions. However, this is not condensed into a standalone theorem, and the numerical experiments lack explicit quantitative measures such as residual norms or convergence rates. In the revision, we will add a dedicated theorem summarizing the bound and include quantitative error analysis (e.g., residual norms and comparison to the predicted decay) in the numerical experiments section. revision: yes

  2. Referee: The weakest assumption—that the perturbation is a general real-valued function with Gevrey decay (no Hankel structure) and that Melnikov conditions hold—is stated but not verified numerically or analytically for the chosen models; § on numerical experiments should include a check that the first and second Melnikov conditions are satisfied for the computed solutions to support the small-divisor elimination.

    Authors: The first and second Melnikov conditions are assumed as part of the standard setup for elliptic lower-dimensional tori to eliminate small divisors, consistent with the KAM framework, and the models are selected with frequencies in the perturbative regime where these are expected to hold. Explicit numerical verification was not included in the current experiments section. We will add a check in the revised numerical experiments (e.g., computing minimal distances to resonances for the obtained solutions) to support the small-divisor control. revision: yes

Circularity Check

0 steps flagged

No significant circularity in derivation chain

full rationale

The paper extends the authors' prior alternating numerical scheme from the 2026 reference but introduces independent analytical steps: simplification of multi-scale analysis via the resolvent identity, a gluing procedure expressing the global inverse linearly in terms of local inverses, and control of Gevrey decay for real-valued (non-polynomial) perturbations lacking Hankel structure. These elements, along with invocation of first and second Melnikov conditions and numerical validation on the Hénon-Heiles and FPU models, do not reduce any claimed result to a fitted parameter, self-definition, or unverified self-citation chain. The derivation remains self-contained with externally demonstrable components.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the system being nearly integrable, the perturbation satisfying Gevrey decay properties as a real-valued function, and the validity of first and second Melnikov conditions to control small divisors; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (2)
  • domain assumption The Hamiltonian system is nearly integrable with a perturbation that is a real-valued function exhibiting Gevrey decay without Hankel structure.
    Invoked to justify the extension of the alternating scheme and the multi-scale analysis.
  • domain assumption The first and second Melnikov conditions hold to eliminate small divisors arising from normal frequencies.
    Required for the construction of elliptic lower-dimensional quasi-periodic solutions.

pith-pipeline@v0.9.0 · 5590 in / 1582 out tokens · 33170 ms · 2026-05-09T16:34:06.843684+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

92 extracted references · 92 canonical work pages

  1. [1]

    1993 , publisher=

    Craig, Walter and Wayne, C Eugene , journal=. 1993 , publisher=

  2. [2]

    Construction of quasi-periodic solutions for

    Bourgain, Jean , journal=. Construction of quasi-periodic solutions for

  3. [3]

    Quasi-periodic solutions of

    Bourgain, Jean , journal=. Quasi-periodic solutions of. 1998 , publisher=

  4. [4]

    , title =

    De Vogelaere, R. , title =. 1956 , number =

  5. [5]

    Ruth, R. D. , title =. IEEE Transactions on Nuclear Science , volume =

  6. [6]

    , title =

    Feng, K. , title =. Proceedings of the Fifth International Symposium on Differential Geometry and Differential Equations, August 1984, Beijing , publisher =

  7. [7]

    , title =

    Feng, K. , title =. Journal of Computational Mathematics , volume =

  8. [8]

    Numerische Mathematik , volume=

    Volume-preserving algorithms for source-free dynamical systems , author=. Numerische Mathematik , volume=. 1995 , publisher=

  9. [9]

    Proceedings of the International Congress of Mathematicians: Rio de Janeiro 2018 , pages=

    Dynamics, numerical analysis, and some geometry , author=. Proceedings of the International Congress of Mathematicians: Rio de Janeiro 2018 , pages=. 2018 , organization=

  10. [10]

    Calvo, M. P. and Sanz-Serna, J. M. , title =. Numerical Analysis 1991 , series =

  11. [11]

    , title =

    Hairer, E. , title =. Annals of Numerical Mathematics , volume =

  12. [12]

    Numerische Mathematik , volume=

    The life-span of backward error analysis for numerical integrators , author=. Numerische Mathematik , volume=. 1997 , publisher=

  13. [13]

    Geometric numerical integration illustrated by the

    Hairer, Ernst and Lubich, Christian and Wanner, Gerhard , journal=. Geometric numerical integration illustrated by the. 2003 , publisher=

  14. [14]

    Numerische Mathematik , volume=

    Symmetric multistep methods over long times , author=. Numerische Mathematik , volume=. 2004 , publisher=

  15. [15]

    The SMAI Journal of computational mathematics , volume=

    Symmetric multistep methods for charged-particle dynamics , author=. The SMAI Journal of computational mathematics , volume=

  16. [16]

    1999 , publisher=

    Shang, Zai-Jiu , journal=. 1999 , publisher=

  17. [17]

    Resonant and

    Shang, Zai-Jiu , journal=. Resonant and. 2000 , publisher=

  18. [18]

    2013 , publisher=

    Mathematical methods of classical mechanics , author=. 2013 , publisher=

  19. [19]

    Feng , title =

    K. Feng , title =. Proceedings of the International Conference on Scientific Computation , year =

  20. [20]

    Yoshida , title =

    H. Yoshida , title =. Celestial Mechanics and Dynamical Astronomy , year =

  21. [21]

    Modulated

    Cohen, David and Hairer, Ernst and Lubich, Christian , journal=. Modulated. 2003 , publisher=

  22. [22]

    Ernst Hairer and Gerhard Wanner , title =

  23. [23]

    SIAM journal on numerical analysis , volume=

    Long-time energy conservation of numerical methods for oscillatory differential equations , author=. SIAM journal on numerical analysis , volume=. 2000 , publisher=

  24. [24]

    V. I. Arnold , title =. Russian Mathematical Surveys , year =

  25. [25]

    Nachrichten der Akademie der Wissenschaften in G\"

    On invariant curves of area-preserving mapping of an annulus , author=. Nachrichten der Akademie der Wissenschaften in G\". 1962 , pages =

  26. [26]

    The imbedding problem for

    Nash, John , journal=. The imbedding problem for. 1956 , publisher=

  27. [27]

    Annali della Scuola Normale Superiore di Pisa

    Moser, Jürgen , title =. Annali della Scuola Normale Superiore di Pisa. Classe di Scienze , volume =

  28. [28]

    Communications in Mathematical Physics , volume =

    Chemin, Jean-Yves and Zhang, Ping , title =. Communications in Mathematical Physics , volume =

  29. [29]

    1987, Geophysical Fluid Dynamics (New York: Springer), doi: 10.1007/978-1-4612-4650-3

    Pedlosky, Joseph , title =. 1987 , edition =. doi:10.1007/978-1-4612-4650-3 , pages =

  30. [30]

    2013 , publisher=

    Matrix Computations , author=. 2013 , publisher=

  31. [31]

    Melnikov, V. K. , title =. Doklady Akademii Nauk SSSR , volume =

  32. [32]

    and Kozlov, Valery V

    Arnold, Vladimir I. and Kozlov, Valery V. and Neishtadt, Anatoly I. , title =. 2006 , volume =

  33. [33]

    Green's Function Estimates for Lattice

    Jean Bourgain , publisher =. Green's Function Estimates for Lattice

  34. [34]

    Invariant tori of

    Yuan, Xiaoping , journal=. Invariant tori of. 1998 , publisher=

  35. [35]

    Lagrange stability for

    Yuan, Xiaoping , journal=. Lagrange stability for. 2000 , publisher=

  36. [36]

    2007 , publisher=

    Approaches to the qualitative theory of ordinary differential equations: dynamical systems and nonlinear oscillations , author=. 2007 , publisher=

  37. [37]

    The applicability of the third integral of motion: some numerical experiments , journal =

    H. The applicability of the third integral of motion: some numerical experiments , journal =

  38. [38]

    Symplectic integration of

    Channell, Paul J and Scovel, Clint , journal=. Symplectic integration of. 1990 , publisher=

  39. [39]

    , TITLE =

    Kanamaru, T. , TITLE =. 2008 , JOURNAL =

  40. [40]

    and Bowen, M

    Witelski, T. and Bowen, M. , TITLE =. 2009 , JOURNAL =

  41. [41]

    1996 , publisher=

    Lectures on entire functions , author=. 1996 , publisher=

  42. [42]

    Analysis of

    Kuksin, Sergej B , volume=. Analysis of. 2000 , publisher=

  43. [43]

    Liu, Jianjun and Yuan, Xiaoping , journal=. A. 2011 , publisher=

  44. [44]

    2021 , publisher=

    Yuan, Xiaoping , journal=. 2021 , publisher=

  45. [46]

    2012 , publisher=

    Matrix analysis , author=. 2012 , publisher=

  46. [47]

    Studies of nonlinear problems

    Fermi, Enrico and Pasta, John and Ulam, Stanislaw , booktitle =. Studies of nonlinear problems. 1965 , note =

  47. [48]

    2006 , publisher=

    Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations , author=. 2006 , publisher=

  48. [49]

    On certain cases of conservation of almost periodic motions with a small change of the

    Melnikov, Viktor Kuz'mich , booktitle=. On certain cases of conservation of almost periodic motions with a small change of the. 1965 , organization=

  49. [50]

    Perturbations of stable invariant tori for

    Eliasson, L H. Perturbations of stable invariant tori for. Annali della Scuola Normale Superiore di Pisa-Classe di Scienze , volume=

  50. [51]

    On elliptic lower dimensional tori in

    P. On elliptic lower dimensional tori in. Mathematische Zeitschrift , volume=. 1989 , publisher=

  51. [52]

    Nonlinearity , volume=

    KAM theory without action-angle variables , author=. Nonlinearity , volume=

  52. [53]

    2016 , publisher=

    The parameterization method for invariant manifolds: from rigorous results to effective computations , author=. 2016 , publisher=

  53. [54]

    Discrete and Continuous Dynamical Systems , volume=

    Computation of whiskered invariant tori and their associated manifolds: New fast algorithms , author=. Discrete and Continuous Dynamical Systems , volume=. 2012 , publisher=

  54. [55]

    A parametrization algorithm to compute lower dimensional elliptic tori in

    Caracciolo, Chiara and Figueras, Jordi-Llu. A parametrization algorithm to compute lower dimensional elliptic tori in. Nonlinearity , volume=. 2025 , publisher=

  55. [56]

    Kolmogorov, A. N. , journal=. On the preservation of conditionally periodic motions under a small change of the. 1954 , note=

  56. [57]

    Arnold, V. I. , journal=. Proof of a theorem of

  57. [58]

    Bourgain, Jean , journal=. On. 1997 , publisher=

  58. [59]

    SIAM Review , volume=

    On the theory of quasiperiodic motions , author=. SIAM Review , volume=. 1966 , publisher=

  59. [60]

    Mathematische Annalen , volume=

    Convergent series expansions for quasi-periodic motions , author=. Mathematische Annalen , volume=. 1967 , publisher=

  60. [61]

    Soviet Mathematics , volume=

    Instability of dynamical systems with several degrees of freedom , author=. Soviet Mathematics , volume=

  61. [62]

    A Lecture on the Classical

    P. A Lecture on the Classical. Proceedings of Symposia in Pure Mathematics , volume=. 2001 , publisher=

  62. [63]

    1995 , publisher =

    Collected Works of Feng Kang , author =. 1995 , publisher =

  63. [64]

    Quasi-periodic solutions of nonlinear Schr

    Yuan, Xiaoping , journal=. Quasi-periodic solutions of nonlinear Schr. 2003 , publisher=

  64. [65]

    Physica D: Nonlinear Phenomena , volume=

    Breathers in nonlinear lattices: Existence, linear stability and quantization , author=. Physica D: Nonlinear Phenomena , volume=. 1997 , publisher=

  65. [66]

    Communications in Mathematical Physics , volume =

    Yuan, Xiaoping , title =. Communications in Mathematical Physics , volume =. 2002 , publisher =

  66. [67]

    KdV & KAM , series =

    Thomas Kappeler and J\". KdV & KAM , series =. 2003 , edition =

  67. [68]

    V. I. Arnold. Proof of a theorem of A . N . K olmogorov on the invariance of quasi-periodic motions under small perturbations of the H amiltonian. Russian Mathematical Surveys, 18 0 (5): 0 9--36, 1963

  68. [69]

    V. I. Arnold. Instability of dynamical systems with several degrees of freedom. Soviet Mathematics, 5: 0 581--585, 1964

  69. [70]

    V. I. Arnold. Mathematical methods of classical mechanics, volume 60. Springer Science & Business Media, 2013

  70. [71]

    S. Aubry. Breathers in nonlinear lattices: Existence, linear stability and quantization. Physica D: Nonlinear Phenomena, 103 0 (1-4): 0 201--250, 1997

  71. [72]

    Bourgain

    J. Bourgain. On M elnikov's persistence problem. Mathematical Research Letters, 4 0 (4): 0 445--458, 1997

  72. [73]

    Bourgain

    J. Bourgain. Green's Function Estimates for Lattice S chr \"o dinger Operators and Applications . Princeton University Press, Princeton, 2005

  73. [74]

    L. H. Eliasson. Perturbations of stable invariant tori for H amiltonian systems. Annali della Scuola Normale Superiore di Pisa-Classe di Scienze, 15 0 (1): 0 115--147, 1988

  74. [75]

    K. Feng. Collected Works of Feng Kang, volume II. Science Press, Beijing, 1995. Edited by Institute of Computational Mathematics and Scientific/Engineering Computing, Academy of Mathematics and Systems Science, Chinese Academy of Sciences; Proofread by Liqun, Cao and Yifa, Tang

  75. [76]

    Fermi, J

    E. Fermi, J. Pasta, and S. Ulam. Studies of nonlinear problems I . In Collected Papers of Enrico Fermi. University of Chicago Press, 1965. Originally published as Los Alamos report LA-1940, 1955

  76. [77]

    Fu and B

    M. Fu and B. Shi. Numerical construction of quasi-periodic solutions beyond symplectic integrators. arXiv preprint arXiv:2602.16275, 2026

  77. [78]

    Gauckler, E

    L. Gauckler, E. Hairer, and C. Lubich. Dynamics, numerical analysis, and some geometry. In Proceedings of the International Congress of Mathematicians: Rio de Janeiro 2018, pages 453--485. World Scientific, 2018

  78. [79]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner. Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, volume 31 of Springer Series in Computational Mathematics. Springer Berlin Heidelberg, second edition, 2006

  79. [80]

    H \'e non and C

    M. H \'e non and C. Heiles. The applicability of the third integral of motion: some numerical experiments. Astronomical Journal, 69: 0 73--79, 1964

  80. [81]

    Kappeler and J

    T. Kappeler and J. P\" o schel. KdV & KAM. Ergebnisse der Mathematik und ihrer Grenzgebiete. 3. Folge / A Series of Modern Surveys in Mathematics. Springer Berlin, Heidelberg, 1 edition, 2003

Showing first 80 references.