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.     

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.     

On robust density fitting in molecules and extended systems

Two theorems for robust density fitting are proved: A. For a given electron density, the best robust density fitting approximation to the Coulomb electron repulsion energy results from the unconstrained Coulomb metric fit, B. For infinite periodic systems, the necessary condition for correct long-range behavior of any robust density fitting approximation of the Coulomb electron repulsion energy is the exact reproduction of the number of electrons.     

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.     

Electron density distribution in potassium bis -(carbonato)cuprate(II)

The electron density distribution in potassiumbis-(carbonato)cuprate(II) has been analyzed using x-ray diffraction data from an earlier structure determination. While the copper-ligand geometry is close to square planar the deformation density near the metal is strongly asymmetric. There are local maxima near the copper atom along the line of the Cu-K vectors. These resemble features found in corresponding regions in normal length metal-metal bonds. The observation is consistent with the long range nature of the Coulomb potential associated with the potassium ion.     

Grid‐based density functional calculations of many‐electron systems

Exploratory variational pseudopotential density functional calculations are performed for the electronic properties of many‐electron systems in the 3D cartesian coordinate grid (CCG). The atom‐centered localized gaussian basis set, electronic density, and the two‐body potentials are set up in the 3D cubic box. The classical Hartree potential is calculated accurately and efficiently through a Fourier convolution technique. As a first step, simple local density functionals of homogeneous electron gas are used for the exchange‐correlation potential, while Hay‐Wadt‐type effective core potentials are employed to eliminate the core electrons. No auxiliary basis set is invoked. Preliminary illustrative calculations on total energies, individual energy components, eigenvalues, potential energy curves, ionization energies, and atomization energies of a set of 12 molecules show excellent agreement with the corresponding reference values of atom‐centered grid as well as the grid‐free calculation. Results for three atoms are also given. Combination of CCG and the convolution procedure used for classical Coulomb potential can provide reasonably accurate and reliable results for many‐electron systems. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

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     

Self-interaction error in density functional theory: a mean-field correction for molecules and large systems

Corrections to the self-interaction error which is rooted in all standard exchange-correlation functionals in the density functional theory (DFT) have become the object of an increasing interest. After an introduction reminding the origin of the self-interaction error in the DFT formalism, and a brief review of the self-interaction free approximations, we present a simple, yet effective, self-consistent method to correct this error. The model is based on an average density self-interaction correction (ADSIC), where both exchange and Coulomb contributions are screened by a fraction of the electron density. The ansatz on which the method is built makes it particularly appealing, due to its simplicity and its favorable scaling with the size of the system. We have tested the ADSIC approach on one of the classical pathological problem for density functional theory: the direct estimation of the ionization potential from orbital eigenvalues. A large set of different chemical systems, ranging from simple atoms to large fullerenes, has been considered as test cases. Our results show that the ADSIC approach provides good numerical values for all the molecular systems, the agreement with the experimental values increasing, due to its average ansatz, with the size (conjugation) of the systems.     

Key thermodynamic characteristics of nucleation on charged and neutral cores of molecular sizes in terms of the gradient density functional theory

The gradient density functional theory and the Carnahan–Starling model formulated for describing the contribution of hard spheres have been used to calculate the profiles of condensate density in small critical droplets formed via homogeneous nucleation, as well as in stable and critical droplets formed via heterogeneous nucleation on solid charged and neutral condensation cores of molecular sizes. The calculations performed for water and argon at different values of condensate chemical potential have yielded the heights of the activation barriers for homoand heterogeneous nucleation as functions of vapor supersaturation at preset system temperatures. The interaction of condensate molecules with a solid core has been described by the resultant potential of molecular attractive forces. In the case of a charged core, the long-range Coulomb potential of electric forces has additionally been taken into account. Dielectric permittivities have been calculated as known functions of the local density of the fluid and temperature. The radius of the equimolecular droplet surface has been chosen as a variable describing the droplet size. Dependences of the chemical potential of condensate molecules in a droplet on its size have been plotted for water and argon with allowance for the action of capillary, electrostatic, and molecular forces. It has been shown that the role of the molecular force potential in heterogeneous nucleation increases with the size of condensation cores.     

On Koopmans' theorem in density functional theory

This paper clarifies why long-range corrected (LC) density functional theory gives orbital energies quantitatively. First, the highest occupied molecular orbital and the lowest unoccupied molecular orbital energies of typical molecules are compared with the minus vertical ionization potentials (IPs) and electron affinities (EAs), respectively. Consequently, only LC exchange functionals are found to give the orbital energies close to the minus IPs and EAs, while other functionals considerably underestimate them. The reproducibility of orbital energies is hardly affected by the difference in the short-range part of LC functionals. Fractional occupation calculations are then carried out to clarify the reason for the accurate orbital energies of LC functionals. As a result, only LC functionals are found to keep the orbital energies almost constant for fractional occupied orbitals. The direct orbital energy dependence on the fractional occupation is expressed by the exchange self-interaction (SI) energy through the potential derivative of the exchange functional plus the Coulomb SI energy. On the basis of this, the exchange SI energies through the potential derivatives are compared with the minus Coulomb SI energy. Consequently, these are revealed to be cancelled out only by LC functionals except for H, He, and Ne atoms.     

Pair density related to one‐electron information for the ground state of spin‐compensated two‐electron systems

The recent study by Joubert on effects of Coulomb repulsions in a many‐electron system has focused attention on an integral identity involving the pair density. This has motivated the derivation presented here of a vectorial differential form related to this integral result. Our differential identity is then illustrated explicitly by using (i) an exact ground‐state wave function for the so‐called Hookean atom having external potential energy (1/2)kr², with k = 1/4, and (ii) Moshinsky's model in which both the interparticle interaction and the external potential are of harmonic type. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Performance of parallel TURBOMOLE for density functional calculations

The parallelization of density functional treatments of molecular electronic energy and first-order gradients is described, and the performance is documented. The quadrature required for exchange correlation terms and the treatment of exact Coulomb interaction scales virtually linearly up to 100 nodes. The RI-J technique to approximate Coulomb interactions (by means of an auxiliary basis set approximation for the electron density) even shows superlinear speedup on distributed memory architectures. The bottleneck is then linear algebra. Demonstrative application examples include molecules with up to 300 atoms and 3000 basis functions that can now be treated in a few hours per geometry optimization cycle in C₁ symmetry. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1746–1757, 1998     

Nonlocal Wigner-like correlation energy density functional: parametrization and tests on two-electron systems

Reparametrization of Wigner's correlation energy density functional yields a very close fit to the correlation energies of the helium isoelectronic sequence. However, a quite different reparametrization is required to obtain an equally close fit to the isoelectronic sequence of Hooke's atom. In an attempt to avoid having to reparametrize the functional for different choices of the one-body potential, we propose a parametrization that depends on global characteristics of the ground-state electron density as quantified by scale-invariant combinations of expectation values of local one-body operators. This should be viewed as an alternative to the density-gradient paradigm, allowing one to introduce the nonlocal dependence of the density functional on the density in a possibly more effective way. Encouraging results are obtained for two-electron systems with one-body potentials of the form r(zeta) with zeta=-12,+12,1, which span the range between the Coulomb potential (zeta=-1) and the Hooke potential (zeta=2).     

A short-range correlation energy density functional with multi-determinantal reference

We introduce a short-range correlation density functional defined with respect to a multi-determinantal reference which is meant to be used in a multi-determinantal extension of the Kohn–Sham scheme of density functional theory based on a long-range/short-range decomposition of the Coulomb electron–electron interaction. We construct the local density approximation for this functional and discuss its performance on the He atom.     

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.     

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.     

Note on the density matrix for an arbitrary number of closed shells in a bare Coulomb field

An explicit relation between the density matrix and its s-state part is analyzed for electron closed shells moving in a bare Coulomb potential. The density matrix has a simple separable form in terms ofr ₁ +r ₂ and |r ₁ –r ₂|. It is demonstrated that for an arbitrary number of closed shells, the off-diagonal dependence is simply polynomial in the |r ₁ –r ₂| coordinate.     

Electron density fitting for the Coulomb problem in relativistic density-functional theory

A density fitting approach for the Coulomb matrix representation within the four-component formulation of relativistic density-functional theory is presented. Our implementation, which uses G-spinor basis sets, shares all the advantages of those found in nonrelativistic quantum chemistry. We show that very accurate Coulomb energies may be obtained using a modest number of scalar auxiliary basis functions for molecules containing heavy atoms. The efficiency of this new implementation is demonstrated in a detailed study of the spectroscopic properties of the gold dimer, and its scaling behavior has been tested by calculations of some closed-shell gold clusters (Au2, Au3+, Au4, Au5+). The algorithm is found to scale as O(N3), just as it does in the nonrelativistic case, and represents a dramatic improvement in efficiency over the conventional approach in the calculation of the Coulomb matrix, with computation times that are reduced to less than 3% for Au2 and up to 1% in the case of Au5+.     

Experimental and theoretical electron density study of estrone

The electron density and the electrostatic potential (ESP) distributions of estrone have been determined using X-ray diffraction analysis and compared with theoretical calculations in the solid and gas phases. X-ray diffraction measurements are performed with a Rigaku Rapid rotating anode diffractometer at 20 K. The electron density in the estrone crystal has been described with the multipole model, which allowed extensive topological analysis and calculation of the ESP. From DFT calculations in the solid state a theoretical X-ray diffraction data set has been produced and treated in the same way as the experimental data. Two sets of single molecule DFT calculations were performed: (a) An electron density distribution was obtained via a single-point calculation with a large basis set at the experimental geometry and subsequently analyzed according to the quantum theory of atoms in molecules (AIM) to obtain the bond and most atomic properties, and (b) another electron density distribution was obtained with a smaller basis set, but at a geometry optimized using the same basis set for the analysis of atomic energies. An interesting locally stabilizing hydrogen-hydrogen bond path linking H(1) and H(11B) is found which represents the first characterization of such bonding in a steroid molecule. AIM delocalization indices were shown to be well correlated to the experimental electron density at the bond critical points through an exponential relationship. The aromaticity of ring A, chemical bonding, the O(1)...O(2) distance necessary for estrogenic activity, and the electrostatic potential features are also discussed.