Wavelength-decomposition-based embedded cluster density approximation for systems with nonlocal electron correlation

Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one-dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short-wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long-wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn-Sham orbitals. This new method is termed λ-ECDA, where λ indicates the wavelength decomposition. The performance of λ-ECDA is examined on a one-dimensional model system: a H₂₄ chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron-electron and electron-ion interactions, and slightly stronger nuclear charges (1.2e ) are assigned to the pseudo-H atoms. Bond stretching energies, RPA correlation potentials, and Kohn-Sham eigenvalues predicted by λ-ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy.     

DFT studies and AIM analysis of AlN-polycycles

The structure and properties of AlN-polycycles were studied by DFT (density functional theory) method. The results of calculations were obtained at B3LYP/6-311G(d, p) level on model species. Topological parameters such as electron density, its Laplacian, kinetic electron energy density, potential electron energy density, and total electron energy density at the ring critical points (RCP) from Bader’s ‘Atoms in molecules’ (AIM) theory were analyzed in detail. These results indicate a good correlation between ρ(3, +1), G(r), H(r), and V(r) averaged values and hardness of AlN-polycycles. The aromaticity of all molecules has been studied by nucleus-independent chemical shift. There is a linear correlation between ΣNICS(0.0)_molecule values and polarizability.     

Correlation factor approach to the correlation energy functional

A global survey of the correlation factor energy functionals and its application to atomic and molecular properties is made. Its performances are compared with those of the density functional theory (DFT) correlation energy functionals, and some interesting conclusions from previous publications are reinforced here; namely, after removing the one-Slater-determinant hypothesis from the Kohn–Sham method, all DFT correlation functionals are able to provide reasonable results in any circumstance, with an additional restriction, for systems having a quasi-degenerate wave function, the DFT correlation functionals must depend explicitly on the on-top density. Acknowledgement.This work has been done under the support of the Spain DGICYT, project n⁰ BQU2001-0883.     

Exact electronic properties in the classically forbidden region of a metal surface

We derive exact properties of the inhomogeneous electron gas in the asymptotic classically forbidden region at a metal–vacuum interface within the framework of local effective potential energy theory. We derive a new expression for the asymptotic structure of the Kohn–Sham density functional theory (KS‐DFT) exchange‐correlation potential energy v_xc(r) in terms of the irreducible electron self‐energy. We also derive the exact asymptotic structure of the orbitals, density, the Dirac density matrix, the kinetic energy density, and KS exchange energy density. We further obtain the exact expression for the Fermi hole and demonstrate its structure in this asymptotic limit. The exchange‐correlation potential energy is derived to be v_xc(z → ∞) = ?α_KS,xc/z, and its exchange and correlation components to be v_x(z → ∞) = ?α_KS,x/z and v_c(z → ∞) = ?α_KS,c/z, respectively. The analytical expressions for the coefficients α_KS,xc and α_KS,x show them to be dependent on the bulk‐metal Wigner–Seitz radius and the barrier height at the surface. The coefficient α_KS,c = 1/4 is determined in the plasmon‐pole approximation and is independent of these metal parameters. Thus, the asymptotic structure of v_xc(z) in the vacuum region is image‐potential‐like but not the commonly accepted one of ?1/4z. Furthermore, this structure depends on the properties of the metal. Additionally, an analysis of these results via quantal density functional theory (Q‐DFT) shows that both the Pauli W_x(z → ∞) and lowest‐order correlation‐kinetic W(z → ∞) components of the exchange potential energy v_x(z → ∞), and the Coulomb W_c(z → ∞) and higher‐order correlation‐kinetic components of the correlation potential energy v_c(z → ∞), all contribute terms of O(1/z) to the structure. Hence correlations attributable to the Pauli exclusion principle, Coulomb repulsion, and correlation‐kinetic effects all contribute to the asymptotic structure of the effective potential energy at a metal surface. The relevance of the results derived to the theory of image states and to KS‐DFT is also discussed. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Application of second-order density functional methods to the calculation of the BeFH potential energy surface

Second-order density functional methods (where the correlation energy depends on the second-order density matrix and on a density functional) are used to introduce the electron correlation in two-configuration direct minimization (TCDM) ab initio electronic energy calculations of three-dimensional potential energy surfaces (PES). We analyzed the behavior of these methods in PES calculations by applying them to the Be + FH reaction. This system was studied also by the usual techniques, allowing a full comparison for the lowest ¹A′ adiabatic state. In particular, we compared the results obtained using the TCDM and multiple reference single and double excitations configuration interaction (MRDCI) methods with the corresponding results obtained using the Colle-Salvetti (CS) and Moscardó-San Fabián (MSF) procedures, within the correlation factor method, using as starting point the TCDM results. We found that the CS and MSF results are in a good overall agreement with the more accurate ab initio results, including the heights of the saddle points and the transition state. © 1997 John Wiley & Sons, Inc.     

Density functional calculation for ozone-ethylene complexes

The density functional method (exchange correlation functional Becke3LYP) is used to calculate the potential energy profiles for the reaction of ozone with ethylene. It is shown that “direct” epoxidation of the C = C double bond demands high activation energy and is unlikely for both themwchemical and photochemical reactions of ozonolysis. Translated fromZhumal Struktumoi Khimii, Vol. 41, No. 2, pp. 240–247,     

Reinvestigation of intramolecular hydrogen bond in malonaldehyde derivatives: An ab initio,AIM and NBO study

The RAHB systems in malonaldehyde and its derivatives at MP2/ 6‐311++G(d,p) level of theory were studied and their intramolecular hydrogen bond energies by using the related rotamers method was obtained. The topological properties of electron density distribution in O? H···O intramolecular hydrogen bond have been analyzed in term of quantum theory of atoms in molecules (QTAIM). Correlations between the H‐bond strength and topological parameters are probed. The results of QTAIM clearly showed that the linear correlation between the electron density distribution at HB critical point and RAHB ring critical point with the corresponding hydrogen bond energies was obtained. Moreover, it was found a linear correlation between the electronic potential energy density, V(r_cp), and hydrogen bond energy which can be used as a simple equation for evaluation of HB energy in complex RAHB systems. Finally, the similar linear treatment between the geometrical parameters, such as O···O or O? H distance, and Lp(O)→σ*_OH charge transfer energy with the intramolecular hydrogen bond energy is observed. © 2010 Wiley Periodicals, Inc., Int J Quantum Chem, 2011     

The effect of electron correlation on the topological and atomic properties of the electron density distributions of molecules

The topological properties of the electron density and the properties of an atom in a molecule are calculated by means of second-order Møller-Plesset perturbation theory (MP2) and compared with the results of configuration interaction calculations (C12) which include all single and double substitutions from the Hartree-Fock reference configuration. A software package for analyzing the effects of electron correlation on the topological properties of the electron density of molecules is described. H₂CO is used to provide a numerical example and to indicate that the number of bond critical points is unaffected by the inclusion of electron correlation. Correlation leads to only a small shift in the positions of bond critical points and a small change in the electron density at bond critical points. It is further shown that the energy of an atom in a molecule can be calculated to an accuracy of 1 kcal/mol and the electron population of an atom to about 0.001e. A statistical method is used to show that the deviation of the MP2 correlation correction relative to the CI2 correlation correction for a variety of atomic properties is about 25%.     

Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Analytical gradients

A full implementation of analytical energy gradients for molecular and periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian‐type orbitals as basis functions. Its key component is a combination of density fitting (DF) approximation and continuous fast multipole method (CFMM) that allows for an efficient calculation of the Coulomb energy gradient. For exchange‐correlation part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097) is extended to energy gradients. Computational efficiency and asymptotic O(N) scaling behavior of the implementation is demonstrated for various molecular and periodic model systems, with the largest unit cell of hematite containing 640 atoms and 19,072 basis functions. The overall computational effort of energy gradient is comparable to that of the Kohn–Sham matrix formation. © 2016 Wiley Periodicals, Inc.     

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.     

NDDO MO Calculations

Siegbahn's potential model as extended by Ellison et al. is used with density matrix elements calculated by the NDDO/2 procedure, to correlate the K-shell binding energy shifts of C, N and O atoms in a few molecules containing only the first-row atoms. The correlation is not superior to that obtained with the CNDO/2 method when only the monopole term is retained in calculating the Madelung potential energy. However, the results are in excellent agreement with experiment when the two-parameters model including the dipole and quadrupole terms is used.     

Additive density functional correlation corrections to single particle theories

The correlation energy functional E_c of the Hartree-Fock density is investigated. It was previously established that E_c produces the exact ground-state energy when added to the Hartree-Fock energy. Except when certain degeneracies occur, it is here shown that E_c is bounded from below when the coordinates of the Hartree-Fock density are scaled uniformly by λ, as λ → ∞. Consequently, approximations to E_c should display this bounded property, which is a second-order perturbation energy. It is also shown that a corresponding result applies to that correlation energy functional, Ë_c, which is to be added to a completed exact exchange-only (in the OPM sense) density functional calculation. Scaling requirements are presented for each order in the perturbation expansion for Ë_c. For instance, the second-order term is dimensionless. © 1997 John Wiley & Sons, Inc.     

Diffusion Monte Carlo study of correlation in the hydrogen molecule

The diffusion Monte Carlo (DMC) method shows that correlation in H₂ produces a set of three spatial changes: (i) an enhancement in the electron density distribution n( r ) in the left and right anti‐binding regions that include separately the immediate vicinity of each of the two nuclei, (ii) a reduction in n( r ) in the binding region intervening between the two nuclei as a counterbalance, and (iii) a concomitant increase in the equilibrium internuclear separation. It is stressed that the correlation energy E^c (= T^c + V^c) for diatomic molecules be defined by the difference in the total energy between the exact and the Hartree–Fock (HF) variational calculations that are performed at individually optimized internuclear separations. It is this definition that makes it possible to involve a significant contribution from a correlation‐induced change in the equilibrium internuclear separation as part of the correlation energy and to relate (i) and (ii) to (iii) in consistency with the electrostatic theorem. The present calculations fulfill the virial theorem to an accuracy of ?V/T = 2.00 for DMC and ?V^HF/T^HF = 2.000 for HF. The present correlation energy E^c = ?0.0408 hartree is not only in good agreement with the most accurate value previously reported, but also can be analyzed into all its components in accordance with the correlational virial theorem 2T^c + V^c = 0. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

Coulomb-Hole–Hartree–Fock functional

We have extended to molecules a density functional previously parametrized for atomic computations. The Coulomb-hole–Hartree–Fock functional, introduced by Clementi in 1963, estimates the dynamical correlation energy by the computations of a Hartree–Fock-type single-determinant wave function, where the Hartree–Fock potential was augmented with an effective potential term, related to a hard Coulomb hole enclosing each electron. The method was later revisited by S. Chakravorty and E. Clementi [Phys. Rev. A 39 , 2290 (1989)], where a Yukawa-type soft Coulomb hole replaced the previous hard hole; atomic correlation energies, computed for atoms with Z = 2 to Z = 54 as well as for a number of excited states, validated the method. In this article, we parametrized a function, which controls the width of the soft Coulomb hole, by fitting the first and second atomic ionization potentials of the atoms with 1 ? Z ? 18. The parametrization has been preliminarily validated by computing the dissociation energy for a number of molecules. A few-determinant version of the Coulomb-hole–Hartree–Fock method, necessary to account for the nondynamic correlation corrections, is briefly discussed. © 1994 John Wiley & Sons, Inc.     

Theoretical Studies of Distyrylbenzene and its Optical Properties

The electronic and geometrical properties of distyrylbenzene (DSB) are investigated by using chemistry theoretical calculation methods. Specifically, the excited state properties are studied by performing ab initio correlation interaction singlet (CIS) and time‐dependent density functional theory; the ground state and Raman activities are computed by density functional theory with the B3LYP method. Eight conformers of distyrylbenzene are found and they are derived from three isomers which are cis, cis‐, cis, trans‐, and trans, trans‐, respectively. The relative energy shows that each isomer of three types is separated with a large energy barrier, but a small energy difference of each conformer is found if they are in the same type. The transition state also shows the barrier between conformers is lower than isomers. The computed excited transition energies using ZINDO/S based on the optimized geometries at a DFT/B3LYP level with 6–31+G show an excellent agreement with experimental absorption spectra.     

Correlation of electron density and bond length to band gap for binary oxides and halides

The best quantity correlated to the electronic energy band gap is found for alkali and alkaline-earth metal oxides and halides with face centered cubic (fcc) structure based on density functional theory and Bader's atom-in-molecule theory. Previous studies show the correlation of the band gap to the ground state electron density at the bond critical point (BCP). Whereas, in quantum mechanics, the gap between the energy levels of one dimensional square well potential is inversely proportional to the square of the width of the well which is the metal–nonmetal chemical bond length in our case. These motivate the proposition of a new quantity Q, the ratio of the density at the BCP to the square of the bond length. Our study reveals that, for the aforementioned materials, the band gap has a strong correlation to Q when they are multiplied by the density at the BCP.     

Thomas–Fermi term as the simplest correction to the Weizsacker term

The nature of the correlation function G(1,2) appearing in the definition of the reduced first order density operator γ′(1,2) = ρ(1)^1/2ρ(2)^1/2G(1,2) is analyzed. It is shown that when G(1,2) is expanded in terms of plane waves in the context of a single-determinant approximation to the wave function, the correction to the Weizsacker term in the kinetic energy density expression is the Thomas–Fermi term.     

Correlation energy,correlated electron density,and exchange‐correlation potential in some spherically confined atoms

We report correlation energies, electron densities, and exchange‐correlation potentials obtained from configuration interaction and density functional calculations on spherically confined He, Be, Be²⁺, and Ne atoms. The variation of the correlation energy with the confinement radius R_c is relatively small for the He, Be²⁺, and Ne systems. Curiously, the Lee–Yang–Parr (LYP) functional works well for weak confinements but fails completely for small R_c. However, in the neutral beryllium atom the CI correlation energy increases markedly with decreasing R_c. This effect is less pronounced at the density‐functional theory level. The LYP functional performs very well for the unconfined Be atom, but fails badly for small R_c. The standard exchange‐correlation potentials exhibit significant deviation from the “exact” potential obtained by inversion of Kohn–Sham equation. The LYP correlation potential behaves erratically at strong confinements. © 2016 Wiley Periodicals, Inc.     

Investigation of correlation between impact sensitivities and bond dissociation energies in some triazole energetic compounds

In this article, density functional theory has been utilized to study on the correlation between impact sensitivities h _50% and the bond dissociation energies (BDEs) of nine triazole energetic explosives. By employing B3LYP and B3P86 method with the 6-311G** basis set, all the molecules have been fully optimized. The BDEs for removal of the NO₂ group in these compounds have also been calculated at the same level. Computed results show that BDEs calculated by B3LYP method are all less than those by B3P86 method. The relationship between the impact sensitivities and the weakest C–NO₂ bond dissociation energy (BDE) values have been investigated. The results indicate a good linear correlation between the impact sensitivity h _50% and the ratio (BDE/E) of the weakest BDE to the total energy E.