NySALT: Nyström-type inference-based schemes adaptive to large time-stepping

Abstract

Large time-stepping is important for efficient long-time simulations of deterministic and stochastic Hamiltonian dynamical systems. Conventional structure-preserving integrators, while being successful for generic systems, have limited tolerance to time step size due to stability and accuracy constraints. We propose to use data to innovate classical integrators so that they can be adaptive to large time-stepping and are tailored to each specific system. In particular, we introduce NySALT, Nyström-type inference-based schemes adaptive to large time-stepping. The NySALT has optimal parameters for each time step learnt from data by minimizing the one-step prediction error. Thus, it is tailored for each time step size and the specific system to achieve optimal performance and tolerate large time-stepping in an adaptive fashion. We prove and numerically verify the convergence of the estimators as data size increases. Furthermore, analysis and numerical tests on the deterministic and stochastic Fermi-Pasta-Ulam (FPU) models show that NySALT enlarges the maximal admissible step size of linear stability, and quadruples the time step size of the Störmer--Verlet and the BAOAB when maintaining similar levels of accuracy.

arXiv comment26 pages, 7 figures

References in corpus (22)

• Learning data driven discretizations for partial differential equations
• A momentum conserving model with anomalous thermal conductivity in low dimension
• Data-driven parameterization of the generalized Langevin equation
• Data-driven non-Markovian closure models
• Symplectic ODE-Net: Learning Hamiltonian Dynamics with Control
• On Learning Hamiltonian Systems from Data
• Convergence of Numerical Time-Averaging and Stationary Measures via Poisson Equations
• Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance
• Effective dynamics using conditional expectations
• Structure Preserving Model Reduction of Parametric Hamiltonian Systems
• A discrete approach to stochastic parametrization and dimensional reduction in nonlinear dynamics
• Long-Run Accuracy of Variational Integrators in the Stochastic Context
• Data-driven model reduction, Wiener projections, and the Koopman-Mori-Zwanzig formalism
• Non-intrusive and structure preserving multiscale integration of stiff ODEs, SDEs and Hamiltonian systems with hidden slow dynamics via flow averaging
• Explicit high-order symplectic integrators for charged particles in general electromagnetic fields
• Crossover from Fermi-Pasta-Ulam to normal diffusive behaviour in heat conduction through open anharmonic lattices
• GRIT: a package for structure-preserving simulations of gravitationally interacting rigid-bodies
• Data-driven model reduction for stochastic Burgers equations
• Stochastic Resonance in the Fermi-Pasta-Ulam Chain
• Accurate and efficient splitting methods for dissipative particle dynamics
• Accurate and Efficient Simulations of Hamiltonian Mechanical Systems with Discontinuous Potentials
• Nonlinearity accelerates the thermalization of the quartic FPUT model with stochastic baths