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 SrTiO3, diamond, and metallic lithium. While SPDFT results for SrTiO3 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.
Basing on our hierarchical equations of motion for time-dependent quantum transport [X. Zheng, G. H. Chen, Y. Mo, S. K. Koo, H. Tian, C. Y. Yam, and Y. J. Yan, J. Chem. Phys. 133, 114101 (2010)], we develop an efficient and accurate numerical algorithm to solve the Liouville-von-Neumann equation. We solve the real-time evolution of the reduced single-electron density matrix at the tight-binding level. Calculations are carried out to simulate the transient current through a linear chain of atoms, with each represented by a single orbital. The self-energy matrix is expanded in terms of multiple Lorentzian functions, and the Fermi distribution function is evaluated via the Padè spectrum decomposition. This Lorentzian-Padè decomposition scheme is employed to simulate the transient current. With sufficient Lorentzian functions used to fit the self-energy matrices, we show that the lead spectral function and the dynamics response can be treated accurately. Compared to the conventional master equation approaches, our method is much more efficient as the computational time scales cubically with the system size and linearly with the simulation time. As a result, the simulations of the transient currents through systems containing up to one hundred of atoms have been carried out. As density functional theory is also an effective one-particle theory, the Lorentzian-Padè decomposition scheme developed here can be generalized for first-principles simulation of realistic systems.
For the calculation of neutral excitations, time-dependent density functional theory (TDDFT) is an exact reformulation of the many-body time-dependent Schrödinger equation, based on knowledge of the density instead of the many-body wavefunction. The density can be determined in an efficient scheme by solving one-particle non-interacting Schrödinger equations—the Kohn–Sham equations. The complication of the problem is hidden in the—unknown—time-dependent exchange and correlation potential that appears in the Kohn–Sham equations and for which it is essential to find good approximations. Many approximations have been suggested and tested for finite systems, where even the very simple adiabatic local-density approximation (ALDA) has often proved to be successful. In the case of solids, ALDA fails to reproduce optical absorption spectra, which are instead well described by solving the Bethe–Salpeter equation of many-body perturbation theory (MBPT). On the other hand, ALDA can lead to excellent results for loss functions (at vanishing and finite momentum transfer). In view of this and thanks to recent successful developments of improved linear-response kernels derived from MBPT, TDDFT is today considered a promising alternative to MBPT for the calculation of electronic spectra, even for solids. After reviewing the fundamentals of TDDFT within linear response, we discuss different approaches and a variety of applications to extended systems.
The problem of quantum dynamics in open systems has gained attention in recent decades and not the least due to the advances made in quantum transport in molecular systems. The main motivation behind quantum transport and molecular electronics is the futuristic goal to be able at some point to replace, or to complement, the silicon-based technology and to make the electronic devices faster. On a fundamental level, one has to deal with time-dependent processes where electron-electron or electron- phonon interactions are of great importance, and they can cause profound quantitative and qualitative changes on the physical and dynamical properties of electronic systems compared to the non-interacting case. Most of the studies of quantum transport have been focused on the steady-state description while neglecting the short-time dynamics. However, the dynamical effects are of great importance since fast- switching processes play a pivotal role in the operation of future devices. We studied the problem of time-dependent electron transport through the Anderson impurity model by using many-body perturbation theory (MBPT) together with Keldysh Green’s functions as well as with time-dependent density functional theory (TDDFT). These methods were compared with numerically exact time-dependent density-matrix renormalization group (tDMRG) method. We found that the many-body perturbation theory results beyond Hartree-Fock approximation were in close agreement with tDMRG results. In addition we studied the possibility of multistablity in the density and current of an interacting nanoscale junction as well as how to reversibly switch between the multiple solutions in time domain. An accurate theoretical treatment of electron correlation even in as simple model as an interacting electron gas at metallic densities still continues to be a challenge; especially description of features in the photoemission spectra due to electron correlations provides a theoretical challenge. The many-body perturbation theory yields a systematic way to study electron-electron (electron-phonon) correlations in various systems. One of the widely used approximations in MBPT is the GW approximation in which the bare interaction line is replaced with screened interaction line in the first order exchange diagram. The GW approximation gives good estimates for the band gap values close to experimental ones but especially the self-consistent GW approximation has a number of deficiencies like washing out of plasmon features and overestimation of bandwidths compared to experiment. One way to improve GW calculations is to include vertex corrections. Unfortunately, the straightforward inclusion of vertex corrections yields negative spectra in some frequency regions. We developed a diagrammatic approach to construct self-energy approximations with positive spectral properties. Our approach consists of expressing a self-energy of response diagram as a product of half-diagrams after which a minimal set of additional diagrams is identified to construct a perfect square. We applied this method to study vertex corrections in a homogeneous electron gas. In addition we analyzed the diagrammatic content of photocurrent with density functional theory. The expression for the photocurrent was obtained as an integral over the Kohn-Sham spectral function renormalized by effective potentials that depend on the exchange correlation kernel of current density functional theory. The expression for the photocurrent gives us the angular dependence of the photocurrent but it does not provide a direct access to the kinetic energy distribution of the photoelectrons.
Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg, Kohn, and Sham [4,5] in the mid-1960s. According to this theory, the fundamental physical information about a many-body system is provided by single-particle densities in a three-dimensional space, which are obtained variationally within a time-independent framework. When compared with other previous formalisms, DFT presents two clear advantages: (1) it is able to treat many-body problems in a sufficiently accurate way and (2) it is computationally simple. This explains why it is one of the most widely used theories to deal with electronic structure-the electronic ground-state energy as a function of the position of the atomic nuclei determines molecular structures and solids, providing at the same time the forces acting on the atomic nuclei when they are not at their equilibrium positions. At present, DFT is used routinely to solve many problems in gas phase and condensed matter. Furthermore, it has made possible the development of accurate molecular dynamics schemes in which the forces are evaluated quantum mechanically ‘‘on the fly.’’ Nonetheless, DFT is a fundamental tool provided the systems studied are relatively large; for small systems, standard methods based on the use of the wave function render quite accurate results [6]. Moreover, it is also worth stressing that all practical applications of DFT rely on essentially uncontrolled approximations [7] (e.g., the local density approximation [4,5], the local spin-density approximation, or generalized gradient approximations [8]), and therefore the validity of DFT is conditioned to its capability to provide fairly good values of the experimental data.
Accurate descriptions of the electronic properties of molecules, solids and other materials are crucial in the development of many modern technologies, and are therefore the subject of much scientific research. In tandem with experimental research, computational simulations play a huge role. Density functional theory (DFT) is the most popular method for performing ground-state electronic structure calculations due to its accuracy in many situations along with its computational efficiency. Its time-dependent extension (TDDFT) allows the modelling of the dynamics of systems of interacting electrons. (TD)DFT is an exact theory in principle, but requires approximations in practice, specifically to the many-body effects of electron exchange and correlation (xc), which are described by some unknown universal functional. Despite much success, it often proves inadequate in many scenarios of increasing interest. Therefore, we aim to contribute to the development of improved approximate density functionals. \n \nIn this thesis, we take the approach of studying small model systems for which we can compute the exact functionals. These systems are designed to exhibit the key features present in realistic quantum systems, such as atoms and molecules, which the exact universal (TD)DFT functional is capable of describing, but are often ill-described by the commonly used approximations. By analysing the nature of these exact functionals, we gain valuable insight into their fundamental properties, together with identifying where the current approximations are lacking, thereby allowing us to suggest how these approximations may be improved.
A linear-scaling first-principles quantum mechanical method is developed to evaluate the optical responses of large molecular systems. Instead of a many-body wave function, the equation of motion is solved for the reduced single-electron density matrix in the time domain. The locality of the reduced single-electron density matrix is utilized to ensure that computational time scales linearly with system size. The two-electron Coulomb integrals are evaluated with the fast multipole method, and the calculation of exchange-correlation quadratures utilizes the locality of an exchange-correlation functional and the integral prescreening technique. As an illustration, the resulting time-dependent density-functional theory is used to calculate the absorption spectra of polyacetylene oligomers and linear alkanes. The linear-scaling of computational time versus the system size is clearly demonstrated.