Following the pioneering investigations of Bader on the topology of molecular electron density, the topology analysis of its sister field viz. molecular electrostatic potential (MESP) was taken up by the authors’ groups. Through these studies, MESP topology emerged as a powerful tool for exploring molecular bonding and reactivity patterns. The MESP topology features are mapped in terms of its critical points (CPs), such as bond critical points (BCPs), while the minima identify electron-rich locations, such as lone pairs and π-bonds. The gradient paths of MESP vividly bring out the atoms-in-molecule picture of neutral molecules and anions. The MESP-based characterization of a molecule in terms of electron-rich and -deficient regions provides a robust prediction about its interaction with other molecules. This leads to a clear picture of molecular aggregation, hydrogen bonding, lone pair–π interactions, π-conjugation, aromaticity and reaction mechanisms. This review summarizes the contributions of the authors’ groups over the last three decades and those of the other active groups towards understanding chemical bonding, molecular recognition, and reactivity through topology analysis of MESP.

Hydrogen, halogen, and dihydrogen bonds in weak, medium and strong regimes (<1 to ∼60 kcal/mol) have been investigated for several intermolecular donor–acceptor (D-A) complexes at ab initio MP4//MP2 method coupled with atoms-in-molecules and molecular electrostatic potential (MESP) approaches. Electron density ρ at bond critical point correlates well with interaction energy (E_nb) for each homogeneous sample of complexes, but its applicability to the entire set of complexes is not satisfactory. Analysis of MESP minimum (V_min) and MESP at the nuclei (V_n) shows that in all D-A complexes, MESP of A becomes more negative and that of D becomes less negative suggesting donation of electrons from D to A leading to electron donor–acceptor (eDA) interaction between A and D. MESP based parameter ΔΔV_n measures donor–acceptor strength of the eDA interactions as it shows a good linear correlation with E_nb for all D-A complexes (R² = 0.976) except the strongly bound bridged structures. The bridged structures are classified as donor–acceptor–donor complexes. MESP provides a clear evidence for hydrogen, halogen, and dihydrogen bond formation and defines them as eDA interactions in which hydrogen acts as electron acceptor in hydrogen and dihydrogen bonds while halogen acts as electron acceptor in halogen bonds.

An electrostatics-based approach has been proposed for probing the weak interactions between lone pair containing molecules and π deficient molecular systems. For electron-rich molecules, the negative minima in molecular electrostatic potential (MESP) topography give the location of electron localization and the MESP value at the minimum (V_min) quantifies the electron-rich character of that region. Interactive behavior of a lone pair bearing molecule with electron deficient π-systems, such as hexafluorobenzene, 1,3,5-trinitrobenzene, 2,4,6-trifluoro-1,3,5-triazine and 1,2,4,5-tetracyanobenzene explored within DFT brings out good correlation of the lone pair–π interaction energy (E_int) with the V_min value of the electron-rich system. Such interaction is found to be portrayed well with the Electrostatic Potential for Intermolecular Complexation (EPIC) model. On the basis of the precise location of MESP minimum, a prediction for the orientation of a lone pair bearing molecule with an electron deficient π-system is possible in the majority of the cases studied.