Solution of time-independent Schrödinger equation by the imaginary time propagation method. (English) Zbl 1110.65096

Summary: Numerical solution of eigenvalues and eigenvectors of large matrices originating from discretization of linear and nonlinear Schrödinger equations using the imaginary time propagation (ITP) method is described. Convergence properties and accuracy of 2nd and 4th order operator-splitting methods for the ITP method are studied using numerical examples. The natural convergence of the method is further accelerated with a new dynamic time step adjustment method. The results show that the ITP method has better scaling with respect to matrix size as compared to the implicitly restarted Lanczos method. An efficient parallel implementation of the ITP method for shared memory computers is also demonstrated.


65M70 Spectral, collocation and related methods for initial value and initial-boundary value problems involving PDEs
81Q05 Closed and approximate solutions to the Schrödinger, Dirac, Klein-Gordon and other equations of quantum mechanics
35Q55 NLS equations (nonlinear Schrödinger equations)
65Y05 Parallel numerical computation
65F15 Numerical computation of eigenvalues and eigenvectors of matrices


Full Text: DOI


[1] Eloranta, J.; Apkarian, V.A., A time dependent density functional treatment of superfluid dynamics: equilibration of the electron bubble in superfluid ^{4}he, J. chem. phys., 117, 10139, (2002)
[2] L. Lehtovaara, J. Eloranta, Small multielectron bubbles in bulk superfluid ^{4}He, in: Proceedings of 24th International Conference on Low Temperature Physics, in press. · Zbl 1046.81022
[3] R.B. Lehoucq, D.C. Sorensen, P.A. Vu, C. Wang, ARPACK: an implementation of an implicitly restarted Arnoldi method for computing some of the eigenvalues and eigenvectors of a large sparse matrix, Dept. Comput. Appl. Math., Rice University, 1996. For further details see ⟨http://www.caam.rice.edu/software/ARPACK⟩.
[4] Auer, J.; Krotscheck, E.; Chin, S.A., A fourth-order real-space algorithm for solving local Schrödinger equations, J. chem. phys., 115, 6841, (2001)
[5] Aichinger, M.; Krotscheck, E., A fast configuration space method for solving local kohn – sham equations, Comput. mater. sci., 34, 188, (2005)
[6] Roy, A.K.; Gupta, N.; Deb, D.M., Time-dependent quantum-mechanical calculation of ground and excited states of anharmonic and double-well oscillators, Phys. rev. A, 65, 012109, (2001)
[7] Lehtovaara, L.; Kiljunen, T.; Eloranta, J., Efficient numerical methods for simulating static and dynamic properties of superfluid helium, J. comp. phys., 194, 78, (2004) · Zbl 1052.82033
[8] Martikainen, J.; Rossi, T.; Toivanen, J., Computation of a few smallest eigenvalues of elliptic operators using fast elliptic solves, Commun. numer. meth. eng., 17, 521, (2001) · Zbl 0988.65098
[9] Isaacson, E.; Keller, H.B., Analysis of numerical methods, (1994), Dover Publications, Inc. New York · Zbl 0168.13101
[10] Press, W.H.; Teukolsky, S.A.; Vetterling, W.T.; Flannery, B.P., Numerical recipes in C: the art of scientific computing, (1999), Cambridge University Press Cambridge, UK
[11] Suzuki, M., ()
[12] Chin, S.A., Phys. lett. A, 226, 344, (1997)
[13] Frigo, M.; Johnson, S.G., FFTW: an adaptive software architecture for the FFT, ICASSP conference Proceedings, 3, 1381, (1998), For more details see
[14] Lehtovaara, L.; Eloranta, J., A 2-level anisotropic electronic system in superfluid ^{4}he, J. low temp. phys., 138, 91, (2005)
[15] Golub, G.H.; Van Loan, C.F., Matrix computations, (1996), John Hopkins University Press Baltimore, MD · Zbl 0865.65009
[16] OpenMP C and C++ application program interface. For further details see http://www.openmp.org.
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.