|
|
| Line 1: |
Line 1: |
| '''Time-dependent density functional theory''' ('''TDDFT''') is a [[quantum mechanical]] theory used in physics and chemistry to investigate the properties and [[Molecular dynamics|dynamics]] of many-body systems in the presence of time-dependent potentials, such as electric or magnetic fields. The effect of such fields on molecules and solids can be studied with TDDFT to extract features like excitation energies, frequency-dependent response properties, and photoabsorption spectra. | | Myrtle Benny is how I'm called and I feel comfortable when individuals use the full title. Since she was 18 she's been operating as a meter reader but home std [http://dateplate.com/dayzfe std testing at home] test she's always wanted her personal company. [http://www.Womenshealthmag.com/sex-and-relationships/sex-lubes Puerto Rico] is exactly where he's been living for years and he will never transfer. It's home std test not at home std testing a typical factor but what she likes performing is base leaping [http://www.hamcass.org/index.php?document_srl=268708&mid=gido this contact form] and now she is attempting to make cash with it.<br><br>My blog - [http://www.smylestream.org/groups/solid-advice-in-relation-to-candida/ at home std testing] |
| | |
| TDDFT is an extension of [[density functional theory]] (DFT), and the conceptual and computational foundations are analogous – to show that the (time-dependent) [[wave function]] is equivalent to the (time-dependent) [[electronic density]], and then to derive the effective potential of a fictitious non-interacting system which returns the same density as any given interacting system. The issue of constructing such a system is more complex for TDDFT, most notably because the time-dependent effective potential at any given instant depends on value of the density at all previous times. Consequently the development of time-dependent approximations for the implementation of TDDFT is behind that of DFT, with applications routinely ignoring this memory requirement.
| |
| | |
| == Overview ==
| |
| | |
| The formal foundation of TDDFT is the '''Runge-Gross''' ('''RG''') '''theorem''' (1984)<ref>{{cite journal|last=Runge|first=Erich|coauthors=Gross, E. K. U.|year=1984|title=Density-Functional Theory for Time-Dependent Systems|journal=Phys. Rev. Lett.|volume=52|issue=12|pages=997–1000|doi=10.1103/PhysRevLett.52.997|bibcode=1984PhRvL..52..997R}}</ref> – the time-dependent analogue of the Hohenberg-Kohn (HK) theorem (1964).<ref>{{cite journal|title=Inhomogeneous electron gas|journal=Phys. Rev.|year=1964|first=P.|last=Hohenberg|coauthors=Kohn, W.|volume=136|issue=3B|pages=B864–B871|doi=10.1103/PhysRev.136.B864|bibcode = 1964PhRv..136..864H }}</ref> The RG theorem shows that, for a given initial wavefunction, there is a unique mapping between the time-dependent external potential of a system and its time-dependent density. This implies that the many-body wavefunction, depending upon 3''N'' variables, is equivalent to the density, which depends upon only 3, and that all properties of a system can thus be determined from knowledge of the density alone. Unlike in DFT, there is no general minimization principle in time-dependent quantum mechanics. Consequently the proof of the RG theorem is more involved than the HK theorem.
| |
| | |
| Given the RG theorem, the next step in developing a computationally useful method is to determine the fictitious non-interacting system which has the same density as the physical (interacting) system of interest. As in DFT, this is called the (time-dependent) Kohn-Sham system. This system is formally found as the [[stationary point]] of an [[Action (physics)|action]] functional defined in the [[Keldysh formalism]].<ref>{{cite journal|last=van Leeuwen|first=Robert|year=1998|title=Causality and Symmetry in Time-Dependent Density-Functional Theory|journal=Phys. Rev. Lett.|volume=80|issue=6|pages=1280–283|doi=10.1103/PhysRevLett.80.1280|bibcode=1998PhRvL..80.1280V}}</ref>
| |
| | |
| The most popular application of TDDFT is in the calculation of the energies of excited states of isolated systems and, less commonly, solids. Such calculations are based on the fact that the linear response function – that is, how the electron density changes when the external potential changes – has poles at the exact excitation energies of a system. Such calculations require, in addition to the exchange-correlation potential, the exchange-correlation kernel – the [[functional derivative]] of the exchange-correlation potential with respect to the density.<ref>{{cite book|last=Casida |first=M. E.|coauthors=C. Jamorski, F. Bohr, J. Guan, and D. R. Salahub|title=Theoretical and Computational Modeling of NLO and Electronic Materials|editor=S. P. Karna and A. T. Yeates|publisher=ACS Press|location=Washington, D.C.|year=1996|url=http://dcm.ujf-grenoble.fr/PERSONNEL/CT/casida/research/chong.ps|page=145–}}</ref><ref>{{cite journal|last=Petersilka|first=M.|coauthors=U. J. Gossmann, and E.K.U. Gross|year=1996|title=Excitation Energies from Time-Dependent Density-Functional Theory|journal=Phys. Rev. Lett.|volume=76|issue=8|pages=1212–1215|doi=10.1103/PhysRevLett.76.1212|bibcode=1996PhRvL..76.1212P}}</ref>
| |
| | |
| ==Formalism==
| |
| ===Runge-Gross theorem===
| |
| {{main|Runge-Gross theorem}}
| |
| | |
| The approach of Runge and Gross considers a single-component system in the presence of a time-dependent [[scalar field]] for which the [[Hamiltonian (quantum mechanics)|Hamiltonian]] takes the form
| |
| | |
| :<math>\hat{H}(t)=\hat{T}+\hat{V}_{\mathrm{ext}}(t)+\hat{W},</math>
| |
| | |
| where ''T'' is the kinetic energy operator, ''W'' the electron-electron interaction, and ''V''<sub>ext</sub>(''t'') the external potential which along with the number of electrons defines the system. Nominally, the external potential contains the electrons' interaction with the nuclei of the system. For non-trivial time-dependence, an additional explicitly time-dependent potential is present which can arise, for example, from a time-dependent electric or magnetic field. The many-body wavefunction evolves according to the [[Schrödinger_equation#Time_dependent_equation|time-dependent Schrödinger equation]] under a single [[initial condition]],
| |
| | |
| :<math>\hat{H}(t)|\Psi(t)\rangle=i\hbar\frac{\partial}{\partial t}|\Psi(t)\rangle,\ \ \ |\Psi(0)\rangle=|\Psi\rangle.</math>
| |
| | |
| Employing the Schrödinger equation as its starting point, the Runge-Gross theorem shows that at any time, the density uniquely determines the external potential. This is done in two steps:
| |
| | |
| #Assuming that the external potential can be expanded in a [[Taylor series]] about a given time, it is shown that two external potentials differing by more than an additive constant generate different [[current density|current densities]].
| |
| #Employing the [[continuity equation]], it is then shown that for finite systems, different current densities correspond to different electron densities.
| |
| | |
| ===Time-dependent Kohn-Sham system===
| |
| | |
| For a given interaction potential, the RG theorem shows that the external potential uniquely determines the density. The Kohn-Sham approaches chooses a non-interacting system (that for which the interaction potential is zero) in which to form the density that is equal to the interacting system. The advantage of doing so lies in the ease in which non-interacting systems can be solved – the wave function of a non-interacting system can be represented as a [[Slater determinant]] of single-particle [[Molecular orbital|orbitals]], each of which are determined by a single [[partial differential equation]] in three variable – and that the kinetic energy of a non-interacting system can be expressed exactly in terms of those orbitals. The problem is thus to determine a potential, denoted as ''v''<sub>s</sub>('''r''',''t'') or ''v''<sub>KS</sub>('''r''',''t''), that determines a non-interacting Hamiltonian, ''H''<sub>s</sub>,
| |
| | |
| :<math>\hat{H}_{s}(t) = \hat{T}+\hat{V}_{s}(t),</math>
| |
| | |
| which in turn determines a determinantal wave function
| |
| | |
| :<math>\hat{H}_{s}(t)|\Phi(t)\rangle=i\frac{\partial}{\partial t}|\Phi(t)\rangle,\ \ \ |\Phi(0)\rangle=|\Phi\rangle,</math>
| |
| | |
| which is constructed in terms of a set of ''N'' orbitals which obey the equation,
| |
| | |
| :<math>\left(-\frac{1}{2}\nabla^{2}+v_{s}(\mathbf{r},t)\right)\phi_{i}(\mathbf{r},t)=i\frac{\partial}{\partial t}\phi_{i}(\mathbf{r},t)\ \ \ \phi_{i}(\mathbf{r},0)=\phi_{i}(\mathbf{r}),</math>
| |
| | |
| and generate a time-dependent density | |
| | |
| :<math>\rho_{s}(\mathbf{r},t)=\sum_{i=1}^{N}|\phi_{i}(\mathbf{r},t)|^{2},</math>
| |
| | |
| such that ''ρ''<sub>s</sub> is equal to the density of the interacting system at all times:
| |
| | |
| :<math>\rho_{s}(\mathbf{r},t)=\rho(\mathbf{r},t).</math>
| |
| | |
| If the potential ''v''<sub>s</sub>('''r''',''t'') can be determined, or at the least well-approximated, then the original Schrödinger equation, a single partial differential equation in 3''N'' variables, has been replaced by ''N'' differential equations in 3 dimensions, each differing only in the initial condition.
| |
| | |
| The problem of determining approximations to the Kohn-Sham potential is challenging. Analogously to DFT, the time-dependent KS potential is decomposed to extract the external potential of the system and the time-dependent Coulomb interaction, ''v''<sub>J</sub>. The remaining component is the exchange-correlation potential:
| |
| | |
| :<math>v_{s}(\mathbf{r},t)=v_{\rm ext}(\mathbf{r},t)+v_{J}(\mathbf{r},t)+v_{\rm xc}(\mathbf{r},t).\,</math>
| |
| | |
| In their seminal paper, Runge and Gross approached the definition of the KS potential through an action-based argument starting from the [[Dirac action]]
| |
| | |
| :<math>A[\Psi]=\int\mathrm{d}t\ \langle\Psi(t)|H-i\frac{\partial}{\partial t}|\Psi(t)\rangle.</math> | |
| | |
| Treated as a functional of the wave function, ''A''[Ψ], variations of the wave function yield the many-body Schrödinger equation as the stationary point. Given the unique mapping between densities and wave function, Runge and Gross then treated the Dirac action as a density functional,
| |
| | |
| :<math>A[\rho]=A[\Psi[\rho]],\,</math>
| |
| | |
| and derived a formal expression for the exchange-correlation component of the action, which determines the exchange-correlation potential by functional differentiation. Later it was observed that an approach based on the Dirac action yields paradoxical conclusions when considering the causality of the response functions it generates.<ref>{{cite book|last=Gross|first=E. K. U. |coauthors=C. A. Ullrich and U. J. Gossman|title=Density Functional Theory|editor=E. K. U. Gross and R. M. Dreizler|publisher=Plenum Press|location=New York|year=1995|series=B|volume=337|isbn=0-387-51993-9}}</ref> The density response function, the functional derivative of the density with respect to the external potential, should be causal: a change in the potential at a given time can not affect the density at earlier times. The response functions from the Dirac action however are symmetric in time so lack the required causal structure. An approach which does not suffer from this issue was later introduced with through an action based on the [[Keldysh formalism]] of complex-time path integration.
| |
| | |
| == Linear response TDDFT == | |
| Linear-response TDDFT can be used if the external perturbation is small in the
| |
| sense that it does not completely destroy the ground-state structure of the system. In this case
| |
| one can analyze the linear response of the system. This is a great advantage as, to first order,
| |
| the variation of the system will depend only on the ground-state wave-function so that we can
| |
| simply use all the properties of DFT.
| |
| | |
| Consider a small time-dependent external perturbation <math>\delta V^{ext}(t)</math>.
| |
| This gives
| |
| :<math>H'(t)=H + \delta V^{ext}(t)</math>
| |
| :<math>H'_{KS}[\rho](t)=H_{KS}[\rho]+\delta V_H[\rho](t)+\delta V_{xc}[\rho](t)+\delta V^{ext}(t)</math>
| |
| and looking at the linear response of the density | |
| :<math>\delta \rho(\mathbf{r}t)= \chi(\mathbf{r}t,\mathbf{r'}t')
| |
| \delta V^{ext}(\mathbf{r'}t')</math>
| |
| :<math>\delta \rho(\mathbf{r}t)=\chi_{KS}(\mathbf{r}t,\mathbf{r'}t')
| |
| \delta V^{eff}[\rho](\mathbf{r'}t')</math>
| |
| where <math>\delta V^{eff}[\rho](t)=\delta V^{ext}(t)+\delta V_H[\rho](t)+\delta V_{xc}[\rho](t)</math>
| |
| Here and in the following it is assumed that primed variables are integrated.
| |
| | |
| Within the linear-response domain, the variation of the Hartree (H) and the
| |
| exchange-correlation (xc) potential to linear order may be expanded with respect to the density variation
| |
| :<math>\delta V_H[\rho](\mathbf{r})=\frac{\delta V_H[\rho]}{\delta\rho}\delta\rho=
| |
| \frac{1}{|\mathbf{r}-\mathbf{r'}|}\delta\rho(\mathbf{r'})</math>
| |
| and
| |
| :<math>\delta V_{xc}[\rho](\mathbf{r})=\frac{\delta V_{xc}[\rho]}{\delta\rho}\delta\rho=
| |
| f_{xc}(\mathbf{r}t,\mathbf{r'}t')\delta\rho(\mathbf{r'})</math>
| |
| Finally, inserting this relation in the response equation for the KS system and comparing
| |
| the resultant equation with the response equation for the physical system yields the Dyson
| |
| equation of TDDFT:
| |
| :<math>\chi(\mathbf{r}_1t_1,\mathbf{r}_2t_2)=\chi_{KS}(\mathbf{r_1}t_1,\mathbf{r}_2t_2)+
| |
| \chi_{KS}(\mathbf{r_1}t_1,\mathbf{r}_2't_2')
| |
| \left(\frac{1}{|\mathbf{r}_2'-\mathbf{r}_1'|}+f_{xc}(\mathbf{r}_2't_2',\mathbf{r}_1't_1')\right)
| |
| \chi(\mathbf{r}_1't_1',\mathbf{r}_2t_2) </math>
| |
| | |
| From this last equation it is possible to derive the excitation energies of the system, as these are simply the poles of the response function.
| |
| | |
| Other linear-response approaches include the Casida formalism (an expansion in electron-hole pairs) and the Sternheimer equation (density-functional perturbation theory).
| |
| | |
| == Key papers ==
| |
| * [http://dx.doi.org/10.1103/PhysRev.136.B864 P. Hohenberg and W. Kohn, Phys. Rev. 136 (1964) B864]
| |
| * [http://dx.doi.org/10.1103/PhysRevLett.52.997 E. Runge and E.K.U. Gross, Phys. Rev. Lett. 52 (1984) 997]
| |
| | |
| == Books on TDDFT ==
| |
| | |
| * M.A.L. Marques, C.A. Ullrich, F. Nogueira, A. Rubio, K. Burke, and E.K.U. Gross (eds.), Time-Dependent Density Functional Theory (Springer-Verlag, 2006). ISBN 978-3-540-35422-2
| |
| | |
| ==TDDFT codes==
| |
| * [http://classic.chem.msu.su/gran/gamess/index.html Firefly]
| |
| * [http://www.msg.ameslab.gov/gamess/ GAMESS-US]
| |
| * [[Gaussian (software)|Gaussian]]
| |
| * [[Amsterdam Density Functional]]
| |
| * [[CP2K]]
| |
| * [[Dalton (program)|Dalton]]
| |
| * [[NWChem]]
| |
| * [[octopus (software)|Octopus]]
| |
| * [[PARSEC]]
| |
| * [[Q-Chem]]
| |
| * [[Spartan (software)|Spartan]]
| |
| * [[TeraChem]]
| |
| * [[TURBOMOLE]]
| |
| * [[YAMBO code]]
| |
| * [http://www.thch.uni-bonn.de/tc/orca/ ORCA]
| |
| * [[Jaguar (software)|Jaguar]]
| |
| * [https://wiki.fysik.dtu.dk/gpaw/ GPAW]
| |
| | |
| ==References==
| |
| {{reflist}}
| |
| | |
| == External links ==
| |
| * [http://tddft.org tddft.org]
| |
| * [http://th.physik.uni-frankfurt.de/~engel/tddft.html Brief introduction of TD-DFT]
| |
| | |
| [[Category:Density functional theory]]
| |
| [[Category:Computational chemistry]]
| |
| [[Category:Computational physics]]
| |
| [[Category:Quantum chemistry]]
| |
| [[Category:Theoretical chemistry]]
| |