An electron density (ED)-based methodology is developed for the automatic identification of intermolecular interactions using pro-molecular density. The expression of the ED gradient in terms of atomic components furnishes the basis for the Independent Gradient Model (IGM). This model leads to a density reference for non interacting atoms/fragments where the atomic densities are added whilst their interaction turns off. Founded on this ED reference function that features an exponential decay also in interference regions, IGM model provides a way to identify and quantify the net ED gradient attenuation due to interactions. Using an intra/inter uncoupling scheme, a descriptor (δginter) is then derived that uniquely defines intermolecular interaction regions. An attractive feature of the IGM methodology is to provide a workflow that automatically generates data composed solely of intermolecular interactions for drawing the corresponding 3D isosurface representations.
We hereby introduce the atomic degree of interaction (DOI), a new concept rooted in the electron density-based independent gradient model (IGM). Capturing any manifestation of electron density sharing around an atom, including covalent and non-covalent situations, this index reflects the attachment strength of an atom to its molecular neighbourhood. It is shown to be very sensitive to the local chemical environment of the atom. No significant correlation could be found between the atomic DOI and various other atomic properties, making this index a specific source of information. However, examining the simple H2 + H reacting system, a strong connection has been established between this electron density-based index and the scalar reaction path curvature, the cornerstone of the benchmark unified reaction valley approach (URVA). We observe that reaction path curvature peaks appear when atoms experience an acceleration phase of electron density sharing during the reaction, detected by peaks of the DOI second derivative either in the forward or reverse direction. This new IGM-DOI tool is only in its early stages, but it opens the way to an atomic-level interpretation of reaction phases. More generally, the IGM-DOI tool may also serve as an atomic probe of electronic structure changes of a molecule under the influence of physicochemical perturbations.
1 January 2014 in Computational and Theoretical Chemistry