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

Self-interaction corrected Kohn-Sham effective potentials using the density-consistent effective potential method

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




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

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn-Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. In this work, we adapt and implement the density-consistent effective potential(DCEP) method of Kohut, Ryabinkin, and Staroverov to obtain SI corrected local effective potentials from the SI corrected Fermi-Lowdin orbitals and density in the FLOSIC scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO-LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations.



قيم البحث

اقرأ أيضاً

We present a method to invert a given density and find the Kohn-Sham (KS) potential in Density Functional Theory (DFT) which shares that density. Our method employs the concept of screening density, which is naturally constrained by the inversion pro cedure and thus ensures the density being inverted leads to a smooth KS potential with correct asymptotic behaviour. We demonstrate the applicability of our method by inverting both local (LDA) and non-local (Hartree-Fock and Coupled Cluster) densities; we also show how the method can be used to mitigate the effects of self-interactions in common DFT potentials with appropriate constraints on the screening density.
One-electron self-interaction and an incorrect asymptotic behavior of the Kohn-Sham exchange-correlation potential are among the most prominent limitations of many present-day density functionals. However, a one-electron self-interaction-free energy does not necessarily lead to the correct long-range potential. This is here shown explicitly for local hybrid functionals. Furthermore, carefully studying the ratio of the von Weizsacker kinetic energy density to the (positive) Kohn-Sham kinetic energy density, $tau_mathrm{W}/tau$, reveals that this ratio, which frequently serves as an iso-orbital indicator and is used to eliminate one-electron self-interaction effects in meta-generalized-gradient approximations and local hybrid functionals, can fail to approach its expected value in the vicinity of orbital nodal planes. This perspective article suggests that the nature and consequences of one-electron self-interaction and some of the strategies for its correction need to be reconsidered.
We extend the range-separated double-hybrid RSH+MP2 method [J. G. Angyan et al., Phys. Rev. A 72, 012510 (2005)], combining long-range HF exchange and MP2 correlation with a short-range density functional, to a fully self-consistent version using the optimized-effective-potential technique in which the orbitals are obtained from a local potential including the long-range HF and MP2 contributions. We test this approach, that we name RS-OEP2, on a set of small closed-shell atoms and molecules. For the commonly used value of the range-separation parameter $mu=0.5$ bohr$^{-1}$, we find that self-consistency does not seem to bring any improvement for total energies, ionization potentials, and electronic affinities. However, contrary to the non-self-consistent RSH+MP2 method, the present RS-OEP2 method gives a LUMO energy which physically corresponds to a neutral excitation energy and gives local exchange-correlation potentials which are reasonably good approximations to the corresponding Kohn-Sham quantities. At a finer scale, we find that RS-OEP2 gives largely inaccurate correlation potentials and correlated densities, which points to the need of further improvement of this type of range-separated double hybrids.
We construct exact Kohn-Sham potentials for the ensemble density-functional theory (EDFT) from the ground and excited states of helium. The exchange-correlation (XC) potential is compared with the quasi-local-density approximation and both single det erminant and symmetry eigenstate ghost-corrected exact exchange approximations. Symmetry eigenstate Hartree-exchange recovers distinctive features of the exact XC potential and is used to calculate the correlation potential. Unlike the exact case, excitation energies calculated from these approximations depend on ensemble weight, and it is shown that only the symmetry eigenstate method produces an ensemble derivative discontinuity. Differences in asymptotic and near-ground-state behavior of exact and approximate XC potentials are discussed in the context of producing accurate optical gaps.
Perdew-Zunger self-interaction correction (PZ-SIC) offers a route to remove self-interaction errors on an orbital-by-orbital basis. A recent formulation of PZ-SIC by Pederson, Ruzsinszky and Perdew proposes restricting the unitary transformation to l ocalized orbitals called Fermi-Lowdin orbitals. This formulation, called the FLOSIC method, simplifies PZ-SIC calculations and was implemented self-consistently using a Jacobi-like (FLOSIC-Jacobi) iteration scheme. In this work we implement the FLOSIC approach using the Krieger-Li-Iafrate (KLI) approximation to the optimized effective potential (OEP). We compare the results of present FLOSIC-KLI approach with FLOSIC-Jacobi scheme for atomic energies, atomization energies, ionization energies, barrier heights, polarizability of chains of hydrogen molecules etc. to validate the FLOSIC-KLI approach. The FLOSIC-KLI approach, which is within the realm of Kohn-Sham theory, predicts smaller energy gaps between frontier orbitals due to the lowering of eigenvalues of the lowest unoccupied orbitals. Results show that atomic energies, atomization energies, ionization energy as an absolute of highest occupied orbital eigenvalue, and polarizability of chains of hydrogen molecules between the two methods agree within 2%. Finally the FLOSIC-KLI approach is used to determine the vertical ionization energies of water clusters.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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