No Arabic abstract
An explicit expression for the quadratic density-response function of a many-electron system is obtained in the framework of the time-dependent density-functional theory, in terms of the linear and quadratic density-response functions of noninteracting Kohn-Sham electrons and functional derivatives of the time-dependent exchange-correlation potential. This is used to evaluate the quadratic stopping power of a homogeneous electron gas for slow ions, which is demonstrated to be equivalent to that obtained up to second order in the ion charge in the framework of a fully nonlinear scattering approach. Numerical calculations are reported, thereby exploring the range of validity of quadratic-response theory.
Linear-response time-dependent density-functional theory (TDDFT) can describe excitonic features in the optical spectra of insulators and semiconductors, using exchange-correlation (xc) kernels behaving as $-1/k^{2}$ to leading order. We show how excitons can be modeled in real-time TDDFT, using an xc vector potential constructed from approximate, long-range corrected xc kernels. We demonstrate for various materials that this real-time approach is consistent with frequency-dependent linear response, gives access to femtosecond exciton dynamics following short-pulse excitations, and can be extended with some caution into the nonlinear regime.
We introduce a new implementation of time-dependent density-functional theory which allows the emph{entire} spectrum of a molecule or extended system to be computed with a numerical effort comparable to that of a emph{single} standard ground-state calculation. This method is particularly well suited for large systems and/or large basis sets, such as plane waves or real-space grids. By using a super-operator formulation of linearized time-dependent density-functional theory, we first represent the dynamical polarizability of an interacting-electron system as an off-diagonal matrix element of the resolvent of the Liouvillian super-operator. One-electron operators and density matrices are treated using a representation borrowed from time-independent density-functional perturbation theory, which permits to avoid the calculation of unoccupied Kohn-Sham orbitals. The resolvent of the Liouvillian is evaluated through a newly developed algorithm based on the non-symmetric Lanczos method. Each step of the Lanczos recursion essentially requires twice as many operations as a single step of the iterative diagonalization of the unperturbed Kohn-Sham Hamiltonian. Suitable extrapolation of the Lanczos coefficients allows for a dramatic reduction of the number of Lanczos steps necessary to obtain well converged spectra, bringing such number down to hundreds (or a few thousands, at worst) in typical plane-wave pseudopotential applications. The resulting numerical workload is only a few times larger than that needed by a ground-state Kohn-Sham calculation for a same system. Our method is demonstrated with the calculation of the spectra of benzene, C$_{60}$ fullerene, and of chlorofyll a.
We analyze possible nonlinear exciton-exciton correlation effects in the optical response of semiconductors by using a time-dependent density-functional theory (TDDFT) approach. For this purpose, we derive the nonlinear (third-order) TDDFT equation for the excitonic polarization. In this equation, the nonlinear time-dependent effects are described by the time-dependent (non-adiabatic) part of the effective exciton-exciton interaction, which depends on the exchange-correlation (XC) kernel. We apply the approach to study the nonlinear optical response of a GaAs quantum well. In particular, we calculate the 2D Fourier spectra of the system and compare it with experimental data. We find that it is necessary to use a non-adiabatic XC kernel to describe excitonic bound states - biexcitons, which are formed due to the retarded TDDFT exciton-exciton interaction.
Using a super-operator formulation of linearized time-dependent density-functional theory, the dynamical polarizability of a system of interacting electrons is given a matrix continued-fraction representation whose coefficients can be obtained from the non-symmetric block-Lanczos method. The resulting algorithm allows for the calculation of the {em full spectrum} of a system with a computational workload which is only a few times larger than that needed for {em static} polarizabilities within time-independent density-functional perturbation theory. The method is demonstrated with the calculation of the spectrum of benzene, and prospects for its application to the large-scale calculation of optical spectra are discussed.
Time-dependent density-functional theory (TDDFT) is a computationally efficient first-principles approach for calculating optical spectra in insulators and semiconductors, including excitonic effects. We show how exciton wave functions can be obtained from TDDFT via the Kohn-Sham transition density matrix, both in the frequency-dependent linear-response regime and in real-time propagation. The method is illustrated using one-dimensional model solids. In particular, we show that our approach provides insight into the formation and dissociation of excitons in real time. This opens the door to time-resolved studies of exciton dynamics in materials by means of real-time TDDFT.