The solution of stationary ODE problems in quantum mechanics by Magnus methods with stepsize control

In solid state physics the solution of the Dirac and Schrodinger equation by operator splitting methods leads to differential equations with oscillating solutions for t he radial direction. For standard time int egrators like Runge–Kutta or multistep methods the stepsize is restricted approximately by the length of the period. In contrast the recently developed Magnus methods allow stepsizes that are substantially larger than one period. They are based on a Lie group approach and incorporate exponential functions and matrix commutators. A stepsize control is implemented and tested. As numerical examples eigenvalue problems for the radial Schrodinger equation and the radial Dirac equation are solved. Further, phase shifts for scattering solutions for hydrogen atoms and copper are computed.  2004 Elsevier B.V. All rights reserved.

[1]  Rajendra Prasad,et al.  Electronic Structure of Disordered Alloys , 1998 .

[2]  T. E. Simos,et al.  Eighth order methods for accurate computations for the Schrödinger equation , 1997 .

[3]  High-algebraic, high-phase lag methods for accurate computations for the elastic- scattering phase shift problem , 1998 .

[4]  A. Iserles,et al.  On the Implementation of the Method of Magnus Series for Linear Differential Equations , 1999 .

[5]  Klarsfeld,et al.  Recursive generation of higher-order terms in the Magnus expansion. , 1989, Physical review. A, General physics.

[6]  Patrick Fischer,et al.  Numerical Solution of the Schrödinger Equation in a Wavelet Basis for Hydrogen-like Atoms , 1998 .

[7]  Ben P. Sommeijer,et al.  Explicit Runge-Kutta (-Nyström) methods with reduced phase errors for computing oscillating solutions , 1987 .

[9]  A. Iserles,et al.  Lie-group methods , 2000, Acta Numerica.

[10]  Fernando Casas,et al.  Improved High Order Integrators Based on the Magnus Expansion , 2000 .

[11]  P. J. Prince,et al.  New Runge-Kutta algorithms for numerical simulation in dynamical astronomy , 1978 .

[12]  A. C. Allison,et al.  The numerical solution of coupled differential equations arising from the Schrödinger equation , 1970 .

[13]  J. Butcher.  The numerical analysis of ordinary differential equations: Runge-Kutta and general linear methods , 1987 .

[14]  A. Ernst,et al.  Theoretical study of magnetic layers of nickel on copper; dead or alive? , 2000 .

[15]  B. R. Johnson.  The renormalized Numerov method applied to calculating bound states of the coupled‐channel Schroedinger equation , 1978 .

[16]  W. Magnus.  On the exponential solution of differential equations for a linear operator , 1954 .

[17]  Ilja Turek,et al.  Electronic Structure of Disordered Alloys, Surfaces and Interfaces , 1996 .

[18]  A. Iserles,et al.  On the solution of linear differential equations in Lie groups , 1999, Philosophical Transactions of the Royal Society of London. Series A: Mathematical, Physical and Engineering Sciences.

[19]  E. Hairer,et al.  Solving Ordinary Differential Equations II: Stiff and Differential-Algebraic Problems , 2010 .

[20]  J. Ros,et al.  High Order Optimized Geometric Integrators for Linear Differential Equations , 2002 .

[21]  A. Messiah.  Quantum Mechanics , 1961 .

[22]  Ernst Hairer,et al.  Solving Ordinary Differential Equations I: Nonstiff Problems , 2009 .

[23]  M. S. Kaschiev,et al.  Finite-element solution of the coupled-channel Schrödinger equation using high-order accuracy approximations , 1995 .

[24]  S. N. Papakostas,et al.  High Phase-Lag-Order Runge-Kutta and Nyström Pairs , 1999, SIAM J. Sci. Comput..

[25]  Marlis Hochbruck,et al.  On Magnus Integrators for Time-Dependent Schrödinger Equations , 2003, SIAM J. Numer. Anal..

[26]  Roy G. Gordon,et al.  New Method for Constructing Wavefunctions for Bound States and Scattering , 1969 .

[27]  Jerry D. Praeger.  Relaxational approach to solving the Schrödinger equation , 2001 .

[28]  Antonios Gonis,et al.  Green functions for ordered and disordered systems , 1992 .

[29]  H. Keller.  Numerical Solution of Two Point Boundary Value Problems , 1976 .