Cheap and easy to implement fourth-order methods for the Schrodinger equation with time-dependent Hamiltonians are ïntroduced. The methods require evaluations of exponentials of simple unidimensional integrals, and can be considered an averaging technique, preserving many of the qualitative propert
Optimal stability polynomials for splitting methods, with application to the time-dependent Schrödinger equation
✍ Scribed by Robert I. McLachlan; Stephen K. Gray
- Publisher
- Elsevier Science
- Year
- 1997
- Tongue
- English
- Weight
- 643 KB
- Volume
- 25
- Category
- Article
- ISSN
- 0168-9274
No coin nor oath required. For personal study only.
✦ Synopsis
We determine optimal stability polynomials p(x) for splitting method solutions of differential equations, building on previous work by L6pez-Marcos, Sanz-Sema and Skeel (1996). The methods have a variety of stage numbers and are up to eighth order. Knowledge of p(x) allows construction of the most stable splitting methods for given complexity. As an illustration, we construct symplectic corrector algorithms (C it<C, where the kernel K is an m-stage splitting method) which approximate the solution of linear Hamiltonian systems. The kernels K that realize the optimal stability polynomials are found for this case. We also discuss the construction of correctors C, and find them for two particularly promising kernels. Numerical calculations for a time-dependent Schr0dinger equation problem confirm the methods' usefulness.
📜 SIMILAR VOLUMES
## Abstract A new approach based upon the Taylor series method is proposed for propagating solutions of the time‐dependent Schrödinger equation. Replacing the spatial derivative of the wave function with finite difference formulas, we derive a recursive formula for the evaluation of Taylor coeffici