Transferability of atomic multipoles in amino acids and peptides for various density partitions

Multipole expansion of electron density distribution is an efficient tool for evaluating the energy of interactions in molecules. In as much as atoms in macromolecules such as peptides are modeled with certain types of atoms derived from small organic molecules, investigating transferability of atomic multipoles for various partitions of molecular electron density is an important issue. In this study, multipole moments up to hexadecapoles for types of atoms present in selected amino acids, as well as di- and tripeptides composed of these amino acids, are computed using three density partitions: Hansen-Coppens aspherical pseudoatoms formalism, Hirshfeld's stockholder partition, and Bader's atoms in molecules theory. Electron density of relevant molecules is derived in a procedure including molecular wave function ab initio calculation for isolated molecules in geometry from X-ray measurements, calculation of theoretical structure factors for molecules put in a pseudocubic cell, and multipole refinement as in crystallographic data processing and computation of multipoles. The results were compared to calculations of multipole moments in AIM and in stockholder density partitions obtained directly from molecular wave functions. The presented comparison does not point unambiguously to any particular influence of multipole refinement on moments obtained from these two partitions. The advantage of stockholder partitioning in terms of transferability of atomic multipoles is affirmed. AIM and pseudoatoms provide slightly less transferable multipoles of lower order. Higher rank of multipole expansion reveals a transferability improvement in the case of AIM and meaningful deterioration for pseudoatoms.