Email updates

Keep up to date with the latest news and content from Chemistry Central Journal and Chemistry Central.

Open Access Research article

Valence atom with bohmian quantum potential: the golden ratio approach

Mihai V Putz

Author Affiliations

Laboratory of Computational and Structural Physical Chemistry, Biology-Chemistry Department, West University of Timişoara, Pestalozzi Street No.16, Timişoara, RO-300115, Romania

Chemistry Central Journal 2012, 6:135  doi:10.1186/1752-153X-6-135

The electronic version of this article is the complete one and can be found online at: http://journal.chemistrycentral.com/content/6/1/135


Received:6 September 2012
Accepted:29 October 2012
Published:12 November 2012

© 2012 Putz; licensee Chemistry Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

The alternative quantum mechanical description of total energy given by Bohmian theory was merged with the concept of the golden ratio and its appearance as the Heisenberg imbalance to provide a new density-based description of the valence atomic state and reactivity charge with the aim of clarifying their features with respect to the so-called DFT ground state and critical charge, respectively.

Results

The results, based on the so-called double variational algorithm for chemical spaces of reactivity, are fundamental and, among other issues regarding chemical bonding, solve the existing paradox of using a cubic parabola to describe a quadratic charge dependency.

Conclusions

Overall, the paper provides a qualitative-quantitative explanation of chemical reactivity based on more than half of an electronic pair in bonding, and provide new, more realistic values for the so-called “universal” electronegativity and chemical hardness of atomic systems engaged in reactivity (analogous to the atoms-in-molecules framework).

Keywords:
Electronegativity; Chemical hardness; Bohmian mechanics; Heisenberg imbalance equation; Slater electronic density

Graphical abstract

Introduction

Recently, the crucial problem regarding whether chemical phenomena are reducible to physical ones has had an increasingly strong impact on the current course of conceptual and theoretical chemistry. For instance, the fact that elements arrange themselves in atomic number (Z) triads in approximately 50% of the periodic system seems to escape custom ordering quantifications [1,2]. The same applies to the following: the fascinating golden ratio (τ) limit for the periodicity of nuclei beyond any physical first-principle constants, which provides specific periodic laws for the chemical realm [3-6]; the fact that atoms have no definite atomic radii in the sense of a quantum operator, and even the Aufbau principle, which, although chemically workable, seems to violate the Pauli Exclusion Principle [7]; at the molecular level, the well-celebrated reaction coordinate, which, although formally defined in the projective energy space, does not constitute a variable to drive optimization in the course of chemical reactions, appearing merely as a consequence of such reactions [8]; the problem of atoms in molecules [9], i.e., how much of the free atoms enter molecules and how much independency the atoms preserve in bonding; and chemical bonding itself, which ultimately appears to be reinterpreted as a special case of bosonic condensation with the aid of bondons – the quantum bosons of chemical bonding, which, without being elementary, imbue chemical compounds with a specific reality [10,11].

In the same context, the specific measure of chemical reactivity, electronegativity (χ), which lacks a definite quantum operator but retains an observable character through its formal identity with the macroscopic chemical potential χ=-μ[12,13], was tasked with carrying quantum information within the entanglement environment of Bohmian mechanics [14-17] and has thus far been identified with the square root of the so-called quantum potential χ = VQ1/2[6].

However, the striking difference between an atom as a physical entity, with an equal number of electrons and protons (thus in equilibrium), and the same atom as a chemical object, with incomplete occupancy in its periphery quantum shells (thus attaining equilibrium by changing accepting or releasing electrons), is closely related to the electronegativity phenomenology in modeling chemical reactivity. Moreover, this difference triggers perhaps the most important debate in conceptual chemistry: the ground vs. valence state definition of an atom.

The difficulty may be immediately revealed by considering the variation in the total energy (of the ground and/or valence state – see below for an explanation of their difference) around the physical equilibrium (neutral atom) attained between the release (by ionization, I) and receipt (through affinity, A) of electrons toward chemical equilibrium (in molecules, chemical bonding). Accordingly, the curve passing through these points apparently only behaves as shown in Figure 1(a), while in all systems (with numerical I and A), the obtained interpolating curve presents a minimum toward accepting electrons (see Figure 1(b)), thus confirming the electronegativity concept as a chemical reality, although with a predicted fractional charge (for example, the critical charge N*) on an atom at chemical equilibrium (i.e., not reducible/comprehensible to/by an ordinary physical description of atoms).

thumbnailFigure 1. The two energy curves (thick lines) for the quantum atom in (a) the apparent or reactive ground state and (b) the shifted or critical ground state.

However, the physical-to-chemical paradox continues in an even more exciting fashion as follows. When, in light of the above discussion, electronegativity is recognized with the two-point limits shown in Figure 1(b), namely [13,18]

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M1">View MathML</a>

(1)

the limits represent tangents to a curve that does not describe chemical equilibrium but an excited state driven by the parabolic form

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M2">View MathML</a>

(2)

which happens to correspond to the celebrated density functional theory (DFT) working energy expression [13,19-21] written in terms of electronegativity and chemical hardness, respectively defined as follows [13,22-25]:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M3">View MathML</a>

(3)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M4','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M4">View MathML</a>

(4)

The point is that curve (2) is not chemically minimized, although it is very often assumed to be in the DFT invoked by the chemical reactivity literature [13,26-29]; however, the curve cannot be considered indicative of a sort of ground state (neither reactive nor critical states of Figure 1). Additionally, by comparing the curves of Figure 1 (a) and (b), the curve of eq. (2) occurs above both the reactive and critical curves of Figure 1; it thus should represent the chemical valence state with which to operate. Therefore, much caution should be taken when working with eq. (2) in assessing the properties of atoms, molecules, atoms in molecules, etc. Nevertheless, this is another case of chemistry not being reducible to physics and should be treated accordingly. It is worth noting that Parr, the “father” of eq. (2) and a true pioneer of conceptual density functional theory [30,31], had tried to solve this dichotomy by taking the “valence as the ground state of an atom in a perturbed environment”. This statement is not entirely valid because perturbation is not variation such that it may be corrected by applying the variational principle to eq. (2), for example. In fact, using such variation should be considered a double variational technique that is necessary to arrive at the celebrated chemical reactivity principles of electronegativity and chemical hardness, as recently shown [32].

The current line of work takes a step forward by employing the double variation of the parabolic energy curve of type (2) to provide the quantum (DFT) valence charge of an atom (say, N**) and to compare it either quantitatively and qualitatively with the chemical critical charge N*. The goal of these efforts is to gain new insight into the valence state and chemical reactivity at the quantum level. To this end, the relation of Bohmian mechanics to the concept of the golden ratio will be essential and will be introduced in the following.

The consequences of the joint consideration of Bohmian mechanics and the golden ratio for the main atomic systems will be explored, and the quantum chemical valence state will be accordingly described alongside the so-called universal electronegativity and chemical hardness, refining the work of Parr and Bartolotti [33] as well as generalizing the previous Bohmian-Boeyens approach [3,4].

Background methods

Two apparently disjoint theories of matter will be employed to characterize the quantum valence of an atom: Bohmian mechanics – furnishing the main equation for total energy – and the fundamental quantum mechanics through the Heisenberg combined with de Broglie principles providing the wave-particle indeterminacy framework in which the golden ration dependency of Z/N naturally appears as quantifying the valence states of atoms considered the “ground state” of the atomic chemical reactivity.

Bohmian mechanics

Because of the need to reduce Copenhagen’s indeterminacy for quantum phenomena, i.e., by associating it the quantum description of “Newtonian” forms of motion, though by preserving probability densities, quantum averages, etc., the so-called “minimalist” quantum theory may be formulated following the Bohm quantum mechanical program as follows.

One begins with the general eikonal wave-function form [14]

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M5','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M5">View MathML</a>

(5)

which represents the mid-way between wave and particle mechanics because it contains both information regarding Hamilton-Jacobi theory and the Wentzel-Kramers-Brillouin (WKB) approximation [34] through the principal phase function S(r,t) while preserving the amplitude relationship with the systems’ quantum density:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M6','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M6">View MathML</a>

(6)

In this framework, the Schrödinger equation,

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M7','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M7">View MathML</a>

(7)

decomposes into real and imaginary parts. The real part can be expressed as follows:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M8','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M8">View MathML</a>

(8)

representing a continuous “fluid” of particles driven by the “guidance” momentum:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M9','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M9">View MathML</a>

(9)

moving under a joint external potential V(r) as well as under the so-called quantum potential influence:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M10','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M10">View MathML</a>

(10)

The consequences are nevertheless huge. For example, this methodology allows for the interpretation of the trajectories orthogonal to constant surfaces, by cancelling the Laplacian of the wave fronts ∇ r2S(r, t) = 0, which are obtained from eqs. (8) and (9) as the quantum equation of motion:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M11','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M11">View MathML</a>

(11)

Equation (11) resembles the classical Newtonian acceleration-force relationship only in a formal way; in fact, it generalizes it: it prescribes acceleration motion even in the absence of an external classical potential. This is essential in explaining why the inter-quark forces increase with the increase in inter-quark distances, no matter how great a separation is considered (a specific quantum effect), due to the presence of a quantum potential that does not fall off with distance as V does. It also nicely explains the observed interference patterns in double-slit experiments in the absence of classical forces. Alike, eq. (11) also appears suited for modeling chemical reactivity for the valence atoms as free particles in a virtually infinite potential environment to characterize their reactive behavior. In this regard, it is worth considering for such atoms the uniform motion by having ∂ p/∂ t = 0 through the time-constant associated wavefront condition and action S(r=cnst.,t)=cnst. (equivalent with Lagrangean constancy), in all given chemical space-points (atomic basins within molecule complex) [35]. This picture is also equivalently to have

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M12','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M12">View MathML</a>

(12)

applied to eq. (8). By doing so, one obtains

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M13','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M13">View MathML</a>

(13)

which can be rearranged as follows:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M14','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M14">View MathML</a>

(14)

such that the total energy of a the valence system is now entirely driven by the quantum potential:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M15','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M15">View MathML</a>

(15)

At this point, one can see that when turning to electronegativity and combining eq. (15) with DFT definition (3), one obtains a generalization of the previous Boeyens formulation [6]:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M16','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M16">View MathML</a>

(16)

which is the variation in the quantum potential with electron exchange under a constant classical or external potential.

However, for a quantum characterization of the valence state, we are interested in how the energy described by eq. (15) varies under a quantum potential (10)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M17','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M17">View MathML</a>

(17)

when the above relations (6) and (10) are substituted into eq. (15).

It is worth noting that although we obtained the total energy (17) in the Bohmian mechanics context, it showcases a clear electronic density dependency, not under a density functional (as DFT would require) but merely as a spatial function, which is a direct reflection of the entanglement behavior of Bohmian theory through the involvement of a quantum potential. However, in most cases, and especially for atomic systems, eq. (17) will yield numerical values under custom density function realizations.

Golden ratio imbalance for valence states of atoms

Atomic stability and periodicity remain major issues in the structural theories of matter; fortunately, they both have been largely solved by wave-particle (W/P) complementarily quantum behavior; phenomenologically, such relationship can be expressed as “WAVEPARTICLE = constant”, while it may be quantized (by Planck’s constant h) in the light of Heisenberg principle as [36]

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M18','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M18">View MathML</a>

(18a)

Remarkably, when fixing the particle’s observable property, say O, while letting wave information to vary, say ΔO, equation (18a) takes the workable form

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M19','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M19">View MathML</a>

(18b)

having as the preeminent realization the Bohr-de Broglie formulation a, leading with the first rationalization of the atomic periodicity [37]. However, when about the atomic chemical reactivity a similar analysis may be provided in terms of the number of electrons to atomic number ratio (N/Z): one may fix the observable (“particle”) character of the reactive atomic system by the ratio itself

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M20','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M20">View MathML</a>

(19a)

while modeling its evolving (“wave”) character by the natural variation of the previous ratio in terms of exchanged electrons respecting the neutral state:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M21','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M21">View MathML</a>

(19b)

When combining eqs. (19a) and (19b) into eq. (18b) on the lowest quantized state (nO=1), the “ground state” of atomic reactivity that is the atom in its valence state so to speak, and within atomic units’ formulation (i.e. by putting h=1, since the actual reactivity quantification involves only numbers with no dimension), one has the so called Heisenberg imbalance equation for valence atoms

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M22','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M22">View MathML</a>

(20a)

that can be rewritten as

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M23','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M23">View MathML</a>

(20b)

Eq. (20b) has the elementary acceptable solution

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M24','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M24">View MathML</a>

(21a)

which establishes, the direct “chemical” connection between the number of electrons and the atomic charge by means of the golden ratio

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M25','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M25">View MathML</a>

(21b)

generalizing the “physical” connection between nuclear (cosmic) synthesis at high pressure and atomic stability in the gas phase (Z=N); one has therefore the actual physical-to-chemical electronic charge – atomic number relationships

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M26','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M26">View MathML</a>

(22)

Worth remarking the results of type (20) and (22), here based on chemical reactivity specialization of Heisenberg type equations (18a) and/or (18b), were previously obtained at the level of neutron-protonic imbalance, inside the atomic nuclei, based on well-founded empirical observations [6]. The present golden ratio appearance is ultimately sustained also by the deviation from the N=Z condition for so-called “quark atoms” (as another way in considering the atoms in a quantum valence state), earlier identified as true matter’s entities responsible for matter’s reactivity at the atomic level [38].

Therefore the atomic structure branching (22) can be regarded as the present golden ratio extension to valence atom and as such employed; actually, its consequences regarding the characterization of the quantum valence states of atoms within the Bohmian quantum potential are the main aims of the present endeavor and will be discussed next.

Atomic implementation and discussion

On Slater density for valence atoms

Density is considered a “goldmine” in current computational and conceptual quantum chemistry due to its link with observable quantities, energy density functionals in particular, as celebrated by DFT [13,20,39,40]. However, to quantitatively approach the chemical phenomenology presented in Figure 1, involving the ionization-to-affinity atomic description, the general Slater [33] density (involving the orbital parameter ξ dependency) will be here employed for the first trial on modeling the combined Bohmian and gold-ratio features of valence atom; it assumes the general (trough still crude) working form:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M27','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M27">View MathML</a>

(23a)

For the reactivity at the valence atomic level, or for some outer shell (n) considered at the atomic frontier, one may assume almost electronic free motion or at least electronic motion under almost vanishing nuclear potential V(r); this way the density (23a), while entering the quantum potential (10) recovers the negative kinetic energy by the virial identity (14). Analytically, since eqs. (6), (10) and (23a), one has ∇ r2ρ1/2 = ξ2ρ1/2 and the actual valence atomic virial realization looks like

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M28','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M28">View MathML</a>

(24)

Equation (24) leaves with the identity:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M29','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M29">View MathML</a>

(25a)

that may be further rewritten with the help of the atomic Bohr-de Broglie relationship (see the note a) to provide the atomic frontier radii shell-dependency

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M30','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M30">View MathML</a>

(25b)

Remarkably, the same result is obtained when employing a far more reach atomic shell structure description, namely when starting with the full atomic radial Schrödinger density [25]

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M31','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M31">View MathML</a>

(26)

and imposing the null-gradient condition [41], ∇rρn(r, ξ) = 0, in accordance with the celebrated Bader condition of electronic flux of atoms-in-molecules [9,42], to yield:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M32','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M32">View MathML</a>

(25c)

The identity between eqs. (25b) and (25c) gives sufficient support to the present Slater density approach eq. (23a) in modeling the valence atoms or the atoms at their frontiers approaching reactivity (i.e. atoms-in-molecules complexes by chemical reactions).

Quantum chemical bonding and reactivity indices

Once convinced by the usefulness of the Slater density form (23a) for the present valence atomic analysis, one will next employ it under the so called Parr-Bartolotti form [33]

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M33','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M33">View MathML</a>

(23b)

such that to obey the N-normalization condition, as required by DFT [43-47],

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M34','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M34">View MathML</a>

(27)

by applying the Slater integral recipe

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M35','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M35">View MathML</a>

(28)

It nevertheless showcases the parametric ξ dependency that can be smeared out by considering the variational procedure

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M36','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M36">View MathML</a>

(29)

upon applying the total atomic energy

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M37','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M37">View MathML</a>

(30)

where the components are individually evaluated within a radial atomic framework with the respective results for [21,48]

● kinetic energy

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M38','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M38">View MathML</a>

(31a)

● nucleus-electronic interaction

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M39','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M39">View MathML</a>

(31b)

● inter-electronic interaction (see also Appendix)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M40','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M40">View MathML</a>

(31c)

With these results, the optimum atomic parameter is quantified by the electronic number as follows:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M41','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M41">View MathML</a>

(32)

which immediately releases the working electronic density

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M42','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M42">View MathML</a>

(33)

Having the completely analytical density in terms of number of reactive electrons as in eq. (33), worth pointing here on the so called sign problem relating with its variation, e.g., its gradient, the gradient of its square root, etc. Although this problem usually arises in density functional theory when specific energy functionals are considered in gradient forms, see for instance ref. [49], there is quite instructive discussing the present behavior and its consequences.

For instance, one can adapt either eqs. (25b) or (25c) through considering the present form (32) for the orbital exponent to be

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M43','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M43">View MathML</a>

(25d)

Here, one combines the frontier and maximum atomic radii with atoms-in-molecules phenomenology, as above indicated, to arrive to the present identification for the number of valence electrons possible to be involved in the same chemical bonding state as being Nbonding in (25d). Accordingly, the Figure 2 reveals interesting features of the present Slater-Parr-Bartolotti atomic density with quantum potential:

● the fact that the (covalent) bond length is proportional to the atomic radii and in inverse correlation with bonding order is well known [50], and this it is also nicely reflected in eq. (25d); however, changing the sign to negative radii as surpassing the threshold 21/5 and fixing in fact the limit Nbonding=4, is consistent with maximum bond order met in Nature; it is also not surprising this self-released limit connects with golden ratio by the golden-spiral optimization of bond-order [51]; more subtle, it connects also with the 4π symmetry of two spherical valence atoms making a chemical bond (Figure 2, inset): such “spinning” reminds of the graviton symmetry [52] (the highest spherical symmetry in Nature, with spin equal 2) and justifies the recent treatments of chemical bonding by means of the quasi-particles known as bondons [10,53], as well as the use of the 4D complex projective geometry in modeling the chemical space as a non-Euclidian one, eventually with a time-space metrics including specific “gravitational effects” describing the bonding [51];

● the “gap” between the atomic systems contributing 2 to 3 electrons to produce chemical bond is about double of the golden ratio, <a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M44','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M44">View MathML</a> ; therefore, this gap marks the passage from the space occupied by a pair of electrons and that required when the third electron is added on the same bonding state: it means that the third electron practically needs one golden measure (τ) to (covalently) share with each of the existing pairing electrons, while increasing the bond order to the level of three; it is therefore a space representation of the Pauli exclusion principle itself, an idea also earlier found in relation with dimensionless representation of a diatomic bonding energy (2τ) at its equilibrium bonding distance τ[54]; when the fourth electron is coming into the previous system, in order the maximum fourth order of bonding to be reach the chemical bonding space is inflating about five times more, yet forbidding further forced incoming electrons into the same space of bonding state as the bonding radius becomes negative in sign.

thumbnailFigure 2. Representation of the bonding length as a function of bonding electrons from valence atoms in molecule(s), based on eq. (25d), while marking the double golden ratio 2τ gap between the bonding lengths of the second and third bonding order, as well as the forbidden chemical bonding region for Nbonding ≥21/5 for the electrons participating in the same bonding. Further connection of chemical bonding and the 4D space to model it is suggested by the inset picture illustrating the 2-fold (4π) spinning symmetry of the adduct atom respecting the bonding direction, after [51].

Having revealed the chemical bonding information carried by the density (33) when considered for combined valence atoms-in-molecules, it is next employed on energetically describing the atomic reactivity as a propensity for allowing electronic exchanging and bonding. As such, it leaves the total quantum (Bohmian) energy in (17) with the compact form

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M45','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M45">View MathML</a>

(34a)

Note that the actual working total energy is not that obtained by replacing the density (33) in eqs. (31a)-(31c) and then in total energy (30) because here the double-variational procedure was considered; that is, the first optimization condition was considered as in eq. (29), and the resulting (optimum) density (33) was then employed in the quantum energy (17), which in turn was obtained by applying the variational eq. (12) to the perceived phase transition in the Bohm eikonal wave-function (5). To emphasize the accuracy of eq. (17) over that of (30) with density (33), when one considers the last case, eq. (30) yields the following non-quadratic form for energy:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M46','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M46">View MathML</a>

(34b)

which is not appropriate for describing the valence state of an atom, as eq. (2) prescribes, despite being similar in form to the Bohmian-based result of eq. (34a). Thus, the previous limitation of the Parr-Bartolotti conclusion [33] and the paradox raised in describing the valence (parabolically) state with the optimized atomic density (33) are here solved by the double (or the orthogonal) variational implementation, as recently proved to be customary for chemical spaces [32]. In the light of this remark one may explain also the sign difference between the “physical” energy (34b) and that obtained for the “chemical” situation (34a): through simple variational procedure for “physical” energy (30) the result (34b) is inherently negative – modeling systems stability in agreement with the upper branch of eq. (22), whereas the double variational algorithm employing optimized density (33) into the Bohmian shaped energy (17) it produces the positive output (34a) associated with activation energy characteristic for chemical reactivity corresponding to the lower branch of eq. (22).

Therefore, to be accurate, one should consider the quantum potential related optimized energy (34a) instead of simply the orbital optimized one of eq. (34b). Therefore, assuming that eq. (34a) appropriately describes the atomic valence state in DFT (see the upper/reactive curve in Figure 1b), the next task is to search for the quantum valence charge for which the valence energy approaches its optimum value (or the “ground state” of the atomic chemical-reactivity, i.e. the previously golden-ratio quantification of the valence atomic state); to this aim, at this point, one can employ the golden ratio relationship (21a) and first rewrite eq. (34a) as

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M47','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M47">View MathML</a>

(35)

which is minimized at the value

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M48','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M48">View MathML</a>

(36)

However, one must again apply the double-variational procedure, now in terms of number of electrons, i.e., reconsidering eq. (36) with the golden ratio at the reactive (chemical) electronic level of eq. (22) such that a second equation is formed

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M49','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M49">View MathML</a>

(37)

with the positive solution

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M50','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M50">View MathML</a>

(38)

This expression avails of the significance of the maximum number of electrons, for a given atom, possibly engaged in a reactive environment by either (or both) accepting or (and) ceding electrons to or from its valence state, see Table 1.

Table 1. Synopsis of the critical charges in the physical ground state (N*) as well as for chemical reactive (valence) state (N**) for atoms of the first four periods of the periodic table of elements, as computed from the minimum point of associated interpolations of ionization and electronic affinities[33]and of eq. (38), respectively

The result of this process is different from the expected physical result (NSTABIL=Z) according to the upper branch of eq. (22), which is higher than the physical one until reaching the carbon system (ZINTRCHANGE=6.8), while continuing below it thereafter (see Figure 3).

thumbnailFigure 3. The comparative shapes of the valence electrons to be engaged in chemical reactivity (continuous curve) computed using eq. (38) based on the combined optimal Bohm total energy (35) with the golden ratio imbalance of eq. (22), respecting the stable physical case (dot-dashed curve), and of their differences (dashed curve); all originate at the 0thatom (the neutron, Z=0).

The above interchange (effective) atomic number through which the chemical (reactive) state is associated with lower charge respecting the physical state may be also be found at the energetic level based on quantum equation (34a), as specialized for the two branches of Figure 3 for the N(Z) dependence. Thus, the chemical (reactive) state takes the analytic form

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M51','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M51">View MathML</a>

(39)

and interchanges with the ground state EQ1(NSTABLE → Z) at the points {3.5,6.8}, as observed also from Figure 4; however, the interchanging point beyond which all chemical atomic systems are more stable in the chemical or reactive state than in the physical ground state is consistently recovered.

thumbnailFigure 4. The same comparative shapes shown in Figure3, here at the level of energy (34a) specialized for the reactive and the stable N(Z) dependencies of Figure3; the various plots successively display increasingly large atomic Z-ranges to better emphasize the chemical vs. physical behavior (see text).

Nevertheless, the energetic analysis also reveals the atomic systems Be, B and C to be situated over the corresponding physical stable states; this may explain why boron and carbon present special chemical phenomenology (e.g., triple electronic bonds and nanosystems with long C-bindings, respectively), which is not entirely explained by ordinary physical atomic paradigms [55-60].

The energetic discourse may be complete with the electronegativity and chemical hardness evaluations by applying the DFT definitions (3) and (4) to physical and chemical energies, respectively. In the first case, expression (34b) is applied to provide the following so-called “universal” forms of Parr and Bartolotti [33]:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M52','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M52">View MathML</a>

(40)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M53','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M53">View MathML</a>

(41)

The result, nevertheless, appears to be an unusually higher increase in chemical hardness than in electronegativity, which certainly cannot be used to model a reactive-engaged tendency because it is more stable (by chemical hardness) than reactive (by electronegativity); it is, however, consistent with the physical stability of the system, provided by the single variational procedure through which eq. (34b) was produced.

Instead, to chemically model reactivity, the double variation procedure is applied and eq. (34a) is substituted into eqs. (3) and (4), though by considering also the double reactive procedure for charge as well, i.e., by considering eq. (38) with the golden ratio information of (22) to respectively yield the results

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M54','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M54">View MathML</a>

(42)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M55','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M55">View MathML</a>

(43)

Remarkably, the actual electronegativity of (42) obtained by the quantum Bohm and golden ratio double procedure yields sensible results similar to those of the single variational approach (40); however, the chemical hardness of (43) is approximately 5-fold lower than its “stable” counterpart (41), affirming therefore the manifestly reactive framework it produces – one described by a quadratic equation (34a) instead of a cubic one (34b).

Charge waves in gauge chemical reactivity

Finally, one considers the chemical reactivity discussion as based on the gauge reaction that equilibrates the chemical bond by symmetrical bond polarities [25]

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M56','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M56">View MathML</a>

(44a)

such that the reactive electrons are varied on the reunited intervals of eq. (1); such analysis was previously employed to fundament systematic electronegativity and chemical hardness definitions by the averaging (through the integration) factor

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M57','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M57">View MathML</a>

(44b)

along the reaction path accounting for the acidic (electron accepting, 0 ≤ N ≤ +1) and basic (electron donating, –1 ≤ N ≤ 0) chemical behaviors.

In this scaled (gauge) context of reactivity, the foregoing discussion is dedicated to investigating the link between the critical ground state charge (N*) and the valence or reactive state (N**). While the first appears as a consequence of naturally fitting the three points in Figure 1 (the ionization, neutral and affinity states), with the effect of biasing the minimum of the energetic curve in Figure 1b with respect to the apparent Parr-DFT curve in Figure 1a, and is thus derived graphically (see Figure 5), the valence charge is based on the combined quantum energy and golden ratio information in eq. (38). Both are reported for the indicated number of atomic systems of the periodic table of elements in Table 1. One notes, for instance, that while the critical ground state charge N* always lies in the range [0.5,1], the valence charge N** may span the interval [0,1]; one may interpret such behavior as being associated with the difference between the fraction ½ and integer “1” in driving the principles of chemical reactivity and the electrophilicity equalization principle in particular, when the “quantum transition” 1/2→1 is required in the energy exchange of chemical systems for it to be valid for both electronegativity and chemical reactivity principles [61]; nevertheless such scaling it is equivalent with above acidic-basic gauge averaging of eq. (44b). This way, the valence charge problem may be extended to the interval [0,2], at its turn seen as a gauge transformation of the chemical reactivity charge domain [−1, +1], where one reencounters the challenging problem of whether the “One electron is less than half what an electron pair is” [62], the response to which is generally complex but may here be approached through the following steps.

thumbnailFigure 5. A graphical interpolation for selected elements of Table1in terms of their ionization, neutral and affinity states, aiming to determine the critical (displaced) charge of the DFT ground state, as prescribed by Figure1b.

First, by employing the data presented in Table 1, one constructs the so-called “continuous” ground and valence charge states by appropriately fitting over the first four periods of elements, here restrained to 10th-order polynomials. This is performed by interpolating every three points of the 32 elements presented in Table 1, although by spanning the atomic number range Z ∈ [1, 53], thus yielding (see also the allied representations of Figure 6):

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M58','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M58">View MathML</a>

(45a)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M59','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M59">View MathML</a>

(45b)

thumbnailFigure 6. The critical ground state and valence charge points for the elements of Table1and their 10th-order continuous interpolations according to eqs. (45a) and (45b).

Equations (45a) and (45b) are then combined into a sort of special charge wave function based on their difference on the golden ratio scale (see Figure 6 for graphical representation)

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M60','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M60">View MathML</a>

(46)

with the peculiar property that its square-integrated form over the Z-range of interpolation gives

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M61','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M61">View MathML</a>

(47)

The result (47) has the following conceptual fundamental quantitative interpretation: the difference between the ground and valence optimum charges is regulated by the golden ratio scale, or in other terms,

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M62','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M62">View MathML</a>

(48)

such that it provides a sort of normalization corrected by the golden ratio value; it also fulfills the interesting relationship:

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M63','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M63">View MathML</a>

(49)

In any case, the present analysis provides the qualitative result that the difference between the critical ground state and optimal valence charges is more than half of an electronic pair, giving rise to the significant notion that chemical reactivity is not necessarily governed by a pair of electrons but governed by no less than half of a pair and is related to the golden ratio (τ > 0.5).

However, fractional values in general and those related to the golden ratio particular, may be interpreted as a consistent manifestation of the quantum mechanical (i.e., wave functional) approach of chemical phenomena, here at the reactivity level. Moreover, the quadratic critical charge function (46), as shown in Figure 7, clearly reveals that a higher contribution to electronic pair chemistry is given by the third period of elements and by the third and fourth transitional elements in particular, a result that nicely agrees with the geometrical interpretation of the chemical bond, particularly the crystal ligand field paradigm of inorganic chemistry [9].

thumbnailFigure 7. The linear and quadratic charge “wave function” of eq. (46).

Also a local analysis of the type of charge that is dominant in atomic stability, i.e., the critical physical ground state or the chemical valence reactive state based on eqs. (45a) and (45b), respectively, may be of considerable utility in refined inorganic chemistry structure-reactivity analysis. To the same extent, it depends on the degree of the polynomials used to interpolate the critical and valence charges over the concerned systems; however, through the present endeavor, we may assert that the analysis should be of the type (48), which in turn remains a sort of integral version of the imbalance equation (20a), in this case for the ground-valence charge gap states of a chemical system.

Conclusions

Aiming to hint at the solution to the current debate regarding the physical vs. chemical definition of an atom and as a special stage of a larger project regarding quantum chemical orthogonal spaces, the present work addresses the challenging problem of defining and characterizing valence states with respect to the ground state within conceptual density functional theory. We are aware of the earlier warnings raised by Parr and Bartolotti and others [18,33,63] regarding the limits of density functional theory and of the total energy of atomic systems combined with a Slater-based working density to provide a quadratic form in terms of system charge, as required by the general theory of chemical reactivity of atoms and molecules in terms of electronegativity and chemical hardness. Fortunately, we discovered that the Bohmian form of the total energy of such atomic systems provides, instead, the correct behavior, although it is only density-function-dependent and not a functional expression. Moreover, this finding was reached through the so-called double variational procedure, which, as emphasized earlier, was likely to reproduce the chemical reactivity principles of electronegativity and chemical hardness in an analytical manner; however, such a double analytical variational approach is consistent with the recent advanced chemical orthogonal spaces approaches of chemical phenomenology [64] as being at least complementary to the physical description of many-electronic systems when they are engaging in reactivity or equilibrium as the atoms-in-molecules Bader theory prescribes [9,42]. With the present Bohmian approach, the total energy is in fact identified with the quantum potential, thus inherently possessing non-locality and appropriate reactivity features, which are manifested even over long distances [10,11,53]; this also generalizes the previous Boeyens electronegativity formulation of electronegativity [5,6] from the direct relationship between a quantum potential and its charge derivative. The double algorithm was also implemented to discriminate the valence from the ground state charges, this time by using the golden ratio imbalance equation as provided by adaption of the Heisenberg type relationship to chemical reactivity for atoms. This corresponds to an analytical unfolding of the physical and chemical imbalance of the electronic charge stability of atomic systems, paralleling the deviation from the equal electron-to-proton occupancy in physical systems toward electron deficiency in the valence states of chemical systems. This dichotomy was implemented by the golden ratio presented in eq. (22). As a consequence, the difference between valence and ground state charge systems is naturally revealed and allows for the explanation of chemical reactivity and bonding in terms of fractional electron pairs, althrough driven by the golden ratio under the so-called physical-to-chemical charge difference wave function and associated normalizations, all of which represent elaborated or integral forms of the basic imbalance atomic equation. The present results are based on 10th-order polynomial fitted over 32 elements from the first 54 elements of the first four periods of periodic table of elements and can be further pursued by performing such systematic interpolations that preserve the golden ratio relationships, as advanced herein; they may also provide a comprehensive picture of how valence electrons may always be projected/equalized/transposed into ground state electrons within the perspective of further modeling chemical reactions when chemical reactivity negotiates the physical molecular stabilization of atoms in molecules.

Endnotes

a For circular orbits, the lowest ones in each atomic shells – including the valence ones, one has ΔOr=2πr, with r the orbital radii thereof, while O=p is the fixed particle’s momentum on that orbit; therefore, when combined into eq. (18b) they provide the celebrated Bohr-de Broglie relationship rp=nħ solving the atomic spectra of Hydrogen atom in principal quantum numbers (n).

Appendix: Semi-classical inter-electronic energy

For the inter-electronic interaction, see Figure 8; in evaluating Vee[ξ] of eq. (31c), the two-electronic density is approximated by the Coulombic two mono-electronic density product, thus neglecting the second-order density matrix effects associated with the exchange-correlation density.

However, for the analytical evaluation of the electron–electron repulsion energy using the density (23b), much care must be taken. For instance, one has to use the electrostatic Gauss theorem, which states that the classical electrostatic potential outside a uniform spherical shell of charge is just what it would be if that charge were localized at the center of the shell and that the potential everywhere inside such a shell is that at the surface, [21,48] see Figure 8. Therefore, the electronic repulsion energy becomes

<a onClick="popup('http://journal.chemistrycentral.com/content/6/1/135/mathml/M64','MathML',630,470);return false;" target="_blank" href="http://journal.chemistrycentral.com/content/6/1/135/mathml/M64">View MathML</a>

(A1)

which recovers the expression presented by eq. (31c), when the Slater integral type of Eq. (28) is also employed. Note that the electron–electron repulsion term was written by also considering the Fermi-Amaldi (N-1)/N factor [13], which ensures the correct self-interaction behavior: when only one electron is considered, the self-interaction energy must be zero, Vee (N→1)→0.

thumbnailFigure 8. Representation of the space regions of the 1st and 2nd electrons, their potential influences and reciprocal interaction [21,48]

Competing interests

The author declares he has no competing interests.

Acknowledgements

This work was supported by CNCS-UEFISCDI agency through the research project TE16/2010-2013 within the PN II-RU-TE-2009-1 framework. Inspiring discussions with Profs. Boeyens (University of Pretoria, South Africa) and von Szentpaly (Stuttgart University, Germany) are kindly acknowledged. Constructive referees’ comments and stimulus is also sincerely thanked. This paper is dedicated to Prof. Robert G. Parr for his pioneering quantum work on atomic valence states.

References

  1. Putz MV: Big chemical ideas in the context: the periodic law and the Scerri’s Periodic Table.

    Int J Chem Model 2011, 3:15-22. OpenURL

  2. Scerri ER: The Periodic Table – Its Story and Its Significance. Oxford-New York: Oxford University Press; 2007. PubMed Abstract OpenURL

  3. Boeyens JCA: Emergent Properties in Bohmian Chemistry. In Quantum Frontiers of Atoms and Molecules. Edited by Putz MV. New York: Nova Publishers Inc; 2011:191. OpenURL

  4. Boeyens JCA: New Theories for Chemistry. Amsterdam: Elsevier; 2005. OpenURL

  5. Boeyens JCA: Chemistry from First Principles. Heidelberg-Berlin: Springer; 2008. OpenURL

  6. Boeyens JCA, Levendis DC: Number Theory and the Periodicity of Matter. Heidelberg-Berlin: Springer; 2008. OpenURL

  7. Kaplan IG: Is the Pauli exclusive principle an independent quantum mechanical postulate?

    Int J Quantum Chem 2002, 89:268-276. Publisher Full Text OpenURL

  8. Scerri ER: Just how ab initio is ab initio quantum chemistry?

    Found Chem 2004, 6:93-116. OpenURL

  9. Bader RFW: Atoms in Molecules - A Quantum Theory. Oxford: Oxford University Press; 1990. OpenURL

  10. Putz MV: The bondons: The quantum particles of the chemical bond.

    Int J Mol Sci 2010, 11:4227-4256. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  11. Putz MV: Quantum Theory: Density, Condensation, and Bonding. Toronto: Apple Academics & CRC Press; 2012. OpenURL

  12. Parr RG, Donnelly RA, Levy M, Palke WE: Electronegativity: the density functional viewpoint.

    J Chem Phys 1978, 68:3801-3808. Publisher Full Text OpenURL

  13. Parr RG, Yang W: Density Functional Theory of Atoms and Molecules. New York: Oxford University Press; 1989. OpenURL

  14. Bohm D: A suggested interpretation of the quantum theory in terms of “hidden” variables. I.

    Phys Rev 1952, 85:166-179. Publisher Full Text OpenURL

  15. Bohm D: A suggested interpretation of the quantum theory in terms of “hidden” variables. II.

    Phys Rev 1952, 85:180-193. Publisher Full Text OpenURL

  16. Bohm D, Vigier JP: Model of the causal interpretation of quantum theory in terms of a fluid with irregular fluctuations.

    Phys Rev 1954, 96:208-216. Publisher Full Text OpenURL

  17. Cushing JT: Quantum Mechanics – Historical Contingency and the Copenhagen Hegemony. Chicago & London: The University of Chicago Press; 1994. PubMed Abstract | Publisher Full Text OpenURL

  18. Von Szentpály L: Modeling the charge dependence of total energy and its relevance to electrophilicity.

    Int J Quant Chem 2000, 76:222-234. Publisher Full Text OpenURL

  19. 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. Publisher Full Text OpenURL

  20. Geerlings P, De Proft F, Langenaeker W: Conceptual density functional theory.

    Chem Rev 2003, 103:1793-1874. PubMed Abstract | Publisher Full Text OpenURL

  21. Putz MV: Contributions within Density Functional Theory with Applications in Chemical Reactivity Theory and Electronegativity. Parkland: Dissertation.com; 2003. OpenURL

  22. Parr RG: Density functional theory.

    Annu Rev Phys Chem 1983, 34:631-656. Publisher Full Text OpenURL

  23. Parr RG, Pearson RG: Absolute hardness: companion parameter to absolute electronegativity.

    J Am Chem Soc 1983, 105:7512-7516. Publisher Full Text OpenURL

  24. Putz MV: Absolute and Chemical Electronegativity and Hardness. New York: Nova Publishers Inc.; 2008. OpenURL

  25. Putz MV: Systematic formulation for electronegativity and hardness and their atomic scales within density functional softness theory.

    Int J Quantum Chem 2006, 106:361-389. Publisher Full Text OpenURL

  26. Chattaraj PK, Parr RG: Density functional theory of chemical hardness.

    Struct Bond 1993, 80:11-25. Publisher Full Text OpenURL

  27. Chattaraj PK, Sengupta S: Popular electronic structure principles in a dynamical context.

    J Phys Chem 1996, 100:16129-16130. OpenURL

  28. Chattaraj PK, Maiti B: HSAB principle applied to the time evolution of chemical reactions.

    J Am Chem Soc 2003, 125:2705-2710. PubMed Abstract | Publisher Full Text OpenURL

  29. Chattaraj PK, Duley S: Electron affinity, electronegativity, and electrophilicity of atoms and ions.

    J Chem Eng Data 2010, 55:1882-1886. Publisher Full Text OpenURL

  30. Ayers PW, Parr RG: Variational principles for describing chemical reactions: reactivity indices based on the external potential.

    J Am Chem Soc 2001, 123:2007-2017. PubMed Abstract | Publisher Full Text OpenURL

  31. Kohn W, Becke AD, Parr RG: Density functional theory of electronic structure.

    J Phys Chem 1996, 100:12974-12980. Publisher Full Text OpenURL

  32. Putz MV: Chemical action concept and principle.

    MATCH Commun Math Comput Chem 2011, 66:35-63. OpenURL

  33. Parr RG, Bartolotti LJ: On the geometric mean principle of electronegativity equalization.

    J Am Chem Soc 1982, 104:3801-3803. Publisher Full Text OpenURL

  34. Kleinert H: Path Integrals in Quantum Mechanics, Statistics, Polymer Physics, and Financial Markets. 3rd edition. Singapore: World Scientific; 2004. OpenURL

  35. Guantes R, Sanz AS, Margalef-Roig J, Miret-Artés S: Atom–surface diffraction: a trajectory description.

    Surf Sci Rep 2004, 53:199-330. Publisher Full Text OpenURL

  36. Putz MV: On Heisenberg uncertainty relationship, its extension, and the quantum issue of wave-particle duality.

    Int J Mol Sci 2010, 11:4124-4139. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Pauling L, Wilson EB: Introduction to Quantum Mechanics with Applications to Chemistry. New York: Dover Publications; 1985. OpenURL

  38. Lackner KS, Zweig G: Introduction to the chemistry of fractionally charged atoms: electronegativity.

    Phys Rev D 1983, 28:1671-1691. Publisher Full Text OpenURL

  39. Hohenberg P, Kohn W: Inhomogeneous electron gas.

    Phys Rev 1964, 136:B864-B871. Publisher Full Text OpenURL

  40. Putz MV: Density functionals of chemical bonding.

    Int J Mol Sci 2008, 9:1050-1095. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  41. Ghosh DC, Biswas R: Theoretical calculation of absolute radii of atoms and ions. Part 1. The atomic radii.

    Int J Mol Sci 2002, 3:87-113. Publisher Full Text OpenURL

  42. Bader RFW: Definition of molecular structure: by choice or by appeal to observation?

    J Phys Chem A 2010, 114:7431-7444. PubMed Abstract | Publisher Full Text OpenURL

  43. Dreizler RM, Gross EKU: Density Functional Theory. Heidelberg: Springer Verlag; 1990. OpenURL

  44. Kryachko ES, Ludena EV: Energy Density Functional Theory of Many Electron Systems. Dordrecht: Kluwer Academic Publishers; 1990. OpenURL

  45. Cramer CJ: Essentials of Computational Chemistry. Chichester: Wiley; 2002. OpenURL

  46. Capelle K: A bird's-eye view of density-functional theory.

    Braz J Phys 2006, 36:1318-1343. Publisher Full Text OpenURL

  47. Jensen F: Introduction to Computational Chemistry. Chichester: John Wiley & Sons; 2007. OpenURL

  48. Parr RG: The Quantum Theory of Molecular Electronic Structure. Reading-Massachusetts: WA Benjamin, Inc.; 1972. OpenURL

  49. Cohen AJ, Mori-Sánchez P, Yang W: Challenges for density functional theory.

    Chem Rev 2012, 112:289-320. PubMed Abstract | Publisher Full Text OpenURL

  50. Petrucci RH, Harwood WS, Herring FG, Madura JD: General Chemistry: Principles & Modern Applications. 9th edition. New Jersey: Pearson Education, Inc.; 2007. PubMed Abstract | Publisher Full Text OpenURL

  51. Boeyens JC, Levendis DC: The structure lacuna.

    Int J Mol Sci 2012, 13:9081-9096. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Hawking S: The Universe in a Nutshell. New York: Bantam Books; 2001. OpenURL

  53. Putz MV, Ori O: Bondonic characterization of extended nanosystems: application to graphene's nanoribbons.

    Chem Phys Lett 2012, 548:95-100. OpenURL

  54. Boeyens JCA: A molecular–structure hypothesis.

    Int J Mol Sci 2010, 11:4267-4284. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  55. March NH: Electron Density Theory of Many-Electron Systems. New York: Academic; 1991. OpenURL

  56. Wentorf RH Jr: Boron: another form.

    Science 1965, 147:49-50. PubMed Abstract | Publisher Full Text OpenURL

  57. Eremets MI, Struzhkin VV, Mao H, Hemley RJ: Superconductivity in Boron.

    Science 2001, 293:272-274. PubMed Abstract | Publisher Full Text OpenURL

  58. van Setten MJ, Uijttewaal MA, de Wijs GA, de Groot RA: Thermodynamic stability of boron: the role of defects and zero point motion.

    J Am Chem Soc 2007, 129:2458-2465. PubMed Abstract | Publisher Full Text OpenURL

  59. Widom M, Mihalkovic M: Symmetry-broken crystal structure of elemental boron at low temperature.

    Phys Rev B 2008, 77:064113. OpenURL

  60. Putz MV (Ed): Carbon Bonding and Structures: Advances in Physics and Chemistry. Dordrecht-London: Springer Verlag; 2011.

    Cataldo F, Milani P (Series Editors): Carbon Materials: Chemistry and Physics, Vol. 5.

    OpenURL

  61. Putz MV, Mingos DMP (Eds): Applications of Density Functional Theory to Chemical Reactivity. Berlin-Heidelberg: Springer Verlag; 2012. [Struct Bond] OpenURL

  62. Ferreira R: Is one Electron less than half what an electron pair is?

    J Chem Phys 1968, 49:2456-2457. Publisher Full Text OpenURL

  63. Bergmann D, Hinze J: Electronegativity and charge distribution.

    Struct Bond 1987, 66:145-190. Publisher Full Text OpenURL

  64. Putz MV: Chemical Orthogonal Spaces. Kragujevac: Kragujevac University Press; 2012.

    [Gutman I (Series Editor): Mathematical Chemistry Monographs, Vol. 14].

    OpenURL