Curly arrows meet electron density transfers in chemical reaction mechanisms: from electron localization function (ELF) analysis to valence-shell electron-pair repulsion (VSEPR) inspired interpretation

Probing the electron density transfers during a chemical reaction can provide important insights, making possible to understand and control chemical reactions. This aim has required extensions of the relationships between the traditional chemical concepts and the quantum mechanical ones. The present work examines the detailed chemical insights that have been generated through 100 years of work worldwide on G. N. Lewis's ground breaking paper on The Atom and the Molecule (Lewis, G. N. The Atom and the Molecule, J. Am. Chem. Soc. 1916, 38, 762–785), with a focus on how the determination of reaction mechanisms can be reached applying the bonding evolution theory (BET), emphasizing how curly arrows meet electron density transfers in chemical reaction mechanisms and how the Lewis structure can be recovered. BET that combines the topological analysis of the electron localization function (ELF) and Thom's catastrophe theory (CT) provides a powerful tool providing insight into molecular mechanisms of chemical rearrangements. In agreement with physical laws and quantum theoretical insights, BET can be considered as an appropriate tool to tackle chemical reactivity with a wide range of possible applications. Likewise, the present approach retrieves the classical curly arrows used to describe the rearrangements of chemical bonds for a given reaction mechanism, providing detailed physical grounds for this type of representation. The ideas underlying the valence-shell-electron pair-repulsion (VSEPR) model applied to non-equilibrium geometries provide simple chemical explanations of density transfers. For a given geometry around a central atom, the arrangement of the electronic domain may comply or not with the VSEPR rules according with the valence shell population of the considered atom. A deformation yields arrangements which are either VSEPR defective (at least a domain is missing to match the VSEPR arrangement corresponding to the geometry of the ligands), VSEPR compliant or pseudo VSEPR when the position of bonding and non-bonding domains are interchanged. VSEPR defective arrangements increase the electrophilic character of the site whereas the VSEPR compliant arrangements anticipate the formation of a new covalent bond. The frequencies of the normal modes which account for the reaction coordinate provide additional information on the succession of the density transfers. This simple model is shown to yield results in very good agreement with those obtained by BET.We wish to thank Professors R. J. Gillespie, Henry H. Rzepa and Patrick Chaquin and L. R. Domingo for stimulating discussions and the referees for their very constructive comment

Reaction of atomic hydrogen with formic acid

Peer reviewe

Describing the molecular mechanism of organic reactions by using topological analysis of electronic localization function

Here, we provide an essay on the analysis of the reaction mechanism at the molecular level; in particular, the evolution of the electron pair, as it is provided by the ELF, is used to decribe the reaction pathway. Then, the reaction mechanism is determined by the topological changes of the ELF gradient field along a series of structural stability domains. From this analysis, concepts such as bond breaking/forming processes, formation/annihilation of lone pairs and other electron pair rearrangements arise naturally along the reaction progress simply in terms of the different ways of pairing up the electrons. To visualize these results some organic reaction mechanisms (the thermal ring aperture of cyclobutene and cyclohexa-1,3-diene) have been selected, indicating both the generality and utility of this type of analysis

How a Quantum Chemical Topology Analysis Enables Prediction of Electron Density Transfers in Chemical Reactions. The Degenerated Cope Rearrangement of Semibullvalene

Recent works on the reaction mechanism for the degenerated Cope rearrangement (DCR) of semibullvalene (SBV) in the ground state prompted us to investigate this complex rearrangement in order to assign experimentally observed contrast features in the simulated electron distribution. We present a joint use of the electron localization function (ELF) and Thom's catastrophe theory (CT) as a powerful tool to analyze the electron density transfers along the DCR. The progress of the reaction is monitored by the structural stability domains of the topology of ELF, while the change between them is controlled by turning points derived from CT. The ELF topological analysis shows that the DCR of SBV corresponds to asynchronous electron density rearrangement taking place in three consecutive stages. We show how the pictures anticipated by drawing Lewis structures of the rearrangement correlate with the experimental data and time-dependent quantum description of the process

The nature of inter- and intramolecular interactions in F2OXe…HX (X= F, Cl, Br, I) complexes

Electronic structure of the XeOF2 molecule and its two complexes with HX (X= F, Cl, Br, I) molecules have been studied in the gas phase using quantum chemical topology methods: topological analysis of electron localization function (ELF), electron density, ρ(r), reduced gradient of electron density |RDG(r)| in real space, and symmetry adapted perturbation theory (SAPT) in the Hilbert space. The wave function has been approximated by the MP2 and DFT methods, using APF-D, B3LYP, M062X, and B2PLYP functionals, with the dispersion correction as proposed by Grimme (GD3). For the Xe-F and Xe=O bonds in the isolated XeOF2 molecule, the bonding ELF-localization basins have not been observed. According to the ELF results, these interactions are not of covalent nature with shared electron density. There are two stable F2OXe…HF complexes. The first one is stabilized by the F-H…F and Xe…F interactions (type I) and the second by the F-H…O hydrogen bond (type II). The SAPT analysis confirms the electrostatic term, Eelst (1) and the induction energy, Eind (2) to be the major contributors to stabilizing both types of complexes.peerReviewe

Following the Molecular Mechanism for the NH3 + LiH → LiNH2 + H2 Chemical Reaction: A Study Based on the Joint Use of the Quantum Theory of Atoms in Molecules (QTAIM) and Noncovalent Interaction (NCI) Index

The molecular mechanism for the NH3 + LiH → LiNH2 + H2 reaction has been elucidated by the combined use of quantum theory of atoms in molecules (QTAIM) and noncovalent interactions (NCI) index. The topology of the electron density, obtained by QTAIM/NCI, is able to identify the evolution of strong and weak interactions, recovering the bonding patterns along the reaction pathway. Thus, the combination of these two techniques is a useful and powerful tool in the study of chemical events, providing new strategies to understand and visualize the molecular mechanisms of chemical rearrangements. Also, for the first time, the topology of the reduced density gradient has been analyzed, taking into account saddle points for the construction of bifurcation trees. This approach has demonstrated the ability of NCI to account for delocalized interactions, very often characteristic of transitions states

New insight into the electronic structure of iron(IV)-oxo porphyrin compound I. A quantum chemical topological analysis

The electronic structure of iron-oxo porphyrin π-cation radical complex Por·+FeIV[DOUBLE BOND]O (S[BOND]H) has been studied for doublet and quartet electronic states by means of two methods of the quantum chemical topology analysis: electron localization function (ELF) η(r) and electron density ρ(r). The formation of this complex leads to essential perturbation of the topological structure of the carbon–carbon bonds in porphyrin moiety. The double C[DOUBLE BOND]C bonds in the pyrrole anion subunits, represented by pair of bonding disynaptic basins Vi=1,2(C,C) in isolated porphyrin, are replaced by single attractor V(C,C)i=1–20 after complexation with the Fe cation. The iron–nitrogen bonds are covalent dative bonds, N→Fe, described by the disynaptic bonding basins V(Fe,N)i=1–4, where electron density is almost formed by the lone pairs of the N atoms. The nature of the iron–oxygen bond predicted by the ELF topological analysis, shows a main contribution of the electrostatic interaction, Feδ+···Oδ−, as long as no attractors between the C(Fe) and C(O) core basins were found, although there are common surfaces between the iron and oxygen basines and coupling between iron and oxygen lone pairs, that could be interpreted as a charge-shift bond. The Fe[BOND]S bond, characterized by the disynaptic bonding basin V(Fe,S), is partially a dative bond with the lone pair donated from sulfur atom. The change of electronic state from the doublet (M = 2) to quartet (M = 4) leads to reorganization of spin polarization, which is observed only for the porphyrin skeleton (−0.43e to 0.50e) and S[BOND]H bond (−0.55e to 0.52e)