Chemical bonding without orbitals
Introduction
The chemical bond [1], [2], [3], [4] is a fundamental concept in chemistry and related sciences. Taking a very simplistic standpoint, the chemical bond somehow serves as glue, binding together the partaking atoms [5], [6]. One part of chemistry is to define and quantify chemical bonding in order to predict molecular stability and how the system will undergo possible changes induced by the environment. Unfortunately, the definition of that chemical bonding is ambiguous. This is due to the fact that atoms in molecules are no physical entities. From a puristic standpoint, the molecule is defined once the number of electrons and the nuclear positions are given. By solving the quantum mechanical equations [7], [5], the statistical distribution of the electrons in the field of the nuclei is determined. In that spirit, the molecule only consists as a whole, namely of all nuclei and all electrons, whereby the latter are usually (but not necessarily) described in terms of orbitals. The equations describe all orbitals at once and their solutions extend over the whole space. This form of wavefunction based quantum mechanics is said to be non-local [7], meaning that changes induced at one position in space might cause noticeable effects far away from the initial position. Such a non-local form of quantum mechanics is not suited for an atomic fragment approach, where those fragments are thought as being subjected to a physical interaction.
Indeed a local from of quantum theory is needed, further allowing for a meaningful separation into atomic fragments. This local form of quantum theory is called density-functional theory (DFT) [8], [9]. In the Hohenberg-Kohn variant[10], the system is described by the electron density itself, rather than by orbitals. This theory is guided by a local causality principle [11], [12], [13], [14], [15], [16], [17], [18] (effective interactions decrease with increasing distance) and, as will be shown in this work, allows for a meaningful separation into atomic fragments, while keeping the initial physical equations unaltered. Although the original ideas of Hohenberg and Kohn dates back more than fifty years, a widespread use of orbital-free DFT has been hampered by the lack of sufficiently accurate functional approximations for the kinetic energy [19]. Whereas considerable progress has been made for one-dimensional systems [20], the treatment of molecular systems or solids requires appropriate three-dimensional functional approximations. First work was based on conventional gradient expansions [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], whereas recent functional development is mainly focused on techniques based on the generalized-gradient-approximation motivated by conjoint arguments [33], [34] or the fulfillment of additional constraints [35], [36], [19]. Information-theory motivated functionals [37], [38], [39], functionals based on response theory [40], [41], and expansions in terms of moment densities [42] have been studied. Recently, a very promising approach based on the Liu-Parr power series expansion [43] has been suggested by Ludeña and coworkers [44], [45], including a detailed investigation concerning the ability of their ansatz to represent the atomic shell structure. Despite much efforts, it still remains difficult to design appropriate functional approximations, which yield properly structured electron densities from variational procedure [46]. This failure can be attributed to an insufficient incorporation of the Pauli exclusion principle [27], [8], [19] in the Pauli kinetic energy and the corresponding Pauli potential [47]. Both terms were intensively studied in the literature [48], [49], [50], [51], [52], [53], [54], [55], [56] as they represent the unknown part of the kinetic energy functional and the corresponding potential, respectively. As a consequence, there is up to now no reliably working orbital-free density functional implementation for the waste variety of chemical substances. There are a few noticeable exceptions treating variational orbital-free density functional approximations [57], [58], [59], [60], [61]. However, widespread applications covering all types of chemical substances without need of recursive adjustment of technical details are still out of reach.
This work offers a new strategy in orbital-free density functional theory based on the reformulation in terms of bifunctionals.
Section snippets
Formulation of the problem
In 1964 Hohenberg and Kohn (HK) rigorously founded density functional theory (DFT) [10]. In their seminal paper they proved the one-to-one correspondence between the energy E of the system and the external potential due to the nuclei (first HK theorem) as well as the minimum principle of that energy functional (second HK theorem) for all electron densities that are associated with some external potential. According to the first HK theorem, the electronic energy of a system , given
Reformulation in terms of bifunctionals
The conventional strategy in the design of energy functionals is to start with an ansatz for the energy in the form of an analytical expression that is an integral over a function f, which in some way contains the electron density (and possibly some other ingredients):From that functional, the potential, being the functional derivative, is obtained analytically, cf. Eq. (3). As mentioned, no ansatz for the kinetic energy has been found that is able to predict chemical
The chemical bond in the atomic fragment approximation
From a chemist’s viewpoint it seems natural to interpret [71], [72], [73] large aggregates as being composed from atomic fragments [74], [75], [76], [77] that somehow hold together. In the following it is shown how to link such a viewpoint with the quantum mechanical equations based on the bifunctional formalism and the atomic fragment approach.
A so-called promolecule [78], [76] is the sum of (usually spherical) atoms centered at the positions of the nuclei for the actual molecule of interest.
Discussion
This work presents an orbital-free implementation of the original Hohenberg-Kohn density functional theory, that is able to predict bonding in molecules. The proposed method is completely parameter-free and does not require analytical ansatzes for the energy functionals. Instead, it is based on the idea that atoms are meaningful pieces of a molecule and thus, a promolecule serves as a good model for the latter. This idea is imposed on the physical equations, originating from density functional
Acknowledgments
The author wishes to thank Dr. M. Kohout for fruitful discussions and substantial encouragement over years. Prof. Dr. M. Ruck is greatly acknowledged for academic support and valuable hints on the manuscript. The Technische Universität Dresden is acknowledged for funding in terms of a Maria-Reiche fellowship.
References (79)
- et al.
Frank discussion o the status of ground-state orbital-free DFT
Adv. Quant. Chem.
(2015) The local potential determining the square root of the ground-state electron density of atoms and molecules from the Schrödinger equation
Phys. Lett. A
(1986)Alternative descriptors of Coulomb systems and their relationship to the kinetic energy
Chem. Phys. Lett.
(2008)- et al.
Finite-temperature orbital-free DFT molecular dynamics: coupling profess and quantum espresso
Comput. Phys. Commun.
(2014) The atom and the molecule
J. Am. Chem. Soc.
(1916)Über Molekülbildung als Frage des Atombaus
Ann. Phys.
(1916)The nature of the chemical bond. IV. The energy of single bonds and the relative electronegativity of atoms
J. Am. Chem. Soc.
(1932)The physical nature of the chemical bond
Rev. Mod. Phys.
(1962)Einführung in die Theoretische Chemie
(2002)- et al.
The Chemical Bond: Fundamental Aspects of Chemical Bonding
(2014)
Modern Quantum Chemistry – Introduction to Advanced Electronic Structure Theory
Density Functional Theory
Density-Functional Theory of Atoms and Molecules
Inhomogeous Electron Gas
Phys. Rev. B
Locality principle in wave mechanics
Proc. Natl. Acad. Sci. USA
Density-matrix electronic structure method with linear system-size scaling
Phys. Rev. B
Density functional and density matrix method scaling linearly with the number of atoms
Phys. Rev. Lett.
Sparsity of the density matrix in Kohn-Sham density functional theory and an assessment of linear system-size scaling methods
Phys. Rev. Lett.
Nearsightedness of electronic matter
Proc. Natl. Acad. Sci.
Nearsightedness of electronic matter seen by a physicist and a chemist
J. Phys. Chem. A
Chemical reactivity descriptors for ambiphilic reagents: dual descriptor, local hypersoftness, and electrostatic potential
J. Phys. Chem. A
Chemical transferability of functional groups follows from the nearsightedness of electronic matter
Proc. Natl. Acad. Sci.
Leading corrections to local approximations
Phys. Rev. B
The calculation of atomic fields
Proc. Cambridge Philos. Soc.
Eine statistische Methode zur Bestimmung einiger Eigenschaften des Atoms und ihre Anwendung auf die Theorie des periodischen Systems der Elemente
Z. Phys.
Zur Theorie der Kernmassen
Z. Phys.
Quantum corrections to the Thomas-Fermi equation
Sov. Phys. JETP
Quantum corrections to the Thomas-Fermi approximation – the Kirzhnits Method
Can. J. Phys.
Sixth-order term of the gradient expansion of the kinetic-energy density functional
Phys. Rev. A
Gradient correction in Thomas-Fermi theory
Phys. Rev. A
Various functionals for the kinetic energy density of an atom or molecule
Phys. Rev. A
Density gradient expansion of the kinetic-energy functional for molecules
Phys. Rev. A
A graph approach to the gradient expansion of density functionals
Int. J. Quant. Chem.
Comparison of kinetic-energy density functionals
Phys. Rev. A
Modified statistical treatment of kinetic energy in the thomas-fermi model
J. Phys. Chem. B
Conjoint gradient correction to the Hartree-Fock kinetic- and exchange-energy density functionals
Phys. Rev. A
Link between the kinetic- and exchange-energy functionals in the generalized gradient approximation
Int. J. Quant. Chem.
Condition on the Kohn-Sham kinetic energy and modern parametrization of the Thomas-Fermi density
J. Chem. Phys.
Progress on new approaches to old ideas: orbital-free density functionals
Cited by (15)
A study of the basis set dependence of the bifunctional expression of the non-interacting kinetic energy for atomic systems
2019, Computational and Theoretical ChemistryCitation Excerpt :This approach requires appropriate density-based approximations for the kinetic energy. Recently, a parameter-free approximation based on a bifunctional expression was introduced [8,9]. In contrast to other known approaches in the literature targeting the energy functional directly [10–41], the proposed bifunctional approach is based on an approximation of the first functional derivative.
Developing orbital-free quantum crystallography: The local potentials and associated partial charge densities
2021, Acta Crystallographica Section B: Structural Science, Crystal Engineering and MaterialsDetailed analysis of deformation potentials with application in orbital-free density functional theory
2021, Acta Crystallographica Section B: Structural Science, Crystal Engineering and MaterialsGeometric potential of the exact electron factorization: Meaning, significance, and application
2023, Physical Review Research