We present an extension of second-principles density functional theory (SPDFT) to perform time-dependent simulations. Our approach calculates the real-time, real-space evolution of the density matrix using the Liouville-von Neumann equation of motion, enabling the determination of optical and transport properties in very large systems, containing tens of thousands of atoms, on modest computational platforms. Unlike other methods, SPDFT is applicable to a wide range of materials, including both metals and insulators. We demonstrate its capabilities by computing the spectra of ${\mathrm{SrTiO}}_{3}$ , diamond, and metallic lithium. While SPDFT results for ${\mathrm{SrTiO}}_{3}$ closely match those obtained from DFT with linear perturbation theory, significant improvements are observed for diamond and metallic lithium. In diamond, the inclusion of electron-electron interactions during the density-matrix evolution yields spectra that more accurately reproduce Bethe-Salpeter equation results than perturbative DFT. In lithium, time-dependent SPDFT captures both interband transitions and the Drude peak, enabling detailed ab initio studies of transport properties beyond usual approximations.We present an extension of second-principles density functional theory (SPDFT) to perform time-dependent simulations. Our approach calculates the real-time, real-space evolution of the density matrix using the Liouville-von Neumann equation of motion, enabling the determination of optical and transport properties in very large systems, containing tens of thousands of atoms, on modest computational platforms. Unlike other methods, SPDFT is applicable to a wide range of materials, including both metals and insulators. We demonstrate its capabilities by computing the spectra of ${\mathrm{SrTiO}}_{3}$, diamond, and metallic lithium. While SPDFT results for ${\mathrm{SrTiO}}_{3}$ closely match those obtained from DFT with linear perturbation theory, significant improvements are observed for diamond and metallic lithium. In diamond, the inclusion of electron-electron interactions during the density-matrix evolution yields spectra that more accurately reproduce Bethe-Salpeter equation results than perturbative DFT. In lithium, time-dependent SPDFT captures both interband transitions and the Drude peak, enabling detailed ab initio studies of transport properties beyond usual approximations.

arXiv:2507.13824Physical review. B./Physical review. BNov 8, 20252 citations

Time-Dependent Density Functional Theory

Takao Tsuneda, Kimihiko Hirao

Jun 19, 20180 citations

Time-dependent density-matrix functional in Liouville space and the optical response of many-electron systems

Vladimir Chernyak, Shaul Mukamel

A time-dependent functional approach is developed for calculating the dynamical optical response and ground-state properties of nonrelativistic fermion systems. The functional (effective action) uses the particle-hole components of the reduced single-electron density matrix as variables and is constructed using the natural orbitals as a reference. It is defined on pairs of trajectories that constitute a Liouville space path in the phase space of the classical analog of the quantum system. A perturbative ``semiclassical'' calculation of the effective Liouville-space action is used to incorporate correlation effects in the linear optical response.A time-dependent functional approach is developed for calculating the dynamical optical response and ground-state properties of nonrelativistic fermion systems. The functional (effective action) uses the particle-hole components of the reduced single-electron density matrix as variables and is constructed using the natural orbitals as a reference. It is defined on pairs of trajectories that constitute a Liouville space path in the phase space of the classical analog of the quantum system. A perturbative ``semiclassical'' calculation of the effective Liouville-space action is used to incorporate correlation effects in the linear optical response.

Physical Review ANov 1, 199513 citations

Non-equilibrium Green's function formalism

Massimiliano Di Ventra

Within the Landauer approach to conduction discussed in the previous chapter, electron interactions have been included only at the mean-field level. This is quite a strong approximation, especially in nanojunctions, where large current densities are common. I therefore want to go beyond this level of description. I can follow two different routes. (1) I can develop a functional theory of quantum correlations via time-dependent effective single-particle equations, as I will do in Chapter 7 where I will use dynamical density-functional theories within the micro-canonical approach to conduction. (2) I can try to solve the time-dependent Schrödinger equation 1.16 – or its mixed-state version, the Liouville-von Neumann equation 1.60 – directly. Written this way, this last proposition seems hopeless. In reality, we can employ a many-body technique, known as the non-equilibrium Green's function formalism (NEGF), also referred to as the Keldysh formalism (Keldysh, 1964; Kadanoff and Baym, 1962), which allows us, at least in principle, to do just that: solve the time-dependent Schrödinger equation for an interacting many-body system exactly, from which one can, in principle, calculate the time-dependent current. This is done by solving equations of motion for specific time-dependent single-particle Green's functions, from which the physical properties of interest, such as the charge and current densities, can be obtained. I have stressed the term “time-dependent” several times, because, as I will show in a moment, the NEGF is “exact” only when one solves the time-dependent Schrödinger equation for a closed quantum system, subject to deterministic perturbations: the system is closed but not necessarily isolated. These perturbations may drive the system far away from its initial state of thermodynamic equilibrium (Keldysh, 1964).Within the Landauer approach to conduction discussed in the previous chapter, electron interactions have been included only at the mean-field level. This is quite a strong approximation, especially in nanojunctions, where large current densities are common. I therefore want to go beyond this level of description. I can follow two different routes. (1) I can develop a functional theory of quantum correlations via time-dependent effective single-particle equations, as I will do in Chapter 7 where I will use dynamical density-functional theories within the micro-canonical approach to conduction. (2) I can try to solve the time-dependent Schrödinger equation 1.16 – or its mixed-state version, the Liouville-von Neumann equation 1.60 – directly. Written this way, this last proposition seems hopeless. In reality, we can employ a many-body technique, known as the non-equilibrium Green's function formalism (NEGF), also referred to as the Keldysh formalism (Keldysh, 1964; Kadanoff and Baym, 1962), which allows us, at least in principle, to do just that: solve the time-dependent Schrödinger equation for an interacting many-body system exactly, from which one can, in principle, calculate the time-dependent current. This is done by solving equations of motion for specific time-dependent single-particle Green's functions, from which the physical properties of interest, such as the charge and current densities, can be obtained. I have stressed the term “time-dependent” several times, because, as I will show in a moment, the NEGF is “exact” only when one solves the time-dependent Schrödinger equation for a closed quantum system, subject to deterministic perturbations: the system is closed but not necessarily isolated. These perturbations may drive the system far away from its initial state of thermodynamic equilibrium (Keldysh, 1964).

Cambridge University Press eBooksAug 7, 20083 citations

Time-dependent density-functional theory

Arno Schindlmayr

Jan 1, 20060 citations

{\mathrm{SrTiO}}_{3}

, diamond, and metallic lithium. While SPDFT results for

{\mathrm{SrTiO}}_{3}

closely match those obtained from DFT with linear perturbation theory, significant improvements are observed for diamond and metallic lithium. In diamond, the inclusion of electron-electron interactions during the density-matrix evolution yields spectra that more accurately reproduce Bethe-Salpeter equation results than perturbative DFT. In lithium, time-dependent SPDFT captures both interband transitions and the Drude peak, enabling detailed ab initio studies of transport properties beyond usual approximations.