Symplectic Runge-Kutta methods for Hamiltonian systems driven by Gaussian rough paths

We consider Hamiltonian systems driven by multi-dimensional Gaussian processes in rough path sense, which include fractional Brownian motions with Hurst parameter H∈(1/4,1/2]. We prove that the phase flow preserves the symplectic structure almost surely and this property could be inherited by symplectic Runge-Kutta methods, which are implicit methods in general. If the vector fields satisfy some smoothness and boundedness conditions, we obtain the pathwise convergence rates of Runge-Kutta methods. When vector fields are linear, we get the solvability of the midpoint scheme for skew symmetric cases, and obtain its pathwise convergence rate. Numerical experiments verify our theoretical analysis.