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 (δg) 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.

Extraction of the chemical interaction signature from local descriptors based on electron density (ED) is still a fruitful field of development in chemical interpretation. In a previous work that used promolecular ED (frozen ED), the new descriptor, δg , was defined. It represents the difference between a virtual upper limit of the ED gradient (∇ρIGM , IGM=independent gradient model) that represents a noninteracting system and the true ED gradient (∇ρ ). It can be seen as a measure of electron sharing brought by ED contragradience. A compelling feature of this model is to provide an automatic workflow that extracts the signature of interactions between selected groups of atoms. As with the noncovalent interaction (NCI) approach, it provides chemists with a visual understanding of the interactions present in chemical systems. ∇ρIGM is achieved simply by using absolute values upon summing the individual gradient contributions that make up the total ED gradient. Hereby, we extend this model to relaxed ED calculated from a wave function. To this end, we formulated gradient-based partitioning (GBP) to assess the contribution of each orbital to the total ED gradient. We highlight these new possibilities across two prototypical examples of organic chemistry: the unconventional hexamethylbenzene dication, with a hexa-coordinated carbon atom, and β-thioaminoacrolein. It will be shown how a bond-by-bond picture can be obtained from a wave function, which opens the way to monitor specific interactions along reaction paths.

The covalent chemical bond is intimately linked to electron sharing between atoms.The recent Independent Gradient Model (IGM) and its δg descriptor provide a way to quantify locally this electron density interpenetration from wavefunction calculations.Each bond has its own IGM-δg pair signature. The present work establishes for the first time a strong link between this bond signature and the physically grounded bond force constant concept. Analyzing a large set of compounds and bonds, the Intrinsic Bond Strength Index (IBSI) emerges from the IGM formulation. Our study shows that the IBSI does not belong to the class of conventional bond orders (like Mulliken, Wiberg, Mayer, delocalization index or ELF), but is rather a new index complementary thereto, related to the bond force. A fundamental outcome of this research is a novel index allowing to range all two-centre chemical bonds by their intrinsic strength in molecular situation. We believe that the IBSI is a powerful and robust tool for interpretation accessible to a wide community of chemists (organic, inorganic chemistry, including transition metal complexes and reaction mechanisms).