Dynamical chaos in the integrable Toda chain induced by time discretization

We use the Toda chain model to demonstrate that numerical simulation of integrable Hamiltonian dynamics using time discretization destroys integrability and induces dynamical chaos. Specifically, we integrate this model with various symplectic integrators parametrized by the time step $\tau$ and measure the Lyapunov time $T_{\Lambda}$ (inverse of the largest Lyapunov exponent $\Lambda$). A key observation is that $T_{\Lambda}$ is finite whenever $\tau$ is finite but diverges when $\tau \rightarrow 0$. We compare the Toda chain results with the nonitegrable Fermi-Pasta-Ulam-Tsingou chain dynamics. In addition, we observe a breakdown of the simulations at times $T_B \gg T_{\Lambda}$ due to certain positions and momenta becoming extremely large (``Not a Number''). This phenomenon originates from the periodic driving introduced by symplectic integrators and we also identify the concrete mechanism of the breakdown in the case of the Toda chain.
Comment: 13 pages, 7 figures