Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

The Journal of Chemical Physics

Propagators for the time-dependent Kohn–Sham equations

Alberto Castro, Miguel A. L. Marques, Ángel Rubio

Published Aug 10, 2004• 600 citations
Publisher

Abstract

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.

DOI
10.1063/1.1774980

Related papers

The time-dependent Kohn–Sham scheme

Carsten A. Ullrich

This chapter introduces the time-dependent Kohn-Sham scheme. The first section shows how an interacting time-dependent many-body system can be replaced by an auxiliary noninteracting system which reproduces the exact time-dependent density. This noninteracting system is called the Kohn-Sham system. It features an effective single-particle potential, consisting of external, Hartree, and exchange-correlation potential. The Kohn-Sham formalism is then generalized to spin-dependent systems. The exchange-correlation potential needs to be approximated in practice. Of particular importance is the adiabatic approximation, which uses any static exchange-correlation potential and evaluates it with time-dependent densities. The last two sections deal with practical issues related to the numerical implementation of the time-dependent Kohn-Sham scheme, namely, the meaning of self-consistency and unitary time propagation. The Crank-Nicholson algorithm is introduced, and important practical concepts such as the predictor-corrector scheme and absorbing boundary conditions are discussed.

Oxford University Press eBooks•Dec 22, 2011•1 citation

Time-dependent potential through an Ansatz for the Kohn-Sham orbitals

Rudolph J Magyar

Given the time-evolution of an electron charge density, the local potential in Kohn-Sham time-dependent density functional theory (KS-TDDFT) can be modeled as a sum of instantaneous and dynamic contributions by assuming a certain form of the time-dependent Kohn-Sham (TD-KS) orbitals. The instantaneous part is obtained numerically using methods from ground-state density functional theory (DFT) and the dynamic part is expressed in terms of a velocity potential that depends on the electron current density. The suggested form of the TD-KS orbitals satisfies several known constraints (orthonormality, N-representability, J-representability), and the domain of validity is shown to depend on the evolution of the instantaneous quantities. Through this decomposition, we can relate time-dependent and ground -state V-representability. The resulting potentials are shown numerically to approximate the exact time-dependent Kohn-Sham potentials for a set of 3 non-singlet two-particle systems (a Kohn-mode, a Coulomb explosion, and a double quantum well) where the exact solutions and reference densities are known or obtained through configuration interaction (CI) approaches.

arXiv:1611.02636•Nov 8, 2016•0 citations

Propagators for Quantum-Classical Models: Commutator-Free Magnus Methods

Adrián Gómez Pueyo, Sergio Blanes, Alberto Castro

We consider the numerical propagation of models that combine both quantum and classical degrees of freedom, usually, electrons and nuclei, respectively. We focus, in our computational examples, on the case in which the quantum electrons are modeled with time-dependent density-functional theory, although the methods discussed below can be used with any other level of theory. Often, for these so-called quantum-classical molecular dynamics models, one uses some propagation technique to deal with the quantum part and a different one for the classical equations. While the resulting procedure may, in principle, be consistent, it can however spoil some of the properties of the methods, such as the accuracy order with respect to the time step or the preservation of the geometrical structure of the equations. Few methods have been developed specifically for hybrid quantum-classical models. We propose using the same method for both the quantum and classical particles, in particular, one family of techniques that proves to be very efficient for the propagation of Schrödinger-like equations: the (quasi)-commutator free Magnus expansions. These have been developed, however, for linear systems, yet our problem is nonlinear: formally, the full quantum-classical system can be rewritten as a nonlinear Schrödinger equation, i.e., one in which the Hamiltonian depends on the system itself. The Magnus expansion algorithms for linear systems require the application of the Hamiltonian at intermediate points in a given propagating interval. For nonlinear systems, this poses a problem as this Hamiltonian is unknown due to its dependence on the state. We approximate it by employing a higher order extrapolation using previous steps as input. The resulting technique can then be regarded as a multistep technique or, alternatively, as a predictor corrector formula.

Journal of Chemical Theory and Computation•Jan 30, 2020•6 citations

Exact time-dependent density-functional potentials for strongly correlated tunneling electrons

M. J. P. Hodgson, J. D. Ramsden, Jacob B. J. Chapman +2 more

By propagating the many-body Schr\"odinger equation, we determine the exact time-dependent Kohn-Sham potential for a system of strongly correlated electrons which undergo field-induced tunneling. Numerous features are entirely absent from the approximations commonly used in time-dependent density-functional theory. The self-interaction correction is strong and time dependent, owing to electron localization, and prominent dynamic spatial potential steps arise from minima in the charge density, as modified by the Coulomb interaction experienced by the partially tunneled electron.

arXiv:1312.2534•Physical Review B•Dec 2, 2013•52 citations

Construction of a Kohn-Sham Iteration Scheme in Time-Dependent Lattice Density-Functional Theory

B.-S. Mehmel

In the present thesis, a rigorous mathematical formulation of time-dependent density-functional theory for lattice systems is derived, serving a formally exact approach towards solving time-dependent many-particle Schrödinger problems. After introducing the necessary mathematical foundations, a one-to-one mapping from external potentials onto electronic densities obtained by solving the Schrödinger equation is introduced. Its properties and the related Banach spaces of scalar potentials and electronic densities are discussed in detail. It allows to precisely formulate the Kohn-Sham Iteration scheme, mapping a given many-particle to an effective single-particle Schrödinger problem, both generating the exact same density. Employing the Banach fixed point theorem, the Kohn-Sham Iteration scheme is proven to be convergent depending on the chosen initial density.

•Jan 1, 2019•0 citations

EXPONENTIAL PROPAGATORS (INTEGRATORS) FOR THE TIME-DEPENDENT SCHRÖDINGER EQUATION

André D. Bandrauk, Huizhong Lu

The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.

Journal of Theoretical and Computational Chemistry•Apr 23, 2013•70 citations
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.
This chapter introduces the time-dependent Kohn-Sham scheme. The first section shows how an interacting time-dependent many-body system can be replaced by an auxiliary noninteracting system which reproduces the exact time-dependent density. This noninteracting system is called the Kohn-Sham system. It features an effective single-particle potential, consisting of external, Hartree, and exchange-correlation potential. The Kohn-Sham formalism is then generalized to spin-dependent systems. The exchange-correlation potential needs to be approximated in practice. Of particular importance is the adiabatic approximation, which uses any static exchange-correlation potential and evaluates it with time-dependent densities. The last two sections deal with practical issues related to the numerical implementation of the time-dependent Kohn-Sham scheme, namely, the meaning of self-consistency and unitary time propagation. The Crank-Nicholson algorithm is introduced, and important practical concepts such as the predictor-corrector scheme and absorbing boundary conditions are discussed.
Given the time-evolution of an electron charge density, the local potential in Kohn-Sham time-dependent density functional theory (KS-TDDFT) can be modeled as a sum of instantaneous and dynamic contributions by assuming a certain form of the time-dependent Kohn-Sham (TD-KS) orbitals. The instantaneous part is obtained numerically using methods from ground-state density functional theory (DFT) and the dynamic part is expressed in terms of a velocity potential that depends on the electron current density. The suggested form of the TD-KS orbitals satisfies several known constraints (orthonormality, N-representability, J-representability), and the domain of validity is shown to depend on the evolution of the instantaneous quantities. Through this decomposition, we can relate time-dependent and ground -state V-representability. The resulting potentials are shown numerically to approximate the exact time-dependent Kohn-Sham potentials for a set of 3 non-singlet two-particle systems (a Kohn-mode, a Coulomb explosion, and a double quantum well) where the exact solutions and reference densities are known or obtained through configuration interaction (CI) approaches.
We consider the numerical propagation of models that combine both quantum and classical degrees of freedom, usually, electrons and nuclei, respectively. We focus, in our computational examples, on the case in which the quantum electrons are modeled with time-dependent density-functional theory, although the methods discussed below can be used with any other level of theory. Often, for these so-called quantum-classical molecular dynamics models, one uses some propagation technique to deal with the quantum part and a different one for the classical equations. While the resulting procedure may, in principle, be consistent, it can however spoil some of the properties of the methods, such as the accuracy order with respect to the time step or the preservation of the geometrical structure of the equations. Few methods have been developed specifically for hybrid quantum-classical models. We propose using the same method for both the quantum and classical particles, in particular, one family of techniques that proves to be very efficient for the propagation of Schrödinger-like equations: the (quasi)-commutator free Magnus expansions. These have been developed, however, for linear systems, yet our problem is nonlinear: formally, the full quantum-classical system can be rewritten as a nonlinear Schrödinger equation, i.e., one in which the Hamiltonian depends on the system itself. The Magnus expansion algorithms for linear systems require the application of the Hamiltonian at intermediate points in a given propagating interval. For nonlinear systems, this poses a problem as this Hamiltonian is unknown due to its dependence on the state. We approximate it by employing a higher order extrapolation using previous steps as input. The resulting technique can then be regarded as a multistep technique or, alternatively, as a predictor corrector formula.
By propagating the many-body Schr\"odinger equation, we determine the exact time-dependent Kohn-Sham potential for a system of strongly correlated electrons which undergo field-induced tunneling. Numerous features are entirely absent from the approximations commonly used in time-dependent density-functional theory. The self-interaction correction is strong and time dependent, owing to electron localization, and prominent dynamic spatial potential steps arise from minima in the charge density, as modified by the Coulomb interaction experienced by the partially tunneled electron.
In the present thesis, a rigorous mathematical formulation of time-dependent density-functional theory for lattice systems is derived, serving a formally exact approach towards solving time-dependent many-particle Schrödinger problems. After introducing the necessary mathematical foundations, a one-to-one mapping from external potentials onto electronic densities obtained by solving the Schrödinger equation is introduced. Its properties and the related Banach spaces of scalar potentials and electronic densities are discussed in detail. It allows to precisely formulate the Kohn-Sham Iteration scheme, mapping a given many-particle to an effective single-particle Schrödinger problem, both generating the exact same density. Employing the Banach fixed point theorem, the Kohn-Sham Iteration scheme is proven to be convergent depending on the chosen initial density.
The time-dependent Schrödinger Equation (TDSE) is a parabolic partial differential equation (PDE) comparable to a diffusion equation but with imaginary time. Due to its first order time derivative, exponential integrators or propagators are natural methods to describe evolution in time of the TDSE, both for time-independent and time-dependent potentials. Two splitting methods based on Fer and/or Magnus expansions allow for developing unitary factorizations of exponentials with different accuracies in the time step △t. The unitary factorization of exponentials to high order accuracy depends on commutators of kinetic energy operators with potentials. Fourth-order accuracy propagators can involve negative or complex time steps, or real time steps only but with gradients of potentials, i.e. forces. Extending the propagators of TDSE's to imaginary time allows to also apply these methods to classical many-body dynamics, and quantum statistical mechanics of molecular systems.