Correlation strength and information entropy

The correlation present in the nondegenerate ground state of an interacting Fermi system is discussed in terms of reduced density matrices and their cumulant expansion. By generalizing a result obtained for the interacting uniform electron gas (correlation induced exchange-hole narrowing), possible measures of the correlation strength in terms of natural occupation numbers (the eigenvalues of the true one-particle density matrix) are introduced. These quantities-the v-order nonidempotency and the information entropy of the natural occupation numbers-result from the correlated many-body wave function and characterize the ground-state correlation in addition to the usual correlation energy. The uniform electron gas serves as a first illustrative example. © 1995 John Wiley & Sons, Inc.     

A zero differential overlap study of chemical binding in ammonia-borane and carbony l-borane

A cndo/2D study of the charge distribution obtained through Mulliken population analysis in the ground state of the title compounds shows that the features of charge distribution found by severalab initio calculations are fairly well reproduced by this method. The one-particle density, the interference density at the mid-point of the bond axis and the kinetic part of the interference energy calculated through the deorthogonalized density matrices over a wide range of intermolecular separation between the donor and the acceptor show that the one-particle density and the interference density steadily grow with decreasing internuclear separation, while the kinetic interference energy starts with negative value at large distance, then decreases and passes through a minima near but above the equilibrium distance and then increases rapidly below it conforming to the characteristic general behaviour of the kinetic component of Morse curve. The orbital pairwise interference density and the corresponding kinetic energy components reveal that the orbitals involved in the covalent binding are σ_2p AO of B and 2S and σ_2p AO of N and C atoms in H₃B-NH₃ and H₃B-CO respectively.     

Direct evaluation of one-electron properties in coupled cluster methods

Summary An analysis of a method for approximate calculations of expectation values for one-electron operators from available coupled cluster amplitudes is presented and illustrated numerically for the polarizability of the Be atom. The one-particle density matrix resulting from the present approach is accurate through the fourth order in the electron correlation perturbation. It has been found that, in order to obtain quantitative agreement between the energy derivative results and the approximate expectation value formalism, the third orderT ₁ T ₂⁽⁰⁾ wave function term must be included into the calculation of the one-particle density matrix. The present method is also considered as a promising tool for calculations of higher-order atomic and molecular properties from high level correlated wave functions.     

The extended Koopmans' theorem Fock operator

By expanding the wave function of a system of N particles in terms of products of functions of one and (N-1) particles, the one-particle, nonlocal operator F?_EKT (extended Koopmans' theorem) is determined. It is shown that although this operator is nonhermitian, its eigenvalues and eigenfunctions represent the ionization energies and occupied orbitals, respectively. The eigenfunctions of F?_EKT are the one-particle functions that enter into the expansion of the wave function of the system as partners of the (N-1)-particle wave functions. The eingenvalues are also one-particle energies that, multipled by the orbital occupancy probalities, enter the expression for the total N-particle energy of the system.     

On the connections between Brueckner–coupled-cluster,density-dependent Hartree–Fock,and density functional theory

The existence of an effective one-particle Hamiltonian in the Brueckner coupled cluster model naturally leads to the definition of an effective interaction G, which is a function of the T₂ amplitudes. Two types of approximations to G are proposed: One is purely phenomenological, while the other is based on approximations to the Brueckner T₂ equation. In both cases, the resulting effective interaction may be viewed as electron-density-dependent. Generalizing Hartree–Fock theory to accommodate density-dependent interactions (DDHF ), a method is obtained that is capable of accounting for correlation effects in an independent particle framework. The heuristic Skyrme force, successfully used in nuclear physics to model nucleon–nucleon interactions, is presented here as an example of an effective electron–electron correlation interaction. Due to the δ-function character of the Skyrme force, it is possible to express the energy in this model by an integral over an energy density, thus formally providing a connection between DDHF and density functional theory for this particular case. An approximation to the Brueckner T₂ equation is also proposed in the coordinate representation. In this model, the density-matrix dependence of T₂ is reduced to a nonlocal electron density dependence by means of an expansion which introduces terms that depend on the gradient of the density. The first term in this expansion amounts to a “local density approximation” to Brueckner coupled cluster theory. © 1995 John Wiley & Sons, Inc.     

Approximate ansatz for the expansion of the spherically averaged wave function in terms of interelectronic separation r12 for the Hookean atom,atomic ions,and the H2 molecule

For the two‐electron Hookean atom, it is first emphasized that, for a specific force constant k = 1/4, the ground‐state wave function has a simple dependence on the interelectronic separation r₁₂, namely, (1 + ½r₁₂)exp(??r). For this two‐electron model, therefore, the study of Rassolov and Chipman on the electron–electron cusp conditions on the spherically averaged wave function for the N electron atomic ions can be generalized to all orders in the interelectronic separation r₁₂. This Hookean model has therefore been used to give some justification for an ansatz for the spherically averaged wave function in atomic ions with N electrons for N ≥ 2. Several approximate two‐electron wave functions satisfying the Rassolov and Chipman conditions were tested and found to give excellent results. Another ansatz has been tested numerically on the ground state of two‐electron atomic ions and the H₂ molecule. Finally, for the Hookean atom a partial differential equation that is essentially for the pair correlation density is given in the Appendix . © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 95: 21–29, 2003     

Radial limit of lithium revisited

Configuration interaction (CI) calculations are carried out for the ground state of lithium using a thoroughly optimized basis set of s-type Slater functions. They establish that the radial limit of the nonrelativistic energy of the ground ²S state of lithium is no higher than −7.448666443E_h. Thus, radial correlation accounts for 35.2% of the total correlation energy. The radial CI wave function predicts a significantly more accurate Fermi contact parameter than the Hartree-Fock wave function. However, the imbalanced treatment of electron correlation in the radial CI wave function leads to an excessively diffuse electron density that is worse than that of the Hartree-Fock wave function. © 1997 John Wiley & Sons, Inc.     

Computation of electronic transition moments: the length versus the velocity representation

Summary We have compared transition moments (TMs) obtained using the length and velocity representations for transitions from the ground state of H₂ to the lowest two¹ _u and two¹ _u ⁺ Rydberg states, theA ¹–X ¹⁺ transition in BH, and theA ¹ _u –X ¹ _g ⁺ transition in C₂. For H₂, the TMs in the length and velocity representations agree well even in cases where the one-particle basis is incomplete and the TM has not converged. For BH and C₂ the TM in the length representation converges rapidly with improvements in the one-particle basis set and is insensitive to inner-shell correlation. In contrast, in the velocity representation convergence with improvements in the one-particle basis is much slower, especially for C₂, and the TMs are significantly changed by inner-shell correlation. Thus the difference between the TMs in the length and velocity representations would not appear to be a viable diagnostic of TM convergence.     

A relationship between atomic correlation energy of neutral atoms and generalized entropy

In a previous article, hereafter named as Paper I, we have showed a relationship between atomic correlation energy of neutral atoms with 2 < Z < 29 and Tsallis entropy. In this article, we generalize this relation showing the link between the atomic correlation energy and a general form of entropy obtained from deformed algebra. The results evidence the role of both q and Δ parameters of the general entropy, in terms of contribution of the long‐range interactions in the correlation energy. The q and Δ values, obtained as best fit of the atomic correlation energies 2 < Z < 29, indicate that this general form of entropy reduces to the Tsallis one, reproducing well the trend of the correlation energy for low Z. Moreover, as a consequence of these values of the parameters, the state atomic wave function is more localized with respect to the wave function calculated in the limit of Shannon entropy. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Statistical mechanical theory of liquid entropy

The multiparticle correlation expansion for the entropy of a classical monatomic liquid is presented. This entropy expresses the physical picture in which there is no free particle motion, but, rather, each atom moves within a cage formed by its neighbors. The liquid expansion, including only pair correlations, gives an excellent account of the experimental entropy of most liquid metals, of liquid argon, and of the hard-sphere liquid. The pair correlation entropy is well approximated by a universal function of temperature. Higher-order correlation entropy, due to n-particle irreducible correlations for n ≥ 3, is significant in only a few liquid metals, and its occurrence suggests the presence of n-body forces. When the liquid theory is applied to the study of melting, we discover the important classification of normal and anomalous melting, according to whether there is not or is a significant change in the electronic structure upon melting, and we discover the universal disordering entropy for melting of a monatomic crystal. Interesting directions for future research are extension to include orientational correlations of molecules, theoretical calculation of the entropy of water, application to the entropy of the amorphous state, and correlational entropy of compressed argon. We clarify the relation among different entropy expansions in the recent literature. © 1994 John Wiley & Sons, Inc.     

Relation between the Slater–Kohn–Sham orbitals generated by the one-body potential V(r) of DFT and the Löwdin natural orbitals for inhomogeneous electron liquids

The Löwdin natural orbitals (NO) are defined as those orbitals which bring the first-order density matrix of a correlated electron assembly into diagonal form. Another one-particle density matrix with the same diagonal elements is the single-particle idempotent Dirac density matrix generated by the one-body potential of density functional theory. Here, we compare the off diagonal form of γ expanded in terms of the Slater–Kohn–Sham (SKS) orbitals generated by V(r) with the NO expansion of Löwdin for general inhomogeneous electron liquids. In particular, the equation of motion of the correlated γ is corrected from that of γ _s , both now containing the one-body potential V(r). To illustrate the theory presented here, we first construct an approximate, albeit accurate, correlated 1DM for the ground state of the He atom and display connections between the resulting NOs and the SKS orbitals. The second example we discuss, but now quite briefly, is that of the inhomogeneous electron liquid in crystalline Si, where the NO expansion is available from the literature.     

Density matrices for correlated Gaussians: Helium and dipositronium

Formulas are derived for the density matrices belonging to an n-particle wave function built on the basis of single-center explicitly correlated Gaussian basis functions. An explicit formula for the first-order density matrix, P(r₁, r′₁), is obtained for computing the probability distribution P(r₁, r₁). Other formulas are derived for matrix elements of the first-order density operator P on a basis of single-particle Gaussian orbitals so that natural orbitals (NOs) can be expressed in such a basis. The method is illustrated for the case of the ground state of the helium atom using the 16-term (geminal) wave function by Singer and Longstaff (E = −2.90233 au) and a set of even-tempered Gaussian orbitals. The resulting natural orbitals compare favorably with natural orbitals from Cl expansions. The method is also applied to our 20 term (trimal) wave function for the ground state of dipositronium (E = −0.51560 au). Analysis is made in this case for pair correlation functions of both the electron-electron and the positron-electron pairs; results include the radial distributions of these pairs and their relative angular momentum. © 1996 John Wiley & Sons, Inc.     

The extended Koopmans' theorem Fock operator and the generalized overlap amplitude one-electron operator

The wave function of a system may be expanded in terms of eigenfunctions of the N −1 electron Hamiltonian times one-particle functions known as generalized overlap amplitudes (GOAS). The one-electron operator whose eigenfunctions are the GOAS is presented, without using an energy-dependent term as in the one-particle Green function or propagator approach. It is shown that this operator and the extended Koopmans' theorem (EKT) one-electron operator are of similar form, but perform complementary roles. The GOA operator begins with one-electron densities and total energies of N −1 electron states to generate the two-matrix and total energy of an N-electron state. The EKT operator begins with the two-matrix of an N-electron state to generate one-electron densities and ionization potentials (or approximations thereto) for N −1 electron states. However, whereas the EKT orbitals must be linearly independent, no such restriction applies to the GOAS. © 1996 John Wiley & Sons, Inc.     

Thermodynamic study on the adsorption and micelle formation of long chain alkyltrimethylammonium chlorides

The surface tension of aqueous solutions of tetradecyl-trimethylammonium chloride (TTAC) and decyltrimethylammonium chloride (DeTAC) were measured as a function of temperature at concentrations below and above the critical micelle concentration under atmospheric pressure. The entropy and energy of adsorption from the monomeric state and from the micellar state and also the entropy and energy of micelle formation for TTAC were evaluated and compared with those of dodecyltrimethyl-ammonium chloride (DTAC). The values of Δ^M _W s and Δ^M _W u for TTAC and DTAC systems show that the micelle formation is driven by the entropy at low temperatures and by the energy at high temperatures. Received: 9 December 1997 Accepted: 4 March 1998     

Reduced hamiltonian orbitals. III. Unitarily invariant decomposition of hermitian operators

A decomposition of an N-particle operator as a sum of N + 1 components is defined such that, in the case of a model system employing a finite one-particle basis set, the decomposition is invariant under unitary transformations of the basis set. Applied to a two-particle Hamiltonian, this decomposition gives rise to the distinction between the independent-particle energy and the coupling energy defined in previous papers. Applied to the reduced density operator for a quantum state, the decomposition corresponds to partitioning the density into irreducible components. This partitioning is illustrated by graphs of electron density for the water molecule.     

N-Representability conditions generated by a one-particle grassmann factor of an antisymmetric function

N-representability conditions for a two-particle density operator implied by positive-semidefiniteness of the projection operator P^N+1(?¹ Λ Ψ^N) are derived and discussed. The operator P^N+1(?¹ Λ Ψ^N) projects onto an (N + 1)-particle antisymmetric function ?¹ Λ Ψ^N, the Grassmann product of a one-particle factor ?¹ and an N-particle factor Ψ^N. The polar subcone ??²_N(g, q) to the set of N-representable two-particle density operators ??²_N which corresponds to these conditions is found. It is shown that its extreme rays belong to two orbits for the action of the unitary group of transformations in one-particle Hilbert space. The facial structure of the convex set ??²_N exposed by elements of ??²_N(g, q) is analyzed. An example of the operator that changes the structure of its bottom eigenspace when the number of fermions N surpasses a certain value is noted. A new approach to the diagonal conditions for N-representability is found. It consists of the decomposition of the N-particle antisymmetric identity operator onto the mutually orthogonal projection operators.     

Basic theory and some applications concerning the ratio of shear viscosity to entropy density in dense fluids

Much interest continues relating to the conjecture from string theory of an inequality satisfied by the ratio η/s in dense fluids, where η is the shear viscosity and s the entropy density. First, we summarise the models which have prompted the proposal of the inequality. Second, we consider a model equation of state which may be appropriate for the dense fluids NH₃ and H₂O, both of which satisfy the inequality.     

Model universal coulomb hole function for many-electron systems

In the framework of the homogeneous electron gas theory we give a model function G_c of the Coulomb hole that can be considered as an approximate universal correlation function for many-electron systems. The function G_c reflects the right asymptotic behavior of the correlation function of an electron gas in high and low density limits and enables one to reproduce experimental correlation energies of a number of atoms of the first and second periods in a local density approximation with the relative error 0.3–4.2%. The estimate of contributions of electrons with parallel or antiparallel spins into correlation energy shows that in the domain of densities typical for atoms of the first and second periods, the Coulomb correlation of electrons with parallel spins is in high extent suppressed by the Fermi correlation.     

Magnetic order and valence fluctuation in a Pu?Ga intermetallic compound studied via a first principles calculation

In order to reveal electronic properties of a plutonium-gallium intermetallic compound (Pu₃Ga), and its potential implication for microscopic mechanisms for effects of Ga doping on the electronic and structural properties, as well as the phase stability of delta-phase Pu Ga alloy, a first principles calculation on the magnetic properties of this system is implemented by using density functional theory (DFT) plus on-site Coulomb repulsion U with nonmagnetic, ferromagnetic, and antiferromagnetic (AFM) orders, while the intermediate correlation effect, which is beyond the scope of pure itinerant and localized electronic model, is investigated by using a many-body technique combining DFT and dynamical mean-field theory considering the dynamical correlation effect due to the incompletely filled Pu 5f orbitals and the relativistic effect by inclusion of spin-orbit coupling (SOC). Our findings show that Pu₃Ga is a bad metal with AFM order, which is in good agreement with the experimental magnetic measurement. SOC further splitting Pu 5f states into j = 5/2 and j = 7/2 manifolds, the former exhibits metallic character, while the latter insulating feature. Occupation analysis establishes that an average occupancy of Pu 5f electrons in Pu₃Ga is n_f = 4.9598, this result together with the spectrum function indicates that 5f electrons in this system might be a localized state with strong valence fluctuation. Additionally, optimization of lattice parameter, density of state, and momentum-resolved electronic spectrum function are also presented.