Propagators for the time-dependent Kohn-Sham equations.

In this paper we address the problem of the numerical integration of the time-dependent Schrodinger equation i partial differential (t)phi=Hphi. In particular, we are concerned with the important case where H is the self-consistent Kohn-Sham Hamiltonian that stems from time-dependent functional theory. As the Kohn-Sham potential depends parametrically on the time-dependent density, H is in general time dependent, even in the absence of an external time-dependent field. The present analysis also holds for the description of the excited state dynamics of a many-electron system under the influence of arbitrary external time-dependent electromagnetic fields. Our discussion is separated in two parts: (i) First, we look at several algorithms to approximate exp(A), where A is a time-independent operator [e.g., A=-iDeltatH(tau) for some given time tau]. In particular, polynomial expansions, projection in Krylov subspaces, and split-operator methods are investigated. (ii) We then discuss different approximations for the time-evolution operator, such as the midpoint and implicit rules, and Magnus expansions. Split-operator techniques can also be modified to approximate the full time-dependent propagator. As the Hamiltonian is time dependent, problem (ii) is not equivalent to (i). All these techniques have been implemented and tested in our computer code OCTOPUS, but can be of general use in other frameworks and implementations.

[1]  藤田 純一,et al.  A.L. Fetter and J.D. Walecka: Quantum Theory of Many-Particle Systems, McGraw-Hill Book Co., New York, 1971, 601頁, 15×23cm, 7,800円. , 1971 .

[2]  Hua Guo,et al.  The Chebyshev propagator for quantum systems , 1999 .

[3]  C. Lubich,et al.  On Krylov Subspace Approximations to the Matrix Exponential Operator , 1997 .

[4]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[5]  M. Feit,et al.  Solution of the Schrödinger equation by a spectral method , 1982 .

[6]  H. Tal-Ezer,et al.  An accurate and efficient scheme for propagating the time dependent Schrödinger equation , 1984 .

[7]  Michele Parrinello,et al.  Integrating the Car–Parrinello equations. I. Basic integration techniques , 1994 .

[8]  A. Fetter,et al.  Quantum Theory of Many-Particle Systems , 1971 .

[9]  Nikos L. Doltsinis,et al.  Quantum Simulations of complex many-body systems:from theory to algorithms , 2002 .

[10]  Marlis Hochbruck,et al.  Exponential Integrators for Large Systems of Differential Equations , 1998, SIAM J. Sci. Comput..

[11]  M. Petersilka,et al.  DENSITY FUNCTIONAL THEORY OF TIME-DEPENDENT PHENOMENA , 1996 .

[12]  Robert van Leeuwen,et al.  KEY CONCEPTS IN TIME-DEPENDENT DENSITY-FUNCTIONAL THEORY , 2001 .

[13]  M. Feit,et al.  Solution of the Schrödinger equation by a spectral method II: Vibrational energy levels of triatomic molecules , 1983 .

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

[15]  Robert E. Wyatt,et al.  Study, extension, and application of Floquet theory for quantum molecular systems in an oscillating field , 1983 .

[16]  L. Reining,et al.  Electronic excitations: density-functional versus many-body Green's-function approaches , 2002 .

[17]  R. Kosloff.  Time-dependent quantum-mechanical methods for molecular dynamics , 1988 .

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

[19]  George F. Bertsch,et al.  Time-dependent local-density approximation in real time , 1996 .

[20]  J. Andrew McCammon,et al.  A comparative study of time dependent quantum mechanical wave packet evolution methods , 1992 .

[21]  S. Mukamel,et al.  Intermolecular forces and nonbonded interactions: Superoperator nonlinear time-dependent density-functional-theory response approach , 2004 .

[22]  R. Nalewajski.  Density Functional Theory II , 1996 .

[23]  Y. Saad.  Analysis of some Krylov subspace approximations to the matrix exponential operator , 1992 .

[24]  Masaru Tsukada,et al.  Efficient method for simulating quantum electron dynamics under the time-dependent Kohn-Sham equation. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[26]  Masuo Suzuki,et al.  General Nonsymmetric Higher-Order Decomposition of Exponential Operators and Symplectic Integrators , 1992 .

[27]  Á. Rubio,et al.  octopus: a first-principles tool for excited electron-ion dynamics. , 2003 .

[28]  D. Chong.  Recent Advances in Density Functional Methods Part III , 2002 .

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

[30]  S. Mukamel.  Principles of Nonlinear Optical Spectroscopy , 1995 .

[31]  T. Park,et al.  Unitary quantum time evolution by iterative Lanczos reduction , 1986 .

[32]  Marlis Hochbruck,et al.  Exponential Integrators for Quantum-Classical Molecular Dynamics , 1999 .

[33]  T. Beck.  Real-space mesh techniques in density-functional theory , 2000, cond-mat/0006239.

[34]  Roi Baer,et al.  A method for ab initio nonlinear electron-density evolution , 2001 .

[35]  Martin Head-Gordon,et al.  Two-body coupled cluster expansions , 2001 .

[36]  André D. Bandrauk,et al.  Exponential split operator methods for solving coupled time-dependent Schrödinger equations , 1993 .

[37]  Masuo Suzuki,et al.  Convergence of unitary and complex decompositions of exponential operators , 1993 .

[38]  Cleve B. Moler,et al.  Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later , 1978, SIAM Rev..

[39]  E. Gross,et al.  Time-dependent density functional theory. , 2004, Annual review of physical chemistry.