Abstract:The previously reported exact potential and multipole moment (EP/MM) method for fast and precise evaluation of the intermolecular electrostatic interaction energies in molecular crystals using the pseudoatom representation of the electron density [Nguyen, Macchi & Volkov (2020), Acta Cryst. A76, 630-651] has been extended to the calculation of the electrostatic potential (ESP), electric field (EF) and electric field gradient (EFG) in an infinite crystal. The presented approach combines an efficient Ewald-type summation (ES) of atomic multipoles up to the hexadecapolar level in direct and reciprocal spaces with corrections for (i) the net polarization of the sample (the `surface term') due to a net dipole moment of the crystallographic unit cell (if present) and (ii) the short-range electron-density penetration effects. The rederived and reported closed-form expressions for all terms in the ES algorithm have been augmented by the expressions for the surface term available in the literature [Stenhammar, Trulsson & Linse (2011), J. Chem. Phys. 134, 224104] and the exact potential expressions reported in a previous study [Volkov, King, Coppens & Farrugia (2006), Acta Cryst. A62, 400-408]. The resulting algorithm, coded using Fortran in the XDPROP module of the software package XD, was tested on several small molecular crystal systems (formamide, benzene, L-dopa, paracetamol, amino acids etc.) and compared with a series of EP/MM-based direct-space summations (DS) performed within a certain number of unit cells generated along both the positive and negative crystallographic directions. The EP/MM-based ES technique allows for a noticeably more precise determination of the EF and EFG and significantly better precision of the evaluated ESP when compared with the DS calculations, even when the latter include contributions from an array of symmetry-equivalent atoms generated within four additional unit cells along each crystallographic direction. In terms of computational performance, the ES/EP/MM method is significantly faster than the DS calculations performed within the extended unit-cell limits but trails the DS calculations within the reduced summation ranges. Nonetheless, the described EP/MM-based ES algorithm is superior to the direct-space summations as it does not require the user to monitor continuously the convergence of the evaluated properties as a function of the summation limits and offers a better precision-performance balance.

Atomic charges derived from semiempirical methods

DDEC6: A Method for Computing Even-Tempered Net Atomic Charges in Periodic and Nonperiodic Materials

Atomic charges for modeling metal-organic frameworks: Why and how

Accurate conformation-dependent molecular electrostatic potentials for high-throughput in silico drug discovery

On the calculation of the electrostatic potential, electric field and electric field gradient from the aspherical pseudoatom model. II. Evaluation of the properties in an infinite crystal

Benchmark Study of Density Functional Theory Methods in Geometry Optimization of Transition Metal-Dinitrogen Complexes

Accuracy of Charge Densities in Electronic Structure Calculations

On the sensitivity of computed partial charges toward basis set and (exchange‐)correlation treatment

Expanding the Scope of Density Derived Electrostatic and Chemical Charge Partitioning to Thousands of Atoms

On the sensitivity of computed partial charges toward basis set and (exchange-)correlation treatment

Development of Nonbonded Models for Metal Cations Using the Electronic Continuum Correction

A Hessian-Based Assessment of Atomic Forces for Training Machine Learning Interatomic Potentials

On combining Thole's induced point dipole model with fixed charge distributions in molecular mechanics force fields

Kernel-based Minimal Distributed Charges: A Conformationally Dependent ESP-Model for Molecular Simulations

Charge equilibration model with shielded long-range Coulomb for reactive molecular dynamics simulations

New-Generation Electron-Propagator Methods for Molecular Electron-Binding Energies

Atomic effective potentials for starting molecular electronic structure calculations

Using molecular similarity to construct accurate semiempirical electron structure theories

Toward a new approach for determination of solute's charge distribution to analyze interatomic electrostatic interactions in quantum mechanical/molecular mechanical simulations

Dipole moments from atomic-number-dependent potentials in analytic density-functional theory

Implementation and Optimization of the Embedded Cluster Reference Interaction Site Model with Atomic Charges