Equidensity orbitals in resultant-information description of electronic states

The equidensity orbitals (EO) are used in the resultant entropic description of molecular states which combines the probability and current contributions in the overall information content. Continuities of the modulus and phase components of electronic wavefunctions are examined, and the Harriman–Zumbach–Maschke (HZM) construction of Slater determinants yielding the prescribed electron density is explored. The conditional probability interpretation of (complex) HZM wavefunctions is formulated, the entropy/information contributions due to the state phase component are summarized, and a nonclassical origin of quantum dynamics of the resultant gradient information, related to average kinetic energy of electrons, is emphasized. The phase equilibria maximizing the resultant-entropy measures are explored, and “thermodynamic” phase minimizing the overall gradient information is determined. It generates finite orbital currents giving rise to the vanishing resultant flow of electrons in the system as a whole. Potential use of atomic and molecular EO bases in electronic structure calculations and interpretations in chemistry is discussed, and illustrative example of Gaussian probability distribution is examined in some detail.


, in the resultant-entropy/information measures S[p, ϕ] = S[p] + S[ϕ] or I[p, ϕ] = I[p] + I[ϕ].
The extrema of such overall global or gradient entropy descriptors determine the phase equilibria in molecules and their fragments [15][16][17][18][19][20]. These generalized IT descriptors allow one to distinguish the information content of states exhibiting the same electron density but differing in their current composition, e.g., the bonded (entangled) and nonbonded (disentangled) states of subsystems in the chemical reactivity theory [21,22]. The entropic phase contributions play an important role in an information description of The following notation is adopted: A denotes a scalar, A is the row or column vector, A represents a square or rectangular matrix, and the dashed symbol Â stands for the quantum-mechanical operator of the physical property A. The logarithm of Shannon's information measure is taken to an arbitrary but fixed base: log = log 2 corresponds to the information content measured in bits (binary digits), while log = ln expresses the amount of information in nats (natural units): 1 nat = 1.44 bits. reactive systems [23][24][25][26][27][28] by precisely specifying the mutually open (bonded) and closed (nonbonded) states of reactants. To paraphrase Prigogine [29], the density distribution establishes the molecular "static" structure of being, while the current pattern establishes the complementary "dynamical" aspect of molecular organization-the structure of becoming: the latter determines the dynamics of the former.
In present analysis, the phase/current feature of electronic states will be emphasized throughout. General quantum states of N electrons yielding the prescribed electron density involve Slater determinants constructed from the (orthonormal) (plane-wave)-type equidensity orbitals (EO) of Harriman, Zumbach and Maschke (HZM) [73,74], each generating the specified molecular probability distribution and exhibiting the density-dependent local phase safeguarding the orbital orthogonality. This construction separates the modulus and phase aspects of multi-electronic states thus offering an attractive framework for studying the role of classical (probability) and nonclassical (current) entropy/ information contributions in chemical processes. We shall examine the construction of orbital phases and electronic currents they generate, providing them an appropriate probabilistic interpretation in terms of conditional probabilities involving electron coordinates. This development will be illustrated for the special case of the Gaussian-type probability distributions.
The continuities of the physical, probability and phase descriptors of electronic states implied by the Schrödinger equation will be summarized, and probability interpretation of the key constructs of this orbital representation will be given. The dynamics of resultant gradient information will be explored, and a nonclassical character of the information production will be emphasized. The overlap extension of the equidensity orbitals [75] will be outlined, and possible applications of EO bases of constituent atoms in electronic structure calculations and interpretations in chemistry will be explored.

Continuities of probability and phase distributions
The evolution of molecular electronic wavefunctions is determined by the Schrödinger equation (SE) of quantum mechanics. It also determines the dynamics of probabilityand current attributes of the abstract (complex) quantum states [26][27][28] and the associated information contributions. It is of interest to examine general implications of SE for the dynamics of these two physical descriptors of the modulus and phase components of molecular wavefunctions. For simplicity, let us first consider a single electron at time t in state � (t)⟩ ≡ � ⟩ , described by the (complex) wavefunction in position representation, where the real functions R(r, t) and ϕ(r, t) stand for its modulus and phase parts, respectively. They determine the electron probability distribution at time t, and the corresponding current density Here, the effective velocity of the probability "fluid," V(r, t) = dr/dt, measures the local current per particle and reflects the state phase-gradient: The state probability density p(r, t) and its flux j(r, t) thus constitute two independent physical degrees-of-freedom of such a general quantum state of an electron: The phase component of molecular states reflects the coherence, current aspect of the system electronic structure.
In molecular scenario, the electron is moving in the external potential v(r) due to the "frozen" nuclei of the familiar Born-Oppenheimer (BO) approximation, thus being described by the Hamiltonian where T (r) denotes its kinetic part. The dynamics of its quantum state ψ(r, t) is determined by SE, (4) V(r, t) = j(r, t)∕p(r, t) = (ℏ∕m)∇ (r, t).
which also generates the associated temporal evolutions of the state probability density and of its phase distribution.
The SE implies the sourceless continuity relation for the particle probability "fluid": The total time derivative dp(r, t)/dt of the electron probability density p(r, t) = p[r(t), t] defines its local production σ p (r, t). It measures the time rate of change in an infinitesimal volume element around r(t) of the probability fluid moving with the local velocity dr/dt = V(r, t), while the partial derivative ∂p(r, t)/∂t refers to the volume element around the fixed point in space. The vanishing total derivative dp/dt of electron probability density thus reads: Using the probability continuity also gives: Therefore, comparing the two preceding equations implies the vanishing divergence of the effective velocity field, determined by the Laplacian of the phase component of molecular quantum states: The probability continuity of Eq. (9) also determines the state modulus dynamics: while the phase dynamics from SE reads: The probability velocity V also drives the phase current The SE then predicts a nonvanishing phase source in the associated continuity equation: Using Eq. (13) finally gives the following production term of a local phase: It follows from these continuity relations that a common velocity of the probability and phase currents, determined by the state phase-gradient, reflects the structure dynamical aspect, while the probability density at the specified time establishes its static facet. The former characterizes the temporal evolution of the latter. In other words, a structure of "becoming" determines the evolution of an instantaneous structure of "being."

Equidensity orbitals
In DFT, one often explores the density-constrained ("entropic") variational principles, e.g., the Levy [7] construction of the universal density functional for the sum of electronic kinetic and repulsion energies. They are also called the "vertical" or "thermodynamic" searches, by analogy to the minimum-energy and maximum-entropy criteria for determining equilibria in ordinary thermodynamics [76]. A related problem of constructing the antisymmetric wave functions of N fermions yielding the prescribed density ρ(r), vital for solving the familiar N-representability problem of DFT, has been solved by Harriman [73], who had used the crucial insights due to Macke [77,78] and Gilbert [79]. Its three-dimensional generalization by Zumbach and Maschke [74] introduces the complete set of Slater determinants build using the density-dependent (orthonormal) equidensity orbitals (EO) of the plane-wave type, each generating the molecular probability distribution the shape factor of the system electron density ρ(r). Here, denotes the (constant) reduced-momentum (wave-number) vector, and Φ q (r) stands for the EO resultant phase, defined by the sum of the orthogonality-phase F q (r) and thermodynamic-phase contribution ϕ(r), common to all occupied EO. The latter has to be eventually determined from the auxiliary maximum resultant-entropy principle, which establishes the system phase equilibrium. The vector function in Eq.
This construction can be thus viewed as transformation of the physical space into itself: for any r = (x, y, z) ∈ ℜ 3 we define the (invertible) density-dependent transformation of r into the new vector f(r) ∈ ℜ 3 , with the familiar plane-wave type orthogonality relation between orbitals. The overall phase Φ q (r) of EO involves the orbital-specific (geometric) internal contribution q·f(r) ≡ F q (r), which enforces the independence of these one-particle states, and "thermodynamic" external phase (r) , common to all EO, The latter has to be determined from a subsidiary variational principle for the extremum of the state resultant entropy [12][13][14][15][16][17][18][19][20][21]36]. By definition, in the HZM construction one sets this transformation between local volume elements as being proportional to molecular probability distribution, df (r) ∝ p(r)dr , with the associated Jacobian determinant These "diagonal" derivatives { f ∕ x } can be expressed in terms of the relevant conditional probabilities [14] involving the position coordinates (see Figs. 1 and 2), defined by the corresponding ratios of the relevant joint probabilities: Here, p(a|b) ≡ p(a, b)/p(b) denotes the conditional probability of event a, given event b, while p(a, b) stands for the joint probability of simultaneously observing the two events. Substituting Eq. (25) into Eq. (24) then indeed generates the postulated density-dependent transformation between the volume elements dr and df: The mutual orthogonality between such EO states then directly follows: In "Appendix 1," the illustrative case of Gaussian probability distribution is examined, where the key concepts of HZM construction have been examined, with the local phase vector f(r) expressed in terms of the error function of probability theory.
To summarize, the Slater determinants build from specific selections of N different equidensity orbitals, then by construction generate the prescribed electron density (r): .
They constitute the complete orthonormal set of N-particle functions capable of representing the molecular ground state of N electrons for the associated electron distribution , the equilibrium density for the system electronic Hamiltonian in the HZM configuration interaction (CI) type expansion: It should be emphasized, however, that such trial functions do not satisfy the phase constraint of Eq. (11), the implication of the probability-continuity relation, since ∇·f(r) ≠ 0. This phase requirement should be ultimately fulfilled by the exact solutions of SE, which reflect the system dynamical conditions.
The EO { q (r) = ⟨r� q ⟩ ≡ ⟨r�q⟩} constitute the complete set of functions capable of representing any quantum state � ⟩ of a single electron or the associated wavefunction in the position representation: A single circle delineates probability of a separate coordinate, e.g., p(x) = ∫ p(r)dydz , the envelope of the two or three overlapping circles represents the probability of the union (sum) of observing the specified (dependent) coordinates, e.g., The probabilities of the joint coordinate events (products), e.g., p(x,y) or p(x,y,z) similarly correspond to the overlap areas (intersections) of two or three circles, respectively. The non-overlapping circles in part (a) represent the independent, exclusive coordinate events They constitute the basis of the Harriman representation of electronic states [14,36]. The coefficient function in the preceding equation, the momentum representation of � ⟩ can be subsequently where the transformed "volume" element Here, (f ) = ⟨f � ⟩ stands for the f representation of the state vector � ⟩ , while q (f ) similarly expresses the state � q ⟩ in terms of the density-dependent spatial variable f. These equations summarize the Fourier transformations between the momentum (�q⟩) and position ( �r⟩ or �f ⟩ ) representations of � ⟩ . Notice that in this vector interpretation, the orthogonality relation of Eq. (27) reads, thus manifesting the completeness of the EO basis yielding the molecular probability density: It is also of interest to express the HZM representation of � ⟩, (f ) = ⟨f � ⟩ , in terms of its standard position or momentum representations: (r) = ⟨r� ⟩ or (k) = ⟨k� ⟩ . Using the familiar Fourier transformation gives: with f (k) = ⟨k�f ⟩ and f (r) = ⟨r�f ⟩ = r(f ) * standing for the momentum and position representations of �f ⟩ , respectively.
T h e a d m i s s i b l e N -e l e c t r o n d e t e r m i n a n t s

Resultant descriptors, phase equilibria and information production
For simplicity, let us first consider a single electron in the specified quantum state � ⟩ of Eq. (1). The average Fisher's [1,2] measure of the gradient information for locality events contained in the electronic probability density p(r) = R(r) 2 [80], called the intrinsic accuracy, is reminiscent of von Weizsäcker's [81] inhomogeneity correction to kinetic energy functional, This classical measure characterizes an effective "narrowness" of the spatial probability distribution, a degree of electronic position determinicity. It represents the complementary measure to Shannon's [3,4] global entropy descriptor, Here, p(z) stands for the probability distribution of the coordinate "z," p(y|z) is the conditional probability of observing coordinate "y" given "z," and p(x|y,z) denotes conditional probability of observing coordinate "x," when one has already simultaneously observed coordinates "y" and "z" which reflects the position uncertainty, i.e., a "spread" of the distribution.
In quantum IT (QIT), these complementary classical descriptors are supplemented by the associated nonclassical complements in the resultant-entropy/information functionals, which combine the probability and phase/current contributions. The intrinsic accuracy concept generalizes naturally into the associated resultant descriptor of the gradient information content in the quantum state � ⟩ itself: This overall Fisher-type measure is defined by the expectation value of the Hermitian information operator Î (r) [61], related to kinetic energy operator T (r) of Eq. (6) . The sign of the latter expresses the fact that nonvanishing current pattern introduces more structural determinicity or order information about the system. This dimensionless measure is seen to also reflect the average kinetic energy of an electron It is of interest to examine the resultant gradient information content of EO. One first observes that q (r). satisfies the pseudo-eigenvalue problem of the momentum operator ̂ (r) = − iħ∇, which directly gives the expectation value of the information operator of Eq. (44): Its first (classical) contribution thus amounts to the Fisher information in probability distribution p(r), while the second (nonclassical) component, due to the orthogonality phase of Eq. (17), amounts to the plane-wave contribution 4q 2 multiplied by the EO-average value of the squared divergence of the phase vector field f(r).
One similarly generalizes the global entropy concept by supplementing the classical (probability) functional with the (negative) nonclassical contribution generated by the state (positive) phase [15][16][17][18][19][20][21]. Its sign implies a decrease in state overall structural indeterminicity due to its nonvanishing current distribution. The resultant global entropy measure then reads while the overall gradient entropy is given by the expectation value The negative character of the phase contribution M[ϕ] in the latter reflects the fact that a presence of a finite current decreases the overall structure-uncertainty (entropy) content in ψ.
The state extrema of these uncertainty descriptors identify the same optimum (nonnegative) equilibrium-phase solution, This entropy equilibrium ("thermodynamic") phase thus generates the associated EO current density reflecting the negative probability-gradient: One also introduces the complex overall entropy concept [11], the quantum expectation value of the non-Hermitian entropy operator ̂ (r) = − 2lnψ(r). In this "vector" measure of the state overall structural uncertainty, the two contributions of Eq. (49) constitute its real and imaginary parts: In the Schrödinger dynamical picture, the time evolution of resultant gradient information, the operator of which does not depend on time explicitly, Î (r) = − 4∇ 2 = σT(r) , results solely from the quantum dynamics of the system state vector � (t)⟩ , (53) j eq. (r) = (ℏ∕m)p(r)∇ eq. (r) = − ℏ∕(2m) ∇p(r).
determined by SE. The time derivative of the resultant gradient information [20], is then generated by the expectation value of the commutator and the integration by parts implies: Hence, the integral production (source) of this information descriptor reads: This derivative is seen to be determined by the current content of molecular electronic state. Therefore, it identically vanishes for the zero-current density everywhere, when the local component of the state phase identically vanishes, thus confirming its nonclassical origin. This qualitative conclusion also follows from the time- while the phase continuity relations [Eqs. (11) and (16)] give the following expression for the nonclassical information source: Therefore, the integral source of resultant gradient information in fact reflects the total time derivative of its nonclassical contribution I [ϕ]. Hence, the associated derivative of the overall gradient entropy of Eq. (50) reads: This result confirms the intuitive expectation that an increase in the state overall structural determinicity (order) information, σ I (t) > 0 implies the associated decrease in the structural indeterminicity (disorder) information (entropy): σ M (t) < 0. The above one-electron development can be straightforwardly generalized into a general N-electron system in the specified electronic state �Ψ(N)⟩ , exhibiting the electron density ρ(r) = Np(r). The N-electron information operator then combines terms due to each particle, and determines the state overall gradient information, proportional to the system average kinetic energy T(N). The relevant separation of the modulus and phase components of N-electron states calls for wavefunctions yielding the specified electron density. For example, in DFT this goal can be effected using the HZM construction, which we have outlined in Sect. 3.

Probability currents and information equilibria
Both parts of the resultant EO phases Φ q in Eq. (17) contribute to the overall probability current generated by the Slater determinant Ψ q (N) describing an electron configuration. The given EO φ q (r) generates the associated current density Hence, the overall current in a trial configuration of N electrons, corresponding to the N-electron operator where denotes the resultant "wave-number" vector of Ψ (N) , independent of the spatial position r, and ∇·f(r) stands for the divergence of its phase vector field f(r). A reference to Eq. (25) and Fig. 2 shows that this divergence corresponds to the probability density of the union (sum) U(r) = x∨y∨z of observing the specified position coordinates, which is represented by the area enclosed by the envelope of three probability circles in the figure [14,36].
Therefore, in the HZM construction of an antisymmetric wavefunction of N electrons the Jacobian of the r → f(r) transformation, expressing a dependence of the volume element df on the spatial volume dr, is proportional to molecular probability function of the product (joint) event involving the spatial coordinates of an electron, while its divergencedetermining the current density-corresponds to the probability of the union of such elementary coordinate events.
The overall current of Eq. (66) consists of the resultant "orthogonality" component j f (r) in HZM construction, determined by the local electron density ρ(r), the average wavenumber vector 〈q〉 and the divergence of the phase vector f(r), supplemented by the overall "thermodynamic" component j ϕ (r) determined by the "thermodynamic" phase-gradient, ∇ϕ(r). The orbital contributions {j l (r)} to the former, can be also regarded as the ∇·f(r)-enhanced plane-wave currents j(q l ; r) of Ψ q (N), with the local enhancement factor reflecting the probability density of Fig. 2.
It is also of interest to determine the information-optimum thermodynamic phase, which minimizes the nonclassical gradient information of EO, for the prescribed molecular probability distribution. Indeed, in such a vertical search, for the fixed (ground-state) probability density p 0 = ρ 0 /N, it is the thermodynamic phase ϕ(r) of EO states {ψ l [p 0 ]}, which is being optimized, since their orthogonality phases are uniquely determined by the electron density ρ 0 : f(r) = f[ρ 0 ; r] and {q l = q l [ρ 0 ]}. The condition of the minimum overall gradient information in configuration Ψ q [p 0 ], for the specified probability distribution p 0 , the shape factor of ρ 0 , then implies the minimum value of the overall nonclassical information, One observes that this information-optimum local phase of EO differs from its entropy-optimum ("equilibrium") analog of Eq. (52). It generates the orbital current in direction of the EO-average wave-number vector 〈q[p 0 ]〉: 108 Page 10 of 15 To summarize, the information-optimum thermodynamic phase of EO is determined by the average "wave-number" vector 〈q〉 in Ψ q (N) = Ψ q 1 ,q 1 ,…,q N (N) . The associated information-"equilibrium" EO thus exhibits the overall phase determined by the deviation of its wave-number vector q l from configuration-average vector 〈q〉: and generates the associated orbital current: Such phase-transformed EO thus gives rise to the vanishing resultant probability current: since The information-optimum equilibria in EO, for the specified electron density, thus ascribe finite orbital currents, which sum up to the vanishing resultant probability current in the system as a whole. Moreover, since in the informationoptimum state of Eq. (78) δq[p] replaces q[p], its resultant gradient information content [see Eq. (48)], exhibits the current contribution proportional to the squared deviation of the EO wave-number vector q from the configuration-average value 〈q〉.

Atomic bases and orbitals reconstructing overlap distributions
The EO states{ q (X) ≡ q X (r) } of constituent atoms X = A, B, … in a molecule, (77) j opt. ≡ j opt. p 0 ; r = (ℏ∕m)p 0 (r)∇ opt. p 0 ; r opt. q l form an admissible (non-orthogonal) molecular basis for electronic structure calculations. They give rise to a nonvanishing (complex) diatomic (X ≠ Y) overlap integral, in the molecular Hilbert space, the element of the associated (Hermitian) EO-overlap matrix: Since the complex AO are characterized by their probability and current distributions, the natural question arises about the chemically favorable combinations of these physical descriptors of AIM subsystems in a molecule. In particular, one would like to address the intriguing question, how to combine the probability currents to facilitate/enhance the chemical bonds between AIM or reactant fragments? The constructive superposition of subsystem probabilities in the bonding MO combinations of the real AO and their destructive superposition pattern in the antibonding MO states have been well established in the familiar MO theory of molecular electronic structure. Hence, the favorable MO combination of the subsystem EO calls for their summation in the probability-bonding combination, while the EO difference, should be associated with the probability-antibonding situation.
The phase/current factors of atomic EO dynamically influence a chemical character of the bonding MO combination describing a given X-Y bond. One recalls that directions of atomic currents in standard EO states, for ϕ = 0, are determined by the orbital wave-number vectors: q X = q and q Y = q′. Therefore, in the probability-bonding combination of Eq. (85), the inward orientation of EO currents, toward the bonding region between the nuclei, should enhance the covalent bond X-Y, while the outward pattern of EO currents, is expected to weaken this bond component. Therefore, the ("static") probability-bonding combination of Eq. (85) can assume both the ("dynamic") current-bonding and current-antibonding characters of Eqs. (87) and (88), respectively. In such EO combinations, the current pattern (the structure of "becoming") determines different time evolutions of the same initial (bonding) probability distribution (the structure of "being"). Similar current-bonding and current-antibonding combinations of the subsystem MO, e.g., substrates in chemical reactions, can be encountered in molecular interactions. Consider two interacting (parallel) benzene rings of reactants, perpendicular to "z" axis passing through their centers. In the bonding combination of complex MO exhibiting finite ring currents, the conrotatory matching of the latter should indeed result in an attraction between the magnetic moments associated with these substrate currents, while their disrotatory pattern should produce a magnetic repulsion.
One can also design the EO bases for the specified probability distributions of molecular fragments, e.g., substrates in a chemical reaction. This application should facilitate an orbital interpretation of molecular interactions. The HZM construction of antisymmetric states can be generalized to generate any probability distribution. In particular, one can design the orthonormal EO of N electrons which yields the prescribed diatomic overlap distribution [75], which represents the geometrical average of atomic probability densities {p Z (r)}, and then construct Slater determinants from alternative selections of different N such orbitals.
Following the development of Sect. 3, this overlap-HZM construction transforms r into with the associated Jacobian relation between volume elements: This can be effected through the overlap-dependent vector field F XY (r) = f[Ω XY , r] of the EO phase component, , It gives the functional Jacobian determinant of Eq. (88): Again, the diagonal derivatives of the preceding equation can be also interpreted as the associated conditional distributions: The orthonormality of such transformed EO conserving the given overlap distribution can be then demonstrated directly, following the substitution of Eq. (27): 108 Page 12 of 15 The standard HZM construction, reproducing the prescribed molecular probability density p(r), can be thus considered as representing a special case of the above overlap development for the equal fragment distributions: The generalized approach, however, allows one to generate the orthonormal basis for any average distribution, e.g., that resulting from the overlap between atomic orbitals, AIM, reactants, etc. Such an EO framework can be thus focused on the valence region of the overlap between the specified pair of AIM, which is mostly responsible for the chemical bond formation between such subsystems. It can be better suited for a more compact, chemical interpretation of the bonding patterns generated by atomic electron distributions.
As an illustrative example, consider an overlap between two spherical-Gaussian densities {G α (A α , R α )} corresponding to subsystems α = (a, b) (see "Appendix 1"), centered in positions {R α }, Their product can be transformed into an effective Gaussian function where ρ ab = (AR a + BR b )/(A + B). It also generates the Gaussian average density, which provides the distribution basis of the associated overlap EO of Eq. (90). As shown in "Appendix 1," the HZM conditional probabilities and local phase-vector components for such an effective spherical-Gaussian distribution are readily available.

Conclusion
In this article, we have first summarized the probability and phase continuities of molecular states, as implied be SE. We have then reexamined the HZM construction of DFT, of antisymmetric states yielding the prescribed electron density, and provided the probability interpretation of its key phase concept. It has been argued that this transformed position (97) p X (r) = p Y (r) = Ω XY (r) = p(r). representation in fact preserves main features of the momentum representation in molecular quantum mechanics. A need for the quantum extension of the familiar classical (probability) measures of the entropy/information content in molecular states, to accommodate the complex wave functions of the quantum-mechanical description, has been stressed, and resultant descriptors combining the probability and phase/ current contributions have been introduced. This quantum generalization of Fisher's gradient information, related to electronic kinetic energy, generates a nonvanishing information source. We have also demonstrated the nonclassical (current) origins of the associated resultant-information production. The present analysis and related treatments of reactivity phenomena [23][24][25][26][27][28] complement the previous classical information approaches to reactive systems, e.g., [82][83][84][85][86].
The orbital flows of electrons have been examined, and the information-optimum "thermodynamic" phases of the information-equilibrium EO states have been determined. It has been argued that the overall current of electron configurations in such information equilibria identically vanishes. The current orientations in the bonding combinations of atomic EO have been examined and the current-bonding or current-antibonding patterns have been identified. Finally, the EO basis reproducing the specified diatomic overlap density has been discussed and an illustrative example of the overlapping Gaussian distributions has been explored.
This article further develops the resultant-information description of molecular states [26][27][28]86]. The electronic state exhibiting a nonvanishing probability current indeed contains more information, compared to the zero-current state of the same electron density, since the former exhibits an additional structural element, and thus also the associated coherence information of the phase/flux pattern, which is missing in the state specification by the probability distribution alone. Following Prigogine [29], one could refer at this point to the now famous picture of a complicated network of vortices on Jupiter. The current structure of "becoming" contains the information about the (subsequent) quantum evolution of the system temporary probability distribution, the structure of "being," thus representing the state of a greater "knowledge" about the system, as reflected by the higher degree of the gradient "order" (determinicity) information, corresponding to lower level of the system "uncertainty," corresponding to smaller gradient "disorder" (indeterminicity) information (entropy) measure.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creat iveco mmons .org/licen ses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
The information-equilibrium thermodynamic phase, for the fixed probability distribution, then results from the extremum of the nonclassical information contribution, which gives the following Euler's equation for the information-optimum phase of EO: where we have recognized the phase continuities of Eq. (11): ∇ 2 F q = ∇ 2 ϕ opt. = 0. It thus follows from the preceding equation that