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.

A modified electronegativity equalization method for fast and accurate calculation of atomic charges in large biological molecules.

Two Improved Electronegativity Equalization Methods for Charge Distribution in Large Scale Non-Uniform System

Σπ Scheme in Atom‐bond Electronegativity Equalization Method

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

Development of an Effective Polarizable Bond Method for Biomolecular Simulation.

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

The ENUF method—Ewald summation based on nonuniform fast Fourier transform: Implementation, parallelization, and application

An optimized charge penetration model for use with the AMOEBA force field.

The ENUF Method -- Ewald Summation based on Non-Uniform Fast Fourier Transform: Implementation, Parallelization, and Application

Accurate Calculation of Relative Binding Free Energies Between Ligands with Different Net Charges

Extended Energy Divide-and-Conquer Method Based on Charge Conservation.

AFNMR: Automated Fragmentation Quantum Mechanical Calculation of NMR Chemical Shifts for Biomolecules

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

Accelerated Quantum Mechanics/Molecular Mechanics Simulations Via Neural Networks Incorporated with Mechanical Embedding Scheme.

Advanced Charge Extended Hückel (CEH) Model and a Consistent Adaptive Minimal Basis Set for the Elements Z=1-103

Developing Differentiable Long-Range Force Field for Proteins with E(3) Neural Network Predicted Asymptotic Parameters

Density-gradient-corrected embedded atom method

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

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

Scalable Indirect Free Energy Method Applied to Divalent Cation-Metalloprotein Binding