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 kJ mol 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.