Atomic properties of selected biomolecules. Part 1. The interpretation of atomic integration errors

Abstract

Reliable atomic properties can be obtained via the theory of "Atoms in Molecules" (AIM) via integration over a finite volume. These integrations are challenging because of the variety and complexity of the shape of the AIM atoms. In practice the integration of a large number of atoms (100-1000, sampled from many molecules) yields integration errors L(Ω) of varying magnitude. We prove that it is impossible to predict the size of an angular Gauss-Legendre grid (outside the β sphere) that guarantees a pre-set error. Hence it is incorrect to assume that a large grid (~23 000 angular grid points) will automatically yield a low L(Ω) value. The erratic relationship between the integration error and the grid size prompts a statistical interpretation of atomic integration, at a purely practical level. More importantly we have investigated the relationship between L(Ω) and seven atomic properties which include volume, energy, and the magnitudes of five electrostatic multipole moments. The electronic population (N(Ω)) and the volume (v(Ω)) of carbon is linearly correlated with L(Ω), enabling the interpolation or extrapolation of N(Ω) and v(Ω). Other properties of carbon and other atoms (N, O, and S) yield low correlation coefficients but occasionally trends can be observed. For example, we find that some properties are systematically underestimated if L(Ω) is negative. This work has led to an estimate of safe error bars of atomic properties for atoms occurring in biological molecules with reasonably sized integration grids. The most stable properties were found to be the energy and the population. Finally, we have observed that the influence of the grid orientation is less if L(Ω) is small, and that population and energy are the least affected.Key words: electron density, topology, atoms in molecules, atomic properties, amino acids.