REVIEW 4 major objections 4 minor 35 references
Multisymplectic structure of nonintegrable Henon-Heiles system
T0 review · 4 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read Henon-Heiles, Kepler, and open Toda admit second invariant symplectic forms.
desk verdict Explicit second symplectic forms for nonintegrable systems are new and checkable, but the paper's central Jacobi verification is omitted and the integrator promise is speculative. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is an invariant bivector field $\tilde{P}$ of type $(2,0)$, found by substituting the polynomial ansatz (1.7) into the invariance equation $L_X P'=0$. The ansatz assumes each bivector entry is quadratic in the momenta with coefficients depending on coordinates, which turns the problem into 60 partial differential equations for 36 unknown functions. For the Henon-Heiles system the solution is $P'=(a_1H+a_2)P+a_3\tilde{P}$, and the supplemental bivector $\tilde{P}$ satisfies $\tilde{P}\,dH=2HX$; after multiplication by $H^{2/3}$ it obeys the Jacobi identity and is compatible with $H^{10/3}P$, so inverting it gives the second symplectic form $\hat{\omega}=H^{-8/3}\tilde{\omega}$. The same mechanism, with different supplemental bivectors, produces the Kepler and Toda results.
What would settle it
Substitute the Henon-Heiles vector field into the computed expressions and check, in a computer algebra system, that $L_X\tilde{P}=0$ and $d(H^{-8/3}\tilde{\omega})=0$ for generic nonzero $a,b$; a single choice of $a,b$ where either identity fails would destroy the claimed second symplectic form.
Extended reading notes
Core claim
The central claim is that the nonintegrable Henon-Heiles system with potential $V=q_1(aq_2^2+bq_1^2)$ possesses a second invariant symplectic form $\hat{\omega}=H^{-8/3}\tilde{\omega}$, obtained from an invariant bivector $\tilde{P}$ that is not proportional to the canonical Poisson bivector $P$. The invariance equation $L_X P'=0$ has the three-parameter solution $P'=(a_1H+a_2)P+a_3\tilde{P}$, and multiplying $\tilde{P}$ by $H^{2/3}$ produces a Poisson bivector satisfying the Jacobi identity, so its inverse is a genuine second symplectic form. The same ansatz yields a nine-parameter family of invariant bivectors for the Kepler problem, giving four Poisson bivectors and three invariant symplectic forms, and a three-parameter family for the open Toda lattice. The periodic Toda lattice, by contrast, admits only the two-parameter canonical family, so the existence of a second symplectic form is not a general consequence of integrability.
Load-bearing premise
The search is restricted to bivectors whose entries are polynomials of degree two in the momenta, so any second invariant symplectic form that is not of that polynomial form would be missed.
Editorial extensions
If this is right
- A multisymplectic integrator for the Henon-Heiles system can in principle preserve both $\omega$ and $\hat{\omega}$, so the resulting discrete flow would keep two independent symplectic structures instead of one.
- For the Kepler problem, the nine-parameter family of invariant bivectors yields four Poisson bivectors and three invariant symplectic forms, giving a richer menu of structures a discretization might preserve.
- The open Toda lattice has a three-parameter family of invariant bivectors, while the periodic Toda lattice has only the canonical two-parameter family, so the existence of a second form is not a general integrability feature.
- The trace of the invariant tensor $N=P'P^{-1}$ reproduces the integrals of motion of the Kepler problem, including the Runge-Lenz vector, which connects the second form to first integrals that the standard symplectic integrator does not preserve.
- Because $\tilde{P}\,dH=2HX$, the second structure is tied to the same flow under a time rescaling, so integrators preserving the second form should also behave well under reparametrized time.
Reading between the lines
- Extension: The failure of the periodic Toda lattice suggests that the presence of two symplectic forms may be tied to open boundary conditions or noncompactness; one could test whether other open versus periodic Toda pairs behave the same way.
- Extension: A practical numerical study comparing a standard symplectic integrator with a hypothetical two-form-preserving integrator on the chaotic Henon-Heiles system would test whether preserving the extra form actually improves long-time statistics; the paper stops before constructing such a scheme.
- Extension: The Kepler family contains complex-valued bivectors in polar coordinates, so a real multisymplectic integrator would need a reality-preserving combination; the paper does not address implementation.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper solves the invariance equation L_X T = 0 for tensor fields of the Hamiltonian vector field with two degrees of freedom. For the Hénon–Heiles potential V = q1(aq2^2 + bq1^2), it presents a three-parameter solution P' = (a1H + a2)P + a3P̃ within the polynomial bivector ansatz (1.7). It then claims that H^{2/3}P̃ is a Poisson bivector and is compatible with H^{10/3}P, so that its inverse gives a second invariant symplectic form. Similar invariant bivectors are given for a family of weight-homogeneous potentials, for the Kepler problem, and for open and periodic G2 Toda lattices; the paper suggests that these structures can be used to construct multisymplectic integrators preserving both symplectic forms.
Significance. If the Jacobi-identity assertions are correct, the paper provides explicit nontrivial invariant Poisson structures for benchmark nonintegrable systems, which is a useful contribution to geometric integration. The construction is not circular: the bivectors are obtained as solutions of an explicit PDE system with a stated ansatz, not fitted to a target result. However, the main weight of the paper rests on algebraic identities that are asserted without derivation, and at least one of the claimed symplectic structures appears to be complex-valued. The result is therefore promising but not yet fully verified.
major comments (4)
- [Section 2, Eqs. (2.9)–(2.10)] The central claim that P̂ = H^{2/3}P̃ is a Poisson bivector is asserted with the words 'verified analytically', but no computation is shown and no auxiliary file or code is provided. Invariance L_XP̃ = 0 and relation (2.11) do not imply the Jacobi identity [[P̂,P̂]] = 0; the Schouten bracket is a separate algebraic condition involving derivatives of both H and P̃. Since the inverse of P̂ is the advertised second symplectic form, this is load-bearing. Please include the full bracket computation, or a verifiable script, for both [[P̂,P̂]] = 0 and the claimed compatibility [[P̂,H^{10/3}P]] = 0, and justify the associated statement Ω = ω² = 4H^{10/3}ω̂². Note also that H^{10/3}P is not automatically a Poisson bivector for nonconstant H, so the compatibility statement needs a precise definition and proof.
- [Section 4, Eq. (4.19) and following paragraph] The 'Poisson bivector' P'_1 contains explicit complex objects, namely the factor i and e^{iϕ}. For a real Hamiltonian system, a complex bivector does not define a real invariant symplectic form unless its real and imaginary parts are separately shown to be Poisson. The paper does not do this, yet it concludes that the Kepler problem has three invariant symplectic forms. Please clarify whether the construction is over the complexification and, if so, state how a real second symplectic form is obtained.
- [Section 5, Proposition 5] For the open Toda lattice the paper proves only that P̃ is invariant under L_X. It does not show that P̃, or any scalar multiple of it, satisfies the Jacobi identity, nor that it is invertible, nor that its inverse is closed. The abstract and conclusion nevertheless claim a second invariant symplectic form for Toda type systems. This gap is load-bearing for that claim; please supply the missing Jacobi and invertibility checks, or restrict the conclusion to invariant bivectors.
- [Section 3, Proposition 2] The proof refers to 'three differential equations, which are omitted for the sake of brevity.' Since the proposition is used to assert preservation of P̃ for the whole weight-homogeneous family, omitting these equations leaves the proof incomplete. They should be included, at least in an appendix or supplementary file.
minor comments (4)
- [Eq. (5.21)] The entry P′^{12}_t = αp1 + βp1 should presumably be αp1 + βp2; as written β is redundant and inconsistent with the potentials in (5.22).
- [Eq. (5.23)] The entries for P̃^{13} and P̃^{24} contain unbalanced parentheses; the display should be corrected.
- [Abstract and Conclusion] The phrase 'nonintegrable Toda type systems' is not clearly identified: the Toda lattices in Section 5.1 are both integrable. Please state which family in Proposition 4 is meant by 'nonintegrable Toda type'.
- [Sections 2 and 3] The restriction of the search to the polynomial ansatz (1.7) is stated in the text but should be repeated in the abstract and conclusion; otherwise 'generic solution' overstates the classification.
Circularity Check
No significant circularity: the invariant bivectors are obtained by direct solution of the invariance equation, and the claimed second symplectic forms are not fitted to any target result.
full rationale
The derivation chain is self-contained. Section 2 introduces the explicit polynomial ansatz (1.7), solves L_X P' = 0, and exhibits the bivector P-tilde in (2.9); the candidate second symplectic form is then obtained by the stated scalar multiplication H^{2/3} and the asserted Jacobi identity (2.10). No parameter is fitted to the quantity being 'predicted', and no target result is fed back into the construction. Propositions 2-5 are established by direct solution of the invariance equation with stated integrability and parameter conditions, so the citations to the author's earlier work [21-23] supply method or terminology, not the load-bearing theorem. The polynomial restriction is explicitly acknowledged in Section 2 as a possible reason for the special weight-homogeneous cases, which is a scope limitation rather than a circular reduction. The unproved statement that H^{2/3}P-tilde satisfies the Jacobi identity is a verification gap and a correctness risk if false, but it is not an input-output equivalence: the Jacobi identity is an independent algebraic condition on a bivector constructed solely from invariance. The Kepler and Toda results likewise follow from solving the invariance and Jacobi equations, not from renaming a known result. Therefore no circular step can be exhibited.
Assumptions & free parameters
assumptions (4)
- standard math Lie derivative calculus and Hamiltonian mechanics, especially L_X ω = 0 and L_X P = 0 for Hamiltonian vector fields.
- domain assumption Invariant bivectors are assumed to be polynomials of degree at most two in momenta, as in ansatz (1.7).
- ad hoc to paper The 'straightforward calculation' proofs are accepted as correct without a shown verification.
- standard math For the Kepler problem, action-angle variables and the Bogoyavlenskij construction are taken as given.
Cite this review
Pith. "Pith review of Multisymplectic structure of nonintegrable Henon-Heiles system." pith.science (2026). https://pith.science/paper/FCWUY34F
@misc{pith2026250203786,
author = {Pith},
title = {Pith review of: Multisymplectic structure of nonintegrable Henon-Heiles system},
year = {2026},
howpublished = {\url{https://pith.science/paper/FCWUY34F}},
note = {Machine review of arXiv:2502.03786}
}
read the original abstract
Multi-symplectic integrators are typically regarded as a discretization of the Hamiltonian partial differential equations. This is due to the fact that, for generic finite-dimensional Hamiltonian systems, there exists only one independent symplectic structure. In this note, the second invariant symplectic form is presented for the nonintegrable Henon-Heiles system, Kepler problem, integrable and non-integrable Toda type systems. This approach facilitates the construction of a multi-symplectic integrator, which effectively preserves both symplectic forms for these benchmark problems.
Reference graph
Works this paper leans on
-
[1]
E., West M., Discrete mechanics and variational integr ators, Acta Numer., 10 (2001), 357-514
Marsden J. E., West M., Discrete mechanics and variational integr ators, Acta Numer., 10 (2001), 357-514
work page 2001
-
[2]
Cambridge University Press, Cam bridge (2005)
Leimkuhler B., Reich S., Simulating Hamiltonian Dynamics, Cambridge Mo nographs on Applied and Computational Mathematics. Cambridge University Press, Cam bridge (2005)
work page 2005
-
[3]
Hairer E., Lubic, C., Wanner, G., Geometric Numerical Integration : Structure-Preserving Algo- rithms for Ordinary Differential Equations, Springer, Berlin (2006)
work page 2006
-
[4]
Betsch P. (ed.): Structure-Preserving Integrators in Nonline ar Structural Dynamics and Flexible Multibody Dynamics, CISM Courses and Lectures, vol. 565. Springe r, Berlin (2016)
work page 2016
-
[5]
Bridges T., Reich S., Multi-symplectic integrators: Numerical sche mes for Hamiltonian PDEs that conserve symplecticity, Phys. Lett. A, 284:4-5 (2001), 184-193
work page 2001
-
[6]
Demoures F., Gay-Balmaz F., Ratiu T.S., Multisymplectic variational in tegrators and space/time symplecticity, Analysis and Applications, 14:3 (2016), 341-391
work page 2016
-
[7]
Hairer E., Lubich C., Wanner G., Geometric numerical integration illus trated by the St¨ ormer- Verlet method, Acta Numerica, 12 (2003), 399-450
work page 2003
-
[8]
Minesaki Y., Nakamura Y., New numerical integrator for the St¨ a ckel system conserving the same number of constants of motion as the degree of freedom, J. Phys . A: Math. Gen. 39 (2006), 9453-9476
work page 2006
Show all 35 references
-
[9]
Elenin G.G., Elenina T.G., Testing of Adaptive Symplectic Conservative Numerical Methods for Solving the Kepler Problem, Comput. Math. and Math. Phys., 58 (201 8), 863-880. 13
-
[10]
Bogoyavlenskij O.I., Theory of tensor invariants of integrable H amiltonian systems. I. Incompatible Poisson structures, Commun. Math. Phys., 180 (1996), 529-586
1996
-
[11]
Bogoyavlenskij O.I., Theory of tensor invariants of integrable H amiltonian systems. II. Theorem on symmetries and its applications, Commun. Math. Phys., 184 (1991 ), 301-365
1991
-
[12]
Bogoyavlenskij O.I., Extended integrability and bi-Hamiltonian sys tems, Comm. Math. Phys., 196 (1998), 19-51
1998
-
[13]
Kozlov V.V., Tensor invariants and integration of differential equ ations, Russ. Math. Surv., 74:1 (2019), 111-140
2019
-
[14]
Jiang K., Ratiu T.S., Zung N.T., Simultaneous local normal forms of d ynamical systems with singular underlying geometric structures, Nonlinearity, 37:10, (20 24), 105013
-
[15]
and Heiles, C., The applicability of the third integral of motion: some numerical experiments, Astronomical Journal, 69 (1964), 73-79
H´ enon, M. and Heiles, C., The applicability of the third integral of motion: some numerical experiments, Astronomical Journal, 69 (1964), 73-79
1964
-
[16]
Escobar-Ruiz A.M., Jim´ enez-Lara L., Ju´ arez-Flores P.M., Montoya-Molina F., Moreno-S´ aenz J., Quiroz-Ju´ arez M.A., Data-driven reconstruction of chaotic dyna mical equations: The H´ enon- Heiles type system, Chaos, Solitons and Fractals, 184 (2024), 115 025
2024
-
[17]
E., Foundations of Mechanics, second e dition, revised and enlarged, with the assistance of Tudor Ratiu and Richard Cushman
Abraham R., Marsden J. E., Foundations of Mechanics, second e dition, revised and enlarged, with the assistance of Tudor Ratiu and Richard Cushman. Benjamin/Cum mings Publishing Co., Inc., Advanced Book Program, Reading, Mass., 1978
1978
-
[18]
On the geometry of multisymple ctic manifolds, J
Cantrijn F., Ibort A., de Le´ on M. On the geometry of multisymple ctic manifolds, J. Aust. Math. Soc. A, 66 (1999), 303-330
1999
-
[19]
Barron T., Shafiee M., Multisymplectic structures induced by sym plectic structures, J. of Geom. Phys., 136 (2019), 1-13
2019
-
[20]
New gener alized Poisson structures, J
de Azc´ arraga J.A., Perelomov A.M., P´ erez Bueno J.C. New gener alized Poisson structures, J. Phys. A, 29:7 (1996), pp. L151-L157
1996
-
[21]
Non- linear Math
Tsiganov A.V., On bi-Integrable Natural Hamiltonian Systems on R iemannian Manifolds, J. Non- linear Math. Phys., 18:2 (2011), 245-268
2011
-
[22]
P hys.A: Math
Tsiganov A.V., On natural Poisson bivectors on the sphere, J. P hys.A: Math. Theor., 44:10 (2011), 105203
2011
-
[23]
Maciejewski A.J., Przybylska M., Tsiganov A.V., On algebraic const ruction of certain integrable and super-integrable systems, Physica D: Nonlinear Phenomena, 2 40:18 (2011), 1426-1448
2011
-
[24]
Kolmogorov A.N., The general theory of dynamical systems and classical mechanics, In: Proc. Intern. Congr. Math., 1, (1954), Amsterdam: North-Holland Pub l. Co., (1957), 315-333
1954
-
[25]
Tsiganov A.V., On maximally superintegrable systems, Regul. Chao t. Dyn., 13:3 (2008), 178-190
2008
-
[26]
Euler L., Probleme un corps ´ etant attir´ e en raison r´ eciproqu e quarr´ ee des distances vers deux points fixes donn´ es, trouver les cas o´ u la courbe d´ ecrite par cecorps sera alg´ ebrique, M´ emoires de l’academie des sciences de Berlin, 16 (1767), 228-249, 1760-176 7
-
[27]
Tsiganov A.V., Superintegrable systems and Riemann-Roch theo rem, J. Math. Phys., 61 (2020), 012701
2020
-
[28]
Jacobi C. G. J., Vorlesungen ¨ uber Dynamik, Gehalten an der Un iversit¨ at zu K¨ onigsberg im Win- tersemester 1842-1843, ed. Alfred Clebsch, Georg Reimer, Berlin , 1866, Revised ed. Gesammelte Werke suppl., Reimer, Berlin, 1884. 14
-
[29]
St¨ ackel P., Ueber die Integration der Hamilton-Jacobi’schen D ifferentialgleichung mittels Sepa- ration der Variabeln, Habilitationsschrift, Halle, 1891
-
[30]
Morbidelli A., Modern Celestial Mechanics: Dynamics in the Solar Sy stem, Taylor Francis, (2002)
2002
-
[31]
Fernandes R. L. , Completely integrable bi-Hamiltonian systems, J. Dyn. Diff. Eq., 6 (1994), 53-69
1994
-
[32]
A., Tsiganov A
Grigoryev Yu. A., Tsiganov A. V., On bi-Hamiltonian formulation of t he perturbed Kepler prob- lem, J. Phys. A: Math. Theor., 48 (2015), 175206
2015
-
[33]
Marmo, G
G. Marmo, G. Vilasi, When do recursion operators generate new conservation laws?, Phys. Lett. B, 277 ( 1992), 137-140
1992
-
[34]
McDaniel A., Smolinsky L., The flow of the G2 periodic Toda lattice, J. Math. Phys. v.38:2 (1997), 926-945
1997
-
[35]
Tsiganov A.V., On tensor invariants for the Euler, Lagrange and Kowalewski cases of rigid body motion, arXiv:2410.10109, (2024). 15
2024 arXiv
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.