Abstract

To use the short iterative Lanczos algorithm to solve the time-dependent Schroedinger equation, one must choose, for a given Lanczos space size, a time step. We compare the derivation of the well-known Lubich and Hochbruck time step from SIAM J. Numer. Anal. 34 (1997) 1911 with the a priori time step we proposed in Mohankumar and Carrington (MC) Comput. Phys. Commun., 181 (2010) 1859 and demonstrate that the MC time step is somewhat larger, i.e., that the MC error bound is tighter. In addition, we use the MC approach to derive an error bound and time step for imaginary time propagation. The error bound we derive is much tighter than the error bound of Stewart and Leyk.

Highlights

  • The short iterative Lanczos (SIL) algorithm [1] for solving the time-dependent Schroedinger equation (TDSE), i ∂ψ( x, t) = Ĥψ( x, t), ∂t (1)is widely used in chemical physics and other fields

  • In typical calculations, ∆t is small enough that 1 − α ∼ 1 and e−my ∼ 1 and in this case the time step ∆t Mohankumar and Carrington (MC) is larger than ∆t L by a factor of (8πm)1/m . ∆t MC is larger than ∆t L whenever α < 0.8

  • Once an error bound has been found, it is straightforward to derive an equation for the best time step to use

Read more

Summary

Introduction

The short iterative Lanczos (SIL) algorithm [1] for solving the time-dependent Schroedinger equation (TDSE) (in atomic units), i. M is small enough that loss of orthogonality of the Lanczos vectors is unimportant. We shall assume that a(t = 0) is the Lanczos starting vector. To use the SIL algorithm, one chooses a value of m, propagates from t = 0 to t = ∆t, and uses a(∆t) as the starting vector and propagates again from t = ∆t to t = 2∆t, etc. Small enough that loss of orthogonality is not a problem, one must choose ∆t. We compare the derivations of the Lubich and MC equations and show that. We derive and test a new equation for ∆t for propagating in imaginary time

Common Starting Point
Lubich Derivation
Mohankumar–Carrington Derivation
Test Calculations for Real-Time Propagation
Imaginary Time Propagation
An Error Bound from a Geometric Series
Test Calculations
Conclusions
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call