Source author record

Eli Kraisler

Eli Kraisler appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

12works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

atoMEC: An open-source average-atom Python code

Average-atom models are an important tool in studying matter under extreme conditions, such as those conditions experienced in planetary cores, brown and white dwarfs, and during inertial confinement fusion. In the right context, average-atom models can yield results with similar accuracy to simulations which require orders of magnitude more computing time, and thus can greatly reduce financial and environmental costs. Unfortunately, due to the wide range of possible models and approximations, and the lack of open-source codes, average-atom models can at times appear inaccessible. In this paper, we present our open-source average-atom code, atoMEC. We explain the aims and structure of atoMEC to illuminate the different stages and options in an average-atom calculation, and to facilitate community contributions. We also discuss the use of various open-source Python packages in atoMEC, which have expedited its development.

preprint2022arXiv

The Geometric Potential of the Exact Electron Factorization: Meaning, significance and application

The theoretical and computational description of materials properties is a task of utmost scientific and technological importance. A first-principles description of electron-electron interactions poses an immense challenge that is usually approached by converting the many-electron problem to an effective one-electron problem. There are different ways to obtain an exact one-electron theory for a many-electron system. An emergent method is the exact electron factorization (EEF) -- one of the branches of the Exact Factorization approach to many-body systems. In the EEF, the Schrödinger equation for one electron, in the environment of all other electrons, is formulated. The influence of the environment is reflected in the potential $v^{\rm H}$, which represents the energy of the environment, and in a potential $v^{\rm G}$, which has a geometrical meaning. In this paper, we focus on $v^{\rm G}$ and study its properties in detail. We investigate the geometric origin of $v^{\rm G}$ as a metric measuring the change of the environment, exemplify how translation and scaling of the state of the environment are reflected in $v^{\rm G}$, and explain its shape for homo- and heteronuclear diatomic model systems. Based on the close connection between the EEF and density functional theory, we also use $v^{\rm G}$ to provide an alternative interpretation to the Pauli potential in orbital-free density functional theory.

preprint2021arXiv

From Kohn-Sham to many-electron energies via step structures in the exchange-correlation potential

Accurately describing excited states within Kohn-Sham (KS) density functional theory (DFT), particularly those which induce ionization and charge transfer, remains a great challenge. Common exchange-correlation (xc) approximations are unreliable for excited states owing, in part, to the absence of a derivative discontinuity in the xc energy ($Δ$), which relates a many-electron energy difference to the corresponding KS energy difference. We demonstrate, analytically and numerically, how the relationship between KS and many-electron energies leads to the step structures observed in the exact xc potential, in four scenarios: electron addition, molecular dissociation, excitation of a finite system, and charge transfer. We further show that steps in the potential can be obtained also with common xc approximations, as simple as the LDA, when addressed from the ensemble perspective. The article therefore highlights how capturing the relationship between KS and many-electron energies with advanced xc approximations is crucial for accurately calculating excitations, as well as the ground-state density and energy of systems which consist of distinct subsystems.

preprint2021arXiv

On the Geometric Potential and the Relationship between the Exact Electron Factorization and Density Functional Theory

There are different ways to obtain an exact one-electron theory for a many-electron system, and the exact electron factorization (EEF) is one of them. In the EEF, the Schrödinger equation for one electron in the environment of other electrons is constructed. The environment provides the potentials that appear in this equation: A scalar potential $v^{\rm H}$ representing the energy of the environment and another scalar potential $v^{\rm G}$ as well as a vector potential that have geometric meaning. By replacing the interacting many-electron system with the non-interacting Kohn-Sham (KS) system, we show how the EEF is related to density functional theory (DFT) and we interpret the Hartree-exchange-correlation potential as well as the Pauli potential in terms of the EEF. In particular, we show that from the EEF viewpoint, the Pauli potential does not represent the difference between a fermionic and a bosonic non-interacting system, but that it corresponds to $v^{\rm G}$ and partly to $v^{\rm H}$ for the (fermionic) KS system. We then study the meaning of $v^{\rm G}$ in detail: Its geometric origin as a metric measuring the change of the environment is presented. Additionally, its behavior for a simple model of a homo- and heteronucler diatomic is investigated and interpreted with the help of a two-state model. In this way, we provide a physical interpretation for the one-electron potentials that appear in the EEF and in DFT.

preprint2015arXiv

Effect of ensemble generalization on the highest-occupied Kohn-Sham eigenvalue

There are several approximations to the exchange-correlation functional in density-functional theory that accurately predict total energy-related properties of many-electron systems, such as binding energies, bond lengths, and crystal structures. Other approximations are designed to describe potential-related processes, such as charge transfer and photoemission. However, the development of a functional which can serve the two purposes simultaneously is a long-standing challenge. Trying to address it, we employ in the current work the ensemble generalization procedure proposed in Phys. Rev. Lett. 110, 126403 (2013). Focusing on the prediction of the ionization potential via the highest occupied Kohn-Sham eigenvalue, we examine a variety of exchange-correlation approximations: the local spin-density approximation, semi-local generalized gradient approximations, and global and local hybrid functionals. Results for a test set of 26 diatomic molecules and single atoms are presented. We find that the aforementioned ensemble generalization systematically improves the prediction of the ionization potential, for various systems and exchange-correlation functionals, without compromising the accuracy of total energy-related properties. We specifically examine hybrid functionals. These depend on a parameter controlling the ratio of semi-local to non-local functional components. The ionization potential obtained with ensemble-generalized functionals is found to depend only weakly on the parameter value, contrary to common experience with non-generalized hybrids, thus eliminating one aspect of the so-called `parameter dilemma' of hybrid functionals.

preprint2015arXiv

Elimination of the asymptotic fractional dissociation problem in Kohn-Sham density functional theory using the ensemble-generalization approach

Many approximations within density-functional theory spuriously predict that a many-electron system can dissociate into fractionally charged fragments. Here, we revisit the case of dissociated diatomic molecules, known to exhibit this problem when studied within standard approaches, including the local spin-density approximation (LSDA). By employing our recently proposed [E. Kraisler and L. Kronik, Phys. Rev. Lett. 110, 126403 (2013)] ensemble-generalization we find that asymptotic fractional dissociation is eliminated in all systems examined, even if the underlying exchange-correlation (xc) is still the LSDA. Furthermore, as a result of the ensemble generalization procedure, the Kohn-Sham potential develops a spatial step between the dissociated atoms, reflecting the emergence of the derivative discontinuity in the xc energy functional. This step, predicted in the past for the exact Kohn-Sham potential and observed in some of its more advanced approximate forms, is a desired feature that prevents any fractional charge transfer between the system's fragments. It is usually believed that simple xc approximations such as the LSDA cannot develop this step. Our findings show, however, that ensemble generalization to fractional electron densities automatically introduces the desired step even to the most simple approximate xc functionals and correctly predicts asymptotic integer dissociation.

preprint2015arXiv

Fundamental gaps with approximate density functionals: the derivative discontinuity revealed from ensemble considerations

The fundamental gap is a central quantity in the electronic structure of matter. Unfortunately, the fundamental gap is not generally equal to the Kohn-Sham gap of density functional theory (DFT), even in principle. The two gaps differ precisely by the derivative discontinuity, namely, an abrupt change in slope of the exchange-correlation (xc) energy as a function of electron number, expected across an integer-electron point. Popular approximate functionals are thought to be devoid of a derivative discontinuity, strongly compromising their performance for prediction of spectroscopic properties. Here we show that, in fact, all exchange-correlation functionals possess a derivative discontinuity, which arises naturally from the application of ensemble considerations within DFT, without any empiricism. This derivative discontinuity can be expressed in closed form using only quantities obtained in the course of a standard DFT calculation of the neutral system. For small, finite systems, addition of this derivative discontinuity indeed results in a greatly improved prediction for the fundamental gap, even when based on the most simple approximate exchange-correlation density functional - the local density approximation (LDA). For solids, the same scheme is exact in principle, but when applied to LDA it results in a vanishing derivative discontinuity correction. This failure is shown to be directly related to the failure of LDA in predicting fundamental gaps from total energy differences in extended systems.

preprint2015arXiv

One-electron self-interaction and the asymptotics of the Kohn-Sham potential: an impaired relation

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 Weizsäcker kinetic energy density to the (positive) Kohn-Sham kinetic energy density, $τ_\mathrm{W}/τ$, 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.

preprint2014arXiv

A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues

We present and test a new approximation for the exchange-correlation (xc) energy of Kohn-Sham density functional theory. It combines exact exchange with a compatible non-local correlation functional. The functional is by construction free of one-electron self-interaction, respects constraints derived from uniform coordinate scaling, and has the correct asymptotic behavior of the xc energy density. It contains one parameter that is not determined ab initio. We investigate whether it is possible to construct a functional that yields accurate binding energies and affords other advantages, specifically Kohn-Sham eigenvalues that reliably reflect ionization potentials. Tests for a set of atoms and small molecules show that within our local-hybrid form accurate binding energies can be achieved by proper optimization of the free parameter in our functional, along with an improvement in dissociation energy curves and in Kohn-Sham eigenvalues. However, the correspondence of the latter to experimental ionization potentials is not yet satisfactory, and if we choose to optimize their prediction, a rather different value of the functional's parameter is obtained. We put this finding in a larger context by discussing similar observations for other functionals and possible directions for further functional development that our findings suggest.

preprint2014arXiv

Higher Ionization Energies of Atoms in Density Functional Theory

Density functional theory (DFT) is an exact alternative formulation of quantum mechanics, in which it is possible to calculate the total energy, the spin and the charge density of many-electron systems in the ground state. In practice, it is necessary to use uncontrolled approximations that can mainly be verified against experimental data. Atoms and ions are simple systems, where the approximations of DFT can be easily tested. We have calculated within DFT the total energies, spin and higher ionization energies of all the ions of elements with 1 $\leq$ Z $\leq$ 29. We find the calculations in close agreement with experiment, with an error of typically less than ca. 1% for 1 $\leq$ Z $\leq$ 29. Surprisingly, the error depends on the electronic configuration of the ion in both local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof general gradient approximation (PBE-GGA) and independent of both self-interaction correction (SIC) and relativistic corrections. Larger errors are found for systems in which the spin-spin correlation is significant, which indicates the possible benefit from an orbital-dependent formulation of the correlation energy functional.

preprint2010arXiv

Ensemble v-representable ab-initio density functional calculation of energy and spin in atoms: atest of exchange-correlation approximations

The total energies and the spin states for atoms and their first ions with Z = 1-86 are calculated within the the local spin-density approximation (LSDA) and the generalized-gradient approximation (GGA) to the exchange-correlation (xc) energy in density-functional theory. Atoms and ions for which the ground-state density is not pure-state v-representable, are treated as ensemble v- representable with fractional occupations of the Kohn-Sham system. A newly developed algorithm which searches over ensemble v-representable densities [E. Kraisler et al., Phys. Rev. A 80, 032115 (2009)] is employed in calculations. It is found that for many atoms the ionization energies obtained with the GGA are only modestly improved with respect to experimental data, as compared to the LSDA. However, even in those groups of atoms where the improvement is systematic, there remains a non-negligible difference with respect to the experiment. The ab-initio electronic configuration in the Kohn-Sham reference system does not always equal the configuration obtained from the spectroscopic term within the independent-electron approximation. It was shown that use of the latter configuration can prevent the energy-minimization process from converging to the global minimum, e.g. in lanthanides. The spin values calculated ab-initio fit the experiment for most atoms and are almost unaffected by the choice of the xc-functional. Among the systems with incorrectly obtained spin there exist some cases (e.g. V, Pt) for which the result is found to be stable with respect to small variations in the xc-approximation. These findings suggest a necessity for a significant modification of the exchange-correlation functional, probably of a non-local nature, to accurately describe such systems. PACS numbers: 31.15.A

preprint2010arXiv

Fractional occupation in Kohn-Sham density-functional theory and the treatment of non-pure-state v-representable densities

In the framework of Kohn-Sham density-functional theory, systems with ground-state densities that are not pure-state v-representable in the non-interacting reference system (PSVR) occur frequently. In the present contribution, a new algorithm, which allows the solution of such systems, is proposed. It is shown that the use of densities which do not correspond to a ground state of their non-interacting reference system is forbidden. As a consequence, the proposed algorithm considers only non-interacting ensemble v-representable densities. The Fe atom, a well-known non-PSVR system, is used as an illustration. Finally, the problem is analyzed within finite temperature density-functional theory, where the physical significance of fractional occupations is exposed and the question of why degenerate states can be unequally occupied is resolved.