E. K. U. Gross

Density-Functional Theory for Time-Dependent Systems

The response of an interacting many-particle system to a time-dependent external field can usually be treated within linear response theory. Due to rapid experimental progress in the field of laser physics, however, ultra-short laser pulses of very high intensity have become available in recent years. The electric field produced in such pulses can reach the strength of the electric field caused by atomic nuclei. If an atomic system is placed in the focus of such a laser pulse one observes a wealth of new phenomena [1] which cannot be explained by traditional perturbation theory. The non-perturbative quantum mechanical description of interacting particles moving in a very strong time-dependent external field therefore has become a prominent problem of theoretical physics. In principle, it requires a full solution of the time-dependent Schrodinger equation for the interacting many-body system, which is an exceedingly difficult task. In view of the success of density functional methods in the treatment of stationary many-body systems and in view of their numerical simplicity, a time-dependent version of density functional theory appears highly desirable, both within and beyond the regime of linear response.

Time-dependent density functional theory: Past, present, and future

Time-dependent density functional theory (TDDFT) is presently enjoying enormous popularity in quantum chemistry, as a useful tool for extracting electronic excited state energies. This article discusses how TDDFT is much broader in scope, and yields predictions for many more properties. We discuss some of the challenges involved in making accurate predictions for these properties.

Fundamentals of time-dependent density functional theory

Part I Theory and Experiment - Why We Need TDDFT.- Part II Basic Theory.- PartIII Advanced Concepts.- Part IV Real-Time Dynamics.- Part V Numerical Aspects.- Part VI TDDFT vs Other Theoretical Techniques.

Reproducibility in density functional theory calculations of solids

A comparison of DFT methods Density functional theory (DFT) is now routinely used for simulating material properties. Many software packages are available, which makes it challenging to know which are the best to use for a specific calculation. Lejaeghere et al. compared the calculated values for the equation of states for 71 elemental crystals from 15 different widely used DFT codes employing 40 different potentials (see the Perspective by Skylaris). Although there were variations in the calculated values, most recent codes and methods converged toward a single value, with errors comparable to those of experiment. Science, this issue p. 10.1126/science.aad3000; see also p. 1394 A survey of recent density functional theory methods shows a convergence to more accurate property calculations. [Also see Perspective by Skylaris] INTRODUCTION The reproducibility of results is one of the underlying principles of science. An observation can only be accepted by the scientific community when it can be confirmed by independent studies. However, reproducibility does not come easily. Recent works have painfully exposed cases where previous conclusions were not upheld. The scrutiny of the scientific community has also turned to research involving computer programs, finding that reproducibility depends more strongly on implementation than commonly thought. These problems are especially relevant for property predictions of crystals and molecules, which hinge on precise computer implementations of the governing equation of quantum physics. RATIONALE This work focuses on density functional theory (DFT), a particularly popular quantum method for both academic and industrial applications. More than 15,000 DFT papers are published each year, and DFT is now increasingly used in an automated fashion to build large databases or apply multiscale techniques with limited human supervision. Therefore, the reproducibility of DFT results underlies the scientific credibility of a substantial fraction of current work in the natural and engineering sciences. A plethora of DFT computer codes are available, many of them differing considerably in their details of implementation, and each yielding a certain “precision” relative to other codes. How is one to decide for more than a few simple cases which code predicts the correct result, and which does not? We devised a procedure to assess the precision of DFT methods and used this to demonstrate reproducibility among many of the most widely used DFT codes. The essential part of this assessment is a pairwise comparison of a wide range of methods with respect to their predictions of the equations of state of the elemental crystals. This effort required the combined expertise of a large group of code developers and expert users. RESULTS We calculated equation-of-state data for four classes of DFT implementations, totaling 40 methods. Most codes agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Even in the case of pseudization approaches, which largely depend on the atomic potentials used, a similar precision can be obtained as when using the full potential. The remaining deviations are due to subtle effects, such as specific numerical implementations or the treatment of relativistic terms. CONCLUSION Our work demonstrates that the precision of DFT implementations can be determined, even in the absence of one absolute reference code. Although this was not the case 5 to 10 years ago, most of the commonly used codes and methods are now found to predict essentially identical results. The established precision of DFT codes not only ensures the reproducibility of DFT predictions but also puts several past and future developments on a firmer footing. Any newly developed methodology can now be tested against the benchmark to verify whether it reaches the same level of precision. New DFT applications can be shown to have used a sufficiently precise method. Moreover, high-precision DFT calculations are essential for developing improvements to DFT methodology, such as new density functionals, which may further increase the predictive power of the simulations. Recent DFT methods yield reproducible results. Whereas older DFT implementations predict different values (red darts), codes have now evolved to mutual agreement (green darts). The scoreboard illustrates the good pairwise agreement of four classes of DFT implementations (horizontal direction) with all-electron results (vertical direction). Each number reflects the average difference between the equations of state for a given pair of methods, with the green-to-red color scheme showing the range from the best to the poorest agreement. The widespread popularity of density functional theory has given rise to an extensive range of dedicated codes for predicting molecular and crystalline properties. However, each code implements the formalism in a different way, raising questions about the reproducibility of such predictions. We report the results of a community-wide effort that compared 15 solid-state codes, using 40 different potentials or basis set types, to assess the quality of the Perdew-Burke-Ernzerhof equations of state for 71 elemental crystals. We conclude that predictions from recent codes and pseudopotentials agree very well, with pairwise differences that are comparable to those between different high-precision experiments. Older methods, however, have less precise agreement. Our benchmark provides a framework for users and developers to document the precision of new applications and methodological improvements.

Time-dependent quantum transport: A practical scheme using density functional theory

We present a computationally tractable scheme of time-dependent transport phenomena within openboundary time-dependent density functional theory. Within this approach all the response properties of a system are determined from the time propagation of the set of “occupied” Kohn-Sham orbitals under the influence of the external bias. This central idea is combined with an open-boundary description of the geometry of the system that is divided into three regions: left/right leads and the device region “real simulation region”. We have derived a general scheme to extract the set of initial states in the device region that will be propagated in time with proper transparent boundary-condition at the device/lead interface. This is possible due to a new modified Crank-Nicholson algorithm that allows an efficient time-propagation of open quantum systems. We illustrate the method in one-dimensional model systems as a first step towards a full first-principles implementation. In particular we show how a stationary current develops in the system independent of the transientcurrent history upon application of the bias. The present work is ideally suited to study ac transport and photon-induced charge-injection. Although the implementation has been done assuming clamped ions, we discuss how it can be extended to include dissipation due to electron-phonon coupling through the combined simulation of the electron-ion dynamics as well as electron-electron correlations.

Exact Coulomb cutoff technique for supercell calculations

Received 23 December 2005; revised manuscript received 31 March 2006; published 26 May 2006 We present a reciprocal space analytical method to cut off the long range interactions in supercell calculations for systems that are infinite and periodic in one or two dimensions, generalizing previous work to treat finite systems. The proposed cutoffs are functions in Fourier space, that are used as a multiplicative factor to screen the bare Coulomb interaction. The functions are analytic everywhere except in a subdomain of the Fourier space that depends on the periodic dimensionality. We show that the divergences that lead to the nonanalytical behavior can be exactly canceled when both the ionic and the Hartree potential are properly screened. This technique is exact, fast, and very easy to implement in already existing supercell codes. To illustrate the performance of the scheme, we apply it to the case of the Coulomb interaction in systems with reduced periodicity as one-dimensional chains and layers. For these test cases, we address the impact of the cutoff on different relevant quantities for ground and excited state properties, namely: the convergence of the ground state properties, the static polarizability of the system, the quasiparticle corrections in the GW scheme, and the binding energy of the excitonic states in the Bethe-Salpeter equation. The results are very promising and easy to implement in all available first-principles codes.

Quantum optimal control theory

The control of quantum dynamics via specially tailored laser pulses is a long-standing goal in physics and chemistry. Partly, this dream has come true, as sophisticated pulse-shaping experiments allow us to coherently control product ratios of chemical reactions. The theoretical design of the laser pulse to transfer an initial state to a given final state can be achieved with the help of quantum optimal control theory (QOCT). This tutorial provides an introduction to QOCT. It shows how the control equations defining such an optimal pulse follow from the variation of a properly defined functional. We explain the most successful schemes to solve these control equations and show how to incorporate additional constraints in the pulse design. The algorithms are then applied to simple quantum systems and the obtained pulses are analysed. Besides the traditional final-time control methods, the tutorial also presents an algorithm and an example to handle time-dependent control targets.

Exact Factorization of the Time-Dependent Electron-Nuclear Wave Function

We present an exact decomposition of the complete wave function for a system of nuclei and electrons evolving in a time-dependent external potential. We derive formally exact equations for the nuclear and electronic wave functions that lead to rigorous definitions of a time-dependent potential energy surface (TDPES) and a time-dependent geometric phase. For the H(2)(+) molecular ion exposed to a laser field, the TDPES proves to be a useful interpretive tool to identify different mechanisms of dissociation.

Excitations in Time-Dependent Density-Functional Theory

An approximate solution to the time-dependent density-functional theory response equations for finite systems is developed, yielding corrections to the single-pole approximation. These explain why allowed Kohn-Sham transition frequencies and oscillator strengths are usually good approximations to the true values, and why sometimes they are not. The approximation yields simple expressions for Görling-Levy perturbation theory results, and a method for estimating expectation values of the unknown exchange-correlation kernel.

Density-functional theory using an optimized exchange-correlation potential

Abstract We have performed self-consistent calculations for first and second row atoms using a variant of density-functional theory, the optimized effective potential method, with an approximation due to Krieger, Li and Iafrate and a correlation-energy functional developed by Colle and Salvetti. The mean absolute deviation of first-row atomic ground-state energies from the exact non-relativistic values is 4.7 mEh in our scheme, as compared to 4.5 mEh in a recent configuration-interaction calculation. The proposed scheme is significantly more accurate than the conventional Kohn-Sham method while the numerical effort involved is about the same as for an ordinary Hartree-Fock calculation.