Do you want to publish a course? Click here

Nonadiabatic time-dependent density-functional theory at the cost of adiabatic local density approximation

183   0   0.0 ( 0 )
 Added by Dmitry R Gulevich
 Publication date 2019
  fields Physics
and research's language is English




Ask ChatGPT about the research

We propose a computationally efficient approach to the nonadiabatic time-dependent density functional theory (TDDFT) which is based on a representation of the frequency-dependent exchange correlation kernel as a response of a set of damped oscillators. The requirements to computational resources needed to implement our approach do not differ from those of the standard real-time TDDFT in the adiabatic local density approximation (ALDA). Thus, our result offers an exciting opportunity to take into account temporal nonlocality and memory effects in calculations with TDDFT in quantum chemistry and solid state physics for unprecedentedly low costs.



rate research

Read More

First-order nonadiabatic coupling matrix elements (fo-NACMEs) are the basic quantities in theoretical descriptions of electronically nonadiabatic processes that are ubiquitous in molecular physics and chemistry. Given the large size of systems of chemical interests, time-dependent density functional theory (TDDFT) is usually the first choice. However, the lack of wave functions in TDDFT renders the formulation of NAC-TDDFT for fo-NACMEs conceptually difficult. The present account aims to analyze the available variants of NAC-TDDFT in a critical but concise manner and meanwhile point out the proper ways for implementation. It can be concluded, from both theoretical and numerical points of view, that the equation of motion-based variant of NAC-TDDFT is the right choice. Possible future developments of this variant are also highlighted.
We study both static and transport properties of model quantum dots, employing density functional theory as well as (numerically) exact methods. For the lattice model under consideration the accuracy of the local-density approximation generally is poor. For weak interaction, however, accurate results are achieved within the optimized effective potential method, while for intermediate interaction strengths a method combining the exact diagonalization of small clusters with density functional theory is very successful. Results obtained from the latter approach yield very good agreement with density matrix renormalization group studies, where the full Hamiltonian consisting of the dot and the attached leads has to be diagonalized. Furthermore we address the question whether static density functional theory is able to predict the exact linear conductance through the dot correctly - with, in general, negative answer.
Reliable and robust convergence to the electronic ground state within density functional theory (DFT) Kohn-Sham (KS) calculations remains a thorny issue in many systems of interest. In such cases, charge sloshing can delay or completely hinder the convergence. Here, we use an approach based on transforming the time-dependent DFT equations to imaginary time, followed by imaginary-time evolution, as a reliable alternative to the self-consistent field (SCF) procedure for determining the KS ground state. We discuss the theoretical and technical aspects of this approach and show that the KS ground state should be expected to be the long-imaginary-time output of the evolution, independent of the exchange-correlation functional or the level of theory used to simulate the system. By maintaining self-consistency between the single-particle wavefunctions and the electronic density throughout the determination of the stationary state, our method avoids the typical difficulties encountered in SCF. To demonstrate dependability of our approach, we apply it to selected systems which struggle to converge with SCF schemes. In addition, through the van Leeuwen theorem, we affirm the physical meaningfulness of imaginary time TDDFT, justifying its use in certain topics of statistical mechanics such as in computing imaginary time path integrals.
Imaginary-time time-dependent Density functional theory (it-TDDFT) has been proposed as an alternative method for obtaining the ground state within density functional theory (DFT) which avoids some of the difficulties with convergence encountered by the self-consistent-field (SCF) iterative method. It-TDDFT was previously applied to clusters of atoms where it was demonstrated to converge in select cases where SCF had difficulty with convergence. In the present work we implement it-TDDFT propagation for {it periodic systems} by modifying the Quantum ESPRESSO package, which uses a plane-wave basis with multiple $boldsymbol{k}$ points, and has the options of non-collinear and DFT+U calculations using ultra-soft or norm-conserving pseudo potentials. We demonstrate that our implementation of it-TDDFT propagation with multiple $boldsymbol{k}$ points is correct for DFT+U non-collinear calculations and for DFT+U calculations with ultra-soft pseudo potentials. Our implementation of it-TDDFT propagation converges to the exact SCF energy (up to the decimal guaranteed by double precision) in all but one case where it converged to a slightly lower value than SCF, suggesting a useful alternative for systems where SCF has difficulty to reach the Kohn-Sham ground state. In addition, we demonstrate that rapid convergence can be achieved if we use adaptive-size imaginary-time-steps for different kinetic-energy plane-waves.
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.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا