Mathematics – Numerical Analysis
Scientific paper
2010-05-12
Mathematics
Numerical Analysis
27 pages, 3 figures, submitted to IMA Journal of Numerical Analysis
Scientific paper
Prior to the recent development of symplectic integrators, the time-stepping operator $\e^{h(A+B)}$ was routinely decomposed into a sum of products of $\e^{h A}$ and $\e^{hB}$ in the study of hyperbolic partial differential equations. In the context of solving Hamiltonian dynamics, we show that such a decomposition give rises to both {\it even} and {\it odd} order Runge-Kutta and Nystr\"om integrators. By use of Suzuki's forward-time derivative operator to enforce the time-ordered exponential, we show that the same decomposition can be used to solve non-autonomous equations. In particular, odd order algorithms are derived on the basis of a highly non-trivial {\it time-asymmetric} kernel. Such an operator approach provides a general and unified basis for understanding structure non-preserving algorithms and is especially useful in deriving very high-order algorithms via {\it analytical} extrapolations. In this work, algorithms up to the 100th order are tested by integrating the ground state wave function of the hydrogen atom. For such a singular Coulomb problem, the multi-product expansion showed uniform convergence and is free of poles usually associated with structure-preserving methods. Other examples are also discussed.
Chin Siu A.
Geiser Jurgen
No associations
LandOfFree
Multi-product operator splitting as a general method of solving autonomous and non-autonomous equations does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Multi-product operator splitting as a general method of solving autonomous and non-autonomous equations, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Multi-product operator splitting as a general method of solving autonomous and non-autonomous equations will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-639871