Darío J. R. Duarte

Topological analysis of aromatic halogen/hydrogen bonds by electron charge density and electrostatic potentials.

AbstractIn this work, the intermolecular distribution of the electronic charge density in the aromatic hydrogen/halogen bonds is studied within the framework of the atoms in molecules (AIM) theory and the molecular electrostatic potentials (MEP) analysis. The study is carried out in nine complexes formed between benzene and simple lineal molecules, where hydrogen, fluorine and chlorine atoms act as bridge atoms. All the results are obtained at MP2 level theory using cc-pVTZ basis set. Attention is focused on topological features observed at the intermolecular region such as bond, ring and cage critical points of the electron density, as well as the bond path, the gradient of the density maps, molecular graphs and interatomic surfaces. The strength of the interaction increases in the following order: F⋅⋅⋅π < Cl⋅⋅⋅π < H⋅⋅⋅π. Our results show that the fluorine atom has the capability to interact with the π−cloud to form an aromatic halogen bond, as long as the donor group is highly electron withdrawing. The Laplacian topology allows us to state that the halogen atoms can act as nucleophiles as well as electrophiles, showing clearly their dual character. FigureA study by electron charge density and electrostatic potentials about aromatic halogen/hydrogen bonds.

Double hole-lump interaction between halogen atoms

In this paper a theoretical study has been carried out to investigate the nature of the unusual halogen-halogen contacts in the complexes R-X···X-R (with R = -H, -Cl, -F and X = Cl, Br, I). AIM, NBO, and MEP analyses have been used to characterize X···X interactions. Formation of the unusual X···X interactions leads to a significant increase of electron charge density in the bonding region between the two halogen atoms. The geometry and stability of these complexes is mainly due to electrostatic interactions lump(X1) → hole(X2) and lump(X2) → hole(X1) [or equivalently [VS,min(X1) → VS,max(X2) and VS,min(X2) → VS,max(X1)] and the charge transfers LP(X1) → σ*(R-X2) and LP(X2) → σ*(R-X1). In other words, these findings suggest that the electrostatic interactions and the charge transfer play a substantial role in determining the optimal geometry of these complexes, as in conventional halogen bonds, even though the dispersion term is the most important attractive term for all the complexes studied here, save one.

Is the decrease of the total electron energy density a covalence indicator in hydrogen and halogen bonds

In this work, halogen bonding (XB) and hydrogen bonding (HB) complexes were studied with the aim of analyzing the variation of the total electronic energy density H(rb) with the interaction strengthening. The calculations were performed at the MP2/6−311++G(2d,2p) level of approximation. To explain the nature of such interactions, the atoms in molecules theory (AIM) in conjunction with reduced variational space self-consistent field (RVS) energy decomposition analysis were carried out. Based on the local virial theorem, an equation to decompose the total electronic energy density H(rb) in two energy densities, (−G(rb)) and 1/4∇2ρ(rb), was derived. These energy densities were linked with the RVS interaction energy components. Through the connection between both decomposition schemes, it was possible to conclude that the decrease in H(rb) with the interaction strengthening observed in the HB as well as the XB complexes, is mainly due to the increase in the attractive electrostatic part of the interaction energy and in lesser extent to the increase in its covalent character, as is commonly considered.

Computational study on the vinyl azide decomposition.

The decomposition mechanism of vinyl azide (CH2CHN3) has been studied by calculations of the electronic structure. In addition, a study based on the topology of the electron charge density distribution and its Laplacian function, within the Quantum Theory of Atoms in Molecules (QTAIM), has been carried out with the aim of comprehending the electron redistribution mechanisms that take place in the formation of vinyl nitrenes. The electronic structure calculations reveal that the decomposition of the s-cis conformer of vinyl azide leads to the formation of ketenimine through a single-step conversion, s-cis-CH2CHN3 → CH2CNH + N2, while the conversion of the s-trans conformer to acetonitrile occurs in two steps, s-trans-CH2CHN3 → cyc-CH2NCH + N2 → CH3CN + N2. The topological analysis of the L(r) function reveals that triplet vinyl nitrene has one lone pair on the valence shell charge concentration (VSCC) of nitrogen and thus could act as a monodentate Lewis base, while singlet vinyl nitrene has two lone pairs on the VSCC of nitrogen and thus could act as a bidentate Lewis base.

Physical meaning of the QTAIM topological parameters in hydrogen bonding.

This work examined the local topological parameters of charge density at the hydrogen bond (H-bond) critical points of a set of substituted formamide cyclic dimers and enolic tautomers. The analysis was performed not only on the total electron density of the hydrogen bonded complexes but also on the intermediate electron density differences derived from the Morokuma energy decomposition scheme. Through the connection between these intermediate electron density differences and the corresponding differences in topological parameters, the meaning of topological parameters variation due to hydrogen bonding (H-bonding) becomes evident. Thus, for example, we show in a plausible way that the potential energy density differences at the H-bond critical point properly describe the electrostatics of H-bonding, and local kinetic energy density differences account for the localization/delocalization degree of the electrons at that point. The results also support the idea that the total electronic energy density differences at the H-bond critical point describe the strength of the interaction rather than its covalent character as is commonly considered.

Adsorption of Alkenes on Acidic Zeolites. Theoretical Study Based on the Electron Charge Density

In the present work, experiments on electron density changes in the adsorption process of alkenes on acidic zeolites, in the framework of atoms in molecules theory (AIM), were carried out. Electron densities were obtained at MP2 and B3LYP levels using a 6-31++G(d,p) basis set. This study explores the energetic and the electron density redistributions associated with O-H...pi interactions. The main purpose of this work is to provide an answer to the following questions: (a) Which and how large are the changes induced on the molecular electron distribution by the formation of adsorbed alkenes? (b) Can a reasonable estimate of the adsorption energy of alkenes on the active site of zeolite be solely calculated from an analysis of the electron densities? We have used topological parameters to determine the strength and nature of the interactions in the active site of the zeolite. All the results derived from the electron density analysis show that the stabilization of the adsorbed alkenes follows the order isobutene > trans-2-butene congruent with 1-butene congruent with propene > ethene, reflecting the order of basicity of C=C bonds, i.e., (C(ter)=C(prim)) > (C(sec)=C(sec)) congruent with (C(prim)=C(sec)) > (C(prim)=C(prim)). In addition, we have found a useful set of topological parameters that are good for estimating the adsorption energy in adsorbed alkenes.

Multicenter (FX)n/NH3 Halogen Bonds (X = Cl, Br and n = 1–5). QTAIM Descriptors of the Strength of the X∙∙∙N Interaction

In the present work an in depth deep electronic study of multicenter XBs (FX)n/NH3 (X = Cl, Br and n = 1–5) is conducted. The ways in which X∙∙∙X lateral contacts affect the electrostatic or covalent nature of the X∙∙∙N interactions are explored at the CCSD(T)/aug-cc-pVTZ level and in the framework of the quantum theory of atoms in molecules (QTAIM). Calculations show that relatively strong XBs have been found with interaction energies lying between −41 and −90 kJ mol−1 for chlorine complexes, and between −56 and −113 kJ mol−1 for bromine complexes. QTAIM parameters reveal that in these complexes: (i) local (kinetics and potential) energy densities measure the ability that the system has to concentrate electron charge density at the intermolecular X∙∙∙N region; (ii) the delocalization indices [δ(A,B)] and the exchange contribution [VEX(X,N)] of the interacting quantum atoms (IQA) scheme, could constitute a quantitative measure of the covalence of these molecular interactions; (iii) both classical electrostatic and quantum exchange show high values, indicating that strong ionic and covalent contributions are not mutually exclusive.

NATURE OF Mδ+···δ+C-Oδ- INTERACTIONS IN METAL CARBONYLS. AN ELECTRONIC STUDY BASED ON THE TOPOLOGY OF THE ELECTRON CHARGE DENSITY DISTRIBUTION AND ITS LAPLACIAN FUNCTION

The nature of the metal-ligand interactions, in the [Ti(CO)6]2-, [V(CO)6]-, [Cr(CO)6], [Mn(CO)6]+, [Fe(CO)6]2+ and [Co(CO)6]3+ complexes has been studied by means of topological analyses of the electron charge density, using the Quantum Theory of Atoms in Molecules (QTAIM) and Electron Function Localization (ELF). The calculations were made using B3LYP method with the 6-311++G(2d,2p) basis set. The results show that the charge transferences (both σ-donation and π-backbonding) and the electrostatic interaction between the lone pair of C atom of the CO molecule and nucleus of the metal species play a key role in stabilizing of these metal complexes. Finally, we have found QTAIM parameters that explaining the nature of the Mδ+•••δ+C-Oδ- interactions in metal carbonyls.