ترغب بنشر مسار تعليمي؟ اضغط هنا

Linear-scaling time-dependent density-functional theory in the linear response formalism

247   0   0.0 ( 0 )
 نشر من قبل Tim Zuehlsdorff
 تاريخ النشر 2013
  مجال البحث فيزياء
والبحث باللغة English




اسأل ChatGPT حول البحث

We present an implementation of time-dependent density-functional theory (TDDFT) in the linear response formalism enabling the calculation of low energy optical absorption spectra for large molecules and nanostructures. The method avoids any explicit reference to canonical representations of either occupied or virtual Kohn-Sham states and thus achieves linear-scaling computational effort with system size. In contrast to conventional localised orbital formulations, where a single set of localised functions is used to span the occupied and unoccupied state manifold, we make use of two sets of in situ optimised localised orbitals, one for the occupied and one for the unoccupied space. This double representation approach avoids known problems of spanning the space of unoccupied Kohn-Sham states with a minimal set of localised orbitals optimised for the occupied space, while the in situ optimisation procedure allows for efficient calculations with a minimal number of functions. The method is applied to a number of medium sized organic molecules and a good agreement with traditional TDDFT methods is observed. Furthermore, linear scaling of computational cost with system size is demonstrated on a system of carbon nanotubes.



قيم البحث

اقرأ أيضاً

Linear-response time-dependent (TD) density-functional theory (DFT) has been implemented in the pseudopotential wavelet-based electronic structure program BigDFT and results are compared against those obtained with the all-electron Gaussian-type orbi tal program deMon2k for the calculation of electronic absorption spectra of N2 using the TD local density approximation (LDA). The two programs give comparable excitation energies and absorption spectra once suitably extensive basis sets are used. Convergence of LDA density orbitals and orbital energies to the basis-set limit is significantly faster for BigDFT than for deMon2k. However the number of virtual orbitals used in TD-DFT calculations is a parameter in BigDFT, while all virtual orbitals are included in TD-DFT calculations in deMon2k. As a reality check, we report the x-ray crystal structure and the measured and calculated absorption spectrum (excitation energies and oscillator strengths) of the small organic molecule N-cyclohexyl-2-(4-methoxyphenyl)imidazo[1,2-a]pyridin-3-amine.
We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energie s and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.
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 excit ation 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.
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 f or 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.
Linear scaling methods for density-functional theory (DFT) simulations are formulated in terms of localised orbitals in real-space, rather than the delocalised eigenstates of conventional approaches. In local-orbital methods, relative to conventional DFT, desirable properties can be lost to some extent, such as the translational invariance of the total energy of a system with respect to small displacements and the smoothness of the potential energy surface. This has repercussions for calculating accurate ionic forces and geometries. In this work we present results from textsc{onetep}, our linear scaling method based on localised orbitals in real-space. The use of psinc functions for the underlying basis set and on-the-fly optimisation of the localised orbitals results in smooth potential energy surfaces that are consistent with ionic forces calculated using the Hellmann-Feynman theorem. This enables accurate geometry optimisation to be performed. Results for surface reconstructions in silicon are presented, along with three example systems demonstrating the performance of a quasi-Newton geometry optimisation algorithm: an organic zwitterion, a point defect in an ionic crystal, and a semiconductor nanostructure.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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