Accurate ab initio density fitting for multiconfigurational self-consistent field methods

Using Cholesky decomposition and density fitting to approximate the electron repulsion integrals, an implementation of the complete active space self-consistent field (CASSCF) method suitable for large-scale applications is presented. Sample calculations on benzene, diaquo-tetra-mu-acetato-dicopper(II), and diuraniumendofullerene demonstrate that the Cholesky and density fitting approximations allow larger basis sets and larger systems to be treated at the CASSCF level of theory with controllable accuracy. While strict error control is an inherent property of the Cholesky approximation, errors arising from the density fitting approach are managed by using a recently proposed class of auxiliary basis sets constructed from Cholesky decomposition of the atomic electron repulsion integrals.     

Normalization and Fermi–Coulomb and Coulomb hole sum rules for approximate wave functions

For approximate wave functions, we prove the theorem that there is a one‐to‐one correspondence between the constraints of normalization and of the Fermi–Coulomb and Coulomb hole charge sum rules at each electron position. This correspondence is surprising in light of the fact that normalization depends on the probability of finding an electron at some position. In contrast, the Fermi–Coulomb hole sum rule depends on the probability of two electrons staying apart because of correlations due to the Pauli exclusion principle and Coulomb repulsion, while the Coulomb hole sum rule depends on Coulomb repulsion. We demonstrate the theorem for the ground state of the He atom by the use of two different approximate wave functions that are functionals rather than functions. The first of these wave function functionals is constructed to satisfy the constraint of normalization, and the second that of the Coulomb hole sum rule for each electron position. Each is then shown to satisfy the other corresponding sum rule. The significance of the theorem for the construction of approximate “exchange‐correlation” and “correlation” energy functionals of density functional theory is also discussed. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

A self-consistent coulomb bath model using density fitting

A self-consistent Coulomb bath model is presented to provide an accurate and efficient way of performing calculations for interfragment electrostatic and polarization interactions. In this method, a condensed-phase system is partitioned into molecular fragment blocks. Each fragment is embedded in the Coulomb bath due to other fragments. Importantly, the present Coulomb bath is represented using a density fitting method in which the electron densities of molecular fragments are fitted using an atom-centered auxiliary basis set of Gaussian type. The Coulomb bath is incorporated into an effective Hamiltonian for each fragment, with which the electron density is optimized through an iterative double self-consistent field (DSCF) procedure to realize the mutual many-body polarization effects. In this work, the accuracy of interfragment interaction energies enumerated using the Coulomb bath is tested, showing a good agreement with the exact results from an energy decomposition analysis. The qualitative features of many-body polarization effects are visualized by electron density difference plots. It is also shown that the present DSCF method can yield fast and robust convergence with near-linear scaling in performance with increase in system size.     

Parallelization of three-center electron repulsion integrals

The parallelization of the three-center electron repulsion integrals arising from the variational fitting of the Coulomb potential is presented. A scheme for dynamical load balancing of the corresponding loop structure is discussed. The implementation in the density functional theory program deMon using the message passing interface is described. The efficiency of the parallelization is analyzed by selected benchmark calculations     

Towards a force field based on density fitting

Total intermolecular interaction energies are determined with a first version of the Gaussian electrostatic model (GEM-0), a force field based on a density fitting approach using s-type Gaussian functions. The total interaction energy is computed in the spirit of the sum of interacting fragment ab initio (SIBFA) force field by separately evaluating each one of its components: electrostatic (Coulomb), exchange repulsion, polarization, and charge transfer intermolecular interaction energies, in order to reproduce reference constrained space orbital variation (CSOV) energy decomposition calculations at the B3LYP/aug-cc-pVTZ level. The use of an auxiliary basis set restricted to spherical Gaussian functions facilitates the rotation of the fitted densities of rigid fragments and enables a fast and accurate density fitting evaluation of Coulomb and exchange-repulsion energy, the latter using the overlap model introduced by Wheatley and Price [Mol. Phys. 69, 50718 (1990)]. The SIBFA energy scheme for polarization and charge transfer has been implemented using the electric fields and electrostatic potentials generated by the fitted densities. GEM-0 has been tested on ten stationary points of the water dimer potential energy surface and on three water clusters (n = 16,20,64). The results show very good agreement with density functional theory calculations, reproducing the individual CSOV energy contributions for a given interaction as well as the B3LYP total interaction energies with errors below kBT at room temperature. Preliminary results for Coulomb and exchange-repulsion energies of metal cation complexes and coupled cluster singles doubles electron densities are discussed.     

Attractive electron–electron interactions within robust local fitting approximations

An analysis of Dunlap's robust fitting approach reveals that the resulting two‐electron integral matrix is not manifestly positive semidefinite when local fitting domains or non‐Coulomb fitting metrics are used. We present a highly local approximate method for evaluating four‐center two‐electron integrals based on the resolution‐of‐the‐identity (RI) approximation and apply it to the construction of the Coulomb and exchange contributions to the Fock matrix. In this pair‐atomic resolution‐of‐the‐identity (PARI) approach, atomic‐orbital (AO) products are expanded in auxiliary functions centered on the two atoms associated with each product. Numerical tests indicate that in 1% or less of all Hartree–Fock and Kohn–Sham calculations, the indefinite integral matrix causes nonconvergence in the self‐consistent‐field iterations. In these cases, the two‐electron contribution to the total energy becomes negative, meaning that the electronic interaction is effectively attractive, and the total energy is dramatically lower than that obtained with exact integrals. In the vast majority of our test cases, however, the indefiniteness does not interfere with convergence. The total energy accuracy is comparable to that of the standard Coulomb‐metric RI method. The speed‐up compared with conventional algorithms is similar to the RI method for Coulomb contributions; exchange contributions are accelerated by a factor of up to eight with a triple‐zeta quality basis set. A positive semidefinite integral matrix is recovered within PARI by introducing local auxiliary basis functions spanning the full AO product space, as may be achieved by using Cholesky‐decomposition techniques. Local completion, however, slows down the algorithm to a level comparable with or below conventional calculations. © 2013 Wiley Periodicals, Inc.     

Intermolecular electrostatic energies using density fitting

A method is presented to calculate the electron-electron and nuclear-electron intermolecular Coulomb interaction energy between two molecules by separately fitting the unperturbed molecular electron density of each monomer. This method is based on the variational Coulomb fitting method which relies on the expansion of the ab initio molecular electron density in site-centered auxiliary basis sets. By expanding the electron density of each monomer in this way the integral expressions for the intermolecular electrostatic calculations are simplified, lowering the operation count as well as the memory usage. Furthermore, this method allows the calculation of intermolecular Coulomb interactions with any level of theory from which a one-electron density matrix can be obtained. Our implementation is initially tested by calculating molecular properties with the density fitting method using three different auxiliary basis sets and comparing them to results obtained from ab initio calculations. These properties include dipoles for a series of molecules, as well as the molecular electrostatic potential and electric field for water. Subsequently, the intermolecular electrostatic energy is tested by calculating ten stationary points on the water dimer potential-energy surface. Results are presented for electron densities obtained at four different levels of theory using two different basis sets, fitted with three auxiliary basis sets. Additionally, a one-dimensional electrostatic energy surface scan is performed for four different systems (H2O dimer, Mg2+-H2O, Cu+-H2O, and n-methyl-formamide dimer). Our results show a very good agreement with ab initio calculations for all properties as well as interaction energies.     

Effect of Coulomb interaction between the electrons on two-electron redox-mediated tunneling

Two-electron non-adiabatic redox-mediated tunneling through a symmetric electrochemical contact with a bridge molecule having one electron energy level participating in tunneling is considered under ambient conditions. It is shown that the current/overpotential dependence in this system can disclose two distinct or overlapping clear-cut maxima depending on the value of the effective Coulomb repulsion energy. This new effect is due to the opening of the channel for tunneling of second electron with the variation of the electrode potential. The system manifests also a rectification effect in the current/bias voltage curve which depends on the value of the effective Coulomb repulsion energy.     

Fundamental importance of the Coulomb hole sum rule to the understanding of the Colle-Salvetti wave function functional

In this paper we consider the general form of the correlated-determinantal wave function functional of Colle and Salvetti (CS) for the He atom. The specific form employed by CS is the basis for the widely used CS correlation energy formula and the Lee-Yang-Parr correlation energy density functional of Kohn-Sham density functional theory. We show the following: (i) The key assumption of CS for the determination of this wave function functional, viz., that the resulting single-particle density matrix and the Hartree-Fock theory Dirac density matrix are the same, is equivalent to the satisfaction of the Coulomb hole sum rule for each electron position. The specific wave function functional derived by CS does not satisfy this sum rule for any electron position. (ii) Application of the theorem on the one-to-one correspondence between the Coulomb hole sum rule for each electron position and the constraint of normalization for approximate wave functions then proves that the wave function derived by CS violates charge conservation. (iii) Finally, employing the general form of the CS wave function functional, the exact satisfaction of the Coulomb hole sum rule at each electron position then leads to a wave function that is normalized. The structure of the resulting approximate Coulomb holes is reasonably accurate, reproducing both the short- and the long-range behavior of the hole for this atom. Thus, the satisfaction of the Coulomb hole sum rule by an approximate wave function is a necessary condition for constructing wave functions in which electron-electron repulsion is represented reasonably accurately.     

Iterative introduction of integral approximations with error bounds in SCF calculations

Another possible application of a previously reported approximation theory for electron repulsion integrals using rigorous error bounds is considered by incorporating the electron density matrix in the approximation scheme. Error bounds are set on the contribution of a given integral to the total energy for a given molecular wave-function; the wave-function is then refined cyclically and additional integrals are computed exactly if necessary until convergence is achieved.     

Poisson-transformed density fitting in relativistic four-component Dirac-Kohn-Sham theory

We present recent developments in the implementation of the density fitting approach for the Coulomb interaction within the four-component formulation of relativistic density functional theory [Belpassi et al., J. Chem. Phys. 124, 124104 (2006)]. In particular, we make use of the Poisson equation to generate suitable auxiliary basis sets and simplify the electron repulsion integrals [Manby and Knowles, Phys. Rev. Lett. 87, 163001 (2001)]. We propose a particularly simple and efficient method for the generation of accurate Poisson auxiliary basis sets, based on already available standard Coulomb fitting sets. Just as is found in the nonrelativistic case, we show that the number of standard auxiliary fitting functions that need to be added to the Poisson-generated functions in order to achieve a fitting accuracy equal or, in some cases, better than that of the standard procedure is remarkably small. The efficiency of the present implementation is demonstrated in a detailed study of the spectroscopic properties and energetics of several gold containing systems, including the Au dimer and the CsAu molecule. The extraction reaction of a H(2)O molecule from a Au(H(2)O)(9) (+) cluster is also calculated as an example of mixed heavy-light-atom molecular systems. The scaling behavior of the algorithm implemented is illustrated for some closed shell gold clusters up to Au(5) (+). The increased sparsity of the Coulomb matrices involved in the Poisson fitting is identified, as are potential computational applications and the use of the Poisson fitting for the relativistic exchange-correlation problem.     

Electron correlation: the many-body problem at the heart of chemistry

The physical interactions among electrons and nuclei, responsible for the chemistry of atoms and molecules, is well described by quantum mechanics and chemistry is therefore fully described by the solutions of the Schr?dinger equation. In all but the simplest systems we must be content with approximate solutions, the principal difficulty being the treatment of the correlation between the motions of the many electrons, arising from their mutual repulsion. This article aims to provide a clear understanding of the physical concept of electron correlation and the modern methods used for its approximation. Using helium as a simple case study and beginning with an uncorrelated orbital picture of electronic motion, we first introduce Fermi correlation, arising from the symmetry requirements of the exact wave function, and then consider the Coulomb correlation arising from the mutual Coulomb repulsion between the electrons. Finally, we briefly discuss the general treatment of electron correlation in modern electronic-structure theory, focussing on the Hartree-Fock and coupled-cluster methods and addressing static and dynamical Coulomb correlation.     

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.     

Grid-free density functional calculations on periodic systems

Density fitting scheme is applied to the exchange part of the Kohn-Sham potential matrix in a grid-free local density approximation for infinite systems with translational periodicity. It is shown that within this approach the computational demands for the exchange part scale in the same way as for the Coulomb part. The efficiency of the scheme is demonstrated on a model infinite polymer chain. For simplicity, the implementation with Dirac-Slater Xalpha exchange functional is presented only. Several choices of auxiliary basis set expansion coefficients were tested with both Coulomb and overlap metric. Their effectiveness is discussed also in terms of robustness and norm preservation.     

Scaling relations,virial theorem,and energy densities for long‐range and short‐range density functionals

Decomposition of the Coulomb electron–electron interaction into a long‐range and a short‐range part is described within the framework of density functional theory, deriving some scaling relations and the corresponding virial theorem. We study the behavior of the local density approximation in the high‐density limit for the long‐range and the short‐range functionals by carrying out a detailed analysis of the correlation energy of a uniform electron gas interacting via a long‐range‐only electron–electron repulsion. Possible definitions of exchange and correlation energy densities are discussed and clarified with some examples. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Electron correlation effects in the Fe dimer

The potential energy surface of the Fe dimer is investigated on the basis of density functional theory in the generalized gradient approximation (GGA). Electron correlation effects are taken into account explicitly within the GGA+U approach. We find a value of 2.20 eV for the Coulomb repulsion parameter U to describe the Fe dimer best, yielding a 9 Sigma(g)- ground state with an interatomic separation of 2.143 A. Agreement of the associated vibrational frequency, binding energy, ionization potential, and electron affinity with experimental data as well as corresponding results calculated within a high-level ab initio approach is improved significantly compared to conventional GGA. The effect of U on calculated geometric and magnetic properties of larger Fe clusters is discussed.     

Phase dilemma in density matrix functional theory

For closed-shell systems, a particular parametrization of coefficients in a configuration interaction (CI) expansion provides a convenient formulation for the search over electronic wave functions constrained by a set of natural orbitals (NOs) and the corresponding occupation numbers that are invoked in every variational construction of the density matrix functional (DMF) V(ee)(Gamma) for the electron--electron repulsion energy. It produces an explicit expression for V(ee) in terms of the Coulomb and exchange integrals over NOs, and an idempotent matrix omega, diagonal elements of which equal the occupation numbers. At the same time, it reveals a very serious bottleneck affecting any rigorous approach to the DMF theory, namely the phase dilemma that stems from the necessity to carry out minimization over a large number of possible combinations of CI coefficient signs. While underscoring its lack of variational nature, a simple approximation for the phase factor products provides a strict derivation for the recently proposed Kollmar-Hess functional.     

Tensor product approximation with optimal rank in quantum chemistry

Tensor product decompositions with optimal separation rank provide an interesting alternative to traditional Gaussian-type basis functions in electronic structure calculations. We discuss various applications for a new compression algorithm, based on the Newton method, which provides for a given tensor the optimal tensor product or so-called best separable approximation for fixed Kronecker rank. In combination with a stable quadrature scheme for the Coulomb interaction, tensor product formats enable an efficient evaluation of Coulomb integrals. This is demonstrated by means of best separable approximations for the electron density and Hartree potential of small molecules, where individual components of the tensor product can be efficiently represented in a wavelet basis. We present a fairly detailed numerical analysis, which provides the basis for further improvements of this novel approach. Our results suggest a broad range of applications within density fitting schemes, which have been recently successfully applied in quantum chemistry.     

Electrochemical transistor based on bridge tunneling contact containing two redox groups

The method of kinetic equations was used to show that a bridge tunneling contact containing two redox groups in the sequential configuration immersed into a solution of electrolyte at the room temperature features pronounced transistor properties at a given set of physical system parameters. Valence electrons of redox groups interact strongly with the classical phonon subsystem of the liquid medium. Debye screening of the electric field in the tunneling gap and Coulomb repulsion between electrons in different redox groups are taken into account. The case of nonadiabatic electron transfer both between redox groups and between electrodes and redox groups is considered in the limit of infinitely high Coulomb repulsion between electrons in a redox group. For sufficiently high absolute values of difference δ between unperturbed energy levels of redox groups, the system features voltammetric characteristics typical for a transistor. The amplification effect appears due to a strong dependence of tunneling current on overpotential. The emphasis is upon the peculiaritiespeculiarities of voltammetric characteristics in the case of asymmetric tunneling contacts.