Symplectic methods based on Pad$\acute{e}$ approximation for some stochastic Hamiltonian systems

TL;DR

This paper develops symplectic numerical schemes based on Padé approximation for stochastic Hamiltonian systems, ensuring structure preservation and specific convergence orders for different system types.

Contribution

Introduces Padé approximation-based symplectic schemes for stochastic Hamiltonian systems, achieving desired order and structure preservation under various conditions.

Findings

Padé approximations $P_{(k,k)}$ preserve symplecticity for linear systems.

Numerical schemes achieve mean-square order $(r+s)/2$ for certain Padé approximations.

Special schemes for additive noise systems attain order $reve r + reve s + 1$ with symplecticity when $reve r=reve s$.

Abstract

In this article, we introduce a kind of numerical schemes, based on Pad$\acute{e}$ approximation, for two stochastic Hamiltonian systems which are treated separately. For the linear stochastic Hamiltonian systems, it is shown that the applied Pad$\acute e$ approximations $P_{(k,k)}$ give numerical solutions that inherit the symplecticity and the proposed numerical schemes based on $P_{(r,s)}$ are of mean-square order $\frac{r+s}{2}$ under appropriate conditions. In case of the special stochastic Hamiltonian systems with additive noises, the numerical method using two kinds of Pad$\acute e$ approximation $P_{(\hat r,\hat s)}$ and $P_{(\check r,\check s)}$ has mean-square order $\check r+\check s+1$ when $\hat r+\hat s=\check r+\check s+2$. Moreover, the numerical solution is symplectic if $\hat r=\hat s$.