Calculation of electrostatic and polarization energies from electron densities

Fast analytical evaluation of intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density. III. Application to crystal structures via the Ewald and direct summation methods

The previously reported exact potential and multipole moment (EP/MM) method for fast and accurate evaluation of the intermolecular electrostatic interaction energies using the pseudoatom representation of the electron density [Volkov, Koritsanszky & Coppens (2004). Chem. Phys. Lett. 391, 170-175; Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536; Nguyen & Volkov (2019). Acta Cryst. A75, 448-464] is extended to the calculation of electrostatic interaction energies in molecular crystals using two newly developed implementations: (i) the Ewald summation (ES), which includes interactions up to the hexadecapolar level and the EP correction to account for short-range electron-density penetration effects, and (ii) the enhanced EP/MM-based direct summation (DS), which at sufficiently large intermolecular separations replaces the atomic multipole moment approximation to the electrostatic energy with that based on the molecular multipole moments. As in the previous study [Nguyen, Kisiel & Volkov (2018). Acta Cryst. A74, 524-536], the EP electron repulsion integral is evaluated analytically using the Löwdin α-function approach. The resulting techniques, incorporated in the XDPROP module of the software package XD2016, have been tested on several small-molecule crystal systems (benzene, L-dopa, paracetamol, amino acids etc.) and the crystal structure of a 181-atom decapeptide molecule (Z = 4) using electron densities constructed via the University at Buffalo Aspherical Pseudoatom Databank [Volkov, Li, Koritsanszky & Coppens (2004). J. Phys. Chem. A, 108, 4283-4300]. Using a 2015 2.8 GHz Intel Xeon E3-1505M v5 computer processor, a 64-bit implementation of the Löwdin α-function and one of the higher optimization levels in the GNU Fortran compiler, the ES method evaluates the electrostatic interaction energy with a numerical precision of at least 10^-5 kJ mol^-1 in under 6 s for any of the tested small-molecule crystal structures, and in 48.5 s for the decapeptide structure. The DS approach is competitive in terms of precision and speed with the ES technique only for crystal structures of small molecules that do not carry a large molecular dipole moment. The electron-density penetration effects, correctly accounted for by the two described methods, contribute 28-64% to the total electrostatic interaction energy in the examined systems, and thus cannot be neglected.

One-particle density matrix polarization susceptibility tensors

The electric field-induced change in the one-electron density has been expressed as a series of the one-particle density matrix susceptibilities interacting with the spatial distribution of the electric field. The analytic approximate expressions are derived at the Hartree-Fock theory, which serves as a basis for the construction of the generalized model that is designed for an arbitrary form of wavefunction and any type of one-particle density matrix. It is shown that it is possible to accurately predict the changes in the one-electron ground-state density of water molecule in a spatially uniform electric field, as well as in spatially non-uniform electric field distribution generated by point charges. When both linear and quadratic terms with respect to the electric field are accounted for, the electric field-induced polarization energies, dipole moments, and quadrupole moments are quantitatively described by the present theory in electric fields ranging from weak to very strong (0.001-0.07 a.u.). It is believed that the proposed model could open new routes in quantum chemistry for fast and efficient calculations of molecular properties in condensed phases.

CheckDen, a program to compute quantum molecular properties on spatial grids

CheckDen, a program to compute quantum molecular properties on a variety of spatial grids is presented. The program reads as unique input wavefunction files written by standard quantum packages and calculates the electron density rho(r), promolecule and density difference function, gradient of rho(r), Laplacian of rho(r), information entropy, electrostatic potential, kinetic energy densities G(r) and K(r), electron localization function (ELF), and localized orbital locator (LOL) function. These properties can be calculated on a wide range of one-, two-, and three-dimensional grids that can be processed by widely used graphics programs to render high-resolution images. CheckDen offers also other options as extracting separate atom contributions to the property computed, converting grid output data into CUBE and OpenDX volumetric data formats, and perform arithmetic combinations with grid files in all the recognized formats.

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.

Universal method for computation of electrostatic potentials

A computational approach to determine electrostatic interaction and gravitational potentials by performing direct numerical integration is presented. The potential is expanded using finite-element functions of arbitrary order. The method does not involve any solutions of systems of linear equations. The potential is instead obtained as a sum of differential contributions. Thus, no boundary conditions for the potential are needed. It is computationally efficient and well suited for parallel computers, since the innermost loops constitute matrix multiplications and the outer ones can be used as parallel indices. Without using prescreening or other computational tricks to speed up the calculation, the algorithm scales as N4/3 where N denotes the grid size.

Determination of noncovalent interaction energies from electronic densities

Starting from the Hellmann-Feynman theorem, an expression is derived for the interaction energy in forming a noncovalently bound complex. It invokes only classical electrostatics, involving the charge distributions of the components as they are in the complex. We propose a method for obtaining these from the total charge distribution. Integration over the electronic densities is carried out by a numerical procedure slightly modified from that of Gavezzotti. We calculate the interaction energies for four molecular dimers at a variety of computational levels. The results are analyzed and compared to the best estimated values available in the literature.