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

Reducing charge delocalization error of density functional theory

118   0   0.0 ( 0 )
 نشر من قبل Jing Kong
 تاريخ النشر 2021
  مجال البحث فيزياء
والبحث باللغة English




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

The charge delocalization error, besides nondynamic correlation, has been a major challenge to density functional theory. Contemporary functionals undershoot the dissociation of symmetric charged dimers A2+, a simple but stringent test, predict a spurious barrier and improperly delocalize charges for charged molecular clusters. We extend a functional designed for nondynamic correlation to treat the charge delocalization error by modifying the nondynamic correlation for parallel spins. The modified functional eliminates those problems and reduces the multielectron self-interaction error. Furthermore, its results are the closest to those of CCSD(T) in the whole range of the dissociation compared with contemporary functionals. It correctly localizes the net positive charge in (CH4)n+ clusters and predicts a nearly constant ionization potential as a result. Testing of the SIE4x4 set shows that the new functional outperforms a wide variety of functionals assessed for this set in the literature. Overall, we show the feasibility of treating charge delocalization together with nondynamic correlation.

قيم البحث

اقرأ أيضاً

138 - Julien Toulouse 2021
A relativistic density-functional theory based on a Fock-space effective quantum-electrodynamics (QED) Hamiltonian using the Coulomb or Coulomb-Breit two-particle interaction is developed. This effective QED theory properly includes the effects of va cuum polarization through the creation of electron-positron pairs but does not include explicitly the photon degrees of freedom. It is thus a more tractable alternative to full QED for atomic and molecular calculations. Using the constrained-search formalism, a Kohn-Sham scheme is formulated in a quite similar way to non-relativistic density-functional theory, and some exact properties of the involved density functionals are studied, namely charge-conjugation symmetry and uniform coordinate scaling. The usual no-pair Kohn-Sham scheme is obtained as a well-defined approximation to this relativistic density-functional theory.
The time-dependent density functional theory (TDDFT) has been broadly used to investigate the excited-state properties of various molecular systems. However, the current TDDFT heavily relies on outcomes from the corresponding ground-state density fun ctional theory (DFT) calculations which may be prone to errors due to the lack of proper treatment in the non-dynamical correlation effects. Recently, thermally-assisted-occupation density functional theory (TAO-DFT) [J.-D. Chai, textit{J. Chem. Phys.} textbf{136}, 154104 (2012)], a DFT with fractional orbital occupations, was proposed, explicitly incorporating the non-dynamical correlation effects in the ground-state calculations with low computational complexity. In this work, we develop time-dependent (TD) TAO-DFT, which is a time-dependent, linear-response theory for excited states within the framework of TAO-DFT. With tests on the excited states of H$_{2}$, the first triplet excited state ($1^3Sigma_u^+$) was described well, with non-imaginary excitation energies. TDTAO-DFT also yields zero singlet-triplet gap in the dissociation limit, for the ground singlet ($1^1Sigma_g^+$) and the first triplet state ($1^3Sigma_u^+$). In addition, as compared to traditional TDDFT, the overall excited-state potential energy surfaces obtained from TDTAO-DFT are generally improved and better agree with results from the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD).
In numerical computations of response properties of electronic systems, the standard model is Kohn-Sham density functional theory (KS-DFT). Here we investigate the mathematical status of the simplest class of excitations in KS-DFT, HOMO-LUMO excitati ons. We show using concentration-compactness arguments that such excitations, i.e. excited states of the Kohn-Sham Hamiltonian, exist for $Z>N$, where $Z$ is the total nuclear charge and $N$ is the number of electrons. The result applies under realistic assumptions on the exchange-correlation functional, which we verify explicitly for the widely used PZ81 and PW92 functionals. By contrast, and somewhat surprisingly, we find using a method of Glaser, Martin, Grosse, and Thirring cite{glaser1976} that in case of the hydrogen and helium atoms, excited states do not exist in the neutral case $Z=N$ when the self-consistent KS ground state density is replaced by a realistic but easier to analyze approximation (in case of hydrogen, the true Schr{o}dinger ground state density). Implications for interpreting minus the HOMO eigenvalue as an approximation to the ionization potential are indicated.
109 - Emmanuel Giner 2020
We extend to strongly correlated molecular systems the recently introduced basis-set incompleteness correction based on density-functional theory (DFT) [E. Giner et al., J. Chem. Phys. 149, 194301 (2018)]. This basis-set correction relies on a mappin g between wave-function calculations in a finite basis set and range-separated DFT (RSDFT) through the definition of an effective non-divergent interaction corresponding to the electron-electron Coulomb interaction projected in the finite basis set. This enables the use of RSDFT-type complementary density functionals to recover the dominant part of the short-range correlation effects missing in this finite basis set. To study both weak and strong correlation regimes we consider the potential energy curves of the H10, N2, O2, and F2 molecules up to the dissociation limit, and we explore various approximations of complementary functionals fulfilling two key properties: spin-multiplet degeneracy (i.e., independence of the energy with respect to the spin projection Sz) and size consistency. Specifically, we investigate the dependence of the functional on different types of on-top pair densities and spin polarizations. The key result of this study is that the explicit dependence on the on-top pair density allows one to completely remove the dependence on any form of spin polarization without any significant loss of accuracy. Quantitatively, we show that the basis-set correction reaches chemical accuracy on atomization energies with triple-zeta quality basis sets for most of the systems studied here. Also, the present basis-set incompleteness correction provides smooth potential energy curves along the whole range of internuclear distances.
Exciton formation leads to J-bands in solid pentacene. Describing these exciton bands represents a challenge for both time-dependent (TD) density-functional theory (DFT) and for its semiempirical analogue, namely for TD density-functional tight bindi ng (DFTB) for three reasons (i) solid pentacene and pentacene aggregates are bound only by van der Waals forces which are notoriously difficult to describe with DFT and DFTB, (ii) the proper description of the long-range coupling between molecules, needed to describe Davydov splitting, is not easy to include in TD-DFT with traditional functionals and in TD-DFTB, and (iii) mixing may occur between local and charge transfer excitons, which may, in turn, require special functionals. We assess how far TD-DFT and TD-DFTB have progressed towards a correct description of this type of exciton by including both a dispersion correction for the ground state and a range-separated hybrid functional for the excited state. Analytic results for parallel-stacked ethylene are derived which go beyond Kashas exciton model in that we are able to make a clear distinction between charge transfer and energy transfer excitons. This is further confirmed when it is shown that range-separated hybrids have a markedly greater effect on charge-transfer excitons than on energy-transfer excitons in the case of parallel-stacked pentacenes. TD-DFT calculations with the CAM-B3LYP functional and TD-lc-DFT calculations lead to negligeable excitonic corrections for the herringbone crystal structure, possibly because of an overcorrection of charge-transfer effects. In this case, TD-DFT calculations with the B3LYP functional or TD-DFTB calculations parameterized to B3LYP give the best results for excitonic corrections for the herringbone crystal structure as judged from comparison with experimental spectra and with Bethe-Salpeter equation calculations from the literature.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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