pith. sign in

arxiv: 1907.10351 · v1 · pith:7FUIBKD7new · submitted 2019-07-24 · 🧮 math.NA · cs.NA

Energy-preserving multi-symplectic Runge-Kutta methods for Hamiltonian wave equations

Chuchu Chen , Jialin Hong , Chol Sim , Kwang Sonwu This is my paper

Pith reviewed 2026-05-24 16:56 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords multi-symplectic methodsRunge-Kutta methodsHamiltonian wave equationsenergy preservationparametric methodsgeometric numerical integrationHamiltonian PDEs
0
0 comments X

The pith

Parametric multi-symplectic Runge-Kutta methods conserve energy in a weaker sense for Hamiltonian wave equations when a suitable parameter is chosen.

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

The paper develops a class of parametric multi-symplectic Runge-Kutta methods for Hamiltonian wave equations. These methods preserve the multi-symplectic structure for any parameter while achieving energy conservation in a weaker sense once a specific parameter is selected. The authors prove the existence of this parameter under the conditions of fixed step sizes and a fixed initial condition. This addresses the known general incompatibility between simultaneous geometric and physical preservation in numerical methods for such equations. Experiments compare the new methods to classical ones and illustrate the energy behavior.

Core claim

A novel class of parametric multi-symplectic Runge-Kutta methods is presented for Hamiltonian wave equations. These methods conserve energy simultaneously in a weaker sense with a suitable parameter. The existence of the parameter enforcing the energy-preserving property is proved under assumptions on the fixed step sizes and the fixed initial condition.

What carries the argument

Parametric multi-symplectic Runge-Kutta method, where a free parameter is adjusted to satisfy an additional discrete energy conservation condition derived from the scheme while the multi-symplectic property holds independently of the parameter value.

If this is right

  • The methods remain multi-symplectic regardless of the parameter value chosen.
  • Energy is conserved in a weaker sense precisely when the parameter satisfies the derived algebraic condition.
  • Existence of the required parameter is guaranteed under the stated assumptions on step sizes and initial data.
  • Numerical tests demonstrate that the energy error remains smaller than for the non-parametric classical multi-symplectic Runge-Kutta method.

Where Pith is reading between the lines

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

  • The parametric construction could be tested on other Hamiltonian PDEs if analogous energy conditions can be written down.
  • Adaptive step-size versions would likely need a parameter recomputed at each step rather than fixed once.
  • The weaker energy property may still bound long-time drift in applications where strict conservation is impractical.
  • Similar parameter tuning might be applied to other structure-preserving discretizations such as finite-difference or collocation schemes.

Load-bearing premise

The existence of the energy-preserving parameter is proved only when step sizes remain fixed and the initial condition is fixed throughout the computation.

What would settle it

Finding fixed step sizes and a fixed initial condition for which no real value of the parameter satisfies the discrete energy equality would show the existence claim fails.

Figures

Figures reproduced from arXiv: 1907.10351 by Chol Sim, Chuchu Chen, Jialin Hong, Kwang Sonwu.

Figure 1
Figure 1. Figure 1: A uniform grid and the unknowns in grid points. [PITH_FULL_IMAGE:figures/full_fig_p019_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Values of α in α-Gauss collocation methods at (x, t)-plane [-50, 50]× [0, 200] with h = 1 and τ = 0.1. Next, considering the periodic boundary condition, in each time level we put to￾gether the above individual nonlinear systems through space axis (including M = 100 grid points), which leads to a nonlinear system with 2100 unknowns. Denoting by X the unknowns, this nonlinear system can be rewritten as F(X)… view at source ↗
Figure 3
Figure 3. Figure 3: Values of α in α-Gauss collocation methods at two fixed moments t=50 and t=150 [PITH_FULL_IMAGE:figures/full_fig_p021_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical results show that in the conservation of the total energy [PITH_FULL_IMAGE:figures/full_fig_p021_5.png] view at source ↗
Figure 4
Figure 4. Figure 4: Time evolution of solution during the time interval t [PITH_FULL_IMAGE:figures/full_fig_p022_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Comparison of numerical errors in the global energy for both multi-symplectic Gauss collocation [PITH_FULL_IMAGE:figures/full_fig_p022_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Comparison of numerical errors in the local energy for both multi-symplectic Gauss collocation [PITH_FULL_IMAGE:figures/full_fig_p023_6.png] view at source ↗
read the original abstract

It is well-known that a numerical method which is at the same time geometric structure-preserving and physical property-preserving cannot exist in general for Hamiltonian partial differential equations. In this paper, we present a novel class of parametric multi-symplectic Runge-Kutta methods for Hamiltonian wave equations, which can also conserve energy simultaneously in a weaker sense with a suitable parameter. The existence of such a parameter, which enforces the energy-preserving property, is proved under certain assumptions on the fixed step sizes and the fixed initial condition. We compare the proposed method with the classical multi-symplectic Runge-Kutta method in numerical experiments, which shows the remarkable energy-preserving property of the proposed method and illustrate the validity of theoretical results.

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 / 1 minor

Summary. The manuscript proposes a novel class of parametric multi-symplectic Runge-Kutta methods for Hamiltonian wave equations. These methods preserve multi-symplecticity while also conserving energy in a weaker sense through the choice of a suitable parameter. The existence of such a parameter is proved under assumptions of fixed step sizes and a fixed initial condition. Numerical experiments compare the proposed methods to classical multi-symplectic Runge-Kutta methods and illustrate the energy-preserving property.

Significance. If the existence result and numerical comparisons hold, the work offers a parametric route to simultaneous geometric and (weaker) physical preservation for Hamiltonian PDEs, where general simultaneous preservation is known to be impossible. The explicit existence proof under the stated assumptions and the reported numerical advantage constitute a concrete contribution to the literature on structure-preserving discretizations.

major comments (2)
  1. [Abstract] Abstract: the existence of the energy-preserving parameter is proved only under the assumptions of fixed step sizes and fixed initial condition. This assumption is load-bearing for the central claim, because the manuscript provides no argument or extension showing that the parameter continues to exist (or that the weaker energy property holds) when step sizes vary or the initial condition changes, both of which are standard in applications of Hamiltonian wave equations.
  2. [Abstract] Abstract: energy conservation is obtained by selecting the parameter whose value is determined to enforce the property. The manuscript should clarify whether this selection is performed by solving an auxiliary equation at each step and whether the resulting scheme remains multi-symplectic for arbitrary choices of the free parameter or only for the energy-enforcing value.
minor comments (1)
  1. The phrase 'in a weaker sense' is used without a precise definition or reference to the specific conserved quantity; adding this definition would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. Below we respond point by point to the major comments.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the existence of the energy-preserving parameter is proved only under the assumptions of fixed step sizes and fixed initial condition. This assumption is load-bearing for the central claim, because the manuscript provides no argument or extension showing that the parameter continues to exist (or that the weaker energy property holds) when step sizes vary or the initial condition changes, both of which are standard in applications of Hamiltonian wave equations.

    Authors: We agree that the existence result is proved only under the stated assumptions of fixed step sizes and fixed initial condition; this is already noted in the abstract and is essential to the proof technique. The manuscript contains no extension or argument for variable step sizes or changing initial conditions. In the revised version we will strengthen the wording in the abstract and add a short remark in the introduction explicitly acknowledging the scope of the result and that extensions to variable steps lie outside the present work. revision: partial

  2. Referee: [Abstract] Abstract: energy conservation is obtained by selecting the parameter whose value is determined to enforce the property. The manuscript should clarify whether this selection is performed by solving an auxiliary equation at each step and whether the resulting scheme remains multi-symplectic for arbitrary choices of the free parameter or only for the energy-enforcing value.

    Authors: The parameter is determined once, under the fixed-step and fixed-initial-condition hypotheses, so that the weaker energy identity holds for the whole trajectory; no auxiliary equation is solved at each step. Multi-symplecticity holds for every value of the free parameter because it is built into the parametric Runge-Kutta tableau; the energy-enforcing choice is merely one specific member of this family. We will insert a clarifying sentence in the revised manuscript stating both facts. revision: yes

Circularity Check

1 steps flagged

Energy preservation obtained by selecting parameter to enforce the property under fixed steps/initial condition

specific steps
  1. fitted input called prediction [Abstract]
    "we present a novel class of parametric multi-symplectic Runge-Kutta methods for Hamiltonian wave equations, which can also conserve energy simultaneously in a weaker sense with a suitable parameter. The existence of such a parameter, which enforces the energy-preserving property, is proved under certain assumptions on the fixed step sizes and the fixed initial condition."

    The energy conservation is not shown to hold for the method independently; instead a parameter is introduced and its value is shown to exist specifically to enforce the property (under restrictive fixed-step/fixed-IC assumptions). The preservation therefore reduces to the definitional choice of that parameter.

full rationale

The paper introduces parametric multi-symplectic RK methods and claims they conserve energy in a weaker sense with a suitable parameter whose existence is proved only under fixed step sizes and fixed initial condition. This reduces the claimed preservation to a quantity enforced by parameter selection rather than an independent derivation, matching the fitted-input-called-prediction pattern. No other circular steps are identifiable from the provided text.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence of a parameter chosen to enforce energy preservation while retaining multi-symplecticity, together with domain assumptions on step sizes and initial data.

free parameters (1)
  • energy-enforcing parameter
    The parameter whose value is selected to make the method energy-preserving; its existence is asserted but no explicit formula or fitting procedure is given in the abstract.
axioms (1)
  • domain assumption fixed step sizes and fixed initial condition allow existence of the parameter
    The existence proof is stated to hold only under these restrictions.

pith-pipeline@v0.9.0 · 5650 in / 1234 out tokens · 18505 ms · 2026-05-24T16:56:51.114883+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

20 extracted references · 20 canonical work pages

  1. [1]

    Z. Ge, J. E. Marsden, Lie-Poisson Hamilton-Jacobi theory and Lie-Poisson inte- grators, Phys. Lett. A 133 (3) (1988) 134–139

    work page 1988

  2. [2]

    C.Kane,J.E.Marsden,M.Ortiz,Symplectic-energy-momentum preserving varia- tional integrators, J. Math. Phys. 40 (7) (1999) 3353–3371

    work page 1999

  3. [3]

    J. C. Simo, N. Tarnow, K. K. Wong, Exact energy-momentum conserving algo- rithms and symplectic schemes for nonlinear dynamics, Comput. Methods Appl. Mech. Engrg. 100 (1) (1992) 63–116

    work page 1992

  4. [4]

    Brugnano, G

    L. Brugnano, G. Gurioli, F. Iavernaro, Analysis of energy and quadratic invariant preserving (EQUIP) methods, J. Comput. Appl. Math. 335 (2018) 51–73

    work page 2018

  5. [5]

    Brugnano, F

    L. Brugnano, F. Iavernaro, D. Trigiante, Energy- and quadratic invariants- preserving integrators based upon Gauss collocation formulae, SIAM J. Numer. Anal. 50 (6) (2012) 2897–2916

    work page 2012

  6. [6]

    D. Wang, A. Xiao, X. Li, Parametric symplectic partitioned Runge-Kutta methods with energy-preserving properties for Hamiltonian systems, Comput. Phys. Com- mun. 184 (2) (2013) 303–310

    work page 2013

  7. [7]

    Y . Gong, J. Cai, Y . Wang, Some new structure-preserving algorithms for general multi-symplectic formulations of Hamiltonian PDEs, J. Comput. Phys. 279 (2014) 80–102

    work page 2014

  8. [8]

    J. Hong, H. Liu, G. Sun, The multi-symplecticity of partitioned Runge-Kutta meth- ods for Hamiltonian PDEs, Math. Comp. 75 (253) (2006) 167–181

    work page 2006

  9. [9]

    J. Hong, Y . Sun, Generating functions of multi-symplectic RK methods via DW Hamilton-Jacobi equations, Numer. Math. 110 (4) (2008) 491–519

    work page 2008

  10. [10]

    Reich, Multi-symplectic Runge-Kutta collocation methods for Hamiltonian wave equations, J

    S. Reich, Multi-symplectic Runge-Kutta collocation methods for Hamiltonian wave equations, J. Comput. Phys. 157 (2) (2000) 473–499. 25

    work page 2000

  11. [11]

    Y . Wang, J. Hong, Multi-symplectic algorithms for Hamiltonian partial differen- tial equations, Commun. Appl. Math. Comput. 27 (2) (2013) 163–230

    work page 2013

  12. [12]

    J. Hong, C. Li, Multi-symplectic Runge-Kutta methods for nonlinear Dirac equa- tions, J. Comput. Phys. 211 (2) (2006) 448–472

    work page 2006

  13. [13]

    J. Hong, X. Liu, C. Li, Multi-symplectic Runge-Kutta-Nystr ¨om methods for non- linear Schr ¨odinger equations with variable coe fficients, J. Comput. Phys. 226 (2) (2007) 1968–1984

    work page 2007

  14. [14]

    J. Hong, S. Jiang, C. Li, Explicit multi-symplectic methods for Klein-Gordon- Schr¨odinger equations, J. Comput. Phys. 228 (9) (2009) 3517–3532

    work page 2009

  15. [15]

    J. Hong, S. Jiang, C. Li, Accuracy of classical conservation laws for Hamiltonian PDEs under Runge-Kutta discretizations, Numer. Math. 112 (1) (2009) 1–23

    work page 2009

  16. [16]

    J. Hong, S. Jiang, C. Li, H. Liu, Explicit multi-symplectic methods for Hamilto- nian wave equations, Commun. Comput. Phys. 2 (4) (2007) 662–683

    work page 2007

  17. [17]

    H. Liu, K. Zhang, Multi-symplectic Runge-Kutta-type methods for Hamiltonian wave equations, IMA J. Numer. Anal. 26 (2) (2006) 252–271

    work page 2006

  18. [18]

    Cohen, E

    D. Cohen, E. Hairer, C. Lubich, Conservation of energy, momentum and actions in numerical discretizations of non-linear wave equations, Numer. Math. 110 (2) (2008) 113–143

    work page 2008

  19. [19]

    Hairer, G

    E. Hairer, G. Wanner, Solving ordinary di fferential equations. II, V ol. 14 of Springer Series in Computational Mathematics, Springer-Verlag, Berlin, 2010, stiff and differential-algebraic problems, Second revised edition

    work page 2010

  20. [20]

    J. Hong, C. Sim, X. Yin, Solvability of concatenated Runge-Kutta equations for second-order nonlinear PDEs, J. Comput. Appl. Math. 245 (2013) 232–241. 26

    work page 2013