Downloads & Free Reading Options - Results

The real-time-propagation formulation of time-dependent density-functional theory (RT-TDDFT) is an efficient method for modeling the optical response of molecules and nanoparticles. Compared to the widely adopted linear-response TDDFT approaches based on, e.g., the Casida equations, RT-TDDFT appears, however, lacking efficient analysis methods. This applies in particular to a decomposition of the response in the basis of the underlying single-electron states. In this work, we overcome this limitation by developing an analysis method for obtaining the Kohn-Sham electron-hole decomposition in RT-TDDFT. We demonstrate the equivalence between the developed method and the Casida approach by a benchmark on small benzene derivatives. Then, we use the method for analyzing the plasmonic response of icosahedral silver nanoparticles up to Ag$_{561}$. Based on the analysis, we conclude that in small nanoparticles individual single-electron transitions can split the plasmon into multiple resonances due to strong single-electron-plasmon coupling whereas in larger nanoparticles a distinct plasmon resonance is formed.

In this paper we employ the hydrodynamic formulation of time-dependent density functional theory to obtain the van der Waal coefficients $C_{6}$ and $C_{8}$ of alkali-metal clusters of various sizes including very large clusters. Such calculation becomes computationally very demanding in the orbital-based Kohn-Sham formalism, but quite simple in the hydrodynamic approach. We show that for interactions between the clusters of same sizes, $C_{6}$ and $C_{8}$ sale as the sixth and the eighth power of the cluster radius rsepectively, and approach the respective classically predicted values for the large size clusters.

Time-dependent density-functional response theory for molecules. Lecture Notes Collection FreeScience.info ID2122 Obtained from http://tddft.org/TDDFT2004/PracticalSessions/papers/casida_1996.ps http://www.freescience.info/go.php?pagename=books&id=2122

Time-dependent density functional theory is widely used to describe excitations of many-fermion systems. In its many applications, 3D coordinate-space representation is used, and infinite-domain calculations are limited to a finite volume represented by a box. For finite quantum systems (atoms, molecules, nuclei), the commonly used periodic or reflecting boundary conditions introduce spurious quantization of the continuum states and artificial reflections from boundary; hence, an incorrect treatment of evaporated particles. These artifacts can be practically cured by introducing absorbing boundary conditions (ABC) through an absorbing potential in a certain boundary region sufficiently far from the described system. But also the calculations of infinite matter (crystal electrons, quantum fluids, neutron star crust) suffer artifacts from a finite computational box. In this regime, twist- averaged boundary conditions (TABC) have been used successfully to diminish the finite-volume effects. In this work, we extend TABC to time-dependent framework and apply it to resolve the box artifacts for finite quantum systems using as test case small- and large-amplitude nuclear vibrations. We demonstrate that by using such a method, one can reduce finite volume effects drastically without adding any additional parameters. While they are almost equivalent in the linear regime, TABC and ABC differ in the nonlinear regime in their treatment of evaporated particles.

An efficient algorithm for time propagation of the time-dependent Kohn-Sham equations is presented. The algorithm is based on dividing the Hamiltonian into small time steps and assuming that it is constant over these steps. This allows for the time-propagating Kohn-Sham wave function to be expanded in the instantaneous eigenstates of the Hamiltonian. The stability and efficiency of the algorithm are tested not just for non-magnetic but also for fully non-collinear magnetic systems. We show that even for delicate properties, like magnetization density, large time-step sizes can be used indicating the stability and efficiency of the algorithm.

Adiabatic approximations in time-dependent density functional theory (TDDFT) will in general yield unphysical time-dependent shifts in the resonance positions of a system driven far from its ground-state. This spurious time-dependence is rationalized in [J. I. Fuks, K. Luo, E. D. Sandoval and N. T. Maitra, Phys. Rev. Lett. {\bf 114}, 183002 (2015)] in terms of the violation of an exact condition by the non-equilibrium exchange-correlation kernel of TDDFT. Here we give details on the derivation and discuss reformulations of the exact condition that apply in special cases. In its most general form, the condition states that when a system is left in an arbitrary state, in the absence of time-dependent external fields nor ionic motion, the TDDFT resonance position for a given transition is independent of the state. Special cases include the invariance of TDDFT resonances computed with respect to any reference interacting stationary state of a fixed potential, and with respect to any choice of appropriate stationary Kohn-Sham reference state. We then present several case studies, including one that utilizes the adiabatically-exact approximation, that illustrate the conditions and the impact of their violation on the accuracy of the ensuing dynamics. In particular, charge-transfer across a long-range molecule is hampered, and we show how adjusting the frequency of a driving field to match the time-dependent shift in the charge-transfer resonance frequency, results in a larger charge transfer over time.

We study time-dependent electron transport through an Anderson model. The electronic interactions on the impurity site are included via the self-energy approximations at Hartree-Fock (HF), second Born (2B), GW, and T-Matrix level as well as within a time-dependent density functional (TDDFT) scheme based on the adiabatic Bethe-Ansatz local density approximation (ABALDA) for the exchange correlation potential. The Anderson model is driven out of equilibrium by applying a bias to the leads and its nonequilibrium dynamics is determined by real-time propagation. The time-dependent currents and densities are compared to benchmark results obtained with the time-dependent density matrix renormalization group (tDMRG) method. Many-body perturbation theory beyond HF gives results in close agreement with tDMRG especially within the 2B approximation. We find that the TDDFT approach with the ABALDA approximation produces accurate results for the densities on the impurity site but overestimates the currents. This problem is found to have its origin in an overestimation of the lead densities which indicates that the exchange correlation potential must attain nonzero values in the leads.

A recent paper (Phys. Rev A. 75, 022513 (2007), arXiv:cond-mat/0602020) challenges exact time-dependent density functional theory (TDDFT) on several grounds. We explain why these criticisms are either irrelevant or incorrect, and that TDDFT is both formally exact and predictive.

Time-Dependent Density Functional Theory (TDDFT) has recently been extended to describe many-body open quantum systems (OQS) evolving under non-unitary dynamics according to a quantum master equation. In the master equation approach, electronic excitation spectra are broadened and shifted due to relaxation and dephasing of the electronic degrees of freedom by the surrounding environment. In this paper, we develop a formulation of TDDFT linear-response theory (LR-TDDFT) for many-body electronic systems evolving under a master equation, yielding broadened excitation spectra. This is done by mapping an interacting open quantum system onto a non-interacting open Kohn-Sham system yielding the correct non-equilibrium density evolution. A pseudo-eigenvalue equation analogous to the Casida equations of usual LR-TDDFT is derived for the Redfield master equation, yielding complex energies and Lamb shifts. As a simple demonstration, we calculate the spectrum of a C$^{2+}$ atom in an optical resonator interacting with a bath of photons. The performance of an adiabatic exchange-correlation kernel is analyzed and a first-order frequency-dependent correction to the bare Kohn-Sham linewidth based on Gorling-Levy perturbation theory is calculated.

Based on our earlier works [Phys. Rev. B 75, 195127 (2007) & J. Chem. Phys. 128, 234703 (2008)], we propose a formally exact and numerically convenient approach to simulate time-dependent quantum transport from first-principles. The proposed approach combines time-dependent density functional theory with quantum dissipation theory, and results in a useful tool for studying transient dynamics of electronic systems. Within the proposed exact theoretical framework, we construct a number of practical schemes for simulating realistic systems such as nanoscopic electronic devices. Computational cost of each scheme is analyzed, with the expected level of accuracy discussed. As a demonstration, a simulation based on the adiabatic wide-band limit approximation scheme is carried out to characterize the transient current response of a carbon nanotube based electronic device under time-dependent external voltages.

The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived starting from time-dependent density functional theory. Effective permittivity and permeability coefficients are obtained.

We explore an asymmetric two-fermion Hubbard dimer to test the accuracy of the adiabatic approximation of time-dependent density functional theory in modelling time-resolved charge transfer. We show that the model shares essential features of a ground state long-range molecule in real-space, and by applying a resonant field we show that the model also reproduces essential traits of the CT dynamics. The simplicity of the model allows us to propagate with an "adiabatically-exact" approximation, i.e. one that uses the exact ground-state exchange-correlation functional, and compare with the exact propagation. This allows us to study the impact of the time-dependent charge-transfer step feature in the exact correlation potential of real molecules on the resulting dynamics. Tuning the parameters of the dimer allows a study both of charge-transfer between open-shell fragments and between closed-shell fragments. We find that the adiabatically-exact functional is unable to properly transfer charge, even in situations where the adiabatically-exact resonance frequency is remarkably close to the exact resonance, and we analyze why.

The density-functional approach to quantum electrodynamics is extending traditional density-functional theory and opens the possibility to describe electron-photon interactions in terms of effective Kohn-Sham potentials. In this work, we numerically construct the exact electron-photon Kohn-Sham potentials for a prototype system which consists of a trapped electron coupled to a quantized electromagnetic mode in an optical high-Q cavity. While the effective current that acts on the photons is known explicitly, the exact effective potential that describes the forces exerted by the photons on the electrons is obtained from a fixed-point inversion scheme. This procedure allows us to uncover important beyond-mean-field features of the effective potential which mark the breakdown of classical light-matter interactions. We observe peak and step structures in the effective potentials, which can be attributed solely to the quantum nature of light, i.e., they are real-space signatures of the photons. Our findings show how the ubiquitous dipole interaction with a classical electromagnetic field has to be modified in real-space in order to take the quantum nature of the electromagnetic field fully into account.

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.

The photoionization cross-section of several atoms (Ar, Xe, Rn, Cs) and ions (Ne-like Ar, H-like and Li-like C) of experimental interest are calculated using the time-dependent response scheme within the framework of local density functional method (DFM). The cross-sections for rare gas atoms calculated using this method agree very well with the experimental data; whereas conventional independent particle model, calculations do not. The polarization effect of the atom brought about by the incident time-varying radiation field is shown to be important in describing observed results. Unlike the independent particle model, this effect effect is treated adequately in the DFM. To study the effect plasma density and temperature on the photoionization cross-section, calculations were also done for the ions mentioned above at various densities and temperatures. For computational simplicity, a simplified model of self-consistent finite temperature DFM was used in which the long range part of the ionic potential was taken as the Debye-screened potential. These calculations were compared with the isolated ion calculations (without any plasma effect). With increasing plasma density, significant shifts of the ionization threshold as well as substantial modifications of photoionization cross-section are obtained. This points out the need for incorporating the effect of surrounding plasma in realistic modeling of atomic properties for dense plasmas.

In this paper the review progress in the density functional theory of time dependent phenomena. Contents: Introduction; fundamentals; exact conditions; linear response; approximated functionals; applications. Lecture Notes Collection FreeScience.info ID1450 Obtained from http://www.physik.fu-berlin.de/~ag-gross/articles/pdf/BG98.pdf http://www.freescience.info/go.php?pagename=books&id=1450

We present a rigorous formulation of the time-dependent density functional theory for interacting lattice electrons strongly coupled to cavity photons. We start with an example of one particle on a Hubbard dimer coupled to a single photonic mode, which is equivalent to the single mode spin-boson model or the quantum Rabi model. For this system we prove that the electron-photon wave function is a unique functional of the electronic density and the expectation value of the photonic coordinate, provided the initial state and the density satisfy a set of well defined conditions. Then we generalize the formalism to many interacting electrons on a lattice coupled to multiple photonic modes and prove the general mapping theorem. We also show that for a system evolving from the ground state of a lattice Hamiltonian any density with a continuous second time derivative is locally $v$-representable.

The application of methods of time-dependent density functional theory (TDDFT) to systems of qubits provided the interesting possibility of simulating an assigned Hamiltonian evolution by means of an auxiliary Hamiltonian having different two-qubit interactions and hence a possibly simpler wave function evolution. In this note we extend these methods to some instances of Lindblad evolution of a spin chain.

The scalar $f_{xc}$ and tensor $\hat{f}_{xc}$ exchange-correlation (xc) kernels are key ingredients of the time-dependent density functional theory and the time-dependent current density functional theory, respectively. We derive a comparatively simple relation between these two kernels under the assumption that the dynamic xc can be considered "weak". We expect our formula to serve as a convenient bridge between the scalar $f_{xc}$ which directly enters many applications and the tensor $\hat{f}_{xc}$ which, due to its locality in space, is much easier to approximate.

Linear response calculations based on the time-dependent density-functional theory are presented. Especially, we report results of the finite amplitude method which we have recently proposed as an alternative and feasible approach to the (quasiparticle-)random-phase approximation. Calculated properties of the giant resonances and low-energy E1 modes are discussed. We found a universal linear correlation between the low-energy E1 strength and the neutron skin thickness.

Recent advances in laser technology allow us to follow electronic motion at its natural time-scale with ultra-fast time resolution, leading the way towards attosecond physics experiments of extreme precision. In this work, we assess the use of tailored pumps in order to enhance (or reduce) some given features of the probe absorption (for example, absorption in the visible range of otherwise transparent samples). This type of manipulation of the system response could be helpful for its full characterization, since it would allow us to visualize transitions that are dark when using unshaped pulses. In order to investigate these possibilities, we perform first a theoretical analysis of the non-equilibrium response function in this context, aided by one simple numerical model of the Hydrogen atom. Then, we proceed to investigate the feasibility of using time-dependent density-functional theory as a means to implement, theoretically, this absorption-optimization idea, for more complex atoms or molecules. We conclude that the proposed idea could in principle be brought to the laboratory: tailored pump pulses can excite systems into light-absorbing states. However, we also highlight the severe numerical and theoretical difficulties posed by the problem: large-scale non-equilibrium quantum dynamics are cumbersome, even with TDDFT, and the shortcomings of state-of-the-art TDDFT functionals may still be serious for these out-of-equilibrium situations.

Motivated by the large interest in the non-equilibrium dynamics of low-dimensional quantum many-body systems, we present a fully-microscopic theoretical and numerical study of the "charge" and "spin" dynamics in a one-dimensional ultracold Fermi gas following a quench. Our approach, which is based on time-dependent current-density-functional theory, is applicable well beyond the linear-response regime and produces both spin-charge separation and spin-drag-induced broadening of the spin packets.

We propose a new method of calculating electronically excited states that combines a density functional theory (DFT) based ground state calculation with a linear response treatment that employs approximations used in the time-dependent density functional based tight binding (TD-DFTB) approach. The new method termed TD-DFT+TB does not rely on the DFTB parametrization and is therefore applicable to systems involving all combinations of elements. We show that the new method yields UV/Vis absorption spectra that are in excellent agreement with computationally much more expensive time-dependent density functional theory (TD-DFT) calculations. Errors in vertical excitation energies are reduced by a factor of two compared to TD-DFTB.

We formulate equations of time-dependent density functional theory (TDDFT) in the co-moving Lagrangian reference frame. The main advantage of the Lagrangian description of many-body dynamics is that in the co-moving frame the current density vanishes, while the density of particles becomes independent of time. Therefore a co-moving observer will see the picture which is very similar to that seen in the equilibrium system from the laboratory frame. It is shown that the most natural set of basic variables in TDDFT includes the Lagrangian coordinate, $\bm\xi$, a symmetric deformation tensor $g_{\mu\nu}$, and a skew-symmetric vorticity tensor, $F_{\mu\nu}$. These three quantities, respectively, describe the translation, deformation, and the rotation of an infinitesimal fluid element. Reformulation of TDDFT in terms of new basic variables resolves the problem of nonlocality and thus allows to regularly derive a local nonadiabatic approximation for exchange correlation (xc) potential. Stationarity of the density in the co-moving frame makes the derivation to a large extent similar to the derivation of the standard static local density approximation. We present a few explicit examples of nonlinear nonadiabatic xc functionals in a form convenient for practical applications.

Floquet formulation of time-dependent density-functional theory is revisited in light of its recent criticism [Maitra and Burke, Chem. Phys. Lett. 359 (2002), 237]. It is shown that Floquet theory is well founded and its criticism has overlooked important points of both the Runge-Gross formalism and Floquet formulation itself. We substantiate our analysis by examples similar to that considered by Maitra and Burke.

We investigate the description of excitonic effects within time-dependent density-functional theory (TDDFT). The exchange-correlation kernel f_xc introduced in TDDFT allows a clear separation of quasiparticle and excitonic effects. Using a diagrammatic representation for f_xc, we express its excitonic part f_xc^Ex in terms of the effective vertex function Lambda. The latter fulfills an integral equation which thereby establishes the exact correspondence between TDDFT and the standard many-body approach based on Bethe-Salpeter equation (BSE).The diagrammatic structure of the kernel in the equation for Lambda suggests the possibility of strong cancellation effects. Should the cancellation take place, already the first-order approximation to f_xc^Ex is sufficient. A potential advantage of TDDFT over the many-body BSE method is thus dependent on the efficiency of the above-quoted cancellation. We explicitly verify this for an analytically solvable two-dimensional two-band model. The calculations confirm that the low-order f_xc^Ex perfectly describes the bound exciton as well as the excitonic effects in the continuous spectrum in a wide range of the electron--hole coupling strength.

We report on a general method for the calculation of the frequency-dependent optical response of clusters based upon time-dependent density functional theory (TDDFT). The implementation is done using explicit propagation in the time domain and a self-consistent program that uses a linear combination of atomic orbitals (LCAO). Our actual calculations employ the SIESTA program, which is designed to be fast and accurate for large clusters. We use the adiabatic local density approximation to account for exchange and correlation effects. Results are presented for the imaginary part of the linear polarizability, Im [\alpha(w)], and the dipole strength function, S(w), of C60 and Na8, compared to previous calculations and to experiment. We also show how to calculate the integrated frequency-dependent second order non-linear polarizability for the case of a step function electric field, \gamma_{step}(w), and present results for C60.

Time-dependent density-functional theory (TDDFT) treats dynamical exchange and correlation (xc) via a single-particle potential, Vxc(r,t), defined as a nonlocal functional of the density n(r',t'). The popular adiabatic local-density approximation (ALDA) for Vxc(r,t) uses only densities at the same space-time point (r,t). To go beyond the ALDA, two local approximations have been proposed based on quantum hydrodynamics and elasticity theory: (a) using the current as basic variable (C-TDDFT) [G. Vignale, C. A. Ullrich, and S. Conti, Phys. Rev. Lett. 79, 4878 (1997)], (b) working in a co-moving Lagrangian reference frame (L-TDDFT) [I. V. Tokatly, Phys. Rev. B 71, 165105 (2005)]. This paper illustrates, compares, and analyzes both non-adiabatic theories for simple time-dependent model densities in the linear and nonlinear regime, for a broad range of time and frequency scales. C- and L-TDDFT are identical in certain limits, but in general exhibit qualitative and quantitative differences in their respective treatment of elastic and dissipative electron dynamics. In situations where the electronic density rapidly undergoes large deformations, it is found that non-adiabatic effects can become significant, causing the ALDA to break down.

The exact exchange potential in time-dependent density-functional theory is defined as an orbital functional through the time-dependent optimized effective potential (TDOEP) method. We numerically solve the TDOEP integral equation for the real-time nonlinear intersubband electron dynamics in a semiconductor quantum well with two occupied subbands. By comparison with adiabatic approximations, it is found that memory effects in the exact exchange potential become significant when the electron dynamics takes place in the vicinity of intersubband resonances.

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 o 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 o 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. Lecture Notes Collection FreeScience.info ID2183 Obtained from http://www.fhi-berlin.mpg.de/th/publications/RPP-70-357-2007.pdf http://www.freescience.info/go.php?pagename=books&id=2183

We identify the key property that the exchange-correlation (XC) kernel of time-dependent density functional theory must have in order to describe long-range charge-transfer excitations. We show that the discontinuity of the XC potential as a function of particle number induces a space -and frequency-dependent discontinuity of the XC kernel which diverges as $r\to\infty$. In a combined donor-acceptor system, the same discontinuity compensates for the vanishing overlap between the acceptor and donor orbitals, thereby yielding a finite correction to the Kohn-Sham eigenvalue differences. This mechanism is illustrated to first order in the Coulomb interaction.

We apply the coupled dynamics of time-dependent density functional theory and Maxwell equations to the interaction of intense laser pulses with crystalline silicon. As a function of electromagnetic field intensity, we see several regions in the response. At the lowest intensities, the pulse is reflected and transmitted in accord with the dielectric response, and the characteristics of the energy deposition is consistent with two-photon absorption. The absorption process begins to deviate from that at laser intensities ~ 10^13 W/cm^2, where the energy deposited is of the order of 1 eV per atom. Changes in the reflectivity are seen as a function of intensity. When it passes a threshold of about 3 \times 1012 W/cm2, there is a small decrease. At higher intensities, above 2 \times 10^13 W/cm^2, the reflectivity increases strongly. This behavior can be understood qualitatively in a model treating the excited electron-hole pairs as a plasma.

I describe the Time-Dependent Superfluid Local Density Approximation, which is an adiabatic extension of the Density Functional Theory to superfluid Fermi systems and their real-time dynamics. This new theoretical framework has been applied to describe a number of phenomena in cold atomic gases and nuclear collective motion: excitation of the Higgs modes in strongly interacting Fermi superfluids, generation of quantized vortices, crossing and reconnection of vortices, excitation of the superflow at velocities above the critical velocity, excitation of quantum shock waves and domain walls in the collisions of superfluid atomic clouds, excitation of collective states in nuclei.

The Runge-Gross [E. Runge, and E. K. U. Gross, Phys. Rev. Lett., 52, 997 (1984)] action functional of time-dependent density-functional theory leads to a well-known causality paradox, i.e., a perturbation of the electronic density in the future affects the response of the system in the present. This paradox is known to be caused by an inconsistent application of the Dirac-Frenkel variational principle. In view of the recent solutions to this problem, the action functional employed by Runge and Gross in their formulation of time-dependent density functional theory is analyzed in the context of the Keldysh contour technique. The time-dependent electronic density, as well as the concept of causality, are extended to the contour. We derive a variational equation that obeys causality and relates the exchange-correlation potential with its kernel, and the functional derivative of the exchange-correlation action functional with respect to the density. It is shown that the adiabatic local-density approximation is a consistent solution of this equation and that the time-dependent optimized potential method can also be derived from it. The formalism presented here can be used to find new approximations methods to the exchange-correlation potential and avoid the causality dilemma.

We discuss an extension of time-dependent density-functional theory by a self-interaction correction (SIC). A strictly variational formulation is given taking care of the necessary constraints. A manageable and transparent propagation scheme using two sets of wavefunctions is proposed and applied to laser excitation with subsequent ionization of a dimer molecule.

To construct an N-representable time-dependent density-functional theory, a generalization to the time domain of the Levy-Lieb (LL) constrained search algorithm is required. That the action is only stationary in the Dirac-Frenkel variational principle eliminates the possibility of basing the search on the action itself. Instead, we use the norm of the partial functional derivative of the action in the Hilbert space of the wave functions in place of the energy of the LL search. The electron densities entering the formalism are $N$-representable, and the resulting universal action functional has a unique stationary point in the density at that corresponding to the solution of the Schr\"{o}dinger equation. The original Runge-Gross (RG) formulation is subsumed within the new formalism. Concerns in the literature about the meaning of the functional derivatives and the internal consistency of the RG formulation are allayed by clarifying the nature of the functional derivatives entering the formalism.

We demonstrate the capabilities of time-dependent density functional theory (TDDFT) for strong-field, short wavelength (soft X-ray) physics, as compared to a formalism based on rate equations. We find that TDDFT provides a very good description of the total and individual ionization yields for Ne and Ar atoms exposed to strong laser pulses. We assess the reliability of different adiabatic density functionals and conclude that an accurate description of long-range interactions by the exchange and correlation potential is crucial for obtaining the correct ionization yield over a wide range of intensities ($10^{13}$ -- $5 \times 10^{15}$ W/cm$^2$). Our TDDFT calculations disentangle the contribution from each ionization channel based on the Kohn-Sham wavefunctions.

We have employed non-adiabatic molecular dynamics based on time-dependent density-functional theory to characterize the scattering behaviour of a proton with the Li$_4$ cluster. This technique assumes a classical approximation for the nuclei, effectively coupled to the quantum electronic system. This time-dependent theoretical framework accounts, by construction, for possible charge transfer and ionization processes, as well as electronic excitations, which may play a role in the non-adiabatic regime. We have varied the incidence angles in order to analyze the possible reaction patterns. The initial proton kinetic energy of 10 eV is sufficiently high to induce non-adiabatic effects. For all the incidence angles considered the proton is scattered away, except in one interesting case in which one of the Lithium atoms captures it, forming a LiH molecule. This theoretical formalism proves to be a powerful, effective and predictive tool for the analysis of non-adiabatic processes at the nanoscale.

This article has been published as a chapter in "Chemical Reactivity Theory: A Density Functional View", ed. P. K. Chattaraj (CRC Press, New York, 2009), ch. 8, p. 105. In it, an overview of the relationship between time-dependent DFT and quantum hydrodynamics is presented, showing the role that Bohmian mechanics can play within the ab-initio methodology as both a numerical and an interpretative tool.

The resonant interaction of laser light with atoms is analyzed from the time-dependent density functional theory perspective using a model Helium atom which can be solved exactly. It is found that in exact-exchange approximation the time-dependent dipole shows Rabi-type oscillations of its amplitude. However, the time-dependent density itself is not well described. These seemingly contradictory findings are analyzed. The Rabi-type oscillations are found to be essentially of classical origin. The incompatibility of time-dependent density functional theory with few-level approximations for the description of resonant dynamics is discussed.

Quantum Optimal Control Theory (QOCT) provides the necessary tools to theoretically design driving fields capable of controlling a quantum system towards a given state or along a prescribed path in Hilbert space. This theory must be complemented with a suitable model for describing the dynamics of the quantum system. Here, we are concerned with many electron systems (atoms, molecules, quantum dots, etc) irradiated with laser pulses. The full solution of the many electron Schr{\"{o}}dinger equation is not feasible in general, and therefore, if we aim to an ab initio description, a suitable choice is time-dependent density-functional theory (TDDFT). In this work, we establish the equations that combine TDDFT with QOCT, and demonstrate their numerical feasibility with examples.

We reformulate and generalize the uniqueness and existence proofs of time-dependent density-functional theory. The central idea is to restate the fundamental one-to-one correspondence between densities and potentials as a global fixed point question for potentials on a given time-interval. We show that the unique fixed point, i.e. the unique potential generating a given density, is reached as the limiting point of an iterative procedure. The one-to-one correspondence between densities and potentials is a straightforward result provided that the response function of the divergence of the internal forces is bounded. The existence, i.e. the v-representability of a density, can be proven as well provided that the operator norms of the response functions of the members of the iterative sequence of potentials have an upper bound. The densities under consideration have second time-derivatives that are required to satisfy a condition slightly weaker than being square-integrable. This approach avoids the usual restrictions of Taylor-expandability in time of the uniqueness theorem by Runge and Gross [Phys.Rev.Lett.52, 997 (1984)] and of the existence theorem by van Leeuwen [Phys.Rev.Lett. 82, 3863 (1999)]. Owing to its generality, the proof not only answers basic questions in density-functional theory but also has potential implications in other fields of physics.

We introduce a new simplified method for computing the electron field emission current in short carbon nanotubes using ab-initio computation in periodic simulation cells. We computed the evolution of the wave functions using Time-Dependent Density Functional Theory, where we have utilized the Crank-Nicholson propagator. We found that in pristine carbon nanotubes, the emitted charge tends to emerge mostly from electrons that are concentrated at the nanotube tip region. The charge beam concentrates into specific channel structures, showing the utility of carbon nanotubes in precision emission applications.

Excitons are electron-hole pairs appearing below the band gap in insulators and semiconductors. They are vital to photovoltaics, but are hard to obtain with time-dependent density-functional theory (TDDFT), since most standard exchange-correlation (xc) functionals lack the proper long-range behavior. Furthermore, optical spectra of bulk solids calculated with TDDFT often lack the required resolution to distinguish discrete, weakly bound excitons from the continuum. We adapt the Casida equation formalism for molecular excitations to periodic solids, which allows us to obtain exciton binding energies directly. We calculate exciton binding energies for both small- and large-gap semiconductors and insulators, study the recently proposed bootstrap xc kernel [S. Sharma et al., Phys. Rev. Lett. 107, 186401 (2011)], and extend the formalism to triplet excitons.

Excitons are electron-hole pairs appearing below the band gap in insulators and semiconductors. They are vital to photovoltaics, but are hard to obtain with time-dependent density-functional theory (TDDFT), since most standard exchange-correlation (xc) functionals lack the proper long-range behavior. Furthermore, optical spectra of bulk solids calculated with TDDFT often lack the required resolution to distinguish discrete, weakly bound excitons from the continuum. We adapt the Casida equation formalism for molecular excitations to periodic solids, which allows us to obtain exciton binding energies directly. We calculate exciton binding energies for both small- and large-gap semiconductors and insulators, study the recently proposed bootstrap xc kernel [S. Sharma et al., Phys. Rev. Lett. 107, 186401 (2011)], and extend the formalism to triplet excitons.

Most applications of time-dependent density-functional theory (TDDFT) use the adiabatic local-density approximation (ALDA) for the dynamical exchange-correlation potential Vxc(r,t). An exact (i.e., nonadiabatic) extension of the ground-state LDA into the dynamical regime leads to a Vxc(r,t) with a memory, which causes the electron dynamics to become dissipative. To illustrate and explain this nonadiabatic behavior, this paper studies the dynamics of two interacting electrons on a two-dimensional quantum strip of finite size, comparing TDDFT within and beyond the ALDA with numerical solutions of the two-electron time-dependent Schroedinger equation. It is shown explicitly how dissipation arises through multiple particle-hole excitations, and how the nonadiabatic extension of the ALDA fails for finite systems, but becomes correct in the thermodynamic limit.

In this work a practical scheme is developed for the first-principles study of time-dependent quantum transport. The basic idea is to combine the transport master-equation with the well-known time-dependent density functional theory. The key ingredients of this paper include: (i) the partitioning-free initial condition and the consideration of the time-dependent bias voltages which base our treatment on the Runge-Gross existence theorem; (ii) the non-Markovian master equation for the reduced (many-body) central system (i.e. the device); and (iii) the construction of Kohn-Sham master equation for the reduced single-particle density matrix, where a number of auxiliary functions are introduced and their equations of motion (EOM) are established based on the technique of spectral decomposition. As a result, starting with a well-defined initial state, the time-dependent transport current can be calculated simultaneously along the propagation of the Kohn-Sham master equation and the EOM of the auxiliary functions.

We show that the time-dependent particle density $n(\vec r,t)$ and the current density ${\vec j}(\vec r,t)$ of a many-particle system that evolves under the action of external scalar and vector potentials $V(\vec r,t)$ and $\vec A(\vec r,t)$ and is initially in the quantum state $|\psi (0)>$, can always be reproduced (under mild assumptions) in another many-particle system, with different two-particle interaction, subjected to external potentials $V'(\vec r,t)$ and $\vec A'(\vec r,t)$, starting from an initial state $|\psi' (0)>$, which yields the same density and current as $|\psi (0)>$. Given the initial state of this other many-particle system, the potentials $V'(\vec r,t)$ and $\vec A'(\vec r,t)$ are uniquely determined up to gauge transformations that do not alter the initial state. As a special case, we obtain a new and simpler proof of the Runge-Gross theorem for time-dependent current density functional theory. This theorem provides a formal basis for the application of time-dependent current density functional theory to transport problems.

We formulate a time-dependent density-matrix functional theory (TDDMFT) approach for higher-order correlation effects like biexcitons in optical processes in solids based on the reduced two-particle density-matrix formalism within the normal orbital representation. A TDDMFT version of the Schr\"odinger equation for biexcitons in terms of one- and two-body reduced density matrices is derived, which leads to finite biexcitonic binding energies already with an adiabatic approximation. Biexcitonic binding energies for several bulk semiconductors are calculated using a contact biexciton model.

The van Leeuwen proof of linear-response time-dependent density functional theory (TDDFT) is generalized to thermal ensembles. This allows generalization to finite temperatures of the Gross-Kohn relation, the exchange-correlation kernel of TDDFT, and fluctuation dissipation theorem for DFT. This produces a natural method for generating new thermal exchange-correlation (XC) approximations.