BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Manoj K. Harbola appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
A general penalty method is presented for the construction of of Kohn-Sham system for given density through Levy's constrained-search. The method uses a functional $S[ρ]$ of one's choice. Different forms of $S[ρ]$ are employed to calculate the kinetic energy and exchange-correlation potential of atoms, jellium spheres, and Hookium and consistency among results obtained from them is shown.
Over the past few years it has been pointed out that direct inversion of accurate but approximate ground state densities leads to Kohn-Sham exchange-correlation (xc) potentials that can differ significantly from the exact xc potential of a given system. On the other hand, the corresponding wavefunction based construction of exchange-correlation potential as done by Baerends et al. and Staroverov et al. obviates such problems and leads to potentials that are very close to the true xc potential. In this paper, we provide an understanding of why the wavefunction based approach gives the exchange-correlation potential accurately. Our understanding is based on the work of Levy, Perdew and Sahni (LPS) who gave an equation for the square root of density (density amplitude) and the expression for the associated effective potential in the terms of the corresponding wavefunction. We show that even with the use of approximate wavefunctions the LPS expression gives accurate effective and exchange-correlation potentials. Based on this we also identify the source of difference between the potentials obtained from a wavefunction and those given by the inversion of the associated density. Finally, we suggest exploring the possibility of obtaining accurate ground-state density from an approximate wavefunction for a system by making use of the LPS effective potential.
We report results from a fast, efficient, and first-principles full-potential N$^{th}$-order muffin-tin orbital (FP-NMTO) method combined with van Leeuwen-Baerends correction to local density exchange-correlation potential. We show that more complete and compact basis set is critical in improving the electronic and structural properties. We exemplify the self-consistent FP-NMTO calculations on group IV and III-V semiconductors. Notably, predicted bandgaps, lattice constants, and bulk moduli are in good agreement with experiments (e.g., we find for Ge $0.86~e$V, $5.57$~Å, $75$~GPa vs. measured $0.74~e$V, $5.66$~Å, $77.2$~GPa). We also showcase its application to the electronic properties of 2-dimensional $h-$BN and $h-$SiC, again finding good agreement with experiments.
We formulate a spin-polarized van Leeuwen and Baerends (vLB) correction to the local density approximation (LDA) exchange potential [Phys. Rev. A 49, 2421 (1994)] that enforces the ionization potential (IP) theorem following Stein et al. [Phys. Rev. Lett. 105, 266802 (2010)]. For electronic-structure problems, the vLB-correction replicates the behavior of exact-exchange potentials, with improved scaling and well-behaved asymptotics, but with the computational cost of semi-local functionals. The vLB+IP corrections produces large improvement in the eigenvalues over that from LDA due to correct asympotic behavior and atomic shell structures, as shown on rare-gas, alkaline-earth, zinc-based oxides, alkali-halides, sulphides, and nitrides. In half-Heusler alloys, this asymptotically-corrected LDA reproduces the spin-polarized properties correctly, including magnetism and half-metallicity. We also considered finite-sized systems [e.g., ringed boron-nitirde (B$_{12}$N$_{12}$) and graphene (C$_{24}$)] to emphasize the wide applicability of the method.
Using an accurate semi-analytic wavefunction for two electron atoms, we construct the external potential for varying strength of electron-electron (e-e) interaction. Using this potential we explicitly calculate the energy of their positive ion and show that the ionization-potential of these systems remains unchanged with respect to the strength parameter for e-e interaction. Furthermore, using total energies of these systems as a function of strength parameter, we provide a new perspective into a variety of hybrid functionals.
Levy-Perdew-Sahni (LPS) provided a proof of the connection between the density and ionization potential of a many-electron system using the equation for the density directly. This equation employs an effective potential which is derived by LPS from the many-electron wavefunction of the system. In this paper, we calculate this potential explicitly by employing an accurate wavefunction for two-electron systems. With this wavefunction connection the LPS equation is then solved self- consistently employing this potential and shown to lead to accurate densities and chemical potentials for these systems.
It is shown for two electron atoms that ground-state wavefunctions of the form \begin{equation} Ψ(\vec{r_{1}}, \vec{r_{2}})=ϕ(\vec{r_{1}})ϕ(\vec{r_{2}})(\cosh ar_{1}+\cosh ar_{2})(1+0.5 r_{12}e^{-b r_{12}}) \end{equation} where $\vec{r_{1}}$ and $\vec{r_{2}}$ are the coordinates of two electrons and $r_{12}=|\vec{r_{1}}-\vec{r_{2}}|$, can be made highly accurate by optimizing $a$, $b$ and $ϕ$. This is done by solving a variationally derived equation for $ϕ$ for a given $a$ and $b$ and finding $a$ and $b$ so that the expectation value of the Hamiltonian is minimum. For the set $\{a, b, ϕ\}$ the values for various quantities obtained from the above wavefunction are compared with those given by $204$-parameter wavefunction of Koga et al.[11] and are found to be matching quite accurately(within ppm) with them.
The correlation energies of various atoms in their excited-states are estimated by modelling the Coulomb hole following the previous work by Chakravorty and Clementi. The parameter in the model is fixed by making the corresponding Coulomb hole to satisfy the exact constraint of charge neutrality.
Solutions of time-independent Schrodinger equation for potentials periodic in space satisfy Bloch theorem. The theorem has been used to obtain solutions of the Schrodinger equation for periodic systems by expanding them in terms of plane waves of appropriate wave-vectors and then diagonalising the resulting Hamiltonian matrix. However, this method can give exact solutions only if an infinite number of plane-waves are used in the wavefunction. Thus the Schrodinger equation has not been solved exactly even for one-dimensional infinitely large periodic systems despite the simplification that the theorem provides. On the other hand, solutions of one-dimensional Schrodinger equation for finite systems are obtained routinely by using standard numerical techniques to integrate the Schrodinger equation directly. In this paper, we extend this approach to infinitely large systems and use Bloch theorem in a novel way to find direct exact numerical solutions of the Schrodinger equation for arbitrary one-dimensional periodic potentials. In addition to providing a different way of solving the Schrodinger equation for such systems, the simplicity of the algorithm renders it a great pedagogical value.
The modified local spin density functional and the related local potential for excited states is tested by employing the ionization potential theorem. The functional is constructed by splitting $k$-space. Since its functional derivative cannot be obtained easily, the corresponding potential is given by analogy to its ground-state counterpart. Further to calculate the highest occupied orbital energy $ε_{max}$ accurately, the potential is corrected for its asymptotic behavior by employing the van Leeuwen and Baerends correction to it. $ε_{max}$ so obtained is then compared with the $Δ$SCF ionization energy calculated using the MLSD functional. It is shown that the two match quite accurately.
We present an approach for accurate calculation of exchange energy and total energy for excited states using time independent density functional formalism. This is done by inclusion of gradient correction into the excited state exchange energy functionals developed by us. We have incorporated Becke and Perdew Wang corrections into our functional and have studied various types of excited states having one and two gaps in occupation of orbitals.
Over the past decade, fundamentals of time independent density functional theory for excited state have been established. However, construction of the corresponding energy functionals for excited states remains a challenging problem. We have developed a method for constructing functionals for excited states by splitting k-space according to the occupation of orbitals. In this paper we first show the accuracy of kinetic energy functional thus obtained. We then perform a response function analysis of the kinetic energy functional proposed by us and show why method of splitting the k-space could be the method of choice for construction of energy functionals for excited states.
We show that our procedure of constructing excited-state energy functionals by splitting k-space, employed so far to obtain exchange energies of excited-states, is quite general. We do so by applying the same method to construct modified Thomas-Fermi kinetic energy functional and its gradient expansion up to the second order for the excited-states. We show that the resulting kinetic energy functional has the same accuracy for the excited-states as the ground-state functionals do for the ground-states.
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.