Next Article in Journal
Betulinic Acid for Cancer Treatment and Prevention
Next Article in Special Issue
Exact and Effective Pair-Wise Potential for Protein-Ligand Interactions Obtained from a Semiempirical Energy Partition
Previous Article in Journal
Targeting Receptor Tyrosine Kinases for Chemoprevention by Green Tea Catechin, EGCG
Previous Article in Special Issue
Bonding in Mercury-Alkali Molecules: Orbital-driven van der Waals Complexes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Density Functionals of Chemical Bonding

Laboratory of Computational and Structural Physical Chemistry, Chemistry Department, West University of Timişoara, Pestalozzi Street No.16, Timişoara, RO-300115, Romania
Int. J. Mol. Sci. 2008, 9(6), 1050-1095; https://doi.org/10.3390/ijms9061050
Submission received: 30 April 2008 / Revised: 9 June 2008 / Accepted: 10 June 2008 / Published: 26 June 2008
(This article belongs to the Special Issue The Chemical Bond and Bonding)

Abstract

:
The behavior of electrons in general many-electronic systems throughout the density functionals of energy is reviewed. The basic physico-chemical concepts of density functional theory are employed to highlight the energy role in chemical structure while its extended influence in electronic localization function helps in chemical bonding understanding. In this context the energy functionals accompanied by electronic localization functions may provide a comprehensive description of the global-local levels electronic structures in general and of chemical bonds in special. Becke-Edgecombe and author’s Markovian electronic localization functions are discussed at atomic, molecular and solid state levels. Then, the analytical survey of the main workable kinetic, exchange, and correlation density functionals within local and gradient density approximations is undertaken. The hierarchy of various energy functionals is formulated by employing both the parabolic and statistical correlation degree of them with the electronegativity and chemical hardness indices by means of quantitative structure-property relationship (QSPR) analysis for basic atomic and molecular systems.

1. Introduction

In Walter Kohn’s lecture, with the occasion of receiving his Nobel Prize in Chemistry [1], back in 1998, for density functional theory (DFT) theory [24], there was formulated a quite provoking assertion affirming that, heuristically, the general eigen-wave-function Ψ (r1,...,rN ) associated to a system of N electrons fails to be a legitimate scientific concept when NN0 ≈ 103.
Nevertheless, this affirmation may be at any time turned in a theorem, eventually as Kohn’s zero DFT theorem, with a proof following the van Vleck prescription of the so called „exponential wall”, leaving with the applicability limits of the conceptually eigen-wave function of multi-electronics systems. However, before proceeding to demonstration, there must be noted that such reality limitation characterized by eigen-wave-function of the multi electronic-systems is not transferable at the quantum mechanics postulates, but providing an alternative quantum scheme, thus paralleling Schrödinger equation, however in a more generally integrated level.
The demonstration of the non-representability of the eigen-function for systems containing more than N0 ≈ 103 electrons involves two aspects: the accuracy of representation by using the eigen-wave function and the possibility of measuring it.
Regarding the accuracy of the Ψ(r1,...,rN) representation there is widely known that it associates with the density of probability smoothly approaching unity, written in the “liberal” form (according to Kohn):
| Ψ * ( r 1 , ... r N ) Ψ ( r 1 , ... , r N ) | = 1 ɛ , ɛ = 10 −2
Now, considering a collection of N’ molecules the total density of probability of this multi-molecular system (and implicit a multi-electronic one) will consequently be:
| ( Ψ * Ψ ) 1 ... ( Ψ * Ψ ) N ' | = ( 1 ɛ ) N ' = exp ( ɛ N ' ) = exp [ 10 2 N ' ] .
For N’=103 molecules in whatever aggregates, e.g. solids, clusters, super-molecules, or biological macro-molecules, the total density of probability will result from (2) as exp[−10]≈5×10−5, meaning that it is determined with much less accuracy compared with the degree of individual eigen-function localization precision in (1). Since each molecular system has at least one electron there follows the threshold limit of N0103 electrons from where the lost in associated wavefunction nature is recorded.
This result, relaying on the exponential form (2), justifies the title of „exponential wall” for the wavefunction limitation.
Then, going to the measurable issue of such eigenfunctions, let’s ask how many bits are necessary for recording its quantum dimension? Assume, again, the working wavefunction Ψ(r1,...,rN) for all the N electrons in a concerned system. The N electrons in system have a total of 3N space variables (in the configuration space); let’s now assume an average of q bits necessary in measuring a single variable from the total of 3N; there results a total information of
B = q 3 N
bits for recording (storing) the total eigen-wave-function of the system. However, a simple evaluation of the dimension (3) shows that for a minimum of q=3 bits/variable and for the above consecrated minimum limit of N=N0≈103 electrons in the system the total yield of necessary bits for recording is about of 101000 order – a truly non-realistic dimension. This can be immediately visualized if one recalls that the total number of baryons (i.e. all fermions and elementary particles of protonic and neutronic type, but not limited only to these) estimated in entirely Universe (summing up all existing atoms and free nuclei in the plasma state, but not only limited to these) gives a result of about 1080 order.
Definitely, the concept of eigen-wave function must be enlarged or modified in such a manner that the quantum description does not be blocked by the exponential wall: from where we can start? Firstly, as was exposed, the eigen-wave function in the configuration space multiplies in an exponential manner the variables accounting for the number and the position of the electrons; thus, the configuration space must be avoided. Then, the density of probability must be reformulated as such the exponential wall for a poly-electronic system be avoided while preserving the dependency of the total number of electrons N.
Fortunately, the above described conceptual project was unfolded in 1963 when Walter Kohn met in Paris (at École Normale Supérieure), during his sabbatical semester, the mate Pierre Hohenberg who was working at the description of the metallic alloys (specially the CuxZn1-x systems) by using quantum traditionally methods of averaging crystalline periodic field. Studies of this type of problems often start from the level of the uniform electronic density referential upon which specific perturbation treatments are applied. From this point Kohn and Hohenberg made two crucial further steps in reformulation of the quantum picture of the matter structure: one referred at the electronic density, and another at the relation between electronic density with the externally applied potential on the electronic system; they were consecrated in the so called Hohenberg-Kohn (KH) theorems [2, 57]. They were the fundaments of new emerging quantum density functional theory that mostly impacted the reformation of the quantum chemistry itself and its foreground principles of structure and transformation.
The present work likes to review some fundamental aspects of density functional theory highlighting on the primer conceptual and computational consequences in electronic localization and chemical reactivity.

2. Primary Density Functional Theory Concepts

2.1. Hohenberg-Kohn theorems

The first Hohenberg-Kohn (HK1) theorem gives space to the concept of electronic density of the system ρ(r) in terms of the extensive relation with the N electrons from the system that it characterizes [8]:
ρ ( r ) d r = N    .
The relation (4) as much simple it could appears stands as the decisive passage from the eigen-wave function level to the level of total electronic density [911]:
ρ ( r ) = N Ψ * ( r , r 2 , ... , r N ) Ψ ( r , r 2 , ... , r N ) d r 2 ... d r N .
Firstly, Eq. (5) satisfies Eq. (4); this can be used also as simple immediate proof of the relation (4) itself. Then, the dependency from the 3N-dimensions of configuration space was reduced at 3 coordinates in the real space, physically measurable.
However, still remains the question: what represents the electronic density of Eq. (5)? Definitely, it neither represents the electronic density in the configuration space nor the density of a single electron, since the N-electronic dependency as multiplication factor of the multiple integral in (5). What remains is that ρ(r) is simple the electronic density (of the whole concerned system) in „r space point. Such simplified interpretation, apparently classics, preserves its quantum roots through the averaging (integral) over the many-electronic eigenfunction Ψ(r1,...,rN) in (5). Alternatively, the explicit non-dependency of density on the wave function is also possible within the quantum statistical approach where the relation with partition function of the system (the global measure of the distribution of energetic states of a system) is mainly considered.
The major consequence of this theorem consists in defining of the total energy of a system as a function of the electronic density function in what is known as the density functional [8, 9, 12]:
E [ ρ ] = F H K [ ρ ] + C A [ ρ ] ,
from where the name of the theory. The terms of energy decomposition in (6) are identified as: the Hohenberg-Kohn density functional
F H K [ ρ ] = T [ ρ ] + V e e [ ρ ]
viewed as the summed electronic kinetic T[ρ] and electronic repulsion Vee [ρ], and the so called chemical action term [1214]:
C A [ ρ ] = ρ ( r ) V ( r ) d r ,
being the only explicit functional of total energy.
Although not entirely known the HK functional has a remarkably property: it is universally, in a sense that both the kinetic and inter-electronic repulsion are independent of the concerned system. The consequence of such universal nature offers the possibility that once it is exactly or approximately knew the HK functional for a given external potential V(r) remain valuable for any other type of potential V’(r) applied on the concerned many-electronic system. Let’s note the fact that V(r) should be not reduced only to the Coulombic type of potentials but is carrying the role of the generic potential applied, that could beg of either an electric, magnetic, nuclear, or even electronic nature as far it is external to the system fixed by the N electrons in the investigated system.
Once “in game” the external applied potential provides the second Hohenberg-Kohn (HK2) theorem. In short, HK2 theorem says that “the external applied potential is determined up to an additive constant by the electronic density of the N-electronic system ground state”.
In mathematical terms, the theorem assures the validity of the variational principle applied to the density functional (6) relation, i.e. [6]
E [ ρ ¯ ] E [ ρ ] δ E [ ρ ] = 0
for every electronic test density ρ̄ around the real density ρ of the ground state.
The proof of variational principle in (9), or, in other words, the one-to-one correspondence between the applied potential and the ground state electronic density, employs the reduction ad absurdum procedure. That is to assume that the ground state electronic density ρ(r) corresponds to two external potentials (V1, V2) fixing two associate Hamiltonians (H1, H2) to which two eigen-total energy (E1, E2) and two eigen-wave functions (Ψ1, Ψ2) are allowed. Now, if eigen-function Ψ1 is considered as the true one for the ground state the variational principle (9) will cast as the inequality:
E 1 [ ρ ] = Ψ 1 * H ^ 1 Ψ 1 d τ < Ψ 2 * H ^ 1 Ψ 2 d τ = Ψ 2 * [ H ^ 2 + ( H ^ 1 H ^ 2 ) ] Ψ 2 d τ
which is further reduced, on universality reasons of the HK functional in (6), to the form:
E 1 [ ρ ] < E 2 [ ρ ] + ρ ( r ) [ V 1 ( r ) V 2 ( r ) ] d r .
On another way, if the eigen-function Ψ2 is assumed as being the one true ground state wavefunction, the analogue inequality springs out as:
E 2 [ ρ ] < E 1 [ ρ ] + ρ ( r ) [ V 2 ( r ) V 1 ( r ) ] d r.
Taken together relations (11) and (12) generate, by direct summation, the evidence of the contradiction [2]:
E 1 [ ρ ] + E 2 [ ρ ] < E 1 [ ρ ] + E 2 [ ρ ] .
The removal of such contradiction could be done in a single way, namely, by abolishing, in a reverse phenomenologically order, the fact that two eigen-functions, two Hamiltonians and respectively, two external potential exist for characterizing the same ground state of a given electronic system. With this statement the HK2 theorem is formally proofed.
Yet, there appears the so called V-representability problem signaling the impossibility of an a priori selection of the external potentials types that are in bi-univocal relation with ground state of an electronic system [1518]. The problem was revealed as very difficult at mathematical level due to the equivocal potential intrinsic behavior that is neither of universal nor of referential independent value. Fortunately, such principial limitation does not affect the general validity of the variational principle (9) regarding the selection of the energy of ground state level from a collection of states with different associated external potentials.
That because, the problem of V-representability can be circumvented by the so called N-contingency features of ground state electronic density assuring that, aside of the N – integrability condition (4), the candidate ground state densities should fulfill the positivity condition (an electronic density could not be negative) [17, 18]:
ρ ( r ) 0 ,    | r | ,
as well as the non-divergent integrability condition on the real domain (in relation with the fact that the kinetic energy of an electronic system could not be infinite – since the light velocity restriction):
| ρ ( r ) 1 / 2 | 2 d r < .
Both (14) and (15) conditions are easy accomplished by every reasonable density, allowing the employment of the variational principle (9) in two steps, according to the so called Levy-Lieb double minimization algorithm [19]: one regarding the intrinsic minimization procedure of the energetic terms respecting all possible eigen-functions folding a trial electronic density followed by the external minimization over all possible trial electronic densities yielding the correct ground state (GS) energy density functional
E G S = min ρ [ min Ψ ρ ( Ψ * ( T + V e e + V ) Ψ d τ ) ] = min ρ [ min Ψ ρ ( Ψ * ( T + V e e ) Ψ d τ ) + ρ ( r ) V ( r ) d r ] min ρ ( F H K [ ρ ] + C A [ ρ ] ) = min ρ ( E [ ρ ] ) .
One of the most important consequences of the HK2 conveys the rewriting of the variational principle (9) in the light of above N-contingency conditions of the trial densities as the working Euler type equation:
δ { E [ ρ ] μ N [ ρ ] } = 0
from where, there follows the Lagrange multiplication factor with the functional definition:
μ = ( δ E [ ρ ] δ ρ ) ρ = ρ ( V )
this way introducing the chemical potential as the fundamental quantity of the theory. At this point, the whole chemistry can spring out since identifying the electronic systems electronegativity with the negative of the density functional chemical potential [9]:
χ = μ .
Up to now, the Hohenberg-Kohn theorems give new conceptual quantum tools for physico-chemical characterization of an electronic sample by means of electronic density and its functionals, the total energy and chemical potential (electronegativity). These positive density functional premises are in next analyzed towards elucidating of the quantum nature of the chemical bond and reactivity [20].

2.2. Optimized energy-electronegativity connection

Back from Paris, in the winter of 1964, Kohn met at the San Diego University of California his new post-doc Lu J. Sham with who propose to extract from HK1 & 2 theorems the equation of total energy of the ground state. In fact, they propose themselves to find the correspondent of the stationary eigen-equation of Schrödinger type, employing the relationship between the electronic density and the wave function.
Their basic idea consists in assuming a so called orbital basic set for the N-electronic system by replacing the integration in the relation (5) with summation over the virtual uni-electronic orbitals ϕ i , i = 1 , N ¯, in accordance with Pauli principle, assuring therefore the HK1 frame with maximal spin/orbital occupancy [21]:
ρ ( r ) = Σ i N n i | φ i ( r ) | 2 ,   0 ni 1 , Σ i ni = N .
Then, the trial total eigen-energy may be rewritten as density functional of eq. (6) nature expanded in the original form [22, 23]:
E [ ρ ] = F H K [ ρ ] + C A [ ρ ] = T [ ρ ] + V e e [ ρ ] + C A [ ρ ] = T s [ ρ ] + j [ ρ ] + { ( T [ ρ ] T s [ ρ ] ) + ( V e e [ ρ ] J [ ρ ] ) } + C A [ ρ ] = Σ i N n i φ i * ( r ) [ 1 2 2 ] φ i ( r ) d r + 1 2 ρ ( r 1 ) ρ ( r 2 ) r 12 d r 1 d r 2 + E x c [ ρ ] + V ( r ) ρ ( r ) d r
where, the contribution of the referential uniform kinetic energy contribution
T s [ ρ ] = Σ i N n i φ i * ( r ) [ 1 2 2 ] φ i ( r ) d r ,
with the inferior index „s” referring to the „spherical” or homogeneous attribute together with the classical energy of Coulombic inter-electronic repulsion
J [ ρ ] = 1 2 ρ ( r 1 ) ρ ( r 2 ) r 12 d r 1 d r 2
were used as the analytical vehicles to elegantly produce the exchange-correlation energy Exc containing exchange (Vee[ρ] – J[ρ]) and correlation (T[ρ] – Ts [ρ]) heuristically introduced terms as the quantum effects of spin anti-symmetry over the classical interelectronic potential and of corrected homogeneous electronic movement, respectively.
Next, the trial density functional energy (21) will be optimized in the light of variational principle (17) as prescribed by the HK2 theorem. The combined result of the HK theorems will eventually furnish the new quantum energy expression of multi-electronic systems beyond the exponential wall of the wave function.
An instructive method for deriving such equation assume the same types of orbitals for the density expansion (20),
ρ ( r ) = N φ * ( r ) φ ( r )
that, without diminishing the general validity of the results, since preserving the N-electronic character of the system, highly simplifies the analytical discourse.
Actually, with the trial density (24) replaced throughout the energy expression in (21) has to undergo the minimization procedure (17) with the practical equivalent integral variant:
δ ( E [ ρ ] μ N [ ρ ] ) δ φ * δ φ * d r = 0 .
Note that, in fact, we chose the variation in the conjugated uni-orbital φ*(r) in (25) providing from (24) the useful differential link:
δ ρ ( r ) = N φ ( r ) δ φ * ( r ) .
Now, unfolding the equation (25) with the help of relations (21) and (24), together with fundamental density functional prescription (4), one firstly gets [24]:
δ δ φ * ( r ) { N 2 φ * ( r ) 2 φ ( r ) d r + J [ ρ ] + E x c [ ρ ] + N V ( r ) φ * ( r ) φ ( r ) d r μ N φ * ( r ) φ ( r ) d r } = 0
By performing the required partial functional derivations respecting the uni-orbital φ*(r) and by taking account of the equivalence (26) in derivatives relating J[ρ] and Exc [ρ] terms, equation (27) takes the further form:
N 2 2 φ ( r ) + N φ ( r ) δ J [ ρ ] δ ρ + N φ ( r ) δ E x c δ ρ + N V ( r ) φ ( r ) μ N φ ( r ) = 0.
After immediate suppressing of the N factor in all the terms and by considering the exchange-correlation potential with the formal definition [4, 9]:
V x c ( r ) = ( δ E x c [ ρ ] δ ρ ( r ) ) V ( r ) ,
equation (28) simplifies as [25, 26]:
[ 1 2 2 + ( V ( r ) + ρ ( r 2 ) | r r 2 | d r 2 + V x c ( r ) ) ] φ ( r ) = μ φ ( r ) .
Moreover, once introducing the so called effective potential:
V e f f ( r ) = V ( r ) + ρ ( r 2 ) | r r 2 | d r 2 + V x c ( r )
the resulted equation recovers the traditional Schrödinger shape:
[ 1 2 2 + V e f f ] φ ( r ) = μ φ ( r ) .
The result (32) is fundamental and equally subtle. Firstly, it was proved that the joined Hohenberg-Kohn theorems are compatible with consecrated quantum mechanical postulates, however, still offering a generalized view of the quantum nature of electronic structures, albeit the electronic density was assumed as the foreground reality. In these conditions, the meaning of functions φ(r) is now unambiguously producing the analytical passage from configuration (3N-D) to real (3D) space for the whole system under consideration. Nevertheless, the debate may still remain because once equation (32) is solved the basic functions φ(r) generating the electronic density (24) and not necessarily the eigen-functions of the original system due to the practical approximations of the exchange and correlation terms appearing in the effective potential (31). This is why the functions φ(r) are used to be called as Kohn-Sham (KS) orbitals; they provide the orbital set solutions of the associate KS equations [3]:
[ 1 2 2 + V e f f ] φ i ( r ) = μ i φ i ( r ) ,   i = 1 , N ¯
once one reconsiders electronic density (24) back with general case (20).
Yet, equations (33), apart of delivering the KS wavefunctions φi(r), associate with another famous physico-chemical figure, the orbital chemical potential μi, which in any moment can be seen as the negative of the orbital electronegativities on the base of the relation (19).
Going now to a summative characterization of the above optimization procedure worth observing that the N-electronic in an arbitrary external V-potential problem is conceptual-computationally solved by means of the following self-consistent algorithm:
  • It starts with a trial electronic density (20) satisfying the N-contingency conditions (14) and (15);
  • With trial density the effective potential (31) containing exchange and correlation is calculated;
  • With computed Veff the equations (33) are solved for ϕ i ( r ) , i = 1 , N ¯;
  • With the set of functions { ϕ i ( r ) } i = 1 , N ¯ the new density (20) is recalculated;
  • The procedure is repeated until the difference between two consecutive densities approaches zero;
  • Once the last condition is achieved one retains the last set { ϕ i ( r ) , μ i = χ i } i = 1 , N ¯;
  • The electronegativity orbital observed contributions are summed up from (33) with the expression:
    - Σ i N χ i = Σ i N n i φ i * ( r ) [ 1 2 2 + V e f f ( r ) ] φ i ( r ) d r = T s [ ρ ] + V e f f ρ ( r ) d r ;
  • Replacing in (34) the uniform kinetic energy, Ts [ρ] from the general relation (21) the density functional of the total energy for the N-electronic system will take the final figure [9, 24]:
    E [ ρ ] = Σ i N χ i 1 2 ρ ( r 1 ) ρ ( r 2 ) | r 12 | d r 1 d r 2 + { E x c [ ρ ] V x c ( r ) ρ ( r ) d r } .
showing that the optimized many-electronic ground state energy is directly related with global or summed over observed or averaged or expected orbital electronegativities. One can observe from (35) that even in the most optimistic case when the last two terms are hopefully canceling each other there still remains a (classical) correction to be added on global electronegativity in total energy. Or, in other terms, electronegativity alone is not enough to better describe the total energy of a many-electronic system, while its correction can be modeled in a global (almost classical) way. Such considerations stressed upon the accepted semiclassical behavior of the chemical systems, at the edge between the full quantum and classical treatments.
However, analytical expressing the total energy requires the use of suitable approximations, whereas for chemical interpretation of bonding the electronic localization information extracted from energy is compulsory. This subject is in next focused followed by a review of the popular energetic density functionals and approximations.

3. Electronic Localization Problem

3.1. From global functional to localization function. Localization in solids

The application of Hohenberg-Kohn theorems consecrate the crucial contributions of the so called “spherical” or homogeneous kinetic and of the exchange-correlation energy terms in a multi-electronic system’s ground state. However, the spherical electronic case corresponds to the non perturbed electronic system for which the Thomas-Fermi (TF) model was already advanced throughout totally ignoring exchange- correlation terms from the total energy shape:
E T F [ ρ ] = T T F [ ρ ] + J [ ρ ] + C A [ ρ ] .
Such a referential picture is most useful in establishing the uniform electronic distribution by indicating the occupation of the all-possible electronic levels in a semiclassical quantum frame (without explicit exchange-correlation involvement). Actually, the Fermi sphere in a momentum space finely defines the total homogeneous kinetic energy as:
τ s ( r ) = p F 2 2 m 0 ,
while the quantum nature of the kinetic energy (37) is covered by involving the quantum (Heisenberg) uncertainty
Δ p x Δ p y Δ p z Δ x Δ y Δ z h 3
in uniform density computation. This suggests that the density of states in the Fermi volume of the impulse pF has to be normalized to the inverse of the cube power of Planck constant h, while the density of electrons is reached by multiplying the density of states with the electron multiplicity 2(1/2)+1=2 for every occupied state. The obtained density-Fermi impulse relationship:
ρ ( r ) d r = 2 4 3 π p F 3 h 3 d r p F ( r ) = ( 3 8 π ) 1 / 3 h ρ ( r ) 1 / 3
allows the Thomas-Fermi kinetic energy unfolding as the density functional [9, 2731]:
T T F [ ρ ] = 3 5 ρ ( r ) τ s ( r ) d r = C T F ρ ( r ) 5 / 3 d r ,   C T F = 3 h 2 10 m 0 ( 3 8 π ) 2 / 3 2.871 [ a . u . ]
with the help of which the total Thomas-Fermi energy functional takes the form:
E T F [ ρ ] = C T F ρ ( r ) 5 / 3 d r + 1 2 ρ ( r 1 ) ρ ( r 2 ) | r 12 | d r 1 d r 2 + V ( r ) ρ ( r ) d r
that can be seen as the first approximation for the density functional total energy (21).
From physical point of view worth noted that the kinetic TF energy exactly corresponds to the total energy of the free electrons in a crystal, V(r)=0 in (41), equivalently with the fact that the electrons are not “feeling” the nuclei, i.e. electrostatic attractions are excluded, being as close each other to avoid reciprocal repelling. Such picture suggests that free electrons are completely non-localized leaving with the condition of complete cancellation of the electronic inter-repulsion; this feature may be putted formally as [32]:
e 2 | r 1 r 2 | λ e 2 | r 1 r 2 | , λ = 0.
However, the model in which the (valence) electrons are completely free and are neither “feeling” the attraction nor the repulsion is certain not properly describing the nature of the chemical bond. In fact, this limitation was also the main objection brought to Thomas-Fermi model and to the atomic or molecular approximation of the homogeneous electronic gas or jellium model in solids. Nevertheless, the lesson is well served because Thomas-Fermi description may be regarded as the “inferior” extreme in quantum known structures while further exchange-correlation effects may be added in a perturbative manner.
The idea of introducing exchange and correlation effects as a perturbation of the homogeneous electronic system could be considered from the interpolation of the energetic terms for 0≤λ≤1 in (41). Parameter λ is defined as a parameter of the electronic coupling, with a slightly (adiabatically) scaling of the perturbation from the homogeneous electronic systems, λ=0, until the maximal inter-electronic interaction, λ=1 (in accordance with Pauli principle). Therefore, the overall interpolation 0 1 [ ] d λ will be spread over the terms which contain the intermediate degree of exchange and correlation interactions; since it accounts for the electronic inter-repulsion while indexing the electronic presence/absence in a given spatial region the degree of electronic localization is in this way furnishes.
The coupling parameter λ will serve as a switcher between the referential Thomas-Fermi uniform case and the full interaction through the density limit:
lim λ 1 ρ λ ( r ) = ρ ( r = r 1 r 2 ) .
Actually, the density (43) has a major role in defining exchange-correlation functionals. To see that let’s firstly consider the conditional electronic density g(r1, r2;λ) indicating that the electronic density in r1 is conditioned by the presence (localization) of another electron (any from the total N in the system) in r2. Mathematically, this is expressed by using the conditional probabilities:
g ( r 1 , r 2 ; λ ) = ρ ( r = r 1 r 2 ) ρ ( r 2 )
fulfilling the Pauli principle by means of the integration rule:
ρ ( r 2 ) g ( r 1 , r 2 ; λ ) d r 2 = ρ ( r 1 r 2 ) d r 2 = 0
saying that the spatial average of the electronic reciprocal constraint vanishes. This behavior opens the possibility in introducing the conditional probability of electronic holes,
h ( r 1 , r 2 ; λ ) = g ( r 1 , r 2 ; λ ) 1 ,
providing the associate integration rule [33]:
ρ ( r 2 ) h ( r 1 , r 2 ; λ ) d r 2 = 1
consecrating a sort of negative normalization of the exchange and correlation density of holes:
ρ x c ( r 1 , r 2 ; λ ) = ρ ( r 2 ) h ( r 1 , r 2 ; λ ) .
Now, once this exchange-correlation hole density is mediated over the coupling factor λ the averaged exchange-correlation density of holes is generated:
ρ x c ¯ ( r 1 , r 2 ) = 0 1 ρ x c ( r 1 , r 2 ; λ ) d λ
allowing the formal writing of exchange-correlation density functional from (29) as a generalized version of the inter-electronic interaction term (23):
E x c [ ρ ] = 1 2 ρ ( r 1 ) ρ x c ¯ ( r 1 , r 2 ) r 12 d r 1 d r 2 = 1 2 ρ ( r 1 ) R x c ¯ ( r 1 , ρ ( r ) ) d r 1
with the help of the introduced radius of the λ-mediated exchange-correlation density of holes:
R x c ¯ 1 ( r 1 , ρ ( r ) ) : = ρ x c ¯ ( r 1 , r 2 , ρ ( r ) ) r 12 d r 2
The radius R x c ¯ ( r ) could be considered as a functional of density (43) with the leading term being defined in the short limit of the distance, i.e. being of the inter-particle average radius order,
4 π 3 r 0 3 = 1 ρ ( r ) r 0 = ( 3 4 π ) 1 / 3 ρ ( r ) 1 / 3 ,
known as the Wigner radius for indexing the volume of a sphere containing (localizing) a single electron (from the total of N) belonging to the density family (43), although, also other quantities accounting for electron localization such as the domain averaged Fermi hole of Ponec [34a] or the electron sharing index (also known as delocalization index) [34b] have been recently proposed (see also the forthcoming discussion).
In these conditions, the inverse radius (51) could be expressed around the inverse of the Wigner radius in a gradient density expansion [32]:
R x c ¯ 1 ( r , ρ ( r ) ) = F 0 [ ρ ( r ) ] + F 21 [ ρ ( r ) ] 2 ρ ( r ) + F 22 [ ρ ( r ) ] i = 1 N ( i ρ ( r ) ) ( i ρ ( r ) ) + ...
while, by considering it back in exchange-correlation energy (50) produces, after the integration by parts, the generalized gradient density functional:
E x c = G 0 [ ρ ( r ) ] ρ ( r ) d r + G 2 [ ρ ( r ) ] ( ρ ( r ) ) 2 d r + [ G 4 [ ρ ( r ) ] ( 2 ρ ( r ) ) 2 + ... ] d r + ...
The restriction to the first term of the series (54) corresponds to the cases where the spatial distance of variation in electronic density highly exceeds the corresponding Wigner radius (52) this way producing the famous local density approximation (LDA) [34c]:
E x c L D A [ ρ ] = e x c [ ρ ( r ) ] ρ ( r ) d r
with exc being the exchange and correlation density per particle, that can be further approximated (see bellow) as [3540]:
e x c [ ρ ] = e x [ ρ ] + e c [ ρ ] = ( 0.458 r 0 ) + ( 0 . 44 r 0 + 7.8 ) [ a . u . ] .
In fact, the LDA stands as the immediate step after TF approximation; it can be extended also for systems with un-pair spins by the so called local spin density approximation (LSDA) [4150]:
E x c [ ρ = ρ + ρ ] = E x c [ ρ ] + E x c [ ρ ] ,
while further inclusion of the gradient terms in (54) establishes the general gradient approximation (GGA).
Worth noting that when undertaken GGA, beside the gradient terms arising in exchange-correlation energy, the gradient correction of the kinetic energy functional has to be as well considered providing terms of which the standard one takes the von Weizsäcker form [51, 52]:
τ W ( r ) = 1 8 | ρ ( r ) | 2 ρ ( r ) .
While more analytical discussions about various approximations and density functionals are bellow presented in a separate chapter, here we would like only to present the practical difference between the local and gradient density approximations for a solid state case.
For instance, Figure 1 presents the band structure and the density of states (DOS) for the R 3 ¯m oxide of the Cobalt transitional metal (CoO) calculated with either LSDA and GGA approximations [53].
Regarding the energy bands there can be noted that, around the Fermi level EF, LSDA approach is less relevant in indicating the energetically gap respecting the GGA computation. The difference is even more drastic in DOS representations for employed approximations in d orbital separation (t2g=a1g+eg’) due to the central ion of cobalt trigonal symmetry coordination. In fact, with LSDA a strong mixing of the orbitals a1g and eg’ is recorded while in the case of GGA-DOS the bands with the symmetry a1g are up and down shifted for the respective down and up spin projections resulting a clear separation from states with eg’ symmetry.
Nonetheless, at the level of bands structure of the solids and crystals an inevitable localization paradox emerges namely, to use the real 3D electronic densities in furnish a localization description in the reciprocal (energy) space.
Recently, it was found a way to avoid the electronic localization paradox through introducing specific electronic localization functions (ELF) in real space. Nevertheless, an ELF should relay on combination of the gradient and homogeneous energetic density functionals, in accordance with Pauli principle, shaping for instance as [54]:
E L F = ( 1 + [ τ s ( r ) τ W ( r ) T T F [ ρ ( r ) ] ] 2 ) 1
by emphasizing the excess of the kinetic energy difference τs(r)- τW(r) „normalized” to the referential kinetic TF homogeneous behavior.
Worth remarking that the localization function (59) acts like a sort of density, with values between 0 and 1 corresponding with maximum delocalization and localization, respectively. This heuristically proposal has the merit to give an analytical reflection of the qualitative valence shell electron pair repulsion (VSEPR) geometric model [56], with the immediate consequence in topological characterization of the chemical bond [57]. In solid state case, the reliability of above ELF in describing the chemical bond in real space is illustrated in Figure 2 for the Li and Sc crystals. Atomic and molecular levels are in next section illustrated with which occasion further ELF characterization and developments are presented.

3.2. Localization in atoms and molecules

The definitions that are currently used in the classification of chemical bonds are often imprecise, as they are derived from approximate theories. Based on the topological analysis local, quantum-mechanical functions related to the Pauli Exclusion Principle may be formulated as “localization attractors” of bonding, non-bonding, and core types. Bonding attractors lie between nuclei core attractors and characterize shared electron interactions. The spatial arrangement of bond attractors allows for an absolute classification of ionic versus covalent bond to be derived from electronic density combined functions [58].
Most modern classifications of the chemical bond are based on Lewis’ theory and rely on molecular-orbital and valence-bond theories with schemes involving linear combination of atomic orbitals (LCAO). However, electron density alone does not easily reveal the consequences of the Pauli Exclusion Principle on bonding nature. While VSEPR theory indicates that the Pauli principle is important for understanding chemical structures, it has been reformulated in terms of maxima of electronic density’s Laplacian −∇2ρ(r) [56]. Next, the exchange-correlation density functional concept was employed to achieve the coordinate-space dynamical correlation in an inhomogeneous electron gas. This way the exchange-correlation energy (50) further re-expresses like [32, 54]
E x c = 1 2 α α ρ ( r 1 ) ρ ( r 2 ) r 12 h X C α α d r 1 d r 2 ,
by accounting for the four α-spin types of interactions through the “hole” functions [40, 59]
h α α ( r 1 , r 2 ) = P 2 α α ( r 1 , r 2 ) ρ α ( r 1 ) ρ α ( r 2 )
where P2(r1,r2) stands for the two-body or pair probability density or correlation probability of the arbitrary electrons „1 and 2”, defined in terms of the N-body wave function Ψ as follows [9]:
P 2 ( r 1 , r 2 ) = N ( N 1 ) ... Ψ * ( r 1 , r 2 , r 3 , ... , r N ) Ψ ( r 1 , r 2 , r 3 , ... , r N ) d r 3 ... d r N   .
Within the density functional theory the electrons of a pair of electrons or a bond can be considered as belonging to an inhomogeneous continuum gas. In analytical terms this was translated as the ELF (59) index as combining the homogeneous and inhomogeneous behaviors of a many-electronic-nuclei system.
Nevertheless, recently, the Markovian analytical shape of an ELF was shown to have the general qualitative form [60]:
E L F = f 1 ( g ( r ) h ( r ) )
with the limiting constrains
lim E L F = { 0 , ρ ( r )     ρ ( r ) 1 , ρ ( r )     ρ ( r )
assuring the fulfillment of the Heisenberg and Pauli principles. To clarify this [60], we make recourse to the Heisenberg principle, comprised in ELF. When density gradient dominates,∇ρ≫ρ then g≫h in (63), and f(∞) should accounted for the infinite error in assigning of momentum, therefore indicating a precisely spatial localization of electrons; thus f(∞) = ∞ and ELF→0. In such, the meaning of ELF is associated with the error in spatial localization of electrons, being zero when the electrons are precisely located. On the contrary, when ρ≫∇ρ then h≫g in (63), and the resulting f(0) indicates the minimum error in defining of momentum and should provide the maximum uncertain of spatial distribution; in such f(0)=1 and ELF→1, where 1 stands here for 100% of coordinate localization error.
In this context, when the inverse of difference in local kinetic terms is involved, the ELF is interpreted as the error in localization of electrons within traps rather than where they have peaks of spatial density, as is frequently misinterpreted in literature [61], albeit recent extensions of ELF have used the correlated (Hartree-Fock) wavefunctions, through the conditional pair probability, however not using the “kinetic energy approach” [6264].
Among various classes of Markovian ELFs the most representative and efficient one was proposed as having the form [60]:
E L F = exp { 3 2 [ g ( r ) h ( r ) ] 2 }
with the components:
g ( r ) = 1 2 i [ φ i ( r ) ] 2 1 8 [ ρ ( r ) ] 2 ρ ( r ) , h ( r ) = 3 10 ( 3 π 2 ) 2 / 3 [ ρ ( r ) ] 5 / 3
being responsible for the gradient (g) and the homogenous (h) density distributions, respectively. In this frame, the ELF information prescribe that as it has values closer to zero as the better electronic localization is providing, according with the limits (64).
Going to a particular application of this scheme the atomic level is firstly presented for the special case of Li atom. The main stages consist in:
  • Choosing the basis of the atomic functions [65]:
    f 1 L i ( r ) = 8.863248 r exp ( 2 . 698 r ) , f 2 L i ( r ) = 0.369721 r 5 / 2 exp ( 0 . 797 r )
such that to fulfill the natural (radial) normalization conditions
0 [ f n L i ( r ) ] 2 d r = 1 , n = 1 , 2
  • Generating the orthonormal orbital eigen-waves, here according with the Gram-Schmidt algorithm among shells and sub-shells:
    φ 1 s L i ( r ) = f 1 L i ( r ) , φ 2 p L i ( r ) = f 2 L i ( r ) , φ 2 s L i ( r ) = C 2 s L i [ φ 2 p L i ( r ) α φ 1 s L i ( r ) ] = 1 . 2226 exp ( 2 . 698 r ) + 0.373222 r 5 / 2 exp ( 0.797 r )
ensuring the additionally constraints:
0 φ 2 s L i ( r ) φ 1 s L i ( r ) d r = 0 ,    0 [ φ 2 s L i ( r ) ] 2 d r = 1 .
  • Generating the working overall electronic density
    ρ L i ( r ) = 2 [ φ 1 s L i ( r ) ] 2 + [ φ 2 s L i ( r ) ] 2
that satisfies the spatial (radial) global N-integral condition (4):
0 ρ L i ( r ) d r = 3
The electronic density (71) is then used for computation of the Markovian ELF (65), while their comparison is in Figure 3 illustrated.
From the Figure 3, there appears that the smooth delocalization of electrons of Li represented by density structure is removed by the electronic localization function by clearly indicating where are the regions where the electronic realm is with less uncertainty detected. This way the ELF indicates merely where the electronic transitions behave like a step-function. In this respect, ELF can be regarded as the complement of electronic density being a better indicator of the regions where the bonding may arise. For instance, in the case of Li atomic structure, the fact that the ELF does not displays localization over the second shell (due to its values approaching unity in this range) indicates a natural tendency for releasing the outermost electron to the (virtual) neighborhood atoms with uncompleted last shell(s) while preserving its delocalization feature across the bond. As such the lithium hydride (LiH) bond is expected to be formed with a certain degree of ionicity in resonance with its covalence: LiH ↔ Li+H.
The reliability of ELF to quantify the local tendency of atoms to form bonds and aggregates can be further exemplified to diatomic molecules, while the particular cases of HF, HCl, HBr and HI structures are considered in Figure 4. In the bonding region, i.e. in the space between the hydrogen and halogen atomic centers in H-X molecules there are represented both the electron densities, computed upon above recipe [65], and the associate Markovian ELFs for the concerned atoms-in-molecules (AIM).
Figure 4 clearly shows that while the crossing of hydrogen and halogen radial densities does not provide the right bonding region in HCl, HBr, and HI cases, the corresponding ELFs cross-lines of AIM finely indicate the frontier of atomic basins in hydracids thus confirming the ELF reliability in identifying chemical bonds and bonding.
One can equally say that in the crossing vicinity of AIM-ELFs the electrons are at the same time completely localized (for bonding with ELFX – ELFH →0) and completely delocalized for atomic systems (with ELFX,H →1), according to above the ELF definition and present signification.
In other words it can be alleged that ELF application on chemical bond helps in identifying the molecular region in which the electrons undergo the transition from the complete delocalization in atoms to complete localization in molecular bonding behavior.
Actually, it also proves that localization issue of ionic and covalent classification of bonds may be solved by a “continuous” quantum reality. Such a feature gives, nevertheless, an in-depth understanding of the quantum nature of the chemical bond by associating the mysterious pairing of electrons issue to an analytical function able to distinguish the narrow regions of molecular space where the Heisenberg and Pauli principles are jointly satisfied through the ELF’s extreme values. Even more such sharp differentiation between 0 and 1 in atomic and molecular ELF values offers the future possibility in quantifying the chemical bond and bonding in the frame of quantum information theory [67].

4. Popular Energetic Density Functionals

Since the terms of total energy are involved in bonding and reactivity states of many-electronic systems, i.e. the kinetic energetic terms in ELF topological analysis or the exchange and correlation density functionals in chemical reactivity in relation with either localization and chemical potential or electronegativity, worth presenting various schemes of quantification and approximation of these functionals for better understanding their role in chemical structure and dynamics.

4.1. Density functionals of kinetic energy

When the electronic density is seen as the diagonal element ρ(r1) = ρ(r1,r1) the kinetic energy may be generally expressed from the Hartree-Fock model, through employing the single determinant ρ(r1,r1), as the quantity [68]:
T [ ρ ] = 1 2 [ r ' 1 2 ρ ( r 1 , r ' 1 ) ] r 1 = r ' 1 d r 1 ;
it may eventually be further written by means of the thermodynamical (or statistical) density functional [69]:
T β = 3 2 ρ ( r ) k B T ( r ) d r = 3 2 ρ ( r ) 1 β ( r ) d r
that supports various specializations depending on the statistical factor particularization β.
For instance, in LDA approximation, the temperature at a point is assumed as a function of the density in that point, β(r) = β(ρ(r)); this may be easily reached out by employing the scaling transformation to be [70]
ρ λ ( r ) = λ 3 ρ ( λ r ) T [ ρ λ ] = λ 2 T [ ρ ] , λ = c t . ,
providing that
β ( r ) = 3 2 C ρ 2 / 3 ( r ) ,
a result that helps in recovering the traditional (Thomas-Fermi) energetic kinetic density functional form
T [ ρ ] = C ρ 5 / 3 ( r ) d r ,
while the indeterminacy remained is smeared out in different approximation frames in which also the exchange energy is evaluated. Note that the kinetic energy is generally foreseen as having an intimate relation with the exchange energy since both are expressed in Hartree-Fock model as determinantal values of ρ(r1,r1), see below.
Actually, the different LDA particular cases are derived by equating the total number of particle N with various realization of the integral
N = 1 2 | ρ ( r 1 , r ' 1 ) | 2 d r 1 d r 1 '
by rewritting it within the inter-particle coordinates frame:
r = 0.5 ( r 1 + r 1 ' ) , s = r 1 r 1 '
as:
N = 1 2 | ρ ( r + s / 2 , r s / 2 ) | 2 d r d s
followed by spherical averaged expression:
N = 2 π ρ 2 ( r ) Γ ( r , s ) d r s 2 d s
with
Γ ( r , s ) = 1 s β ( r ) +
The option in choosing the Γ(r, s) series (82) so that to converge in the sense of charge particle integral (81) fixes the possible cases to be considered [68]:
  • the Gaussian resummation uses:
    Γ ( r , s ) Γ G ( r , s ) = exp ( s 2 β ( r ) )
  • the trigonometric (uniform gas) approximation looks like:
    Γ ( r , s ) Γ T ( r , s ) = 9 ( sin t t cos t ) 2 t 6 , t = s 5 β ( r )
In each of (83) and (84) cases the LDA-β function (76) is firstly replaced; then, the particle integral (81) is solved to give the constant C and then the respective kinetic energy density functional of (77) type is delivered; the results are [68]:
  • in Gaussian resummation:
    T G L D A = 3 π 2 5 / 3 ρ 5 / 3 ( r ) d r ,
  • whereas in trigonometric approximation
    T T F L D A = 3 10 ( 3 π 2 ) 2 / 3 ρ 5 / 3 ( r ) d r
one recovers the Thomas-Fermi formula type that closely resembles the original TF (40) formulation.
In next one will consider the non-local functionals; this can be achieved through the gradient expansion in the case of slowly varying densities – that is assuming the expansion [52]:
T = d r [ τ ( ρ ) + τ ( ρ ) ] = d r Σ m = 0 [ τ 2 m ( ρ ) + τ 2 m ( ρ ) ] = d r Σ m = 0 τ 2 m ( ρ ) = d r τ ( ρ )
The first two terms of the series respectively covers: the Thomas Fermi typical functional for the homogeneous gas
τ 0 ( ρ ) = 3 10 ( 6 π 2 ) 2 / 3 ρ 5 / 3
and the Weizsäcker related first gradient correction:
τ 2 ( ρ ) = 1 9 τ W ( ρ ) = 1 72 | ρ | 2 ρ .
They both correctly behave in asymptotic limits:
τ ( ρ ) = { τ 0 ( ρ ) = τ 2 ( ρ ) ρ ( far    from    nucleus ) 9 τ 2 ( ρ ) = τ W ( ρ ) = 1 8 | ρ | 2 ρ ρ ( close    to    nucleus )
However, an interesting resummation of the kinetic density functional gradient expansion series (87) may be formulated in terms of the Padé-approximant model [71]:
τ ( ρ ) = τ o ( ρ ) P 4 , 3 ( x )
with
P 4 , 3 ( x ) = 1 + 0.95 x + a 2 x 2 + a 3 x 3 + 9 b 3 x 4 1 0.05 x + b 2 x 2 + b 3 x 3
and where the x-variable is given by
x = τ 2 ( ρ ) τ 0 ( ρ ) = 5 108 1 ( 6 π 2 ) 2 / 3 | ρ | 2 ρ 8 / 3 ,
while the parameters a2, a3, b2, and b3 are determined by fitting them to reproduce Hartree-Fock kinetic energies of He, Ne, Ar, and Kr atoms, respectively [72]. Note that Padé function (92) may be regarded as a sort of generalized ELF susceptible to be further used in bonding characterizations.

4.2. Density functionals of exchange energy

Starting from the Hartree-Fock framework of exchange energy definition in terms of density matrix [73],
K [ ρ ] = - 1 4 | ρ ( r 1 , r ' 1 ) | 2 | r 1 r ' 1 | d r 1 d r 1 ' ,
within the same consideration as before, we get that the spherical averaged exchange density functional
K = π ρ 2 ( r ) Γ ( r , s ) d r s d s
takes the particular forms [68]:
  • in Gaussian resummation:
    K G L D A = 1 2 1 / 3 ρ 4 / 3 ( r ) d r ;
  • and in trigonometric approximation (recovering the Dirac formula):
    K G L D A = 3 4 ( 3 π ) 1 / 3 ρ 4 / 3 ( r ) d r .
Alternatively, by paralleling the kinetic density functional previous developments the gradient expansion for the exchange energy may be regarded as the density dependent series [74]:
K = Σ n = 0 K 2 n ( ρ ) = d r Σ n = 0 k 2 n ( ρ ) = d r k ( ρ )
while the first term reproduces the Dirac LDA term [75, 76]:
k 0 ( ρ ) = 3 2 ( 3 4 π ) 1 / 3 ρ 4 / 3
and the second term contains the density gradient correction, with the Becke proposed approximation [77]:
k 2 ( ρ ) = b | ρ | 2 ρ 4 / 3 ( 1 + d | ρ | 2 ρ 8 / 3 ) a
where the parameters b and d are determined by fitting the k0+k2 exchange energy to reproduce Hartree-Fock counterpart energy of He, Ne, Ar, and Kr atoms, and where for the a exponent either 1.0 or 4/5 value furnishes excellent results. However, worth noting that when analyzing the asymptotic exchange energy behavior, we get in small gradient limit [77]:
k ( ρ ) ρ k 0 ( ρ ) 7 432 π ( 6 π 2 ) 1 / 3 | ρ | 2 ρ 4 / 3
whereas the adequate large-gradient limit is obtained by considering an arbitrary damping function as multiplying the short-range behavior of the exchange-hole density, with the result:
k ( ρ ) ρ c ρ 4 / 5 | ρ | 2 / 5
where the constant c depends of the damping function choice.
Next, the Padé-resummation model of the exchange energy prescribes the compact form [74]:
k ( ρ ) = 10 9 k 0 ( ρ ) P 4 , 3 ( x )
with the same Padé-function (92) as previously involved when dealing with the kinetic functional resummation. Note that when x=0, one directly obtains the Ghosh-Parr functional [78]:
k ( ρ ) = 10 9 k 0 ( ρ ) .
Moreover, the asymptotic behavior of Padé exchange functional (103) leaves with the convergent limits:
k ( ρ ) = { 10 9 ( k 0 + 15 17 7 432 π ( 6 π 2 ) 1 / 3 | ρ | 2 ρ 4 / 3 ) x 0 ( SMALL    GRADIENTS ) 12 π ρ 2 | ρ | 2 x ( LARGE    GRADIENTS )
Once again, note that when particularizing small or large gradients and fixing asymptotic long or short range behavior, we are recovering the various cases of bonding modeled by the electronic localization recipe as provided by ELF’s limits (64).
Another interesting approach of exchange energy in the gradient expansion framework was given by Bartolotti through the two-component density functional [79]:
K [ ρ ] = C ( N ) ρ ( r ) 4 / 3 d r + D ( N ) r 2 | ρ | 2 ρ 2 / 3 d r ,
where the N-dependency is assumed to behave like:
C ( N ) = C 1 + C 2 N 2 / 3 , D ( N ) = D 2 N 2 / 3
while the introduced parameters C1, C2, and D2 were fond with the exact values [8082]:
C 1 = 3 4 π 1 / 3 , C 2 = 3 4 π 1 / 3 [ 1 ( 3 π 2 ) 1 / 3 ] , D 2 = π 1 / 3 729 .
Worth observing that the exchange Bartolotti functional (106) has some important phenomenological features: it scales like potential energy, fulfills the non-locality behavior through the powers of the electron and powers of the gradient of the density, while the atomic cusp condition is preserved [83].
However, density functional exchange-energy approximation with correct asymptotic (long range) behavior, i.e. satisfying the limits for the density
lim r ρ σ = exp ( a σ r )
and for the Coulomb potential of the exchange charge, or Fermi hole density at the reference point r
lim r U X σ = 1 r , σ = α ( o r ) , β ( o r ) s p i n s t a t e s
in the total exchange energy
K [ ρ ] = 1 2 σ ρ σ U X σ d r   ,
was given by Becke via employing the so called semiempirical (SE) modified gradient-corrected functional [77]:
K S E = K 0 β σ ρ σ 4 / 3 x σ 2 ( r ) 1 + γ x σ 2 ( r ) d r , K 0 = d r k 0 [ ρ ( r ) ] ,     x σ ( r ) = | ρ σ ( r ) | ρ σ 4 / 3 ( r )
to the working single-parameter dependent one [84]:
K B 88 = K 0 β σ ρ σ 4 / 3 ( r ) x σ 2 ( r ) 1 + 6 β x σ ( r ) sinh 1 x σ ( r ) d r
where the value β = 00042[a.u.] was found as the best fit among the noble gases (He to Rn atoms) exchange energies; the constant aσ is related to the ionization potential of the system.
Still, having different exchange approximation energetic functionals as possible worth explaining from where such ambiguity eventually comes. To clarify this, it helps in rewriting the starting exchange energy (94) under the formally exact form [85]:
K [ ρ ] = σ ρ σ ( r ) k [ ρ σ ( r ) ] g [ x σ ( r ) ] d r
Where the typical components are identified as:
k [ ρ ] = A X ρ 1 / 3 , A X = 3 2 ( 3 4 π ) 1 / 3 ,
while the gradient containing correction g(x) is to be determined.
Firstly, one can notice that a sufficiency condition for the two exchange integrals (111) and (114) to be equal is that their integrands, or the exchange potentials, to be equal; this provides the leading gradient correction:
g 0 ( x ) = 1 2 U X ( r ( x ) ) k [ ρ ( r ( x ) ) ]
with r(x) following from x(r) by (not unique) inversion.
Unfortunately, the above “integrity” condition for exchange integrals to be equal is not also necessary, since any additional gradient correction
g ( x ) = g 0 ( x ) + Δ g ( x )
fulfills the same constraint if it is chosen so that
ρ 4 / 3 ( r ) Δ g ( x ( r ) ) d r = 0
or, with the general form:
Δ g ( x ) = f ( x ) ρ 4 / 3 ( r ) f ( x ( r ) ) d r ρ 4 / 3 ( r ) d r
being f(x) an arbitrary function.
Nonetheless, if, for instance, the function f(x) is specialized so that
f ( x ) = g 0 ( x )
the gradient correcting function (117) becomes:
g ( x ) = 1 2 A X ρ ( r ) U X ( r ) d r ρ 4 / 3 ( r ) d r α X
recovering the Slater’s famous Xα method for exchange energy evaluations [86, 87]:
K [ ρ ] = α X A x ρ 4 / 3 ( r ) d r .
Nevertheless, the different values of the multiplication factor αX in (122) can explain the various forms of exchange energy coefficients and forms above. Moreover, following this conceptual line the above Becke’88 functional (113) can be further rearranged in a so called Xα-Becke88 form [88]:
K X B 88 = α X B σ ρ σ 4 / 3 ( r ) [ 2 1 / 3 + x σ 2 ( r ) 1 + 6 β X B x σ ( r ) sinh 1 x σ ( r ) ] d r
where the parameters αXB and βXB are to be determined, as usually, throughout atomic fitting; it may lead with a new workable valuable density functional in exchange family.

4.3. Density functionals of correlation energy

The first and immediate definition of energy correlation may be given by the difference between the exact and Hartree-Fock (HF) total energy of a poly-electronic system [89]:
E c [ ρ ] = E [ ρ ] E H F [ ρ ]
Instead, in density functional theory the correlation energy can be seen as the gain of the kinetic and electron repulsion energy between the full interacting (λ = 1) and non-interacting (λ = 0) states of the electronic systems [90]:
E c λ [ ρ ] = ψ λ | ( T ^ + λ V ^ e e ) | ψ λ ψ λ = 0 | ( T ^ + λ V ^ e e ) | ψ λ = 0 .
In this context, taking the variation of the correlation energy (125) respecting the coupling parameter λ [91, 92],
λ E c λ [ ρ ] λ = E c λ [ ρ ] + ρ ( r ) r δ E c λ [ ρ ] δ ρ ( r ) d r ,
by employing it through the functional differentiation with respecting the electronic density,
λ V c λ [ ρ ] λ V c λ [ ρ ] = r V c λ + ρ ( r 1 ) r 1 1 δ 2 E c λ [ ρ ] δ ρ ( r ) δ ρ ( r 1 ) d r 1 ,
one obtains the equation to be solved for correlation potential V c λ = δ E c λ [ ρ ] / δ ρ; then the correlation energy is yielded by back integration:
E c λ [ ρ ] = V c λ ( r , [ ρ ] ) ρ ( r ) d r
from where the full correlation energy is reached out by finally setting λ = 1.
When restricting to atomic systems, i.e. assuming spherical symmetry, and neglecting the last term of the correlation potential equation above, believed to be small [90], the equation to be solved simply becomes:
λ V c λ [ ρ ] λ V c λ [ ρ ] = r V c λ
that can really be solved out with the solution:
V c λ = A p λ p + 1 r p
with the integration constants Ap and p.
However, since the equation (129) is a homogeneous differential one, the linear combination of solutions gives a solution as well. This way, the general form of correlation potential looks like:
V c λ = p A p λ p + 1 r p .
This procedure can be then iterated by taking further derivative of (127) with respect to the density, solving the obtained equation until the second order correction over above first order solution (131),
V c λ = p 1 A p 1 λ p 1 + 1 r p 1 + p 2 A p 2 λ 2 p 2 + 1 r p 2 r p 2 ρ .
By mathematical induction, when going to higher orders the K-truncated solution is iteratively founded as:
V c λ = p k = 1 K A p k λ p k + 1 r p r p ρ k 1
producing the λ-related correlation functional:
E c λ [ ρ ] = p k = 1 K 1 k A p k λ p k + 1 r p ρ k
and the associate full correlation energy functional (λ=1) expression:
E c [ ρ ] = p k = 1 K 1 k A p k r p ρ k .
As an observation, the correlation energy (135) supports also the immediate not spherically (molecular) generalization [90]:
E c [ ρ ] = l m n k = 1 K 1 k A l m n k x 1 x m x n ρ k .
Nevertheless, for atomic systems, the simplest specialization of the relation (135) involves the simplest density moments 〈ρ〉 = N and 〈rρ〉 that gives:
E c [ ρ ] = A c 0 N + A c 1 r ρ .
Unfortunately, universal atomic values for the correlation constants Ac0 and Ac1 in (137) are not possible; they have to be related with the atomic number Z that on its turn can be seen as functional of density as well. Therefore, with the settings
A c 0 = C c 0 ln Z , A c l = C c 1 Z
the fitting of (137) with the HF related correlation energy (124) reveals the atomic-working correlation energy with the form [90]:
E c = 0.16569 N ln Z + 0.0000401 Z r ρ
The last formula is circumvented to the high-density total correlation density approaches rooting at their turn on the Thomas-Fermi atomic theory. Very interesting, the relation (139) may be seen as an atomic reflection of the (solid state) high density regime ( rs <1) given by Perdew et al. [75, 80, 9396]:
E c P Z [ ρ ] = d r ρ ( r ) ( 0.048 0.0116 r s + 0.0311 ln r s + 0.0020 r s ln r s )
in terms of the dimensionless ratio
r s = r 0 a 0
between the Wigner-Seitz radius r0 = (3/4πρ)1/3 , see relation (52), and the first Bohr radius a0 = ħ2 / me2.
Instead, within the low density regime ( rs ≥ 1) the first approximation for correlation energy goes back to the Wigner jellium model of electronic fluid in solids thus providing the LDA form [97, 98]:
E c W L D A [ ρ ] = ε c [ ρ ( r ) ] ρ ( r ) d r ,
where
ε c [ ρ ( r ) ] = 0.44 7.8 + r s
is the correlation energy per particle of the homogeneous electron gas with density ρ, see relations (55) and (56).
However, extended parameterization of the local correlation energy may be unfolded since considering the fit with an LSDA (ρ and ρ) analytical expression by Vosko, Wilk and Nusair (VWN) [99],
E c V W N [ ρ , ρ ] = ε c [ ρ ( r ) , ρ ( r ) ] ρ ( r ) d r
while further density functional gradient corrected Perdew (GCP) expansion will look like:
E c G C P [ ρ , ρ ] = d r ε c [ ρ ( r ) , ρ ( r ) ] ρ ( r ) + d r B [ ρ ( r ) , ρ ( r ) ] | ρ ( r ) | 2 +
where the Perdew recommendation for the gradient integrant has the form [100]:
B c P [ ρ ( r ) , ρ ( r ) ] = B c [ ρ ] exp ( b [ ρ ] f | ρ | ρ 7 / 6 ) d ( m )
with
B c [ ρ ] = ρ 4 / 3 C [ ρ ]
being the electron gas expression for the coefficient of the gradient expansion. The normalization in (146) is to the spin degeneracy:
d ( m ) = 2 1 / 3 [ ( 1 + m 2 ) 5 / 3 + ( 1 m 2 ) 5 / 3 ] 1 / 2 , m = ρ ρ ρ ,   ρ = ρ + ρ ,
while the exponent containing functional
b [ ρ ] = ( 9 π ) 1 / 6 C [ ρ ] c [ ρ ]
is written as the ratio of the asymptotic long-range density behavior to the current one, and is controlled by the cut-off f exponential parameter taking various values depending of the fitting procedures it subscribes (0.17 for closed shells atoms and 0.11 for Ne particular system [101, 102]).
More specifically, we list bellow some nonlocal correlation density functionals in the low density (gradient corrections over LDA) regime:
  • the Rasolt and Geldar paramagnetic case (ρ = ρ = ρ/2) is covered by correlation energy [101, 103]:
    E c R G [ ρ ] = c 1 + c 2 + c 3 r s + c 4 r s 2 1 + c 5 r s + c 6 r s 2 + c 7 r s 3
with c1=1.667·10−3, c1=2.568·10−3, c3=2.3266·10−2, c4=7.389·10−6, c5=8.723, c6=0.472, c7=7.389·10−2 (in atomic units).
  • The gradient corrected correlation functional reads as [102, 104]:
    E c G C = d r ε c [ ρ , ρ ] ρ ( r ) + d r B c p [ ρ , ρ ] c [ ρ ] = 2 π / 4 ( 6 π 2 ) 4 / 3 , f = 0.17 | ρ ( r ) | 2 + + 9 π 4 ( 6 π 2 ) 4 / 3 ( 0.17 ) 2 d r ( | ρ | 2 ρ 4 / 3 + | ρ | 2 ρ 4 / 3 )
  • The Lee, Yang, and Parr (LYP) functional within Colle-Salvetti approximation unfolds like [105]:
    E c L Y P = a c b c d r γ ( r ) ξ ( r ) ( Σ σ ρ σ ( r ) Σ i | φ i σ ( r ) | 2 1 4 Σ σ ρ σ ( r ) Δ ρ σ ( r ) 1 4 | ρ ( r ) | 2 + 1 4 ρ ( r ) Δ ρ ( r ) ) a c d r γ ( r ) η ( r ) ρ ( r )
where
γ ( r ) = 4 ρ ( r ) ρ ( r ) ρ ( r ) 2 , η ( r ) = 1 + d c ρ ( r ) 1 / 3 ,   ξ ( r ) = ρ ( r ) 5 / 3 η ( r ) exp [ c c ρ ( r ) 1 / 3 ]
and the constants: ac=0.04918, bc=0.132, cc=0.2533, dc=0.349.
  • The open-shell (OS) case provides the functional [98]:
    E c O S = d r a s ρ ( r ) + b s | ρ ( r ) | ρ ( r ) 1 / 3 c s + d s ( | ρ | ρ 4 / 3 + | ρ | ρ 4 / 3 ) + r s 1 ζ 2
with the spin-dependency regulated by the factor ζ = (ρ − ρ)/(ρ + ρ) approaching zero for closed-shell case, while the specific coefficients are determined through a scaled-minimization procedure yielding the values: as=−0.74860, bs=−0.06001, cs=3.60073, ds=0.900000.
  • Finally, Perdew and Zunger (PZ) recommend the working functional [106]:
    E c P Z 0 [ ρ ] = d r ρ ( r ) α p 1 + β 1 p r s + β 2 p r s
with the numerical values for the fitting parameters founded as: αp=−0.1423, β1p=1.0529, β2p=0.3334.

4.4. Density functionals of exchange-correlation energy

Another approach in questing exchange and correlation density functionals consists in finding them both at once in what was defined as exchange-correlation density functional (29). In this regard, following the Lee and Parr approach [107], the simplest starting point is to rewrite the inter-electronic interaction potential
V e e ρ 2 ( r 1 , r 2 ) r 12 d r 1 d r 2
and the classical (Coulombic) repulsion
J = 1 2 ρ ( r 1 ) ρ ( r 2 ) r 12 d r 1 d r 2 ,
appeared in the formal exchange energy (VeeJ) in (21), by performing the previously introduced coordinate transformation (79), followed by integration of the averaged pair and coupled densities (denoted with over-bars) over the angular components of s:
V e e = 4 π d r s d s ρ ¯ 2 ( r , s ) ,
J = 2 π d r s d s ρ ( r + s / 2 ) ρ ( r s / 2 ) ¯ .
Now, the second order density matrix in (158) can be expressed as
ρ ¯ 2 ( r , s ) = 1 2 ρ ( r + s / 2 ) ρ ( r s / 2 ) ¯ [ 1 + F 1 ( r , s ) ]
with the help of the introduced function F1 (r, s) carrying the form
F 1 ( r , s ) = exp [ α ( r ) s ] 1 + α ( r ) { 1 + [ α ( r ) s ] 2 F 2 ( r , s ) }
so that the cusp condition for ρ̄2 (r, s)
ln ρ ¯ 2 ( r , s ) s | s = 0 = 1
to be satisfied for a well behaved function of a Taylor series expansion type
F 2 ( r , s ) = k = 0 a k ( r ) [ α ( r ) s ] k
when α(r) stands for a suitable function of r as well, see bellow.
On the other side, the average ρ ( r + s / 2 ) ρ ( r s / 2 ) ¯ in (159) and (160) supports a Taylor expansion [108]:
ρ ( r + s / 2 ) ρ ( r s / 2 ) ¯ = ρ 2 ( r ) [ 1 2 τ w ( r ) 3 ρ ( r ) s 2 + ]
with
τ w ( r ) = 1 8 | ρ ( r ) | 2 ρ ( r ) 1 8 2 ρ ( r )
being the Parr modified kinetic energy of Weizsäcker type [9].
Inserting relations (158)-(165) in (VeeJ) difference it is eventually converted from the “genuine” exchange meaning into practical exchange-correlation energy characterized by the density functional form:
E x c = 2 π d r s d s   ρ ( r + s / 2 ) ρ ( r s / 2 ) ¯ F 1 ( r , s ) = 2 π d r ρ 2 ( r ) 1 + α ( r ) d s s exp [ α ( r ) s ] { 1 2 τ w ( r ) 3 ρ ( r ) s 2 + ... } { 1 + [ α ( r ) s ] 2 Σ k = 0 a k ( r ) s ] k }
Making use of the two possible multiplication of the series in (166), i.e. either by retaining the α(r) containing function only or by including also the density gradient terms in the first curled brackets, thus retaining also the term containing τw (r) function, the so called I-xc or II-xc type functionals are respectively obtained.
Now, laying aside other variants and choosing the simple (however meaningfully) density dependency
α ( r ) = k ρ 1 / 3 ( r ) ,    k = constant
the provided exchange-correlation functionals are generally shaped as [107]:
E s c I = 1 k 2 d r ρ 4 / 3 ( r ) A s c ( r ) 1 + k ρ 1 / 3 ( r ) , E x c I I = 1 k 2 d r ρ 4 / 3 ( r ) 1 + k ρ 1 / 3 ( r ) [ B x c ( r ) + 2 3 τ w ( r ) ρ 5 / 3 ( r ) C x c ( r ) ] .
These functionals are formally exact for any κ albeit the resumed functions Axc(r), Bxc(r), and Cxc(r) are determined for each particular specialization.
Going now to the specific models, let’s explore the type I of exchange-correlation functionals (168). Firstly, they can further undergo simplification since the reasonable (atomic) assumption according which
k ρ 1 / 3 ( r ) 1 , r .
Within this frame the best provided model is of Xα-Padé approximation type, containing N-dependency [107]:
E x c I ( X α ) = a 0 X α 1 + a 1 X α / N 1 + a 2 X α / N ρ 4 / 3 ( r ) d r
with a0=0.7475, a1=17.1903, and a2=14.1936 (atomic units).
When the condition (169) for κ is abolished the Wigner-like model results, again, having the best approximant exchange-correlation model as the Padé form [107]:
E x c I ( W i g ) = a 0 W i g 1 + a 1 W i g / N 1 + a 2 W i g / N ρ 4 / 3 ( r ) 1 + k I ( W i g ) ρ 1 / 3 ( r ) d r
with a0Wig=0.76799, a1Wig=17.5943, a2Wig=14.8893, and κI (Wig)=4.115·10−3 (atomic units).
Turning to the II-type of exchange-correlation functionals, the small density condition (169) delivers the gradient corrected Xα model, taking its best fitting form as the N-dependent Padé approximant [107]:
E x c I I ( X α ) = b 0 X α 1 + b 1 X α / N 1 + b 2 X α / N ρ 4 / 3 ( r ) d r c 0 X α ρ 1 / 3 ( r ) τ w ( r ) d r
with b0=0.7615, b1=1.6034, b2=2.1437, and c2=6.151·10−2 (atomic units), while when laying outside the (169) condition the gradient corrected Wigner-like best model is proved to be without involving the N-dependency [107]:
E s c I I ( W i g ) = b 0 W i g ρ 4 / 3 ( r ) 1 + k I I ( W i g ) ρ 1 / 3 ( r ) d r c 0 W i g ρ 1 / 3 ( r ) τ w ( r ) 1 + k I I ( W i g ) ρ 1 / 3 ( r ) d r
with b0Wig=0.80569, c0Wig=3.0124·10−3, and κII(Wig)=4.0743·10−3 (atomic units).
Still, a Padé approximant for the gradient-corrected Wigner-type exchange-correlation functional exists and it was firstly formulated by Rasolt and Geldar [103] with the working form [109]:
E x c R G = E x c L D A ( o r X α ) + B x c R G [ ρ ( r ) ] | ρ ( r ) | 2 ρ 1 / 3 ( r ) d r
with BxcRG given with the Padé form:
B x c R G [ ρ ( r ) ] = 1 × 10 3 c 1 R G 1 + c 2 R G r s + c 3 R G r s 2 1 + c 4 R G r s + c 5 R G r s 2 + c 6 R G r s 3
having the fitted coefficients c1RG=2.568, c2RG=9.0599, c3RG=2.877·10−3, c4RG=8.723, c5RG=0.472, and c3RG=7.389·10−2 (atomic units). Some studies also consider the nonlocal correction in (174) premultiplied by the 10/7 factor which was found as appropriate procedure for atomic systems.
Finally, worth noting the Tozer and Handy general form for exchange-correlation functionals viewed as a sum of products of powers of density and gradients [110]:
E x c T H = F x c ( ρ , ρ , ζ ) d r
with
F x c = a b c d ω a b c d R a S b X c Y d = a b c d ω a b c d f a b c d ( r )
where
R a = ρ a + ρ a , S b = m 2 b see equation (148) for  m d e f i n i t i o n , X c = ζ c + ζ c 2 ρ 4 c / 3 , Y d = ( ζ 2 + ζ 2 - 2 ζ ρ 8 / 3 ) d
and
ζ = | ρ | ,   ζ = | ρ | ,   ζ = ρ ρ ,   ρ = ρ + ρ .
The coefficients ωabcd of (177) are determined through minimization procedure involving the associated exchange-correlation potentials V x c ( ) a b c d ( r ) = δ f a b c d ( r ) / δ ρ ( ) ( r ) in above (176) functional. The results would depend upon the training set of atoms and molecules but presents the advantage of incorporating the potential information in a non-vanishing asymptotical manner, with a semi-empirical value. Moreover, its exact asymptotic exchange-correlation potential equals chemical hardness [9, 12, 14, 20] for open-shell being less than that for closed shell systems, thus having the merit of including chemical hardness as an intrinsic aspect of energetic approach, a somewhat absent aspect from conventional functionals so far.
However, since electronegativity and chemical hardness closely relate with chemical bonding, their relation with the total energy and component functionals is in next at both conceptual and applied levels explored.

5. Testing (χ, η) Quadratic Dependency Among Several Energetic Density Functionals

5.1. Proof of the E=E(χ,η) quadratic dependency

Employing the Kohn-Sham equation (33) for the (reactivity) equalized chemical potential (minus electronegativity) eigenvalues, among atomic or molecular one-electronic orbitals
μ i = μ j = = μ ,
in quantum mechanically sense, we immediately get:
μ = φ ( r ) [ 1 2 2 + V e f f ] φ ( r ) d r φ ( r ) φ ( r ) d r = φ ( r ) H ^ φ ( r ) d r φ ( r ) φ ( r ) d r = E N
For proper characterization of the chemical systems the quantum (Ehrenfest) version of the fundamental Newtonian law [111], linking the observed force with the minus of the potential gradient F = −∇V, is here written for chemical reactions modeled throughout the charge transfer along the reaction path as:
F μ = d d N ( μ ) .
Now, combining the last two equations projected on the reaction path we successively obtain:
F μ = d d N ( E N ) = E N 2 1 N ( d E d N ) V .
At this point, taking for the derivative in (183) the finite correspondence in what regarding the chemical potential formal (absolute) definition (181),
( d E d N ) V = μ = E N ,
the chemical potential energy equation is unfolded as:
E μ = μ N + F μ N 2 .
The remaining issue is to clarify the chemical potential related force meaning in above equation. In this respect, by considering the electronegativity-chemical potential relationship (19), the associate electronegativity energy equation becomes:
E χ = χ N + F χ N 2
while the involved electronegativity related force is seen as the reactive force of the chemical potential (182),
F χ = F μ
recovering the chemical hardness index:
F χ = d d N ( χ ) = η .
With these considerations, the chemical reactivity energy equation (185) is finally displayed as
E = χ N + η N 2
equally supporting the differential (statistical, thermodynamical) variant:
d E = χ d N + η ( d N ) 2
based on the above finite-differential (184) equivalence between the total (E, N) or exchanged (dE, dN) energy and charge, respectively.
What there was actually proofed is that the quadratic chemical reactivity equations for total energy in both finite and differential fashions may be derived employing the Kohn-Sham (as Schrödinger reminiscence) equation for chemical potential eigenvalue combined with the chemical quantum version of the Ehrenfest theorem involving the force concept and its active-reactive peculiar property for the chemical potential and electronegativity, respectively. No particular assumptions were considered being all above arguments only on first principles grounded. Thus, the exposed demonstration is of general value cutting much discussion in the last decades on the viability of the second order truncation in the total energy expansion in terms of chemical reactivity indices, viz. electronegativity and chemical hardness concepts.
However, a computational test for this behavior is in next addressed.

5.2. Atomic and molecular analysis of the energetic quadratic bilinear (χ, η) dependency

The general bivariate equation linking the energy (various) functionals with electronegativity and chemical hardness, either for atomic and molecular systems, works out with form:
E = a + b χ + c η
with coefficients a, b, and c being determined throughout consecrated statistical analysis methods [112]. The result of such correlation will lead with two kind of information:
  • the degree of correlation itself between the employed energy functional and the couple of electronegativity-chemical hardness structural indices; this is measured by the standard correlation factor [113]:
    r = 1 Σ i ( y i I N P U T y i F I T ) 2 Σ i ( y i I N P U T y i I N P U T ¯ ) 2
  • the degree of parabolic dependency by checking whether the chemical hardness coefficient (c) is the square of the electronegativity coefficient (b) thus giving the opportunity of introducing the so-called sigma-pi reactivity index
    σ π = s i g n ( b ) c b 2 parabolic    E = E ( χ , η ) 1
At this point worth nothing that as σπ → −1 as better the energy fulfils the correlation shape with the established parabolic equation (189). The b and c signs open further discussion among the allowed combinations for the correlation equation (191). Firstly, let’s observe that the signs of b and c give information about the signs of electronegativity and chemical hardness, respectively. Then, while considering the finite-difference approximations to the electronegativity and chemical hardness definitions, in a Koopmans theorem environment [114], they take the working forms [9, 14, 115, 116]:
χ ( E N N ) V ( r ) E N + 1 E N 1 2 = ( E N 1 E N ) + ( E N E N + 1 ) 2 = I P + E A 2 H O M O + L U M O 2
η = ( 2 E N N 2 ) V ( r ) E N + 1 2 E N + E N 1 = I P E A L U M O H O M O
in terms of lowest unoccupied molecular orbital (LUMO) and highest occupied molecular orbital (HOMO) energies, although similar relations (in terms of ionization potential IP and electron affinity EA) hold for atomic systems as well. Nevertheless, the expressions (194) and (195) help in deciding that only the (−, +) and (+, +) combinations for coefficients (b, c) are allowed, based on the fact that only maximum hardness values (i.e. positive values of c coefficient) are in accordance with the reactivity criteria of maximum hardness principle [117119]; in other words, the correlation combination carrying negative values for chemical hardness coefficient in (191) are less physically probably since that would imply that EA>0 and LUMO<0, for atoms and molecules, respectively.
Finally, the sigma-pi index (193) can be used in defining another reactivity index, namely the efficiency index:
Π = | σ π | p a r a b o l i c E = E ( χ , η ) 1
measuring the power of electronic exchange characterizing the chemical reactivity which is maximized by the exact parabolic dependence of the energetic functional respecting the electronegativity-chemical hardness jointly influence.
With these, the Tables IIII and IVV display the tested energetic functionals against various electronegativity and chemical hardness scales for selected atoms and molecules, while the Tables VI and VII show their quantitative (χ & η) structural– (energy functionals) property relationships (QSPR), respectively.
At a glance, Tables VI and VII reveals that the parabolic energetic dependency (189) on the electronegativity and chemical hardness indices is not exactly recovered since all sign combinations in the coefficients c & b as well as poor correlation itself with the (χ, η) indices is revealed, i.e. widely deviating from the optimum (−,+) and maximum correlation coefficient of (192).
However, for atomic analysis, Table VI gives us that the finite-difference scale of electronegativity and chemical hardness provides invariable right (−, +) sign combination in coefficients (b, c), excepting kinetic energy functionals cases, however revealing both scarce reactive efficiency (196) and bilinear correlation (192) in almost all treated cases. Notably, the higher reactivity efficiency on studied atomic system was furnished by the correlation energy merely than the total one, yet with a low correlation factor.
The situation is somehow changed for molecular systems, see Table VII, where, albeit not better records are noted in the sense of the maximum reactive efficiency (196), again excepting the correlation functional case rather than the total energy, a higher correlation factor in semiclassical treatment of the electronegativity and chemical hardness scales is obtained; such behavior was recently validated for hard and soft acids and bases reactions as well [24]. At the same time, the kinetic energy functionals do not fit at all with standard (−,+) correlation scheme for signs of (c, b) coefficients in (191).
Finally, in Figure 5 the atomic and molecular parabolic reactivity efficiencies (196) of various energetic functionals are plotted as charts with the widths emphasizing the bilinear electronegativity & chemical hardness correlation values of (192) factor as abstracted from Tables VI and VII. The comparative behavior reveals many interesting features to be taking into account for further conceptual and computational studies on quantum chemistry:
  • the correlation energy appears to provide acceptable parabolic shapes in both atomic and molecular cases, with better bilinear regression for molecular analysis, while strongly depending on the electronegativity and chemical hardness atomic models and scales;
  • the kinetic energy, while displaying poor parabolic shape at atomic level behaves with negative chemical hardness in molecular systems, probably due the positive contribution in bonding that compete with stabilization (localization) of the electrons within internuclear basin;
  • exchange and exchange-correlation functionals reveal similar reactive (parabolic) efficiency as well as close bivariate regression correlation factors for both atomic and molecular cases, leaving with the impression that the exchange contribution is dominant in exchange-correlation functionals since cancelling somehow the behavior of the correlation part of the functional.
  • overall, the total energy, although with correlation factors in the range of its components’ regressions does not fit with parabolic reactive theoretical prescription (189), at least for present employed set of atoms and molecules.
However, all revealed parabolic features have appeared on the statistical correlation basis and not upon an individual atomic or molecular analysis. At this point we can refer to studies [126] that clearly illustrate the individual parabolic behavior of the total energy against the total number of electrons. On the other side, the systematic focus on atomic and molecular large class of systems should provide more or less the same phenomenological results.
The present dichotomy, while confronting the conceptual and computational behavior of the energetic functionals respecting the parabolic bilinear expansion in terms of electronegativity and chemical hardness, gives the indication that the derivation of the energetic functionals (and of the total energy in special) starting from electronegativity and chemical hardness functionals [14, 123], may eventually conclude with better energy correlation both as parabolic shape and statistical regression with fundamental (χ, η) couple of reactivity indices. Otherwise, we will be forced to admit that the density functionals do not properly model the reactivity through a parabolic fit in total or exchanged electrons of the involved systems at least in a statistical sense.

6. Conclusions

Why another quantum theory of the matter? Why the Density Functional Theory? Its pioneers and promoters, Walter Kohn, Axel Becke and Robert G. Parr – just to name a few, affirm that modeling electronic density features in 3D real space produces much benefits in both conceptually understanding of bonding as well in electronic localization with a bonding impact for the concerned N-many-electronic systems [127129].
The present paper aimed to unitarily presents the quantum fundaments of DFT as well to prospect directions in electronic bonding and reactivity characterization by using the density functionals and localization functions. Such theoretical tools allow the rationalization of the electronic structures while offering an analytical understanding and prediction of the experimental data along of the modeled reality.
In DFT, the density and their combinations in the density functionals of the total energy plays a primordial role. It fulfills the N-contingency, assures the total energy minimization, influences the different levels of approximation, i.e. local density or gradient density frameworks, controls the bonding through electronic localization functions, and decides upon reactivity through the electronic exchange relating the electronegativity and chemical hardness indices. Studying the electronic density properties of electronic systems guaranties the universal treatment of whatever systems, no matter how rich in electrons is, from atoms, to molecules, to solids and aggregates.
Going to link the inner structural information with the manifested chemical reactivity the various energetic density functionals of atoms and molecules were presented, emphasizing on different levels of approximation and of quantum containing information. Nevertheless, their values were assumed merely as properties of particular systems when combined with associate reactive peculiarities quantified by electronegativity and chemical hardness. A theoretical link, based on first quantum mechanical principles, was also provided at a level of chemical potential-chemical force couple for the electronegativity and chemical hardness actions, respectively.
Actually, the parabolic shape of the energetic dependence on electronegativity and chemical hardness ( E = a + bχ + cη) is theoretically proofed with the help of potential-force physical picture of interaction and dynamics, while the computational cases were treated throughout the introduced sigma-pi index σπ = sign(b)c / b2 or parabolic reactive efficiency ∑Π = |σπ|.
The numerical results of such correlation on selected, however limited, set of atomic and molecular systems, were intrigued since there was found out that only the correlation functionals provide appropriate parabolic and statistical behavior among considered atoms and molecules, while all other total energy components and the total energy itself posse both surprisingly low parabolic and correlation factors with the electronegativity and chemical hardness structural parameters. Nevertheless, at atomic level the finite-difference approximation for electronegativity and chemical hardness scales was found as most suitable in providing a correct (−, +) sign combination in parabolic (b, c) coefficients, whereas molecular systems comes out within such agreement when semiclassical scale was adopted. However, the conceptual-computational parabolic reactive dichotomy of energy functionals against the reactivity indices may be eventually in future avoided by generating energy functionals directly from electronegativity and chemical hardness functionals [130]. Such paradox is to be abolished when comprehensive quantum chemical knowledge of structure and reactivity will be unfolded - step by step, from close to closer.

Acknowledgements

The kind support of the Molecular Diversity Preservation International Foundation – MDPI for publishing this work and for recognizing the seminal role of “Chemical Bond and Bonding” theoretical and computational forum in actual Chemistry, as well the involvement of the Romanian National Council of Scientific Research in Universities – CNCSIS in promoting excellence exploratory research programs are hearty saluted.

References

  1. Kohn’s Nobel Prize: http://nobelprize.org/nobel_prizes/chemistry/laureates/1998/index.html.
  2. Hohenberg, P; Kohn, W. Inhomogeneous Electronic Gas. Phys. Rev. 1964, 136, 864–871. [Google Scholar]
  3. Kohn, W; Sham, LJ. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, 1133–1138. [Google Scholar]
  4. Kohn, W; Becke, AD; Parr, RG. Density Functional Theory of Electronic Structure. J. Phys. Chem. 1996, 100, 12974–12980. [Google Scholar]
  5. Keller, J. On the Formulation of the Hohenberg-Kohn-Sham Theory. Int. J. Quantum Chem. 1986, 20, 767–768. [Google Scholar]
  6. Ernzerhof, M. Density-Functional Theory as an Example for the Construction of Stationarity Principles. Phys. Rev. A 1994, 49, 76–79. [Google Scholar]
  7. DePristo, AE. Hohenberg-Kohn Density-Functional-Theory as an Implicit Poisson equation for Density Changes from Summed Fragment Densities. Phys. Rev. A 1996, 54, 3863–3869. [Google Scholar]
  8. Bamzai, AS; Deb, BM. The Role of Single-Particle Density in Chemistry. Rev. Mod. Phys. 1981, 53, 95–126. [Google Scholar]
  9. Parr, RG; Yang, W. Density Functional Theory of Atoms and Molecules; Oxford University Press: New York, 1989. [Google Scholar]
  10. Dreizler, RM; Gross, EKU. Density Functional Theory; Springer Verlag: Heidelberg, 1990. [Google Scholar]
  11. March, NH. Electron Density Theory of Many-Electron Systems; Academic Press: New York, 1991. [Google Scholar]
  12. Putz, MV. Contributions within Density Functional Theory with Applications to Chemical Reactivity Theory and Electronegativity; Disertation. Com; Parkland: Florida, 2003; http://www.dissertation.com/book.php?method=ISBN&book=1581121849.
  13. Putz, MV; Chiriac, A; Mracec, M. Foundations for a Theory of the Chemical Field. II. The Chemical Action. Rev. Roum. Chimie 2001, 46, 1175–1181. [Google Scholar]
  14. Putz, MV. Unifying Absolute and Chemical Electronegativity and Hardness Density Functional Formulations through the Chemical Action Concept. In “Progress in Quantum Chemistry Research”; Hoffman, Erik O, Ed.; Nova Science Publishers Inc.: New York, 2007; pp. 59–121. [Google Scholar]
  15. Chen, J; Stott, MJ. V-Representability for Systems of a Few Fermions. Phys. Rev. A 1991, 44, 2809–2814. [Google Scholar]
  16. Chen, J; Stott, MJ. V-Representability for Systems with Low Degeneracy. Phys. Rev. A 1991, 44, 2816–2821. [Google Scholar]
  17. Kryachko, ES; Ludena, EV. Formulation of N-and V-Representable Density-Functional Theory. I. Ground States. Phys. Rev. A 1991, 43, 2179–2192. [Google Scholar]
  18. Kryachko, ES; Ludena, EV. Formulation of N-and V-Representable Density-Functional Theory. II. Spin-Dependent Systems. Phys. Rev. A 1991, 43, 2194–2198. [Google Scholar]
  19. Levy, M; Perdew, J. The Constrained Search Formulation of Density Functional Theory. In ”Density Functional Methods in Physics”; Dreizler, RM, da Providencia, J, Eds.; NATO ASI Series B: Physics; Volume 123, Plenum Press: New York, 1985; pp. 11–31. [Google Scholar]
  20. Putz, MV. Can Quantum-Mechanical Description of Chemical Bond Be Considered Complete? In ”Quantum Chemistry Research Trends”; Mikas, P. Kaisas, Ed.; Nova Science Publishers Inc.: New York, 2007; Expert Commentary. [Google Scholar]
  21. Janak, JF. Proof that ∂E / ∂ni = εi in Density Functional Theory. Phys. Rev. B 1978, 18, 7165–7168. [Google Scholar]
  22. Moscardo, F; San-Fabian, E. Density-Functional Formalism and the Two-Body Problem. Phys. Rev. A 1991, 44, 1549–1553. [Google Scholar]
  23. Neal, HL. Density Functional Theory of One-Dimension Two-Particle Systems. Am. J. Phys. 1998, 66, 512–516. [Google Scholar]
  24. Putz, MV; Chiriac, A. Quantum Perspectives on the Nature of the Chemical Bond. In “Advances in Quantum Chemical Bonding Structures”; Putz, MV, Ed.; Transworld Research Network: Kerala, India, 2008; pp. 1–43. [Google Scholar]
  25. Flores, JA; Keller, J. Differential Equations for the Square Root of the Electronic Density in Symmetry-Constrained Density-Functional Theory. Phys. Rev. A 1992, 45, 6259–6262. [Google Scholar]
  26. Nagy, A. Kohn-Sham Equations for Multiplets. Phys. Rev. A 1998, 57, 1672–1677. [Google Scholar]
  27. Garcia-Gonzales, P; Alvarellos, JE; Chacon, E. Kinetic-Energy Density Functional: Atoms and Shell structure. Phys. Rev. A 1996, 54, 1897–1905. [Google Scholar]
  28. Chan, GKL; Handy, NC. Kinetic-Energy Systems, Density Scaling, and Homogeneity Relations in Density Functional Theory. Phys. Rev. A 1999, 59, 2670–2679. [Google Scholar]
  29. Bartolotti, LJ; Acharya, PK. On the Functional Derivative of the Kinetic Energy Density Functional. J. Chem. Phys. 1982, 77, 4576–4585. [Google Scholar]
  30. Dawson, KA; March, NH. Slater Sum in One Dimension: Explicit Kinetic Energy Functional. Phy. Lett. 1984, 106A, 158–160. [Google Scholar]
  31. Baltin, R. The Three-Dimensional Kinetic Energy Density Functional Compatible with the Exact Differential Equation for its Associated Tensor. J. Chem. Phys. 1987, 86, 947–952. [Google Scholar]
  32. Becke, AD. Correlation Energy of an Inhomogeneous Electron Gas: A Coordinate-Space. J. Chem. Phys. 1988, 88, 1053–1062. [Google Scholar]
  33. Ayers, PW; Levy, M. Sum Rules for Exchange and Correlation Potentials. J. Chem. Phys. 2001, 115, 4438–4443. [Google Scholar](b)Koch, W; Holthausen, MC. A Chemist's Guide to Density Functional Theory; Wiley-VCH: Weinheim, 2000. [Google Scholar]
  34. Ponec, R; Cooper, DL. Anatomy of bond formation. Bond length dependence of the extent of electron sharing in chemical bonds from the analysis of domain-averaged Fermi holes. Faraday Discuss. 2007, 135, 31–42. [Google Scholar]Matito, E; Solà, M; Salvador, P; Duran, M. Electron sharing indexes at the correlated level. Application to aromaticity calculations. Faraday Discuss. 2007, 135, 325–345. [Google Scholar]March, NH; Howard, IA; Van Doren, VE. Recent Progress in Constructing Nonlocal Energy Density Functionals. Int. J. Quantum Chem. 2003, 92, 192–204. [Google Scholar]
  35. Zhao, Q; Morrison, RC; Parr, RG. From Electron Densities to Kohn-Sham Kinetic Energies, Orbital Energies, Exchange-Correlation Potentials, and Exchange-Correlation Energies. Phys. Rev. A 1994, 50, 2138–2142. [Google Scholar]
  36. Gritsenko, OV; Schipper, PRT; Baerends, EJ. Ensuring Proper Short-Range and Asymptotic Behavior of the Exchange-Correlation Kohn-Sham Potential by Modeling with a Statistical Average of Different Orbital Model Potential. Int. J. Quantum Chem. 2000, 76, 407–419. [Google Scholar]
  37. Zhao, Q; Parr, RG. Local Exchange-Correlation Functional: Numerical Test for Atoms and Ions. Phys. Rev. A 1992, 46, R5320–R5323. [Google Scholar]
  38. Lam, KC; Cruz, FG; Burke, K. Viral Exchange-Correlation Energy Density in Hooke’s Atom. Int. J. Quantum Chem. 1998, 69, 533–540. [Google Scholar]
  39. Gaspar, R; Nagy, A. Local-Density-Functional Approximation for Exchange-Correlation Potential. Application of the Self-Consistent and Statistical Exchange-Correlation Parameters to the Calculation of the Electron Binding Energies. Theor. Chim Acta 1987, 72, 393–401. [Google Scholar]
  40. Levy, M. Density-Functional Exchange Correlation Through Coordinate Scaling in Adiabatic Connection and Correlation Hole. Phys. Rev. A 1991, 43, 4637–4645. [Google Scholar]
  41. Guo, Y; Whitehead, MA. Generalized Local-Spin-Density-Functional Theory. Phys. Rev. A 1991, 43, 95–108. [Google Scholar]
  42. Dunlap, BI. Symmetry and Spin Density Functional Theory. Chem. Phys. 1988, 125, 89–97. [Google Scholar]
  43. Harrison, JG. Electron Affinities in the Self-Interactions-Corrected Local Spin Density Approximation. J. Chem. Phys. 1987, 86, 2849–2853. [Google Scholar]
  44. Becke, AD. Completely Numerical Calculations on Diatomic Molecules in the Local-Density Approximation. Phys. Rev. A 1986, 33, 2786–2788. [Google Scholar]
  45. Manoli, SD; Whitehead, MA. Generalized Exchange Local-Spin-Density-Functional Theory: One-Electron Energies and Eigenvalues. Collection Czechoslovak Chem. Commun. 1988, 53, 2279–2307. [Google Scholar]
  46. Filippeti, A. Electron Affinity in Density-Functional Theory in the Local Spin-Density Approximation. Phys. Rev. A 1998, 57, 914–919. [Google Scholar]
  47. Guo, Y; Whitehead, MA. Application of Generalized Exchange Local-Spin-Density-Functional Theory: Electronegativity, Hardness, Ionization Potential, and Electron Affinity. Phys. Rev. A 1989, 39, 2317–2323. [Google Scholar]
  48. Dunlap, BI; Andzelm, J. Second Derivatives of the Local-Density-Functional Total Energy when the Local Potential is Fitted. Phys. Rev. A 1992, 45, 81–86. [Google Scholar]
  49. Dunlap, BI; Andzelm, J; Mintmire, JW. Local-Density-Functional Total Energy Gradients in the Linear Combination of Gaussian-type Orbitals Method. Phys. Rev. A 1990, 42, 6354–6358. [Google Scholar]
  50. Liu, S. Local-Density Approximation, Hierarchy of Equations, Functional Expansion, and Adiabatic Connection in Current-Density-Functional Theory. Phys. Rev. A 1996, 54, 1328–1336. [Google Scholar]
  51. Romera, E; Dehesa, JS. Weizsäcker Energy of Many-Electron Systems. Phys. Rev. A 1994, 50, 256–266. [Google Scholar]
  52. Murphy, DR. Sixth-Order Term of the Gradient Expansion of the Kinetic-Energy Density Functional. Phys. Rev. A 1981, 24, 1682–1688. [Google Scholar]
  53. Dufek, P; Blaha, P; Sliwko, V; Schwarz, K. Generalized-Gradient-Approximation Description of Band Splittings in Transition-Metal Oxides and Fluorides. Phys. Rev. B 1994, 49, 10170–10175. [Google Scholar]
  54. Becke, AD; Edgecombe, KE. A Simple Measure of Electron Localization in Atomic and Molecular Systems. J. Chem. Phys. 1990, 92, 5397–5403. [Google Scholar]
  55. Silvi, B; Gatti, C. Direct Space Representation of the Metallic Bond. J. Phys. Chem. A 2000, 104, 947–953. [Google Scholar]
  56. Bader, RFW; Gillespie, RJ; MacDougall, PJ. A Physical Basis for the VSEPR Model of Molecular Geometry. J . Am. Chem. Soc. 1988, 110, 7329–7336. [Google Scholar]
  57. Silvi, B; Savin, A. Classification of the Chemical Bonds Based on Topological Analysis of Electron Localization Functions. Nature 1994, 371, 683–686. [Google Scholar]
  58. Cioslowski, J. Nonnuclear Attractors in the LI2 Molecule. J. Phys. Chem. 1990, 94, 5496–5498. [Google Scholar]
  59. Liu, S; Ayers, PW; Parr, RG. Alternative Definition of Exchange-Correlation Charge in Density Functional Theory. J. Chem. Phys. 1999, 111, 6197–6203. [Google Scholar]
  60. Putz, MV. Markovian Approach of the Electron Localization Functions. Int. J. Quantum Chem. 2005, 105, 1–11. [Google Scholar]
  61. Santos, JC; Tiznado, W; Contreras, R; Fuentealba, P. Sigma-pi Separation of the Electron Localization Function and Aromaticy. J. Chem. Phys. 2000, 120, 1670–1673. [Google Scholar]Scemama, A; Chaquin, P; Caffarel, M. Electron Pair Localization Function: A Practical Tool to Visualize Electron Localization in Molecules from Quantum Monte Carlo Data. J. Chem. Phys. 2000, 121, 1725–1735. [Google Scholar]Soncini, A; Lazzeretti, P. Nuclear Spin-Spin Coupling Density Functions and the Fermi Hole. J. Chem. Phys. 2003, 119, 1343–1349. [Google Scholar]Silvi, B. The Spin-Pair Compositions as Local Indicators of the Nature of the Bonding. J. Phys. Chem. 2003, 107, 3081–3085. [Google Scholar]
  62. Matito, E; Silvi, B; Duran, M; Solà, M. Electron Localization Function at the Correlated Level. J. Chem. Phys. 2006, 125, 024301. [Google Scholar]
  63. Kohout, M; Pernal, K; Wagner, FR; Grin, Y. Electron Localizability Indicator for Correlated Wavefunctions. I. Parallel-Spin Pairs. Theor. Chem. Acc. 2004, 112, 453–459. [Google Scholar]
  64. Jensen, F. On the Accuracy of Numerical Hartree–Fock Energies. Theor. Chem. Acc. 2005, 113, 187–190. [Google Scholar]
  65. Gombás, P; Szondy, T. Solutions of the Simplified Self-Consistent Field for All Atoms of the Periodic System of elements from Z=2 to Z=92; Adam Hilger Ltd.: London, 1970. [Google Scholar]
  66. Putz, MV; Ionaşcu, C; Hulesch, O. Global and Local Indicators of Chemical Bonding. “Proceedings of the VIII-th International Symposium Young People and Multidisciplinary Research”, Timişoara; Welding Publishing House, Association of Multidisciplinary Research of the West Zone of Romania: Timişoara, Romania, 11–12 Mai 2006; pp. 632–641, (ISBN-10 973-8359-39-2, ISBN-13 978-8359-39-0). [Google Scholar]
  67. Antolin, J; Zarzo, A; Angulo, JC. Maximum Entropy Analysis of the Electron-Pair Density in Many-Electron Systems. Phys. Rev. A 1994, 50, 240–246. [Google Scholar]
  68. Lee, C; Parr, RG. Gaussian and Other Approximations to the First-Order Density Matrix of Electronic System, and the Derivation of Various Local-Density-Functional-Theories. Phys. Rev. A 1987, 35, 2377–2383. [Google Scholar]
  69. Zhao, Q; Parr, RG. Quantities Ts[n] and Tc[n] in Density-Functional Theory. Phys. Rev. A 1992, 46, 2337–2342. [Google Scholar]
  70. Ou-Yang, H; Levy, M. Nonuniform Coordinate Scaling Requirements in Density-Functional Theory. Phys. Rev. A 1990, 42, 155–159. [Google Scholar]
  71. DePristo, AE; Kress, JD. Kinetic-Energy Functionals via Padé Approximations. Phys. Rev. A 1987, 35, 438–441. [Google Scholar]
  72. Liberman, DA; Albritton, JR; Wilson, BG; Alley, WE. Self-Consistent-Field Calculations of Atoms and Ions Using a Modified Local-Density Approximation. Phys. Rev. A 1994, 50, 171–176. [Google Scholar]
  73. Levy, M; Ernzerhof, M; Gorling, A. Exact Local Exchange Potential from Fock Equations at Vanishing Coupling Constant, and δTc / δn from Wave-Function Calculations at Full Coupling Constant. Phys. Rev. A 1996, 53, 3963–3973. [Google Scholar]
  74. Cedillo, A; Robles, J; Gazquez, JL. New Nonlocal Exchange-Energy Functional from a Kinetic-Energy-Density Padé-Approximant Model. Phys. Rev. A 1988, 38, 1697–1701. [Google Scholar]
  75. Perdew, JP; Yue, W. Accurate and Simple Density Functional for the Electronic Exchange Energy: Generalized Gradient Approximation. Phys. Rev. B 1986, 33, 8800–8802. [Google Scholar]
  76. Manoli, SD; Whitehead, MA. Generalized-Exchange Local-Spin-Density-Functional Theory: Calculation and Results for Non-Self-Interaction-Corrected and Self-Interaction-Corrected Theories. Phys. Rev. A 1988, 38, 3187–3199. [Google Scholar]
  77. Becke, AD. Density Functional Calculations of Molecular Bond Energies. J. Chem. Phys. 1986, 84, 4524–4529. [Google Scholar]
  78. Ghosh, SK; Parr, RG. Phase-Space Approach to the Exchange Energy Functional of Density-Functional Theory. Phys. Rev. A 1986, 34, 785–791. [Google Scholar]
  79. Bartolotti, LJ. A new Gradient Expansion of the Exchange Energy to Be Used in Density Functional Calculations on Atoms. J. Chem. Phys. 1982, 76, 6057–6059. [Google Scholar]
  80. Perdew, JP; Chevary, JA; Vosko, SH; Jackson, KA; Pederson, MR; Singh, DJ; Fiolhais, C. Atoms, Molecules, Solids, and Surfaces: Applications of the Generalized Gradient Approximation for Exchange and Correlation. Phys. Rev. B 1992, 46, 6671–6687. [Google Scholar]
  81. Wang, Y; Perdew, JP; Chevary, JA; Macdonald, LD; Vosko, SH. Exchange Potentials in Density-Functional Theory. Phys. Rev. A 1990, 41, 78–85. [Google Scholar]
  82. Alonso, JA; Girifalco, LA. Nonlocal Approximation to the Exchange Potential and Kinetic Energy of an Inhomogeneous Electron Gas. Phys. Rev. B 1978, 17, 3735–3743. [Google Scholar]
  83. Levy, M; Gorling, A. Density-Functional Exchange Identity from Coordinate Scaling. Phys. Rev. A 1996, 53, 3140–3150. [Google Scholar]
  84. Becke, AD. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar]
  85. Taut, M. Generalized Gradient Correction for Exchange: Deduction from the Oscillator Model. Phys. Rev. A 1996, 53, 3143–3150. [Google Scholar]
  86. Slater, JC. A Simplification of the Hartree-Fock Method. Phys. Rev. 1951, 81, 385–390. [Google Scholar]
  87. Slater, JC; Johnson, KH. Self-Consistent-Field Xα Cluster Method for Polyatomic Molecules and Solids. Phys. Rev. B 1972, 5, 844–853. [Google Scholar]
  88. Lee, C; Zhou, Z. Exchange-Energy Density Functional: Reparametrization of Becke’s Formula and Derivation of Second-Order Gradient Correction. Phys. Rev. A 1991, 44, 1536–1539. [Google Scholar]
  89. Senatore, G; March, NH. Recent Progress in the Field of Electron Correlation. Rev. Mod. Phys. 1994, 66, 445–479. [Google Scholar]
  90. Liu, S; Nagy, A; Parr, RG. Expansion of the Density-Functional Energy Components Ec and Tc in Terms of Moments of the Electron Density. Phys. Rev. A 1999, 59, 1131–1134. [Google Scholar]
  91. Ou-Yang, H; Levy, M. Theorem for Functional Derivatives in Density-Functional Theory. Phys. Rev. A 1991, 44, 54–58. [Google Scholar]
  92. Nagy, A; Liu, S; Parr, RG. Density-Functional Formulas for Atomic Electronic Energy Components in Terms of Moments of the Electron Density. Phys. Rev. A 1999, 59, 3349–3354. [Google Scholar]
  93. Perdew, JP. Density-Functional Approximation for the Correlation Energy of the Inhomogeneous Electron Gas. Phys. Rev. B 1986, 33, 8822–8824, with Erratum, Phys. Rev. B 1986, 34, 7406. [Google Scholar]
  94. Wang, Y; Perdew, JP. Spin Scaling of the Electron-Gas Correlation Energy in the High-Density Limit. Phys. Rev. B 1989, 43, 8911–8916. [Google Scholar]
  95. Seidl, M; Perdew, JP; Levy, M. Strictly Correlated Electrons in Density-Functional Theory. Phys. Rev. A 1999, 59, 51–54. [Google Scholar]
  96. Perdew, JP; Burke, K; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar]Zhang, Y; Yang, W. Comment on: Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1998, 80, 890. [Google Scholar](c)Perdew, Burke, and Ernzerhof Replay, idem, 891.
  97. Perdew, JP; Ernzerhof, M; Zupan, A; Burke, K. Nonlocality of the Density Functional for Exchange and Correlation: Physical Origins and Chemical Consequences. J. Chem. Phys. 1998, 108, 1522–1531. [Google Scholar]
  98. Wilson, LC; Levy, M. Nonlocal Wigner-Like Correlation-Energy Density Functional Through Coordinate Scaling. Phys. Rev. B 1990, 41, 12930–12932. [Google Scholar]
  99. Vosko, SJ; Wilk, L; Nusair, M. Accurate Spin-Dependent Electron Liquid Correlation Energies for Local Spin Density Calculations: A Critical Analysis. Can. J. Phys. 1980, 58, 1200–1211. [Google Scholar]
  100. Perdew, JP. Density-Functional Approximation for the Correlation Energy of the Inhomogeneous Electron Gas. Phys. Rev. B 1986, 33, 8822–8824. [Google Scholar]
  101. Savin, A; Stoll, H; Preuss, H. An Application of Correlation Energy Density Functionals to Atoms and Molecules. Theor Chim Acta 1986, 70, 407–419. [Google Scholar]
  102. Savin, A; Preuss, H; Stoll, H. Non-Local Effects on Atomic and Molecular Correlation Energies Studies with a Gradient-Corrected Density Functional. In “Density Matrices and Density Functionals”; Erhahy, R, Smith, VH, Eds.; Reidel Publishing Company, 1987; pp. 457–465. [Google Scholar]
  103. Rasolt, M; Geldart, DJW. Exchange and Correlation Energy in a Nonuniform Fermion Fluid. Phys. Rev. B 1986, 34, 1325–1328. [Google Scholar]
  104. Savin, A; Wedig, U; Preuss, H; Stoll, H. Molecular Correlation Energies Obtained with a Nonlocal Density Functional. Phys. Rev. Lett. 1984, 53, 2087–2089. [Google Scholar]
  105. Lee, C; Yang, W; Parr, RG. Development of the Colle-Salvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar]
  106. Perdew, JP; Zunger, A. Self-Interaction Correction to Density-Functional Approximations for Many-Electron System. Phys. Rev. B 1981, 23, 5048–5079. [Google Scholar]
  107. Lee, C; Parr, RG. Exchange-Correlation Functional for Atoms and Molecules. Phys. Rev. A 1990, 42, 193–199. [Google Scholar]
  108. Berkowitz, M. Exponential Approximation for the Density Matrix and the Wigner’s Distribution. Chem Phys. Lett. 1986, 129, 486–488. [Google Scholar]
  109. Lee, H; Bartolotti, LJ. Exchange and Exchange-Correlation Functionals Based on the Gradient Correction of the Electron Gas. Phys. Rev. A 1991, 44, 1540–1542. [Google Scholar]
  110. Tozer, DJ; Handy, NC. The Development of New Exchange-Correlation Functionals. J. Chem. Phys. 1998, 108, 2545–2555. [Google Scholar]
  111. Mandl, F. Quantum Mechanics; John Wiley & Sons: Chichester, 1992. [Google Scholar]
  112. Miller, JN; Miller, JC. Statistics and Chemometrics for Analytical Chemistry, fourth edition; Prentice Hall: Harlow, 2000. [Google Scholar]
  113. Lyons, L. A Practical Guide to Data Analysis for Physical Science Students; Cambridge University Press: Cambridge, 1993. [Google Scholar]
  114. Koopmans, T. Uber die Zuordnung von Wellen Funktionen und Eigenwerter zu den Einzelnen Elektronen Eines Atom. Physica 1934, 1, 104–113. [Google Scholar]
  115. Parr, RG; Donnelly, RA; Levy, M; Palke, WE. Electronegativity: the Density Functional Viewpoint. J. Chem. Phys. 1978, 68, 3801–3807. [Google Scholar]
  116. Parr, RG; Pearson, RG. Absolute Hardness: Companion Parameter to Absolute Electronegativity. J. Am. Chem. Soc. 1983, 105, 7512–7516. [Google Scholar]
  117. Chattaraj, PK; Liu, GH; Parr, RG. The Maximum Hardness Principle in the Gyftpoulos-Hatsopoulos Three-Level Model for an Atomic or Molecular Species and Its Positive and Negative Ions. Chem. Phys. Lett. 1995, 237, 171–176. [Google Scholar]
  118. Ayers, PW; Parr, RG. Variational Principles for Describing Chemical Reactions: the Fukui Function and Chemical Hardness Revisited. J. Am. Chem. Soc. 2000, 122, 2010–2018. [Google Scholar]
  119. Putz, MV. Maximum Hardness Index of Quantum Acid-Base Bonding. MATCH Commun. Math. Comput. Chem. accepted.
  120. Grabo, T; Gross, EKU. Density-Functional Theory Using an Optimized Exchange-Correlation Potential. Chem. Phys. Lett. 1995, 240, 141–150. [Google Scholar]
  121. Mulliken, RS. A New Electroaffinity Scale: Together with Data on Valence States and An Ionization Potential and Electron Affinities”. J. Chem. Phys. 1934, 2, 782–793. [Google Scholar]
  122. Lackner, KS; Zweig, G. Introduction to the Chemistry of Fractionally Charged Atoms: Electronegativity. Phys. Rev. D. 1983, 28, 1671–1691. [Google Scholar]
  123. Putz, MV. Systematic Formulation for Electronegativity and Hardness and Their Atomic Scales within Density Functional Softness Theory. Int. J. Quantum Chem. 2006, 106, 361–386. [Google Scholar]
  124. Putz, MV. Semiclassical Electronegativity and Chemical Hardness. J. Theor. Comp. Chem. 2007, 6, 33–47. [Google Scholar]
  125. Bratsch, SG. A Group Electronegativity Method with Pauling Units. J. Chem. Educ. 1985, 62, 101–103. [Google Scholar]
  126. Putz, MV; Russo, N; Sicilia, E. On The Applicability of The HSAB Principle through The Use of Improved Computational Schemes for Chemical Hardness Evaluation. J. Comput. Chem. 2004, 25, 994–1003. [Google Scholar]
  127. Gross, EKU; Oliveira, LN; Kohn, W. Density-Functional Theory for Ensembles of Fractionally Occupied States.I. Basic Formalism. Phys. Rev. A 1988, 37, 2809–2820. [Google Scholar]
  128. Oliviera, LN; Gross, EKU; Kohn, W. Density-Functional Theory for Ensembles of Fractionally Occupied States. II. Applications to the He Atom. Phys. Rev. A 1988, 37, 2821–2833. [Google Scholar]
  129. Nesbet, RK. Fractional Occupation Numbers in Density-Functional Theory. Phys. Rev. A 1997, 56, 2665–2668. [Google Scholar]
  130. Putz, MV. Absolute and Chemical Electronegativity and Hardness; Nova Publishers Inc.: New York, USA, 2008; https://www.novapublishers.com/catalog/product_info.php?products_id=7678.
Figure 1. Left: the anti-ferromagnetic structure CoO; right: the band structure and the density of state (DOS) in LSDA and GGA approximations, respectively [53]. The upper and down arrows are associated with the spin orbital projections.
Figure 1. Left: the anti-ferromagnetic structure CoO; right: the band structure and the density of state (DOS) in LSDA and GGA approximations, respectively [53]. The upper and down arrows are associated with the spin orbital projections.
Ijms 09 01050f1
Figure 2. The localization domains for Li (left) and Sc (right) crystals based on the electronic localization function ELF (59) [55].
Figure 2. The localization domains for Li (left) and Sc (right) crystals based on the electronic localization function ELF (59) [55].
Ijms 09 01050f2
Figure 3. Comparison of the radial density given by (71) with the electron localization function (65) with components (66) for the simplified self-consistent approximation (67)–(70) for Li atomic structure [66].
Figure 3. Comparison of the radial density given by (71) with the electron localization function (65) with components (66) for the simplified self-consistent approximation (67)–(70) for Li atomic structure [66].
Ijms 09 01050f3
Figure 4. Comparative analysis of the charge density contours, electronic localization functions (ELFs), and radial densities for the H (dashed lines), F, Cl, Br, and I (full lines) atoms in molecular combinations HF, HCl, HBr, and HI, respectively [24].
Figure 4. Comparative analysis of the charge density contours, electronic localization functions (ELFs), and radial densities for the H (dashed lines), F, Cl, Br, and I (full lines) atoms in molecular combinations HF, HCl, HBr, and HI, respectively [24].
Ijms 09 01050f4
Figure 5. Charts of energetic functionals showing their reactive efficiency (196) as the height of bars for the parabolic expansion and of their statistical correlation (192) as the width of bars, in terms of electronegativity and chemical hardness, resuming atomic analysis of Tables IIII - in (a), and molecular one from Tables IV and V - in (b), by employing representative values of Tables VI and VII (highlighted on grey background), respectively.
Figure 5. Charts of energetic functionals showing their reactive efficiency (196) as the height of bars for the parabolic expansion and of their statistical correlation (192) as the width of bars, in terms of electronegativity and chemical hardness, resuming atomic analysis of Tables IIII - in (a), and molecular one from Tables IV and V - in (b), by employing representative values of Tables VI and VII (highlighted on grey background), respectively.
Ijms 09 01050f5
Table I. Atomic kinetic, exchange, and correlation, energies (in hartrees) from various schemes of computations. The exact values are computed with Hartree-Fock densities.
Table I. Atomic kinetic, exchange, and correlation, energies (in hartrees) from various schemes of computations. The exact values are computed with Hartree-Fock densities.
AtomsKinetic energyExchange energyCorrelation energy
TexactT0T0+T2TPadéKexactK0KB88EcexactEc(139)Ec
He2.861682.560542.878502.87639−1.0260−0.884−1.025−0.0425−0.0215−0.0681
Li7.432736.700627.505047.44941−1.7812−1.538−1.775−0.0454−0.0486−0.0815
Be14.573013.128614.646614.4223−2.6669−2.312−2.658−0.0945−0.0820−0.1192
B24.529122.072024.522824.2089−3.7438−3.272−3.728−0.1247−0.1197−0.1625
C37.688634.014437.598837.2533−5.0444−4.459−5.032−0.1566−0.1609−0.2091
N54.400949.477154.385254.0643−6.5971−5.893−6.589−0.1850−0.2050−0.2567
O74.809467.896574.357374.1625−8.1752−7.342−8.169−0.2579−0.2512−0.3035
F99.409390.459898.642998.6959−10.003−9.052−10.02−0.332−0.2996−0.3510
Ne128.547117.761127.829128.221−12.108−11.03−12.14−0.390−0.3498−0.3987
Na161.859148.809161.093161.718−14.017−12.79−14.03−0.398−0.3892−0.4137
Mg199.614184.017198.749199.578−15.994−14.61−16.00−0.443−0.4351−0.4491
Al241.877223.443240.868242.008−18.069−16.53−18.06−0.480−0.4809−0.4863
Si288.854267.315287.659289.139−20.280−18.59−20.27−0.521−0.5308−0.5308
Cl459.482426.865457.321460.117−27.512−25.35−27.49−0.714−0.6901−0.6710
Ar526.817490.017524.289527.617−30.185−27.86−30.15−0.787−0.7459−0.7190
: from Ref. [71] and references therein;
: from Ref. [107] and references therein;
: from Ref. [84];
: from eq. (139) and Ref. [90];
: from fitting equation Ec=−0.04682N + 0.005753<ρ2/3>−0.00096<ρ1/3>, see Ref. [90] and references therein
Table II. Atomic exchange-correlation and total energies (in hartrees) from various schemes of computations. The exact values are computed with Hartree-Fock densities.
Table II. Atomic exchange-correlation and total energies (in hartrees) from various schemes of computations. The exact values are computed with Hartree-Fock densities.
AtomsExchange-Correlation energyTotal energy
ExcexactExcI(Xα)ExcII(Xα)ExcI(Wig)ExcII(Wig)EtotexactEtotxc(RG)Etotxc(LDA)EtotBLYPEtotPW91
He−1.0685−1.0604−1.0566−1.0633−1.0654−2.9042−3.0317−2.8601−2.9071−2.9000
Li−1.8266−1.8048−1.8134−1.8093−1.8108−7.4781−7.6473−7.3704−7.4827−7.4742
Be−2.7614−2.7260−2.7522−2.7325−2.7342−14.6675−14.8911−14.4966−14.6615−14.6479
B−3.8685−3.8126−3.8415−3.8215−3.8177−24.6538−24.9158−24.4097−24.6458−24.6299
C−5.2010−5.1127−5.1338−5.1248−5.1121−37.8163−38.1305−37.5095−37.8430−37.8265
N−6.7821−6.6400−6.6440−6.6558−6.6321−54.4812−54.8681−54.1287−54.5932−54.5787
O−8.4331−8.3599−8.3405−8.3796−8.3450−75.0271−75.4597−74.5979−75.0786−75.0543
F−10.325−10.327−10.277−10.350−10.305−99.741−100.235−99.247−99.7581−99.7316
Ne−12.498−12.551−12.466−12.579−12.524−128.937−129.522−128.403−128.9730−128.9466
Na−14.415−14.462−14.382−14.488−14.445−162.257−162.862−161.624−162.293−162.265
Mg−16.437−16.482−16.424−16.504−16.484−200.058−200.705−199.340−200.093−200.060
Al−18.549−18.566−18.542−18.583−18.593−242.357−243.028−241.533−242.380−242.350
Si−20.801−20.774−20.791−20.784−20.830−289.356−290.063−288.435−289.388−289.363
Cl−28.226−28.115−28.272−28.092−28.281−460.196−461.005−458.963−460.165−460.147
Ar−30.972−30.827−31.037−30.789−31.035−527.605−528.452−526.267−527.551−527.539
: from Ref. [107];
: from Eqs. (170)(173) and Ref. [107];
: from Ref. [109];
: from Ref. [120]
Table III. Values (in hartrees) of the structural indices electronegativity (χ), chemical hardness (η), in finite-difference [121, 122], density functional [123], and semiclassical [124] modes of computations for atoms of Tables I and II.
Table III. Values (in hartrees) of the structural indices electronegativity (χ), chemical hardness (η), in finite-difference [121, 122], density functional [123], and semiclassical [124] modes of computations for atoms of Tables I and II.
LevelFinite-DifferenceFunctionalSemiclassical
AtomsχFDηFDχDFTηDFTχSCηSC
He0.450940.458661.211321.661890.570380.2172
Li0.110990.161340.151050.087840.004120.00334
Be0.126060.217940.442480.445790.008930.0047
B0.156560.149211.153621.341050.015260.00588
C0.229330.183392.763322.96950.022790.00684
N0.256160.278945.795664.913630.031390.0076
O0.278940.2256610.65055.916940.040720.00816
F0.382210.2598316.91294.377070.050610.00849
Ne0.393610.4013223.7119−0.087470.060790.00864
Na0.102900.106210.231530.187430.000110.00005
Mg0.095550.183390.498710.531420.000180.00007
Al0.118340.103271.046311.198820.000260.00008
Si0.172000.126062.108052.307240.000360.00009
Cl0.305770.1712011.57667.76920.000740.00012
Ar0.282990.2980617.88319.088570.000880.00013
Table IV. Molecular kinetic, exchange, correlation, exchange-correlation, and total energies (in hartrees) from various schemes of computations. The exact values are computed with HF densities.
Table IV. Molecular kinetic, exchange, correlation, exchange-correlation, and total energies (in hartrees) from various schemes of computations. The exact values are computed with HF densities.
MoleculesKinetic
Exchange
Correlation
Exch.-corr.
Total energy
T0T0+T2KexactKPBEEcVWNEcGCPExcexactExcPBEEtotBLYPEtotTH
H21.1401.125−0.657−0.648−95·10−3−47·10−3−0.698−0.691−1.169−1.178
LiH7.9788.003−2.125−2.105−219·10−3−93·10−3−2.212−2.188−8.068−8.070
CH440.05040.141−6.576−6.536−593·10−3−328·10−3−6.883−6.836−40.502−40.515
H2O76.15075.477−8.910−8.917−664·10−3−365·10−3−9.292−9.241−76.448−76.433
HF100.13799.242−10.378−10.385−704·10−3−380·10−3−10.779−10.720−100.48−100.455
N2109.115108.242−13.094−13.128−945·10−3−506·10−3−13.665−13.580−109.559−109.54
O2149.843148.369−16.290−16.358−1110·10−3−599·10−3−16.958−16.887−150.384−150.337
F2198.892196.729−19.872−19.951−1302·10−3−697·10−3−20.661−20.564−199.599−199.533
: from Ref. [97];
: from Ref. [101];
: from Ref. [110]
Table V. Values (in hartrees) of the structural indices electronegativity (χ), chemical hardness (η), compute by means of the group method [14, 125] within the finite-difference [121, 122], density functional [123], and semiclassical [124] modes of computations for molecules of Table IV.
Table V. Values (in hartrees) of the structural indices electronegativity (χ), chemical hardness (η), compute by means of the group method [14, 125] within the finite-difference [121, 122], density functional [123], and semiclassical [124] modes of computations for molecules of Table IV.
Level Molecule sFinite-Difference
Functional
Semiclassical
χFDηFDχDFTηDFTχSCηSC
H20.263870.23700.263840.237040.263870.23705
LiH0.156260.1920.192120.128180.008110.006596
CH40.256160.22390.322160.290510.084680.03064
H2O0.268710.23310.390970.348590.093350.0229
HF0.30770.24790.519640.449740.084930.01639
N20.256160.27895.795664.913630.031390.00761
O20.278940.225710.65055.916940.040720.00816
F20.368980.259816.91294.377070.050610.00849
Table VI. Coefficients in bilinear correlation of the energies of Tables I and II against the electronegativity and chemical hardness of Table III within experimental finite difference (FD), density functional theory (DFT), and semiclassical (SC) modes of computations, respectively. The deviation from the parabolic expansion E=a+bχ+cη in terms of the index σπ =sign(b) c/b2 as well as the correlation factor of the QSPR model (r) are both calculated. The color code indicates: the (+, +) for magenta background and (−, +) for green background sign combinations of (b, c) in (191), respectively, while the grey background is used in highlighting the representative (σπ, r) couple for each energy type of employed functionals.
Table VI. Coefficients in bilinear correlation of the energies of Tables I and II against the electronegativity and chemical hardness of Table III within experimental finite difference (FD), density functional theory (DFT), and semiclassical (SC) modes of computations, respectively. The deviation from the parabolic expansion E=a+bχ+cη in terms of the index σπ =sign(b) c/b2 as well as the correlation factor of the QSPR model (r) are both calculated. The color code indicates: the (+, +) for magenta background and (−, +) for green background sign combinations of (b, c) in (191), respectively, while the grey background is used in highlighting the representative (σπ, r) couple for each energy type of employed functionals.
Method of
QSPR results
Method of
QSPR results
Energ y(χ, η)abcσπrEnergy(χ, η)abcσπr
Texact
FD194.59741.55−951.38−0.00170.33Excexact
FD−15.30−44.8060.89−0.03030.37
DFT51.732.0931.517.20.62DFT−6.56−0.19−1.5340.840.58
SC
186.38
−2315.7
5147.4
−0.001
0.35
SC
−13.76
73.07
−128.0
−0.024
0.38
T0
FD179.67687.26−881.14−0.00190.33ExcI(Xα)
FD−15.28−44.5660.69−0.03060.37
DFT46.991.9129.418.070.62DFT−6.56−0.2−1.5139.110.57
SC
172.42
−2185.3
4874.1
−0.001
0.35
SC
−13.72
72.04
−125.5
−0.024
0.38
T0+T2
FD193.83736.94−946.29−0.00170.33ExcII(Xα)
FD−15.28−44.8160.92−0.03030.37
DFT51.592.0731.367.330.62DFT−6.52−0.19−1.5340.810.57
SC
185.59
−2312.1
5142
−0.001
0.35
SC
−13.75
74.72
−132.4
−0.024
0.38
TPadé
FD194.61741.17−951.43−0.00170.33ExcI(Wig)
FD−15.29−44.5360.67−0.03060.37
DFT51.542.0831.587.30.62DFT−6.58−0.2−1.5038.780.57
SC
186.42
−2334.4
5197.4
−0.001
0.35
SC
−13.72
71.4
−123.8
−0.024
0.38
Kexact
FD−14.91−43.6259.38−0.03120.37ExcII(Wig)
FD−15.30−44.8060.95−0.03040.37
DFT−6.37−0.19−1.4943.170.57DFT−6.54−0.2−1.5239.980.57
SC
−13.4
72.74
−128.78
−0.0243
0.38
SC
−13.76
74.1
−130.7
−0.024
0.38
K0
FD−13.59−40.2954.69−0.03370.37Etotexact
FD−194.98−742.7952.91−0.00170.33
DFT−5.72−0.17−1.3946.960.58DFT−51.91−2.10−31.547.150.62
SC
−12.24
68.77
−123.42
−0.026
0.38
SC
−186.73
2316
−5148
−0.001
0.35
KB88
FD−14.89−43.6159.33−0.03120.37Etotxc(RG)
FD−195.55−743.9954.5−0.00170.33
DFT−6.37−0.19−1.4942.350.57DFT−52.27−2.11−31.567.10.62
SC
−13.39
71.88
−126.54
−0.0245
0.38
SC
−187.24
2315.
−5143.
−0.001
0.35
Ecexact
FD−0.39−1.201.52−1.0550.38Etotxc(LDA)
FD−194.26−740.7950.13−0.00170.33
DFT−0.19−0.0075−0.03595.980.59DFT−51.59−2.09−31.477.180.62
SC
−0.362
0.232
0.997
18.528
0.38
SC
−186.1
2313.
−5142.
−0.001
0.35
Ec(139)
FD−0.40−1.131.55−1.2070.39EtotBLYP
FD−194.99−742.7952.86−0.00170.33
DFT−0.19−0.0056−0.031081.30.57DFT−51.94−2.1−31.537.150.62
SC
−0.356
0.851
−0.563
−0.778
0.41
SC
−186.74
2315.
−5145
−0.001
0.35
EcFD−0.41−1.101.44−1.18080.40EtotPW91FD−194.97−742.7952.81−0.00170.33
DFT−0.22−0.0063−0.03735.50.59DFT−51.92−2.1−31.547.150.62
SC−0.374−0.2242.113−42.240.42SC−186.72315.−5145.−0.0010.35
,
: from associate energetic atomic values of Table I
Table VII. Coefficients in bilinear correlation of the energies of Table IV against the respective electronegativity and chemical hardness of the Table V. The computational meaning of the output values and the color code are the same as those of Table VI caption.
Table VII. Coefficients in bilinear correlation of the energies of Table IV against the respective electronegativity and chemical hardness of the Table V. The computational meaning of the output values and the color code are the same as those of Table VI caption.
Method of
QSPR results
Method of
QSPR results
Energy yχ & ηabcσπrEnergy yχ & ηabcσπr
T0
FD−192.71821.61238.590.000350.77T0+T2
FD−190.58811.48237.960.00040.77
DFT40.4468.1254.5010.0680.89DFT40.1998.0124.4970.070.89
SC
82.14
542.16
−977.9
−0.0033
0.56
SC
81.4
537.74
−969.47
−0.003
0.56
Kexact
FD18.18−70.1−38.020.00770.74KPBE
FD18.324−70.48−38.250.00770.74
DFT−5.234−0.646−0.8051.930.89DFT−5.219−0.65−0.8111.920.89
SC
−9.78
−48.27
94.96
−0.04
0.61
SC
−9.81
−48.15
95.
−0.04
0.61
EcVWN
FD1.039−4.06−2.740.1670.71EcGCP
FD0.59−2.27−1.4960.2890.72
DFT−0.423−0.035−0.0647.830.87DFT−0.225−0.019−0.03394.080.86
SC
−0.71
−3.13
6.22
−0.63
0.64
SC
−0.36
−2.1
3.7
−0.88
0.64
Excexact
FD18.95−72.63−40.10.00760.74ExcPBE
FD18.85−72.45−39.620.00750.74
DFT−5.457−0.668−0.8461.8970.89DFT−5.421−0.666−0.841.8980.89
SC−10.19−50.1698.71−0.0390.61SC−10.13−49.9798.28−0.040.61
EtotBLYPFD
193.12
−824.17
−238.96
0.00035
0.77
EtotTH
FD
193.06
−823.8
−239.04
0.0004
0.77
DFT−40.67−8.149−4.5140.0680.89DFT−40.67−8.145−4.5140.0680.89
SC−82.48−544.35981.7−0.0030.56SC−82.46−544.3981.5−0.0030.56

Share and Cite

MDPI and ACS Style

Putz, M.V. Density Functionals of Chemical Bonding. Int. J. Mol. Sci. 2008, 9, 1050-1095. https://doi.org/10.3390/ijms9061050

AMA Style

Putz MV. Density Functionals of Chemical Bonding. International Journal of Molecular Sciences. 2008; 9(6):1050-1095. https://doi.org/10.3390/ijms9061050

Chicago/Turabian Style

Putz, Mihai V. 2008. "Density Functionals of Chemical Bonding" International Journal of Molecular Sciences 9, no. 6: 1050-1095. https://doi.org/10.3390/ijms9061050

APA Style

Putz, M. V. (2008). Density Functionals of Chemical Bonding. International Journal of Molecular Sciences, 9(6), 1050-1095. https://doi.org/10.3390/ijms9061050

Article Metrics

Back to TopTop