Chemical structure and bonding are key features and concepts in chemical systems which are used in deriving structure-property relationships, and hence in predicting physical and chemical properties of compounds. Even though the contemporary high standards in determination, using both theoretical methods and experimental techniques, questions of chemical bonds as well as their evolution along a reaction pathway are still highly controversial. This paper presents a working methodology to determine the structure and chemical reactivity based on the quantum chemical topology analysis.QTAIM and ELF frameworks, based on the topological analysis of the electron density and the electron localization function, respectively, have been used. We have selected two examples studied by the present approach, to show its potential: (i) QTAIM study on the α -Ag 2 WO 4 , for the simulation of Ag nucleation and formation on α -Ag 2 WO 4 provoked in this crystal by the electron-beam irradiation. (ii) An ELF and Thom´s catastrophe theory study for the reaction pathway associated with the decomposition of stable planar hypercoordinate carbon species, CN 3 Mg 3 + .3Chemistry is the science of substances: their structure, their properties, and the reactions that change them into other substances, as Linus These procedures are hugely successful and they can be considered as an energeticsdriven approach to computational theoretical and chemistry. Then, this research field constitutes a central application of quantum chemistry to the study of stationary points of PESs [53] and the pathways connecting them. The shape of an adiabatic PES is conventionally regarded as a function of the nuclear skeleton, ignoring the wealth of information that can be provided by the total electronic charge density distribution ρ(r).The driving force for the nuclear motion is the potential, which arises from In this context, the electronic structure of a molecule is described in real space terms using topological analysis. Beyond the well-known approach of the QTAIM, which relies on the properties of the electron density ρ(r) when atoms interact, topological analysis of different scalar functions can be employed, such as the electron localization function (ELF) method [90][91][92][93]. Originally, ELF can be developed based on the conditional pair density [90] or can be derived from the kinetic energy density [94].In the latter context, ELF is seen as the scaled difference between the positive kinetic energy density and the von Weizsäcker term [95], whereby the scaling function is the 8 kinetic energy density of the homogeneous electron gas (Thomas-Fermi term) [96,97]. The pictures of a molecule as drawn by the QTAIM and ELF analysis are complementary. Thus, the QTAIM analysis is based on the topology of the electron density ρ(r), whereas the ELF analysis is based on the topology of the electron localization function using the conditional probability for the same spin pairs which is closely related to the local excess of kinetic energy due to the Paul...