Symplectic Error of Implicit Symplectic Integrators: A Qualitative Structural Analysis
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911611672854528 |
|---|---|
| author | Gajdoš, Matěj Brichta, Ondřej Kučera, Václav |
| author_facet | Gajdoš, Matěj Brichta, Ondřej Kučera, Václav |
| contents | We study how inexact nonlinear solvers lead to a loss of exact symplecticity in the Symplectic Euler (SE) and Stormer-Verlet (SV) schemes when applied to general nonseparable Hamiltonian systems. These schemes are implicit and require nonlinear solvers in practice. Here, we consider a fixed number $M$ of fixed-point iterations (FPI). While SE is exactly symplectic under exact solves, a finite $M$ gives only pseudo-symplecticity. Compared to previous results, we provide a more qualitative, block-wise characterization of the induced pseudo-symplecticity by analyzing the resulting perturbations to the matrix of symplectic structure $J$. We prove that the perturbed matrix $\tilde{J}$ is skew-symmetric, that one diagonal block vanishes identically (depending on the SE variant), and that the remaining blocks are $O(h^{M+1})$ perturbations of their counterparts in $J$, with time step $h$. A quadratic Hamiltonian example shows these bounds are sharp. Extending to compositions, we quantify how SV inherits distinct decay orders across different blocks of the symplectic defect. As a corollary, we show that the perturbation of volume preservation in phase space arises solely from the off-diagonal blocks of $\tilde{J}$, and we bound the induced energy error along trajectories. Numerical experiments on a tokamak magnetic-field Hamiltonian, where q-implicit SE is fully nonlinear (requiring FPI) but p-implicit SE is linearly implicit, confirm the sharpness of the theory and highlight the gap to the exactly symplectic counterpart. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2604_19272 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Symplectic Error of Implicit Symplectic Integrators: A Qualitative Structural Analysis Gajdoš, Matěj Brichta, Ondřej Kučera, Václav Numerical Analysis 65P10 We study how inexact nonlinear solvers lead to a loss of exact symplecticity in the Symplectic Euler (SE) and Stormer-Verlet (SV) schemes when applied to general nonseparable Hamiltonian systems. These schemes are implicit and require nonlinear solvers in practice. Here, we consider a fixed number $M$ of fixed-point iterations (FPI). While SE is exactly symplectic under exact solves, a finite $M$ gives only pseudo-symplecticity. Compared to previous results, we provide a more qualitative, block-wise characterization of the induced pseudo-symplecticity by analyzing the resulting perturbations to the matrix of symplectic structure $J$. We prove that the perturbed matrix $\tilde{J}$ is skew-symmetric, that one diagonal block vanishes identically (depending on the SE variant), and that the remaining blocks are $O(h^{M+1})$ perturbations of their counterparts in $J$, with time step $h$. A quadratic Hamiltonian example shows these bounds are sharp. Extending to compositions, we quantify how SV inherits distinct decay orders across different blocks of the symplectic defect. As a corollary, we show that the perturbation of volume preservation in phase space arises solely from the off-diagonal blocks of $\tilde{J}$, and we bound the induced energy error along trajectories. Numerical experiments on a tokamak magnetic-field Hamiltonian, where q-implicit SE is fully nonlinear (requiring FPI) but p-implicit SE is linearly implicit, confirm the sharpness of the theory and highlight the gap to the exactly symplectic counterpart. |
| title | Symplectic Error of Implicit Symplectic Integrators: A Qualitative Structural Analysis |
| topic | Numerical Analysis 65P10 |
| url | https://arxiv.org/abs/2604.19272 |