Abstract pool for the NCI workshop
<<< Topological Approaches to Intermolecular Interactions workshop main page
Download the Book of Abstracts
Plenary lectures
Paul Popelier
Manchester Institute of Biotechnology (MIB) and School of Chemistry University of Manchester, Great Britain
Quantum Chemical Topology: on Bonds and Force Fields Watch the talk
Quantum Chemical Topology (QCT)[1,2] is an umbrella method that includes QTAIM[3,4] as a special case. The central idea of QCT is to partition through a gradient vector field, and apply the language and insights of dynamical systems. This talk has two distinct parts (that could be unified in future work).
The first part[5] discusses how to draw a molecule from a molecular wave function. The spatial distribution of atoms in a molecule in the form of chemical graphs is obtained for a set of molecules, using their corresponding domain-averaged exchange-correlation energies (Vxc). Conveniently, such energies are transferable (for 1,n interactions in saturated linear hydrocarbons) and can provide an accurate estimation of the covalent-like contribution between pairs of given interacting topological atoms A and B.
The second part focuses on the electrostatic interaction in a novel topological force field for biomolecular modeling. Topological atoms are boxes with a particular shape and a finite volume. If the coordinates change then the shapes of the atoms change too, as well as their multipole moments. This complex relationship is captured by a machine learning technique called kriging. Here I will explore how these ideas6 can be used to enhance the realism of the electrostatic energy, [7,8] and put polarisation and charge transfer on the same footing, without having a polarisation catastrophe.
[1] Popelier, P. L. A.; Brémond, É. A. G. Int.J.Quant.Chem. 2009, 109, 2542.
[2] Popelier, P. L. A. In Structure and Bonding. Intermolecular Forces and Clusters, Ed, D.J.Wales; Springer: Heidelberg, Germany, 2005; Vol. 115, p 1.
[3] Bader, R. F. W. Atoms in Molecules. A Quantum Theory.; Oxford Univ. Press, Great Britain, 1990.
[4] Popelier, P. L. A. Atoms in Molecules. An Introduction.; Pearson Education: London, Great Britain, 2000.
[5] Garcia-Revilla, M.; Francisco, E.; Popelier, P.L.A. ; Martin-Pendas, A. M. ChemPhysChem 2013, 14, 1211.
[6] Popelier, P. L. A. AIP Conf.Proc. 2012, 1456, 261.
[7] Mills, M. J. L.; Popelier, P. L. A. Theor.Chem.Acc. 2012, 131, 1137.
[8] Kandathil, S. M.; Fletcher, T. L.; Yuan, Y.; Knowles, J.; Popelier, P. L. A. J.Comput.Chem. 2013, in press, DOI: 10.1002/jcc.23333.
Frank De Proft
Eenheid Algemene Chemie (ALGC), Vrije Universiteit Brussel (VUB), Pleinlaan 2, B-1050 Brussels, Belgium
Combined use of DFT based Reactivity Indices and the Non-covalent Index in the Study of Intermolecular Interactions Watch the talk
Conceptual Density Functional Theory (sometimes also called DFT based reactivity theory or Chemical DFT) has proven to be an ideal framework for the introduction of chemical reactivity descriptors [1]. These indices are defined as response functions of the energy E of the system with respect to either the number of electrons N, the external potential v(r) or both. These definitions have afforded their non-empirical calculation and applications in many fields of chemistry have been performed, often combined with principles such as the electronegativity equalization principle of Sanderson and Pearson’s hard and soft acids and bases (HSAB) and maximum hardness principles. In this talk, we focus on recent studies in which a combined use of these reactivity indices and the Non-Covalent Interaction (NCI) Index [2] was carried out in order to scrutinize intermolecular interactions. In a first part, we focus on halogen bonding. Halogen bonds between the trifluoromethyl halides CF3Cl, CF3Br and CF3I, and dimethyl ether, dimethyl sulfide, trimethylamine and trimethyl phosphine were investigated using the HSAB concept with conceptual DFT reactivity indices, the Ziegler-Rauk type energy decomposition analysis [3], the Natural Orbital for Chemical Valence framework (NOCV) [4], and the NCI index. It is found that the relative importance of electrostatic and orbital (charge transfer) interactions in these halogen bonded complexes varies as a function of both the donor and acceptor molecules [5]. Hard and soft interactions were distinguished and characterized by atomic charges, electrophilicity and local softness indices. Dual descriptor plots indicate an orbital σ-hole on the halogen similar to the electrostatic σ-hole manifested in the molecular electrostatic potential. The characteristic signal found in the reduced density gradient versus electron density diagram corresponds to the non-covalent interaction between contact atoms in the NCI plots, which is the manifestation of halogen bonding within the NCI theory. The unexpected C–X bond strengthening observed in several cases was rationalized within the MO framework.
In a second part of the talk, we present recent results on the use of the above mentioned quantities in the study of metal-metal interactions [6].
[1] (a) R. G. Parr and W. Yang, Density Functional Theory of Atoms and Molecules, Oxford University Press, New York, 1989. (b) R. G. Parr and W. Yang, Ann. Rev. Phys. Chem. 46, 701 (1995). (c) H. Chermette, J. Comput. Chem. 20, 129 (1999). (d) P. Geerlings, F. De Proft and W. Langenaeker, Chem. Rev. 103, 1793 (2003). (e) P. W. Ayers, J. S. M. Anderson and L. J. Bartolotti, Int. J. Quant. Chem. 101, 520 (2005).
[2] (a) E. R. Johnson, S. Keinan, P. Mori-Sanchez, J. Contreras-Garcia, A. J. Cohen and W. T. Yang, Journal of the American Chemical Society 132, 6498 (2010). (b) J. Contreras-Garcia, E. R. Johnson, S. Keinan, R. Chaudret, J. P. Piquemal, D. N. Beratan and W. T. Yang, J. Chem. Theor. Comput. 7, 625 (2011).
[3] (a) F. M. Bickelhaupt and E. J. Baerends, Reviews in Computational Chemistry, 15, 1 (2000); (b) T. Ziegler and A. Rauk, Theoretica Chimica Acta 46, 1 (1977).
[4] M. P. Mitoraj, A. Michalak and T. Ziegler, J. Chem. Theor. Comput. 5, 962 (2009).
[5] B. Pinter, N. Nagels, W. A. Herrebout and F. De Proft, Chem. Eur. J. 19, 518 (2013).
[6] B. Pinter, L. Broeckaert, J. Turek, A. Růžička and F. De Proft, in preparation.
Weitao Yang
Department of Chemistry, Duke University, Durham, N.C. 27708, U.S.A.
Non-Covalent Interaction Analysis in Fluctuating Environments and Exchange-Correlation Energies from Response Properties Watch the talk
To apply the NCI analysis to fluctuating environments as in solution phase, we developed a new Averaged NonCovalent Interaction (i.e., aNCI) index along with a fluctuation index to characterize magnitude of interactions and fluctuations. We applied aNCI for various systems including solute-solvent and ligand-protein noncovalent interactions. For water and benzene molecules in aqueous solution, solvation structures and the specific hydrogen bond patterns were visualized clearly. For the Cl-+CH3Cl SN2 reaction in aqueous solution, charge reorganization influences over solvation structure along SN2 reaction were revealed. For ligand-protein systems, aNCI can recover several key fluctuating hydrogen bond patterns that have potential applications for drug design. Therefore, aNCI, as a complementary approach to the original NCI method, can extract and visualize noncovalent interactions from thermal noise in fluctuating environments.
In the second part of the talk, we will present an adiabatic connection to formulate the exchange-correlation energy in terms of response properties. This formulation of the exchange-correlation energy opens new channels for density functional approximations based on the many-body perturbation theory. We illustrate the potential of such approaches with an approximation based on the Random Phase Approximation. This resulting method has many highly desirable properties. It has minimal delocalization error with a nearly linear energy behavior for systems with fractional charges, describes van der Waals interactions similarly and thermodynamic properties significantly better than the conventional RPA, and eliminates static correlation error for single bond systems. Most significantly, it is the first known functional with an explicit and closed-form dependence on the occupied and unoccupied orbitals that captures the energy derivative discontinuity in strongly correlated systems.
Collaborators: Pan Wu, Robin Chaudret, Xiangqian Hu, Helen van Aggelen and Yang Yang
Manuel Yáñez
Departamento de Química. Facultad de Ciencias. Módulo 13. Universidad Autónoma de Madrid. Campus de Excelencia UAM-CSIC. Cantoblanco. 28049-Madrid. Spain.
Playing around with Beryllium Bonds Watch the talk
Beryllium bonds share many characteristics with conventional hydrogen bonds (HBs) [1], and as in HBs the proton donor may become significant distorter, the formation of Be bonds are also accompanied by a dramatic distortion of the beryllium derivative. This deformation plays a crucial role as far as the relative stability trends are concerned since it amounts to about a 30% of the total interaction energy, but more importantly significantly enhances the electron acceptor capacity of Be containing Lewis acid [2]. When these two closed-shell kinds of interactions are present in the same system, like in the complexes between BeX2 derivatives and imidazole dimer, clear cooperativity effects are observed, which result in a mutual reinforcement of both non-covalent interactions. This permits to use the Be bonds to modulate or tune the strength of inter- and intra-molecular HBs [3]. Indeed, for more complex clusters, as those involving water trimers, both cooperative and anti-cooperative effects are found [4]. More importantly Be bonds lead to a significant change in the intrinsic properties of the Lewis base participating in the bond, in particular on its intrinsic acidity which is dramatically enhanced [5]. The possibility of building up ditopic systems, in which the acid site is a BeX group opens interesting possibilities to design new polymeric materials.
[1] M. Yáñez., P. Sanz, O. Mó, I. Alkorta and J. Elguero, J. Chem Theor. Comput. 2009, 5, 2763-2771.
[2] A. Martín-Sómer, A. M. Lamsabhi, O. Mó and M. Yáñez, Comput. Theor. Chem. 2012, 998, 74-79.
[3] O. Mó, M. Yáñez, I. Alkorta and J. Elguero, J. Chem. Theory Comput. 2012, 8, 2293-2300.
[4] L. Albrecht, R. J. Boyd, O. Mo and M. Yanez, Phys. Chem. Chem. Phys. 2012, 14, 14540-14547.
[5] O. Mó, M. Yáñez, I. Alkorta, J. Elguero, J. Mol. Mod. 2013. dx.doi.org/10.1007/s00894-012-1682-y
Slawomir J. Grabowski
Faculty of Chemistry, University of the Basque Country UPV/EHU, and Donostia International Physics Center (DIPC),P.K. 1072, 20080 Donostia, Spain IKERBASQUE, Basque Foundation for Science, 48011 Bilbao, Spain
Non-covalent interactions: characteristics and mechanisms of formation - the topological approach Watch the talk
Numerous non-covalent interactions are characterized by the electron charge transfer from the Lewis base unit to the Lewis acid [1]. This is connected with the other processes reflected by the change of geometrical, energetic and topological parameters. For example, different characteristics of the hydrogen bond and various criteria of the existence of this interaction were discussed in the literature [2]. One can mention the topological criteria of Koch and Popelier [3,4]. On the other hand, the hydrogen bond mechanism was discussed in terms of NBO method [5]. Very recently it was found that the hydrogen bond, the halogen bond and other non-covalent interactions are steered by the same processes [6]. This is reflected by the same changes of parameters. For example, the A-H...B hydrogen bond formation is connected with the increase of the positive charge of H-atom and the decrease of its volume. The same changes are observed for X-halogen atom in the A-X...B halogen bond. Various similarities and differences between numerous non-covalent interactions may be discussed.
References
[1] Lipkowski, P.; Grabowski, S. J.; Leszczynski, J. J. Phys. Chem. A 2006, 110, 10296–10302.
[2] Grabowski, S.J. Chem.Rev. 2011, 11, 2597-2625.
[3] Koch, U.; Popelier, P.L.A. J.Phys.Chem. 1995, 99, 9747-9754.
[4] Popelier, P. Atoms in Molecules. An Introduction, Prentice Hall, Pearson Education Limited 2000.
[5] Alabugin, I.V.; Manoharan, M.; Peabody, S.; Weinhold, F. J.Am.Chem.Soc. 2003, 125, 5973-5987.
[6] Grabowski, S.J. Phys.Chem.Chem.Phys. accepted
A. Martín Pendás
Universidad de Oviedo. Spain
Decay rate of delocalization indices: towards a real space image of the insulating or metallic character of a material Watch the talk
Real space theories of the chemical bond have come of age, and are now grouped together under the umbrella name of Quantum Chemical Topology (QCT) [1], as originally proposed by Paul Popelier. Most of what is known in chemical bonding has been successfully reformulated using QCT techniques, that provide a wealth of quantities which are invariant under general orbital transformations. However, there is still no clear indication of what type of real space quantities, if any, might be used to discriminate between metals and insulators, in spite of several initial proposals, like those of Gatti and Silvi, or Contreras et. al. Recently, Baranov and Kohout have proposed that localization or delocalization indices might hold the clue [4]. However, their unrestricted Kohn-Sham approach is not devoid of possible criticisms. Here we address the problem from a molecular point of view, examining the analytic behavior of localization and delocalization indices in analytically solvable models and pointing towards the rate of decay of these quantities as a possible measure of long-range (metallic-like) or short-range (insulating-like) electron delocalization in materials [5].
[1] Paul L. A. Popelier, Struct. Bond. 2005, 115, 1.
[2] B. Silvi, C. Gatti, J. Phys. Chem. A, 2000, 104.
[3] J. Contreras-García, A. Martín Pendás, J. M. Recio, B. Silvi , J. Chem. Theory Comput. , 2009, 5, 164.
[4] A. Baranov, M. Kohout, Acta Cryst. A 2011, 67, C115.
[5] A. Martín Pendás, in preparation.
Oral contributions
Bertrand Castro,a Robin Chaudret,b Gino Ricci,c Gerhard Kretzschmar,d Volker Kraft,d Michael Kurz,d Kai Rossen,d and Odile Eisenstein,a
Analysis of weak interactions for asymmetric synthesis: the case of the selective hydrogenation of artemisinic acid by diazene
a Institut Charles Gerhardt, UMR 5253 CNRS Université Montpellier 2, France
b IFP Energie Nouvelles, Rueil-Malmaison, France
c Sanofi Chimie, 45 chemin de Meteline 04200 Sisteron, France
d Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, Building G 838, 65926 Frankfurt am Main, Germany
Artemesinate has been found very effective against malaria. Unfortunately its extraction from a Chinese medical herb Artemisia annua lead to poor and highly variable yields. Consequently, an efficient synthetic route to artemisinin is needed. While several total syntheses of artemisinin were known, none of them was suitable for the required large scale, low cost preparation. Sanofi in collaboration with Amyris Biotechnologies Inc have developed a very efficient synthesis of artemisine.[1] This reaction involves a selective hydrogenation reaction of artemesinic acid followed oxygenation by singlet O2. A very efficient and low cost large scale synthesis of dehydro artemesinic acid was obtained by using diazene (HN=NH), since the diastereosiomeric ratio of the hydrogenation was over 99%.[2,3]
The selectivity of the hydrogenation step is very surprising since no steric effects between the two reactants are present. In addition, hydrogenation by diazene occurs usually with very low energy barrier. This lead to another paradox since reactions with low energy barriers are not supposed to be very selective.
A computational study of this reaction was carried out with DFT calculations. The calculations reproduce almost quantitatively the experimental diastereoisomeric ratio. Several transition states leading to either the major or the minor products, associated to various conformations of the substrate, were identified and a single transition state is responsible for the production of most of the major product. Remarkably, the conformation of the substrate at this preferred transition state is that preferred for the isolated substrate. NCI and NCI-ELF are used to determine the network of attractive interactions that favor this unique conformation, which cannot be understood on the basis of steric effects. The key interaction exists in the reactant and not in the product, which rationalizes why this reaction, which has a reactant-like transition state, is so selective. 2D-NMR studies are underway at SANOFI to analyze the conformational preference of the reactant.
[1] a) http://www.youtube.com/watch?v=7xh8lFHEQDw . b) D. K. Ro, E. M. Paradise, M. Ouellet, K. J. Fisher, K. L. Newman, J. M. Ndungu, K. A. Ho, R. A. Eachus, T. S. Ham, J. Kirby, M. C. Chang, S. T. Withers, Y. Shiba, R. Sarpong, J. D. Keasling, Nature 2006, 440, 940. c) V. Hale, J. D. Keasling, N. Renninger, T. T. Diagana, Am. J. Trop. Med. Hyg. 2007, 77, 198. d) H. Tsuruta, C. J. Paddon, D. Eng, J. R. Lenihan, T. Horning, L. C. Anthony, R. Regentin, J. D. Keasling, N. S. Renninger, J. D. Newman, PloS ONE 2009, 4, e4489. J. Dhainaut, A. Dlubala, R. Guevel, A. Medard, G. Oddon, N. Raymond, J. Turconi, PTC W0 2011/026365 A1, International Publication Date 10.03.2011.
[2] M. P. Feth, K. Rossen, A. Burgard, Org. Proc. Res. & Development 2013, 17, 282–293.
[3] K. Kraft, G. Kretzschmar, K. Rossen, patent PTC WO 2011/030223 A2.
Piotr Kasprzycki, Piotr Ciąćka, Piotr Fita, Czesław Radzewicz
Institute of Experimental Physics, Faculty of Physics, University of Warsaw, ul. Hoza 69, 00-681 Warsaw, Poland
The correlation between the proton transfer time and electron density in the inner cavity of porphycenes
Numerical solutions of the Schrödinger equation can nowadays be obtained for reasonably large organic and organometallic molecules with useful accuracy using approximate ab initio or DFT methods. These approximate solutions are also available for porphycene which is an important chromophore [1]. Porphycene and its derivatives show useful properties for material sciences and are also considered excellent photosensitive media for applications in photodynamical therapy of cancer [2]. Several papers have discussed the energy structure and dynamic aspects of hydrogen transfer in this systems [4, 5]. We report on the results of NCI analysis - theoretical tool based on Bader approach [6] - performed for symmetrically and asymmetrically substituted porphycene derivatives. We have carried out DFT calculations at the level B3LYP/6-311(d,p) of theory. These calculations were followed by NCI analysis. Calculations showed a correlation between the value of electron density gradient localized between nitrogen and hydrogen atoms in the inner cavity of a molecule and the time of the proton transfer in the examined systems. Proton transfer times were obtained in a femtosecond transient absorption anisotropy experiment, which is a polarisation-sensitive variant of an ultrafast pump-probe technique. Double hydrogen transfer leads to the change of transition dipole moment directions - transition moment of chemically identical tautomerization product forms an angle of α = (72±2)◦ with that of initially excited form. This change is manifested as variation in transmitted light intensity measured in two orthogonal polarizations. Under appropriate conditions, rise/decay times of anisotropy reflect the kinetics of double hydrogen transfer. We show that quantum theory of atoms in molecules (QTAIM) based on Bader approach is useful not only for detection of covalent bond but also for hydrogen bonds and weak van der Waals interaction. NCI analysis provides the possibility to compare strength of hydrogen bonds in porphycene derivatives, which allows one to predict the proton transfer time for new molecules. Additionally, these approach have helped us to design structure of a derivative with the proton transfer time of the order of tens of picoseconds. Furthermore, substitution in 2,7 and in 10,19 position of porphycene can induce differentiation of hydrogen bonds. The effect is especially evident in the case of 10,19-di-methylporphycene (2MPc) and 2,7-di-t-buthylporphycene (DTBP). This behavior may be the result of asymmetric distribution of electron density in the center of the molecule. This fact can be significant for understanding the mechanism of double hydrogen transfer between inner nitrogen atoms.
References
[1] D. Sanchez-Garcia, J. L. Sessler, Chem. Soc. Rev., 37, 215, (2008).
[2] J. C. Stockert, M. Canete, A. Juarranz, A. Villanueva, R. W. Horobin, J. I. Borrell, J. Teixido, S. Nonell, Curr. Med. Chem., 14, 997, (2007) .
[4] Ł. Walewski, J. Waluk, B. Lesyng, J. Phys. Chem. A., 114, 2313–2318, (2010).
[5] M. Gil, J. Waluk J. Am. Chem. Soc. 5, 129, (2007).
[6] E. R. Johnson, S. Keinan, , J. Contreras-García, R. Chaudret, J-P Piquemal, D. Beratant, W. Yang, J. Chem. Theo. and Comp., 25, (2011).
Vincent Tognetti and Laurent Joubert
Normandy Univ., COBRA UMR 6014 & FR 3038, Mont Saint Aignan, 76821 cedex, France
QTAIM/DFT descriptors for the description of “weak” interactions
We present here part of our recent work, which aims to design new molecular and atomic descriptors by building a bridge between Density Functional Theory (DFT) and Bader’s Atomsin- molecules (AIM) theory,1 two frameworks that share the same primary ingredient: the electron density. A local point of view will first be presented, dealing with the so-called bond critical points (BCPs). New BCP properties will be investigated,2,3,4,5 based on reduced density gradients and on the variations rates of kinetic, exchange and correlation density energies. Interestingly, these quantities reveal particularly suited to the description of weak interactions, like hydrogen or agostic bonds. Then, we will discuss,6 using Pendás’ Interacting Quantum Atoms scheme,7 to what extent the existence or the absence of such BCPs can be rationalized in terms of integrated interatomic energies, enlightening the fundamental role of exchange channels.8 We will consequently present a fast and accurate protocol for the evaluation of the exchange interaction between two atoms,9 making another promising step in the pathway linking DFT and AIM.
References
[1] R. F. W. Bader, Atoms in Molecules: A Quantum Theory, Oxford University Press (1994).
[1] V. Tognetti, L. Joubert, P. Cortona and C. Adamo, J. Phys. Chem. A 113 12322-12327 (2009).
[3] V. Tognetti, L. Joubert and C. Adamo, J. Chem. Phys. 132 211101 (2010).
[4] V. Tognetti and L. Joubert, J. Phys. Chem. A 115 5505-5515 (2011).
[5] V. Tognetti, L. Joubert, R. Raucoules, T. De Bruin and C. Adamo, J. Phys. Chem. A 116 5472-5479 (2012).
[6] V. Tognetti and L. Joubert, J. Chem. Phys. 138 024102 (2013).
[7] M. A Blanco, A. Martín Pendás and E. Francisco, J. Chem. Theory Comput. 1 1096-1109 (2005).
[8] A. Martín Pendás, E. Francisco, M. A. Blanco and C. Gatti, Chem. Eur. J. 13 9362-9371 (2007).
[9] V. Tognetti and L. Joubert, in preparation.
Verònica Postils, Eduard Matito, Marc Garcia-Borràs, Miquel Solà and Josep M. Luis
Institut de Química Computacional i Catàlisi and Departament de Química, Campus de Montilivi s/n, University of Girona, Girona, 17071, Spain
Characterization and Identification of Electrides
The electrides[1-2] are intriguing chemical species with an electron not formally assigned to any atom. This situation is, however, completely different to that given in a metal where the electrons are delocalized between positively charged metal ions. The electron in an electride acts as a formal anion, which is bonded to positively charged species in the molecule. This particular feature of electrides prompts very particular chemical and physical properties: they are powerful reducing reagents, exhibit exalted electric linear and non-linear optical properties as well as a particular magnetic behavior.
In this work, we analyze the electronic structure and the identification of several electride structures by means of the Quantum Theory of Atoms in Molecules (QTAIM) and the Electron Localization Function (ELF) [3]. Our results show that these tools make possible the classification of candidate species as electrides or not. It was already proved that one could distinguish the electride behavior in insulating high-pressure forms of alkali metals from ELF analysis [4] but now we show that with QTAIM and ELF is possible to characterize the electride behavior in all sorts of molecules.
[1] Dye, J. L. Science, 2003, 301, 607-608.
[2] Dye, J. L. Science, 1990, 247, 663
[3] Matito, E.; Solà, M. Coord. Chem. Rev., 2009, 253, 647-665
[4] Marqués et al., Phys Rev. Lett. 2009, 103, 115501.
Mercedes Alonso, Paul Geerlings and Frank De Proft
Department of General Chemistry (ALGC), Free University of Brussels (VUB), Pleinlaan 2, Brussels, 1050, Belgium
Fine-tuning of the conformation of expanded porphyrins using conceptual DFT descriptors and non-covalent index
Density functional theory (DFT) is a powerful tool to examine a broad variety of structures and reactions involving increasingly larger molecules. Moreover, it is a well-suited theory for the introduction of chemical concepts, also known as reactivity indices [1]. Recently, we have shown that DFT calculations in combination with NCI analysis [2] are very useful in aiding the design of viable Möbius aromatic systems and molecular switches of expanded porphyrins [3-5]. These macrocycles can adopt a variety of intriguing structures, with Hückel, Möbius and twisted-Hückel topologies, which can be interconverted under certain conditions. The vast structural diversity exhibited by expanded porphyrins has led to diverse applications such as non-linear optical materials and molecular sensors, although the control of molecular topology is crucial for most of applications.
Conformational preferences and interconversion pathways in penta-, hexa- and heptaphyrins have been investigated using different hybrid functionals. The stability of Hückel and Möbius conformers depends on intramolecular hydrogen bonding, ring strain, aromaticity and steric effects of the meso-substituents. We proposed a set of descriptors to quantify independently the contribution of the different factors [3]. Two torsional parameters are used to quantify ring strain. Aromaticity has been quantified using energetic, structural, magnetic and reactivity criteria. Importantly, the relative hardness (a conceptual DFT descriptor) is found to be a good descriptor of the macrocyclic aromaticity in Hückel and Möbius porphyrinoids, whereas the dual descriptor is highly useful in predicting reactivity. Steric effects of the substituents are analysed in terms of the non-covalent interaction index (NCI). NCI is capable of distinguishing intramolecular hydrogen bonds, π- π interactions, CH···π and repulsive steric clashes, based on the peaks that appear in the reduced density gradient at low densities. While the relative stabilities of [4n + 2] expanded porphyrins are mainly governed by ring strain, intramolecular hydrogen-bonding and steric effects controls the stability of [4n] congeners. Appealingly, conformational control of pentaphyrins can be achieved by changing the type and number of meso-substituents [5].
References
[1] P. Geerlings and F. De Proft, and W. Langenaeker, Chem. Rev. 103, 1793 (2003).
[2] E. R. Johnson, S. Keinan, P. Mori-Sánchez, J. Contreras-García, A. J. Cohen and W. Yang J. Am. Chem. Soc. 132 6498-6506 (2010).
[3] M. Alonso, P. Geerlings, and F. De Proft, Chem. Eur. J. 18 10916-10928 (2012).
[4] M. Alonso, P. Geerlings, and F. De Proft, Chem. Eur. J. 19 1617-1628 (2013).
[5] M. Alonso, P. Geerlings, and F. De Proft, J. Org. Chem. accepted. DOI: 10.1021/jo4003823.
Hélène GERARD
Université Pierre et Marie Curie - Paris 06
Interactions within heterobimetallic complexes
Main group metal cations (we will report about lithium, aluminium and zinc), when associated to anions (including carbanions) in non-aqueous solvent form a variety of structures going from separated ion pairs to fully polymerized ladders. Among these arrangement, a four-membered ring alternating cations and anions have been shown of major interest in developping enantioselective catalytic processes when lithium is used as a sole cation.[1]
The report in litterature of X-Ray structures exhibiting a similar core but with Zn or Al for one of the cation partner [2] has open the way to the transposition of this process to other organometallic species, and thus of a large variety of alternative selectivities. In this connection, it would be tempting to decribe, on theoretical ground, the electronic properties of these hetero-bimetallic species in order to help design of new experimental processes.
Various methods for the theoretical analysis of the electronic structure of these species have been used. In good accordance with the experimental properties of these species, they yield a very contrasted picture depending not only of the system examined, but also of the modeling tool at stake. Bridging the gap between these different models of molecular interactions brings essential information to understand the complexity of the solution behaviors.
[1] “Catalytic enantioselective nucleophilic addition of organolithium derivatives: pitfalls and opportunities”, A. Harrison-Marchand, H. Gérard, J. Maddaluno, New J. Chem., 2012, 36, 2441–2446 and references therein.
[2] a. “Structure and Reaction Pathway of TMP-Zincate: Amido Base or Alkyl Base?”, M. Uchiyama, Y. Matsumoto, D. Nobuto, T. Furuyama, K. Yamaguchi, K. Morokuma*, J. Am. Chem. Soc., 2006, 128, 8748-8750. b. “An Aluminum Ate Base: Its Design, Structure, Function, and Reaction Mechanism” H. Naka, M. Uchiyama, Y. Matsumoto, A. E. H. Wheatley, M. McPartlin, J. V. Morey, Y. Kondo, J. Am. Chem. Sc. 2007, 129, 1921-1930.
Sílvia Simon, Laia Guillaumes
Institut de Química Computacional i Catàlisis (IQCC) and Departament de Química, Universitat de Girona, Campus de Montilivi, E-17071 Girona, Catalonia (Spain)
Delocalization Indices for Non-Covalent Interaction
Delocalization indices (DI) have been demonstrated to help of electronic structure of molecules, and more concretely in the analysis of HB nature.[1,2] Previous results pointed out that an electron-pair analysis in the framework of the AIM theory could be used to characterize hydrogen bonding in terms of the number of electrons shared between the atoms involved in the HB. As these indices are related to the polarity of the bond, they can be used to classify hydrogen bonds from weak to strong ones. Not only HB but also non-conventional hydrogen bond, eg. Dihydrogen bonded complexes, has been characterized by using DI.[3] Furthermore, these indices were also used in the study of Resonance Assisted Intramolecular Hydrogen Bond (RAHB). In that case, for a series of o-hydroxyaryl ketones, an intramolecular HB was formed giving an extra ring formation, finding out a good relation between DIs (between proton and proton acceptor atoms), and the stabilization energy when forming the HB.[4-6]
Different partition schemes can be defined to identify the atomic contributions to one property like DI. One of the partitions that has been gained more interest is to decompose the 3D space into atomic domains. Two main groups can be distinguished within this scheme, those with sharp boundaries between atomic domains (e.g. AIM, Voronio Cells, ELF) and those which assign a non-negative atomic weight functions to define the atomic basin in all the 3D space. The comparison of these 3D schemes will be also one of the points to discuss in the present work.
(1) X. Fradera, J. Poater, S. Simon, M. Duran, M. Sola, Theoretical Chemistry Accounts 108 (2002) 214.
(2) J. Poater, X. Fradera, M. Sola, M. Duran, S. Simon, Chemical Physics Letters 369 (2003) 248.
(3) D. Hugas, L. Guillaumes, M. Duran, S. Simon, Computational and Theoretical Chemistry 998 (2012) 113.
(4) M Palusiak, S. Simon, M. Sola, Journal of Organic Chemistry 71 (2006) 5241.
(5) M. Palusiak, S. Simon, M. Sola, Journal of Organic Chemistry 74 (2009) 2059.
(6) M. Palusiak, S. Simon, M. Sola, Chemical Physics 342 (2007) 43.
Arturo Espinosa
Departamento de Química Orgánica, Facultad de Química, Campus de Espinardo, Universidad de Murcia, 30100 Murcia, SPAIN.
Intra- and intermolecular NCIs balance in biscarbazole-trisurea (BC3U) receptors and their pyrophosphate complexes
As a part of recent studies on selective (optical, fluorescent or electrochemical) sensing of ions or neutral molecules, a number of derivatives of the novel biscarbazole-trisurea (BC3U) family of receptors have been prepared [1]. The prototypical peripherycally bis-p-nitrophenyl substituted term (Figure 1, left) shows a deeply coiled conformation in DMSO solution, according to the most stable computed [2] geometry obtained with the ORCA electronic structure package [3]. Addition of the biologically highly significant pyrophosphate anion selectively promotes remarkable experimental changes (optical, NMR shifts) that agree with the proposed computed (same level) structure (Figure 1, right). Complex formation requires a favourable energy balance between cleavage of several intramolecular NCIs in the receptor, involving both H-bonding and aromatic stacking (of either π,π- and T-types) and formation of new intermolecular (supramolecular) NCIs of mainly H-bonding type. All NCIs are displayed by means of colour-coded RDG (reduced density gradient) isosurfaces generated with the NCIplot program [4] and conveniently quantified by several typical bond-strength parameters.
Figure 1: Example of computed structures of a BC3U receptor and its pyrophosphate complex.
References
[1] G. Sánchez, A. Espinosa, D. Curiel, A. Tárraga, P. Molina, manuscript in preparation.
[2] COSMO(DMSO)/B3LYP-D3/def2-TZVP level of theory.
[3] ORCA An ab initio, DFT and semiempirical SCF-MO package. Written by F. Neese, Max Planck Institute for Bioinorganic Chemistry, D-45470 Mülheim/Ruhr, 2012. Version 2.9.1. Website: http://www.mpibac.mpg.de/bac/logins/neese/description.php. F. Neese, "The ORCA program system", WIREs Comput Mol Sci 2012, 2, 73–78.
[4] J. Contreras-García, E. R. Johnson, S. Keinan, R. Chaudret, J-P. Piquemal, D. N. Beratan, W. Yang. J. Chem. Theory Comput. 2011, 7, 625-632.
Kasper Mackeprang
Department of Chemistry, University of Copenhagen
Weak intramolecular OH···π hydrogen bonding in methallyl and allyl carbinol
The weak intramolecular hydrogen bond in methallyl carbinol and allyl carbinol has been identified by intracavity laser photoacoustic spectroscopy as an observed redshift of the third OH-stretching overtone transition of the most abundant conformer of each molecule. Furthermore, the vibrational overtone spectrum indicates that the intramolecular hydrogen bond in methallyl carbinol is stronger than that in allyl carbinol as the redshift corresponding to the former molecule is larger. We have used the non-covalent interaction index to characterize and determine the strength of the intramolecular hydrogen bond in both molecules. It was confirmed that the strength of the hydrogen bond in methallyl carbinol was slightly larger than that in allyl carbinol.
Marc Garcia-Borràs,a Sílvia Osuna,a Marcel Swart,a,b Josep M. Luis,a and Miquel Solàa
a Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Campus Montilivi, 17071 Girona, Spain
b Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, 08010 Barcelona, Spain
e-mail: marc.garcia@udg.edu
Intermolecular Dispersion Interaction Effects on the Diels-Alder Endohedral Metallofullerene Monoadduct stabilities. The Cp/Cp*-La@C2v-C82 example.
In the present work, we have studied the regioselectivity of the Diels-Alder (DA) addition of cyclopentadiene (Cp) and 1,2,3,4,5-pentamethylcyclopentadiene (Cp*) on La@C2v-C82 EMF at the B3LYP-(D)/TZP//BLYP-D/DZP level of theory [4]. Experimentally, the major product of the Cp* addition was characterized by X-ray crystallography [5]. By contrast, the product of the Cp DA reaction was still not known [6]. Moreover, both products present very different experimental stabilities being the Cp* adduct thermally much more stable than the Cp one, which rapidly decomposes into reactants at room conditions [5]. This different behavior suggests that the two products differ from each other.
Our results show that, surprisingly, the regioselectivity of both reactions mentioned above is equivalent, the same C-C bond of the EMF being the most favored attacked during the DA reaction. And, moreover, their thermal stability is practically the same with similar retro-reaction barriers. The different product stability is caused by the formation of the initial reactant complex due to the intermolecular dispersion interactions. In the Cp* case, in contrast to the Cp case, the formation of the adduct is largely favored,, being the intermolecular dispersion energy term much larger for the Cp*-EMF adduct than for the Cp-EMF system. This implies a rapid dissociation of the latter, and on the contrary a dynamic equilibrium between the initial complex and the final product for the Cp* diene. Thus, we have found that the key role played by the initial reactant complex finally determines the thermal stability of the final adducts in this particular cycloaddition reaction.
[1] S. Osuna, M. Swart, M. Solà, Phys. Chem. Chem. Phys. 2011, 13, 3561.
[2] M. Garcia-Borràs, S. Osuna, J. M. Luis, M. Swart, M. Solà, Chem. Eur. J. 2012, 18, 7141.
[3] M. Garcia-Borràs, A. Romero-Rivera, S. Osuna, J. M. Luis, M. Swart, M. Solà, S. J. Chem. Theory Comput. 2012, 8, 1671.
[4] M. Garcia-Borràs, M. Swart, J. M. Luis, M. Solà, Chem. Eur. J. 2013, 19, 4468.
[5] Y. Maeda, S. Sato, K. Inada, H. Nikawa, M. Yamada, N. Mizorogi, T. Hasegawa, T. Tsuchiya, T. Akasaka, T. Kato, Z. Slanina, S. Nagase, Chem. Eur. J. 2010, 16, 2193-2197.
[6] Y. Maeda, J. Miyashita, T. Hasegawa, T. Wakahara, T. Tsuchiya, T. Nakahodo, T. Akasaka, N. Mizorogi, K. Kobayashi, S. Nagase, T. Kato, N. Ban, H. Nakajima, Y. Watanabe, J. Am. Chem. Soc. 2005, 127, 12190-12191.
Eduard Matito, Ferran Feixas, Juan M. Barroso, Miquel Solà and Jesus M. Ugalde
Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Campus de Montilivi s/n, 17071 Girona, Catalonia, Spain
Multicenter Bond Index: a versatile tool to characterize electron delocalization and aromaticity.
Multicenter indices (MCI) [1-10] are used in a number of situations such as the analysis of conjugation and hyperconjugation effects [7], to identify agostic bonds, to account for electron distributions in molecules [4-5] or to study aromaticity in both organic [9] and all-metal compounds [10-11]. The formula of MCI involves the intuitive n-order central moment of the electron population [12-13], which depends on the n-order RDM. The calculation of the n-RDM is overwhelmingly expensive for correlated wavefunctions and, therefore, there are so few studies of MCI available in the literature. In this talk we review the applications of MCI and suggest a new approximation for the 3-RDM that provides more accurate MCI values than the approximations suggested thus far in the literature [15-16].
[1] Giambiagi, M. et al. A. P. Phys. Chem. Chem. Phys. 2 (2000) 3381.
[2] Bultinck, P et al. J. Phys. Org. Chem. 18 (2005) 706.
[3] Cioslowski, J. et al. J. Phys. Chem. A 111 (2007) 6521.
[4] Pendás, A. M. et al. Phys. Chem. Chem. Phys. 9 (2007) 1087.
[5] Pendás, A. M. et al. J. Chem. Phys. 127 (2007) 144103.
[6] Ponec, R. et al. J. Phys. Chem. A 101 (1997) 1738.
[7] Feixas, F. et al. J. Phys. Chem. A 115 (2011) 13104.
[8] Bultinck, P. Faraday Discuss. 135 (2007) 347.
[9] Feixas, F. et al. J. Comput. Chem. 29 (2008) 1543.
[10] Feixas, F. et al. Theor. Chem. Acc. 128 (2011) 419.
[11] Feixas, F. et al. J. Chem. Theory Comput. 6 (2010) 1118.
[12] Giambiagi, M. et al. Struct. Chem. 1 (1990) 423.
[13] Francisco, E. et al. J. Chem. Phys 126 (2007) 094102.
[14] Feixas, F. et al. in preparation.
[15] Colmenero, F. et al. Int. J. Quantum Chem. 51 (1994) 369.
[16] Mazziotti, D. A. Phys. Rev. A 57 (1998) 4219.
Alberto Otero-de-la-Roza and Erin Johnson
Chemistry and Chemical Biology, School of Natural Sciences, University of California, Merced, 5200 North Lake Road, Merced, California 95343, USA
Non-Covalent Interactions using the Exchange-Hole Dipole Moment
Conventional density functionals do not include the physics of London dispersion, which is necessary to model intermolecular interactions between non-polar species. In this talk, we present some applications of the exchange-hole dipole moment (XDM) dispersion functional[1,2], which we have implemented for production calculations of molecules [3,4] and solids[5,6]. The applications include: i) accurate prediction of sublimation enthalpies and crystal geometries, ii) enantiomeric excess amplification from lattice energies, iii) modelling of electride crystals, and iv) friction on graphite and carbonaceous materials.
[1] A. D. Becke and E. R. Johnson, J. Chem. Phys. 122, 154104 (2005)
[2] A. D. Becke and E. R. Johnson, J. Chem. Phys. 127, 154108 (2007)
[3] A. Otero-de-la Roza and E. R. Johnson, J. Chem. Phys. 138, 054103 (2013)
[4] A. Otero-de-la Roza and E. R. Johnson, J. Chem. Phys. 138, 204109 (2013)
[5] A. Otero-de-la-Roza and E. R. Johnson, J. Chem. Phys. 136, 174109 (2012)
[6] A. Otero-de-la-Roza and E. R. Johnson, J. Chem. Phys. 137, 054103 (2012)
Christine Lepetit,aJordi Poater,b Miquel Solà,bJulia Contreras-Garcia,cYves Canac,aDavit Zargarian,d Remi Chauvin a
a Laboratoire de Chimie de Coordination du CNRS, BP 44099, 205 route de Narbonne, 31077 Toulouse Cedex 4. France
b Institut de Quimica Computacional and Departament de Quimica, Universitat de Girona, Campus de Montilivi, 17071 Girona, Catalonia, Spain
c Laboratoire de Chimie Théorique, (UPMC) CC137, 4, place Jussieu 75252 Paris cedex 05. France
dDépartement de Chimie, Université de Montréal, Québec H3C 3J7. Canada
- christine.lepetit@lcc-toulouse.fr
Theoretical studies of the metal-phenylene interaction in a P(CH)P pincer rhodium(I) complex
Amidiniophosphines are cationic electron-poor ligands in transition metal complexes with potential original catalytic properties [1]. In the pincer ligand series, based on the meta-phenylene bridge, the first unsymmetrical nickel(II) complex with a phosphinite and an imidazoliophosphine extremities was recently reported [2]. However switching to rhodium(I) and to two cationic imidazoliophosphine extremities, a P(CH)P pincer complex 2 was obtained instead of the targeted PCP pincer complex 1 [3]. On the basis of DFT studies, an alternative pathway to the usually invoked C-H oxidative addition, is proposed to explain the formation of 2 and its reluctance to yield the PCP pincer 1 [3-4].
Experimental and calculated spectra of 2 are in good agreement with a weak rhodium-phenylene interaction. The signature of a Rh-h1-C1 bonding is found from ELF (electron localization function) topological analysis. Further comparative studies using molecular orbital analysis and other theoretical tools such as NCI (non covalent interactions) analysis and MCI (multicenter delocalization indices) [5] allow for refining the characterization of this particular metal-aryl bonding. To the best of our knowledge, it has been little discussed in the literature up to now, although it may be inserted as a missing entry in the well-known series of agostic and anagostic interactions [6].
References
[1] Y. Canac, C. Maaliki, I. Abdellah, R. Chauvin New. J. Chem. 2012, 36, 17.
[2] B. Vabre, Y. Canac, C. Duhayon, R. Chauvin, D. Zargarian Chem. Commun. 2012, 48, 10446.
[3] C. Barthes, C. Lepetit, Y. Canac, C. Duhayon, D. Zargarian, R. Chauvin Inorg. Chem. 2013, 52, 48.
[4] B. Vabre, M. L. Lambert, A. Petit, D. H. Ess, D. Zargarian Organometallics 2012, 31, 6041.
[5] a) P. Bultinck, R. Ponec, S. Van Damme J. Phys. Org. Chem. 2005, 18, 706. b) I. Mayer, P. Salvador Chem. Phys. Lett. 2004, 383, 368.
[6] M. Brookhart, M. L. H. Green, G. Parkin PNAS, 2007, 104, 6908.
Piotr de Silvaa,b, Jacek Korchowieca, Tomasz A. Wesolowskib
aK. Gumiński Department of Theoretical Chemistry, Faculty of Chemistry, Kraków Jagiellonian University, R. Ingardena 3, 30-060 Kraków, Poland
bUniversité de Genève, Département de Chimie Physique 30, quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
Probing Aromaticity and Non-Covalent Interactions with Single Exponential Decay Detector (SEDD)
Single Exponential Decay Detector (SEDD) [1] is an orbital-free descriptor of bonding patterns in molecular systems. Calculation of SEDD requires only the density ans its first and second derivatives
SEDD is able to distinguish such localized elements of electronic structure like bonds, lone pairs and core electrons. In this contribution the recent applications of SEDD are outlined. This includes investigating electron delocalization in aromatic system [2] as well as non-covalent interactions.
[1] P. de Silva, J. Korchowiec, T. A. Wesolowski, ChemPhysChem 2012, 13, 3462-3465.
[2] P. de Silva, J. Korchowiec, N. Ram J.S., T. A.Wesolowski, Chimia (in press).
Ferran Feixas, Steffan Lindert, Meekyum Kim, William Sinko, J. Andrew McCammon
Department of Pharmacology, University of California, San Diego, 9500 Gilman Drive, La Jolla, 92039, California, USA'
The role of cis-prenyl diphosphate synthases in bacterial cell-wall biosynthesis: non-covalent interactions and its implications for drug discovery
The rise in resistance to antibiotics is becoming a major problem. Thus, there is a need for new drugs, with new chemical structures, that target alternative enzymes. In this work, we focus on the cis-farnesyl-diphosphate synthase (cis-FPPS) and cis-decaprenyl-diphosphate synthase (cis-DPPS), two enzymes found in Mycobacterium Tuberculosis that work in series to synthesize decaprenyl-diphosphate, which is essential for the bacterial cell-wall biosynthesis. Interestingly, both enzymes present similar 3D-structures and catalytic mechanism. However, the length of the products synthesized by the two enzymes is remarkably different, cis-FPPS produces a short 15-carbon product called farnesyl-diphosphate, while cis-DPPS generates a larger 50-carbon decaprenyl-diphosphate.
The first aim of this work is to perform MD simulations and docking studies to analyze the influence of protein flexibility on the design of cis-prenyl-diphosphate synthase inhibitors. To this end, we make use of the Relaxed Complex Scheme (RCS), a type of ensemble docking, which relies on previously determined conformations from MD simulations to perform docking studies to select the most promising compounds that will be tested experimentally [1]. This strategy has been successfully applied to design promising inhibitors in related systems such as undecaprenyl diphosphate synthases (UPPS)[2]. The second aim is to study in deep detail the interaction between the receptor and some of the experimentally active inhibitors obtained in the first step. In particular, we focus our attention on the analysis of the non-covalent interactions on the hydrophobic cavity of the active site of cis-FPPS and DPPS enzymes. We will see how a deeper understanding of intermolecular drug-receptor non-covalent interactions can help us on the drug discovery process.
[1] W. Sinko, S. Lindert, J. A. McCammon, Chem. Bio. Drug. Des. 2013, 41.
[2] W. Zhu, Y. Zhang, W. Sinko, M. E. Hensler, J. Olson, K. J. Molohon, S. Lindert, R. Cao, K. Li, K. Wang, Y. Wang, Y. Liu, A. Sankovsky, C.A.F. Oliveira, D. A. Mitchell, V. Nizet, J.A. McCammon, E. Oldfield, PNAS, 2013, 123
Robin Chaudret, Natacha Gillet, Pan Wu, Julia Contreras-Garcia, Jean-Philip Piquemal, Weitao Yang
IFP Energies nouvelles 1 et 4 Av de Bois-Préau 92852 Rueil-Malmaison CEDEX
Non covalent interaction within reactive and dynamic systems : a presentation of ELF/NCI and aNCI
Weak interactions play a predominant role in various problems such as ligand docking, catalytic stereoselectivity or formation of supramolecular structures but also in initiating chemical reactions. These interactions can be very different in nature or energy but they all share a same characteristic: the difficulty for a chemist to localize them. Indeed their energies are really small and therefore their geometries or any other experimental properties are subject to strong variations. It is then particularly difficult to characterize them energetically or spatially. For the theoretical chemist, the only way to identify them is to use very precise (and computationally expensive) methods which can be incompatible with the size of a realistic systems, the exploration of the different possible conformations or the research of a reaction profile. The recently developed NonCovalent Interaction (NCI) index [1-3] allows a characterization and an estimation of these weak interactions. This approach is based on the study of the reduced density gradient as a function of the density and allows a visualization of the non-covalent interactions. The latter are colored in function of their strength. The first part of the talk will briefly describe the NCI analysis. This new methods appear therefore perfectly adapted to visualize the early steps of a reaction (when the interaction is still quite weak) but when the interaction gets stronger, the NCI validity decreases. In such a range, however, other approaches such as the Electron Localization Function [4-6] (ELF) have shown great success. The coupling of the two methods should allow an understanding of every steps of the reaction, from the weak interaction up to the bond formation. The second part of the talk will describe the ELF-NCI [7] approach and demonstrate its interest on several reaction mechanisms.
Finally, the importance of the environment effect on a reaction is now largely recognized and is mainly studied through QM/MM calculations. If NCI can localize the weak interactions, it is, however, strongly dependent on the input structure, which makes it unsuitable for analyzing fluctuating systems such as solute/solvent interactions or protein/ligand binding. In the last part of the talk we will present the new average Non Covalent Interaction (aNCI) approach [8] especially designed to treat such problems. The theory and several examples including benzene solvation, solvent effect along a reaction mechanism and ligand docking will be presented.
[1] Johnson, E. R.; Keinan, S.; Journal of the American Chemical Society 2010, 132, 6498-6506.
[2] Contreras-García, J.; Johnson, E. R.; Journal of Chemical Theory and Computation 2011.
[3] Contreras-García, J.; Yang, W.; Johnson, E. R. The Journal of Physical Chemistry A 2011, 115, 12983-12990.
[4] Becke, A. D.; Edgecombe, K. E. Journal of Chemical Physics 1990, 92, 5397-5403.
[5] Silvi, B.; Savin, A. Nature 1994, 371, 683-686.
[6] Chaudret, R.; Cisneros, G. A.; Chemistry-A European journal 2011, 17, 2833-2837.
[7] Gillet, N., Chaudret, R.; J. Chem. Theory Comput., 2012, 8 3993–3997.
[8] Wu, P., Chaudret, R.; J. Chem. Theory Comput., Article ASAP; DOI: 10.1021/ct4001087.
Gabriele Saleh, Carlo Gatti, Leonardo Lo Presti
Department of Chemistry, Università degli Studi di Milano, via C.Golgi 19, 20159 Milano Italy
Exploring non-covalent interactions in molecular crystals through the analysis of their experimental electron density distribution
The so-called ‘NCI descriptor’ introduced by Johnson et al. [1], besides allowing an easy-to-catch pictorial visualization of Non-Covalent Interactions (NCI), presents the big advantage of being based on electron density (ED) distribution, allowing therefore to investigate NCI theoretically and experimentally on the same ground. Indeed, ED is an observable and measurable quantity (from X-ray diffraction experiments) which, at the same time, can be straightforwardly obtained from the wavefunction. Quite recently, we have applied the 'NCI descriptor' for the first time to experimentally derived ED [2]. We have also compared the experimental results with the ones obtained from fully periodical DFT calculations, to assess the reliability of experimentally derived NCI picture. Moreover, we have contrasted the picture obtained from 'NCI descriptor' with the one obtained from Quantum Theory of Atoms in Molecules [3] to show how the two approaches may fruitfully complement each other. In a separated work [4] we carried out a systematic comparison between the NCI picture obtained from promolecular ED (i.e. the ED obtained by superposing the atomic densities) and the one from experimentally derived ED, to explore to what extent the promolecular ED can be used to investigate NCI, as proposed by Contreras et al. [5] . The results relative to the application of NCI descriptor to experimental, theoretical and promolecular ED of several systems will be presented. In addition, preliminary results obtained by mapping the energy density on the Reduced Density Gradient isosurfaces will be shown.
[1] E. R. Johnson, S. Keinan, P. Mori-Sanchez, J. Contreras-García, A. J.Cohen, W. Yang, J. Am. Chem. Soc. 2010, 132, 6498 – 6506
[2] G. Saleh, C. Gatti, L. Lo Presti, J. Contreras-Garcia, Chem. Eur. J. 2012, 18, 15523–15536
[3] R. F. W. Bader, Atoms In Molecules: A Quantum Theory, Clarendon Press, Oxford, 1990
[4] G. Saleh, C. Gatti, L. Lo Presti, Comput. Theor. Chem. 2012,998,148–163
[5] J. Contreras-García, E.R. Johnson,S. Keinan, R.Chaudret , J-P. Piquemal, D.N. Beratan, W. Yang, J. Chemical Theory and Computation, 2011, 7, 625-632
Monica Calatayud
UPMC Univ Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005,Paris, France
CNRS, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France Institut Universitaire de France
Reactivity indices and bonding in V2O5 based materials
Vanadium oxide is stable as gas-phase cluster, single-crystal or as supported catalyst. In this talk the similarities and differences in chemical reactivity and bonding in these compounds will be described by means of conceptual DFT methods (molecular electrostatic potential, Fukui function) and bonding theories (QTAIM, ELF).
[1] A systematic density functional theory study of VxOy
+ and VxOy (x=1-4, y=1-10). M. Calatayud, J. Andrés, A. Beltran; J. Phys. Chem. A 105 (2001) 9760
[2] An AIM and ELF study on the interaction between O2 and VxOy + / VxOy (x=1-2, y=1-5) systems. M. Calatayud, S. Berski, A. Beltran and J. Andrés; Theor. Chem. Acc. 108 (2002) 12
[3] A theoretical study on the structure, energetics and bonding of VOx + and VOx (x=1-4) systems. M. Calatayud, B. Silvi, A. Beltran, J. Andres; Chem. Phys. Lett. 333 (2001) 493
[4] The hierarchy of localization basins: a tool for the understanding of chemical bonding exemplified by the analysis of the VOx and VOx+ (x=1-4) systems. M. Calatayud, J. Andrés, A. Beltran, B. Silvi; Theor. Chem. Acc. 105 (2001) 299
[5] Understanding the reactivity and catalytic properties of gas-phase, crystal and supported V2O5 systems using Density Functional Theory based Reactivity Indices. M. Calatayud, F. Tielens, F. De Proft Chem. Phys. Lett. 456 (2008) 59
Kati Finzel
Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Straße 48, 01187 Dresden, Germany
email: kati.finzel@cpfs.mpg.de
Indicators based on electron density inhomogeneity
The approach of the ω-restricted space partitioning [1] (originally developed for the derivation of the electron localizability indicator, ELI [2]) was used to generate the charge sampling functionals Cp [3]. Whereas the space partitioning is controlled by the number of pairs in case of ELI-D, the control function is the electron density inhomogeneity in case of Cp. The electron density inhomogeneity was defined in terms of a distance measure between the electron density and its average value within a microcell. The distance depends on the inhomogeneity parameter p, which allows to choose the best parameter p upon certain criteria.
Assuming that one can approximate the volume occupied by a pair (a two-particle property) by the volume enclosing a certain amount of electron density inhomogeneity (a one-particle property), a least square fit searching for the inhomogeneity parameter popt for which Cp resembles best ELI-D was done for the atoms Li to Xe. The values of this optimal parameter popt where found to concentrate around p=0.6 [4]. As the optimal value is nearly independent of the atom type, the resulting functional C0.6 was applied to molecular systems and solids. Whereas the topologies of C0.6 and ELI-D are quite similar for atoms, both indicators show different behavior in the bonding region of molecular systems and solids.
[1] A. M. Pendas, M. Kohout, M. A. Blanco and E. Francisco. Beyond Standard Charge Density Topological Analysis. In C. Gatti and P. Macchi, editors, Modern Charge Density Analyses. Springer, London, 2010.
[2] M. Kohout. A Measure of Electron Localizability. Int. J. Quantum Chem. 97:651-658,2002.
[3] K. Wagner and M. Kohout. Atomic shell structure based on inhomogeneity measures of the electron density. Theor. Chem. Acc., 128:39-46,2011.
[4] K. Finzel, Yu. Grin, M. Kohout. Chemical bonding descriptors based on electron density inhomogeneity measure: a comparison with ELI-D, Theor. Chem. Acc. 131, 2, 1-8, Seq. No.:1106, 2012.
Benoit de Courcy
UPMC Univ Paris 06, UMR 7616, Laboratoire de Chimie Théorique, F-75005,Paris, France
CNRS, UMR 7616, Laboratoire de Chimie Théorique, F-75005, Paris, France Institut Universitaire de France
Unravelling Non Covalent Interactions within Flexible Biomolecules: from electron density topology to gas phase spectroscopy
NCI will be applied to document the presence of intra-and inter-molecular interactions in flexible molecules of biological interest.
NCI will be applied to a series of stable conformations of isolated molecules as an interpretative technique to decipher the different physical interactions at play in systems like the one illustrated in Figure 1. Examples are chosen among neutral molecular systems exhibiting a large diversity of interactions, in particular inter- and intra-molecular hydrogen bonds, and for which an unambiguous structural determination was achieved using state-of-the-art conformer-specific vibrational spectroscopy. The interactions documented range from weak intra-molecular H-bonds in simple amino-alcohols to complex bonding patterns observed in model peptides as well as in chiral bimolecular systems. NCI successfully (Figure 2) allows spectroscopists to provide an elegant visualization of the interactions that potentially impact their spectroscopic vibrational probes, namely the OH and NH stretching motions.
In this contribution I will illustrate the power and the conditions of use of the NCI technique, with the aim to provide an easy tool for all chemists, experimentalists and theoreticians, for visualizing and characterizing interactions in systems involving complex H-bonding patterns.
Dylan Jayatilaka
University of Western Australia, 35 Stirling Hightway, Nedlands, Perth, WA 6009 Australia
Some thoughts about efficient implementation of accurate intermolecular interaction energies
Accurate intermolecular interaction energies require distributed multipoles and polarizabilities. However the monopole polarizabilities are difficult to deal with. In addition, monopole charges are by themselves long-ranged; a neutral molecule at long range has no monopole term. It is proposed to use a charge density partition for neutral species into neutral atoms. Thus only dipole and higher terms are required. Furthermore, field dependent partitioning allows monopole polarizabilities to be eliminated. Equations for this are derived and presented. Finally, it is pointed out that the so-called exchange-repulsion term arises largely from the electrostatic penetration term, which may be calculated exactly from an electrostatic model, as pointed out long ago by Spackman. Are we then ready for accurate ab-initio derived force fields?
Poster session
Agnieszka Roztoczyńska, Justyna Kozłowska, Paweł Lipkowski, Wojciech Bartkowiak
Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 53-370, Wrocław, Poland
Influence of spatial confinement on the cooperativity in the hydrogen-bonded clusters
Cooperativity is a well-known phenomenon and a steering factor of many chemical and physical properties of molecular matter [1-5]. In the present contribution the influence of chemical compression on the hydrogen bond cooperativity has been investigated. For that purpose the cooperativity effects in the isolated and confined (HF)n and (HCN)n clusters were considered. Two models of spatial confinement have been applied in our study. The first one is the cylindrical model harmonic potential, which certainly reflects the general aspects of the confinement effect. The second, less simplified model based on the supermolecular approach: molecular cages within which one can distinguish helium and carbon nanotubes. In order to get an insight into the nature of the cooperative effects in confined spaces a detailed analysis of the energetics and the topological features of electron density according to the concept of the “Quantum Theory of Atoms in Molecules” (QTAIM) proposed by Bader [6,7] have been performed. The properties in question have been analyzed using the ONIOM (M06-2X/6-311++G(2df,2pd):M06-2X/6-31G(d)) method. Moreover, in the case of the harmonic potential calculations, the second-order Møller-Plesset perturbation theory (MP2) has been applied as well. As a part of the study, the comparison of the data obtained within different models of orbital compression have been performed in order to establish the correspondence between the chemical environment and its approximate representations.
[1] A. Korpfen, J. Phys. Chem. 100, 13474 (1996).
[2] T. Kar, S. Scheiner, J. Phys. Chem. A 108, 91619168 (2004).
[3] M. Ziółkowski, S. J. Grabowski, J. Leszczyński, J. Phys. Chem. A 110, 65146521 (2006).
[4] H.J. Song, H.M. Xiao, H.S. Dong, J. Chem. Phys. 125, 074308 (2006).
[5] S. J. Grabowski, Theor. Chem. Acc. 132, 1347 (2013).
[6] R. F. W. Bader, Atoms In Molecules, A Quantum Theory; Oxford University Press: Oxford, U.K., 1990.
[7] R. F. W. Bader, Chem. Rev. 91, 893 (1991).
Paweł Lipkowski, Justyna Kozłowska, Agnieszka Roztoczyńska, Wojciech Bartkowiak
Theoretical Chemistry Group, Institute of Physical and Theoretical Chemistry, Wrocław University of Technology, Wybrzeże Wyspiańskiego 27, 50-370 Wrocław, Poland
Behavior of model hydrogen bonded complexes under the influence of external pressure
Despite of the numerous theoretical and experimental investigations concerning the hydrogen bond phenomenon as well as to the confinement effect there is still a small number of reports on the behavior of the H-bonded systems in the presence of external pressure [1-4]. In order to render the influence of pressure one can apply a model confining potentials which allows to describe the pure valence repulsion contribution. In the present study we have investigated the influence of external pressure, represented by the cylindrical harmonic potential V_conf=1/2(x^2+y^2)ω^2, on the selected linear hydrogen bonded complexes. The model dimeric systems: HF…HF, HCN…HCN, HCN…HCCH have been chosen as a case study. All calculations have been performed employing B3LYP, M06-2X and MP2 methods together with the 6-311++G(2df,2pd) basis set. Moreover, the varied external pressure has been controlled by the oscillator strength (ω=0.0-0.8). The studied complexes have been oriented along the z-axis of the Cartesian coordinate system and fully optimized in the presence of the two-dimensional harmonic oscillator potential. A detailed analysis of changes in topological parameters of hydrogen bonds resulting from the orbital compression has been performed according to the concept of the “Quantum Theory of Atoms in Molecules” (QTAIM) proposed by Bader [5,6]. Furthermore, an energetic analysis has shown a different manner of the interaction energy in the studied H-bonded complexes.
[1] W. Wang, D. Wang, Y. Zhang, B. Ji, A. Tion, J. Chem. Phys. 134 (2011) 054317.
[2] M.C. Gordillo, J. Marti, Chem. Phys. Lett. 329 (2000) 341.
[3] M. Jabłoński, M. Solà, J. Phys. Chem. A 114 (2010) 10253.
[4] G. Miño, R. Contreras, Chem. Phys. Lett. 486 (2010) 119.
[5] R. F. W. Bader, Atoms In Molecules, A Quantum Theory; Oxford University Press: Oxford, U.K., 1990.
[6] R.F. Bader, Chem. Rev. 91 (1991) 893.
Daniela Josa
Centro Singular de Investigación en Química Biolóxica e Materiais Moleculares (CIQUS), Universidade de Santiago de Compostela, Rúa Jenaro de la Fuente, s/n, Santiago de Compostela, 15782, Spain
A DFT-D study of stacking interactions between buckybowls and fullerenes
The first strong evidence of π···π interactions between corannulene and C60 was published in 2007, when Sygula and co-workers synthesized the molecular tweezers made up two units of corannulene that can trap one fullerene (Figure 1) [1]. Since then, the interest by concave-convex π···π interactions has been revived. The modification of the tweezers by a different functionalization of buckybowls, which constitute the main part of tweezers, could improve their efficiency and selectivity. Therefore, the aim of this work is to carry out a detailed study of the effects that can enhance the stacking interactions between buckybowls and fullerenes to achieve our ultimate goal that is predict how to modify the molecular tweezers to improve their efficiency and selectivity. All complexes studied were optimized at the B97-D/TZVP level using resolution of identity approximation (RI) implemented in TURBOMOLE 5.10 program suite [2]. Counterpoise corrections were applied to all reported interaction energies [3].
References
[1] A. Sygula, F.R. Fronczek, R. Sygula, P. W. Rabideau, M. M. Olmstead. J. Am. Chem. Soc. 2007, 129, 3842-3843.
[2] R. Ahlrichs, M. Bär, M. Häser, H. Horn, C. Kölmel. Chem. Phys. Lett. 1989, 162, 165-169.
[3] S.F. Boys, F. Bernardi. Mol Phys. 1970, 19, 553-566.
Olga Syzgantseva, Vincent Tognetti, Laurent Joubert
UMR CNRS 6014 & FR 3038, Université de Rouen et INSA de Rouen
Analysis of halogen bonding using topological indexes and IQA energy decomposition scheme
Halogen bonding is a directional interaction occurring between a halogen and an electronegative atom of a Lewis base. It is namely responsible for ligand binding and molecular folding in biological systems and is of importance in crystal engineering [1,2]. While predominantly electrostatic character of this interaction seems to be established [3], the influence of chemical environment on its strength is not completely understood. In the present study we use topological indexes and an interacting quantum atoms (IQA) energy decomposition scheme [4] to explain how primary and secondary interactions affect the strength of a halogen bond. For this purpose, we use a series of model complexes described within DFT applying a dispersion correction. IQA energy decomposition scheme allows to extract pairwise interatomic contributions from total intermolecular interaction energy between halogen-bonded complex moieties and to estimate the proportion between Coulomb and exchange contributions [5]. Besides, the Source Function [6] analysis complements the description of the impact of chemical surroundings on the formation of halogen bond.
[1] P. Auffinger, F.A. Hays, E. Westhof, S. Ho, Proc. Nat. Acad. Sci, 2004, 101, 16789 – 16794
[2] P. Metrangolo, F. Meyer, T. Pilati, G. Resnati, G. Terraneo, Angew. Chem. Int. Ed., 2008, 47, 6114 – 6127
[3] P. Politzer, K.E. Riley, F.A. Bulat, J.S. Murray, Comp. Theor. Chem., 2012, 998, 2–8
[4] M. A. Blanco, A. Martín Pendás, E. Francisco, J. Chem. Theory Comput., 2005, 1, 1096 – 1109
[5] O.A. Syzgantseva, V. Tognetti, L. Joubert, in preparation
[6] R.F.W. Bader, C.Gatti, Chem. Phys. Lett., 1998, 287, 233 – 238
Rodrigo J. Alvarez-Mendeza**, Jorge M. del Campob, José L. Gázquezc and Alberto Velaa
a Departamento de Química, Cinvestav, Av. IPN 2508, Colonia San Pedro Zacatenco 07360, México, D.F. 07360 México
b Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México, D. F. 04510 México
c Departamento de Química, Universidad Autónoma Metropolitana-Iztapalapa, Av. San Rafael Atlixco 186, México, D. F. 09340 México
Chemical interactions and dimensionless scalar fields
An analysis capable of identifying both strong and weak chemical interactions is presented. It is based on a generalization of the analysis presented by Zupan et. al [1,2] that consist in evaluating the integral
where 
is the diagonal element of the k-nt reduced density matrix and 
is the Heaviside function. 
represents the numbers of electrons in k-nt arrangement localized in the region constrained by the 
scalar field isosurface 
. For k=1, the derivative with respect to the φ yields the correspondent density of states and, at the same time, its bounds allow to localize the values that are chemically relevant for the scalar field.
Many examples to ilustrate the behavior of this analysis are presented using the exchange and correlation dimensionless gradients and the normalized electron localization function [3].
References
[1] A. Zupan, J.P. Perdew, K. Burke and M. Causa. Int. J. Quantum Chem. 61, 835 (1997)
[2] A. Zupan, K. Burke, M. Ernzerhof y J.P. Perdew, J. Chem. Phys. 106, 10184 (1997)
[3] Analysis of the generalized gradient approximation for the exchange energy. J. L. Gázquez , J. M. del Campo, S. B. Trickey, R. J. Alvarez-Mendez, A.Vela. Chapter in "Concepts and Methods in Modern Theoretical Chemistry, Vol. 1: electronic structure an reactivity". In honor of Professor B.M. Deb; S.K. Ghosh and P.K. Chattaraj eds. Taylor & Francis / CRC Press. (2013)
- Presenting autor
Rodrigo A. Cormanich,a,b Luke Crawford,a Michael Bühl,a Cláudio F. Tormenab and Roberto Rittnerb
a EastChem School of Chemistry, University of St Andrews, UK
b Chemistry Institute, University of Campinas, Brazil.
- email: rac22@st-andrews.ac.uk
Conformational analysis study of Ac-gly-NHMe, CF3-C(O)-gly-NHMe and Ac-gly-NH(Me)2 through theoretical calculations and NMR spectroscopy
In any attempt to understand protein folding processes, the characterization of the conformational preferences of short peptides constitutes an unavoidable first step. Indeed, there is a great interest in recognizing short peptide equilibrium geometries and to extrapolate the results obtained in this way to polypeptides and proteins [1]. However, free amino acid compounds and short peptide models (Ac-R-NHMe; R=amino acid) conformational geometry preferences have been commonly attributed to intramolecular hydrogen bonding (IHB) formation [2]. In particular, Ac-gly-NHMe has been the target of several studies [3], but despite steric and hyperconjugative effects, which are ubiquitous and influence the energies and geometries of even simpler systems [4], C5 and C7 Ac-gly-NHMe conformers stabilities have been explained solely by the possibility of a five- and a seven-membered N-H...O IHB formation, respectively. In this way, the present study aims to investigate the interactions that rule the conformational isomerism of Ac-gly-NHMe, CF3-C(O)-gly-NHMe and Ac-gly-N(Me)2 compounds through experimental (obtainment of 3JHH and 1JNH spin-spin coupling constants values) and theoretical calculations, using DFT funcionals with dispersion corrections (DFT-D and DFT-D3) and in the framework of the quantum theory of atoms in molecules (QTAIM) and natural bond orbitals (NBO). The results obtained by these methods indicate that Ac-gly-NHMe C5 conformer should not form an IHB, while conformer C7 forms a strong and stable N-H...O seven-membered IHB (Scheme 1). Indeed, Ac-gly-NHMe conformational preferences are ruled by IHB, but the interplay of hyperconjugative and steric interactions with IHB show to be important as well. Also, the CF3-C(O)-gly-NHMe C7 conformer has a weaker IHB force in comparison with the Ac-gly-NHMe C7 conformer, which increases the C5 conformer population and indicate that C7 stability is strongly seven-membered IHB force dependent. In addition, the Ac-gly-N(Me)2 compound cannot form a seven-membered N-H...O IHB and only the C5 conformer is observed for this compound. Thus, our analysis indicates that Ac-gly-NHMe conformational preferences are ruled by IHB in the gas phase, but that steric and hyperconjugative interactions have also an important contribution and cannot be ignored. QM/MM calculations are currently being performed in order to understand the behavior of these compounds in solution.
| Ac-gly-N(Me)2 | |||||
| C5 | C7 | C5 | C7 | C5 | |
|
|
|
|
| |
| %P | 12% | 88% | 44% | 55% | 100% |
Scheme 1: C5 and C7 Ac-gly-NHMe, CF3-C(O)-gly-NHMe and Ac-gly-N(Me)2 conformers molecular graphs obtained from the QTAIM calculations. Conformational populations (%P) and the eletronic density (ρ) values at the intramolecular hydrogen bonding BCP are indicated.
[1] M. F. Masman, R. F. Murphy, R. D., Enriz, A. M. Rodriguez J. Phys. Chem. 126 10682 (2007).
[2] R. A. Cormanich, L. C. Ducati, R. Rittner J. Mol. Struct. 1014, 12 (2012).
[3] H. Fujitani; A. Matsuura; S. Sakai; H. Sato; Y. Tanida J. Chem. Theory Comput. 5 1155–1165 (2009).
[4] R. A. Cormanich, Freitas M. P., J. Org. Chem.74, 8384 (2009).
I. González-Velosoa, D. Josaa, J. Rodríguez-Oteroa, E. M. Cabaleiro-Lagob
aCentro de investigación en Química Biolóxica e Materiais Moleculares, CIQUS, Universidade de Santiago de Compostela Rúa Jenaro de la Fuente s/n, Santiago de Compostela, 15782, SPAIN bDepartamento de Química Física, Universidade de Santiago de Compostela Facultade de Ciencias, Avda. de Alfonso X o Sabio s/n, Lugo, 27002, SPAIN
Theoretical study of the interaction between planar and curved compounds: perfluoro-ortho-phenylenmercury and corannulene
Recently, Filatov et al.[1] have performed an experimental study of the interaction of perfluoro-ortho-phenylenemercury [Hg3] with corannulene and monoindenocorannulene. The first compound is planar; however the other two are curved bowls. Therefore, these systems allow studying the mutual structural influences between surfaces of planar and nonplanar molecules upon their interaction. This kind of interactions has generated considerable interest in materials chemistry. According to the experimental study, significant geometry adjustments of both interacting partners were found in the solid-state complexes. So, Hg3 adopts highly bent configurations, which are substantially deviating from planarity, whereas flattening of the bowl-shaped polyarenes is rather more modest. corannulene Hg3
In the present work a theoretical study was performed at the DFT and ab initio levels to test the experimental fact, finding that corannulene flattened slightly while [Hg3] curved more sharply. It also checked that the union of these atoms is due basically to dispersion energy between molecules. For that reason the DFT functionals have been corrected by an empirical dispersion term as designed by Grimme (DFT-D)[2]. Finally it was found that the substituents on the [Hg3] also affect the interaction energy between these monomers, as well as the deformation between the molecules.
[1] A. S. Filatov, E. A. Jackson, L. T. Scott, M. A. Petrukhina, Angew. Chem. Int. Ed. 2009, 48, 8473.
[2] S. Grimme, S. Ehrlich, L. Goerik, J. Comput. Chem., 2011, 32, 1456.
Roberto Álvarez Botoa, Mónica Calatayudb, Julia Contreras Garcíab and A. Martín Pendása
aUniversidad de Oviedo. Spain
bLaboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris, France
The role of packing in the platinum-octylamine system: a NCI analysis
The binding of long alkyl molecules to metal surfaces has been widely investigated due to its many technological applications such as films, electronic junctions, coating or tribology [1]. The intermolecular interactions play a key role not only in the physicochemical properties of these systems, but also in their stabilization. These interactions come from the vicinity of the neighboring molecules forming a two-dimensional layer. Despite the amount of investigations concernig metal-supported self assembled monolayers (SAMs), little is known about the dispersion forces present in these 2D arrangements. In the present work the platinum (100)- octylamine(NH2C8H17)has been chosen as model system to show the role of packing. A bonding analysis based on NCI index [2] was applied to both, the octylamanine dimer and the 2D arrangement to ilustrate the nature of the intermolecular interactions. Additionally, a detailed study on the structural parameters of the alkyl-system is reported to gain understanding the factors affecting the stability of the system. All the calculations were carried out at DFT theory level combined with Grimme D2 approach.
[1] Caroline Salzemann; Christophe Petit. Influence of Hydrogen on the Morphology of Platinium and Palladium Nanocrystals. Langmuir, 2012, 28, 4835-4841.
[2] E. R. Johnson, S. Keinan, P. Mori-Sánchez, J. Contreras-García, A. J. Cohen and W. Yang J. Am. Chem. Soc. 132, 6498-6506 (2010).
M. Merced Montero-Campillo,Al Mokhtar Lamsabhi, Otilia Mó and Manuel Yáñez
Departamento de Química, Campus de Cantoblanco, Universidad Autónoma de Madrid, 28049 Madrid SPAIN-e-mail: mm.montero (at) uam.es
Modulating Properties By Beryllium Bonds: The Squaric Acid Case
Non-covalent interactions usually involve closed-shell species, as is the case of hydrogen bonds, beryllium bonds, halogen bonds or Van der Waals complexes. As an electropositive atom with available empty orbitals that can accept electrons, beryllium atom can behave similarly to hydrogen, establishing the so called beryllium bonds when interacting with electron donors. These weak forces may play a very important role in the properties of a given molecular system and, sometimes, a cooperative effect between them is possible when more than one is present. We have recently studied some cases of this cooperativity between hydrogen and beryllium bonds. Chemical properties as acidity can significantly change when the system is perturbed by one or more weak interactions of this kind. It is also possible to modulate self-association in some molecular systems through the formation of beryllium bonds, as we recently shown for squaric acid dimers.
Squaric acid molecule presents a beautiful crystalline assembly formed by planar layers in which the acid units interact through very weak hydrogen bonds. Our group investigated the effect of complexation with beryllium hydride in this molecule and its consequences on its self-assembly capacity. We also studied sulfur and selenium derivatives of squaric acid, as well as the acidity of the whole squaric acid series, analyzing how this property changes on-going through complexation with beryllium hydrides. G4 method was used for the electronic structure calculations. Atoms in Molecules theory (AIM), the electron localization function (ELF) and the Non-Covalent Interaction (NCI) index have been used in this study.
References
[1] Montero-Campillo, M. M; Lamsabhi, A. M.; Mó, O.; Yáñez, M. J. Mol. Mod. 2012 DOI: 10.1007/s00894-012-1603-0
[2] Montero-Campillo, M. M.; Lamsabhi, A.M.; Mó, O.; Yáñez, M. Squaric Acid Derivatives and BeH2 Complexes: The Role of Beryllium Bonds On Spontaneous H2 Lost and Acidity Enhancement. In preparation.
Martine Castellà-Ventura and Emile Kassab
Laboratoire de Chimie Théorique, CNRS-UMR7616, Université Pierre et Marie Curie, Paris, France
DFT study of the non-covalent interation effects on the energetics of 4,4'-bipyridine adsorption in H-ZSM-5 zeolite
The effects of the zeolite framework on the adsorption of bidentate 4,4’-bipyridine (44BPY) ligand in the straight channel of H-ZSM-5 (Z) has been investigated by DFT calculations using M06-2X functional to account for dispersive van der Waals interactions [1]. The straight channel of H-ZSM-5 is simulated by a quantum cluster constituted of 32 tetrahedral units (T). This cluster has two Al atoms located at positions sufficiently distant from each other allowing the bidentate 44BPY to interact with the two Brønsted acid sites of the zeolite. The two pyridyl rings of 44BPY have distinct zeolite environment: only one pyridyl ring (PY1) is completely surrounded by the channel wall, while the second pyridyl ring (PY2) is located in a region opened towards the zigzag channel. Thus the confinement effects exerted from the zeolite framework on the two pyridyl rings are different.
The minimum energy pathways of the double proton transfer from H-ZSM-5 to 44BPY ligand were calculated. The potential energy surfaces are characterized by two minima corresponding to the monodentate 44BPYH+/Z- and bidentate 44BPYH22+/Z2- ion pair complexes. No energy minimum is found for neutral hydrogen bonding structure. Two monodentate ion pair complexes differing by the protonated pyridyl ring may be formed, which leads then to the formation of the bidentate 44BPYH22+/32T2- complex via two distinct pathways (Figure). The relative stability between all adsorption complexes involved in this double proton transfer does not exceed 5 kcal mol-1. The energetics of the double proton transfer reaction is very sensitive to the confinement effects. Our results clearly show that the adsorption energy of 44BPY is mainly due to the effects of the zeolite framework executed through dispersive van der Waals interactions. These interactions decrease the relative stability of mono and bidentate complexes and stabilize the transition state connecting them.
[1] Y. Akacem, M. Castellà-Ventura and E. Kassab, J. Phys. Chem. A 2012, 116, 1261-1271
Yannick Geboes
Department of Chemistry: Research group Cryospectroscopy, University of Antwerp (UA), Groenenborgerlaan 171, B-2020 Antwerp, Belgium.
Halogen Bonding versus Hydrogen Bonding: Battle of the Weak Noncovalent Interactions
During the past decade, a growing attention has been given towards the noncovalent interaction in which a covalently bonded halogen atom interacts with an electron rich site (e.g. a lone pair in a Lewis base), the so-called halogen bond.[1] The name halogen bond has been chosen in order to stress the many similarities this interaction has with the better known hydrogen bonds, [2] such as its directionality and strength (between 5 kJ/mol and 180 kJ/mol for halogen bonds [3] ), which enables halogen bonds to compete with weak hydrogen bonds.
Even though many applications for halogen bonds already exist, experimental data allowing the rationalisation of their exact nature still remains scarce. Moreover, almost no experimental information on the structural competition between C-X∙∙∙Y halogen and C-H∙∙∙Y hydrogen bonding is available. [4] In our current research, we aim to rectify this by collecting and investigating experimental data on halogen and hydrogen bonded complexes formed by a series of combined halogen/hydrogen donors with F, O, N, S and P containing acceptor molecules, using infrared and Raman measurements of cryogenic solutions.
The experimental study of complexes formed by these noncovalent interactions by means of cryospectroscopy in liquid noble gases offers valuable insights into the formation of these complexes, since this measuring technique yields the opportunity to measure isolated complexes at equilibrium conditions, thus enabling the determination of thermodynamic properties of these complexes, such as the complexation enthalpy. [5-6] Furthermore, the experiments are supported by ab initio and Conceptual Density Functional Theory calculations, Monte Carlo simulations and Molecular Dynamics calculations.
During this poster session, results of the competition between the halogen and hydrogen bonded complex of CHF2I and trimethylamine will be presented.
[1] P. Politzer, J.S. Murray, ChemPhysChem 2013, 14, 278.
[2] A.C. Legon, Struct. Bond. 2008, 126, 17.
[3] P. Metrangolo, H. Neukirch, T. Pilati and G. Resnati, Acc. Chem. Res. 2005, 38, 386-395.
[4] (a) T. Di Paolo, C. Sandorfly, Can. J. Chem. 1974, 52, 3612. (b) C.B. Aakeröy, N.C. Schultheiss, A. Rajbanshi, J. Desper, C. Moore, Crystal Growth & Design 2009, 9, 432.
[5] B.J. van der Veken, J. Phys. Chem. 1996, 100, 17436-17438.
[6] D. Hauchecorne, N. Nagels, B.J. van der Veken, W.A. Herrebout, Phys. Chem. Chem. Phys. 2012, 14, 681.
E. Rezabala, T. Schäfera,b
aNanoBioSeparations Group, POLYMAT (UPV/EHU), Donostia-San Sebastian, Spain
bIKERBASQUE, Basque Foundation for Science, Bilbao, Spain
A first principle approach to solvation by lmidazolium based ionic liquids
Ionic liquids (ILs) bear highly versatile and powerful solvation properties, providing ideal media for chemical reactions and catalysis, liquid-liquid extraction and gas separation. The easy tuneability of their properties constitutes an ideal framework for their task oriented adjustment. In order to fully benefit from this great potential, and design ILs, predict their physical properties or even choose suitable ones for a specific reaction, it is necessary to provide a link between the fundamental properties of the system (such as electronic and molecular structure) and its specific (macromolecular) physical and chemical properties.
The physical properties of the ILs stem from the particularly strong Coulomb (long-range) and dispersion (short-range) interactions between the ion pairs. The equilibrium between these two forces results in rather flat potential energy surfaces and provides properties that fundamentally diverge from those of regular salts and solvents. Besides, the capability to form a wide range of strong intermolecular interactions boosts the solubility of many compounds in ILs.
As a first step towards the comprehensive understanding of IL solvation processes, systematic ab initio computational studies of single, isolated ion pairs was previously carried out. A direct relationship between the gas phase electronic structure and the experimentally observed behavior was found, and the nature of the interactions between the methylimidazolium based ILs and a set of different solutes was described in a comprehensive way. This fundamental knowledge set the grounds for predicting the behaviour of different ILs and solvated systems, and shed light into the solvation of bigger systems as DNA, nanostructures, or enzymes in ILs.[1]
In the present work, we extend our model to several ion pairs, with the aim of understanding the microsolvation of the solute and the role of the bulk on the solvation process. As a compromise between the real system and the affordable size for quantum chemical calculations, we have chosen a small, prototypical ionic liquid, [MMIM][Cl], and small solutes, namely H2O and CO2, for their intrinsic interest and fundamentally distinct solvation patterns observed previously. Ab initio MD calculations are carried out to study the structural characteristics of the systems, and the interaction between the solute and the surrounding molecules is analized by energy decomposition analysis (EDA) and topological non covalent interactions (NCI). This provides a complete and dynamic picture of the solvation process in ILs. In fact, we observe that the effect of the surrounding environment differ remarkably in both scenarios, and the patterns observed provide highly valuable information, applicable to more complex solutes.
[1] J. Phys. Chem. B, 2013, 117 (2), pp 553–562
Ana Martín-Sómer, Manuel Yáñez
Departamento de Química. Facultad de Ciencias. Módulo 13. Universidad Autónoma de Madrid. Campus de Excelencia UAM-CSIC. Cantoblanco. 28049-Madrid. Spain.
Acidity enhancement triggered by GaH3 association to phosphines
The complexes formed by the interaction between a series of phosphines R-PH2 (R=H, CH3, c-C3H5, C6H5, C6H5-CH2) and GaH3 have been investigated by means of high-level DFT calculations. These complexes behave as much stronger acids in the gas phase than the isolated phosphines. This dramatic acidity enhancement, which can be as high as 163 kJ mol-1, results from a much greater stabilization of the anionic deprotonated species with respect to the neutral one, upon GaH3 association. This effect depends quantitatively on the nature of the substituent R and is smaller for R=C6H5 because of the conjugation of the P lone pair with the aromatic system.
The analysis of the electron rearrangements of the phosphine and the galane moieties indicates that the P-Ga bond is significantly stronger in the anion than in the corresponding neutral phosphine-galane complex, because the deprotonated phosphine is a much better electron donor than the neutral phosphine. The delocalization of the newly created P lone pair with the P-Ga bonding density as well as the dispersion of the extra electron density into the GaH3 group also contributes to the increased stability of the deprotonated species.
M. M. San Franciscoa, B. Braïdab, A. Savinb, and A. Martín Pendása
aUniversidad de Oviedo. Spain
bLaboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris, France
Maximum Probability Domains and its application on some systems
Among the several methods that belong to the group of interpretative tools which are based on the partition of the real space upon the analysis of the wavefunction, the Maximum Probability Domains (MPD) [1,2] method provides a new landscape that allows to recover the classical Lewis picture of molecules and may open new views in terms of chemical bonding. Basically, a MPD domain is a region of the three dimensional space which maximizes the probability of finding a given number of electrons in a system. In this present work, we introduce the basics concepts of the MPD method and its application in some molecular systems is carried out to investigate the chemical bond. As well, the effect of electronic correlation is examined due to the application of Quantum Monte Carlo techniques.
[1] O. Mafra Lopez Jr., B. Braïda, M. Causà, A. Savin, Progress in Theoretical Chemistry and Physics vol 22, 173. ed Hoggan,
Springer UK, London (2011), 175-186.
[2] A. Scemama, M. Caffarel, A. Savin, J. Comp, Chem. 2007, 28, 442-454.
Daniel Menendez-Crespo a, Alberto Otero-de-la-Roza b, and Victor Luaña a
a Departamento de Quimica Fisica y Analitica, Universidad de Oviedo, E-33006 Oviedo, Spain
b Chemistry and Chemical Biology, School of Natural Sciences, University of California, Merced, 5200 North Lake Road, Merced, California 95343, USA
Calcite and aragonite: the effect of thermodynamic conditions on the chemical bonding of minerals
The calcite and aragonite phases of CaCO3 are explored using electronic structure calculations (with the quantum espresso package [1]) and the quasi harmonic thermal model [2,3]. The topological analysis of the electron density is carried out with the newest version of the critic2 code [4]. This is used to examine the pressure and thermal dependence of the basin integrated atomic properties: volume, charge, kinetic energy. This example portrays the capabilities in the new critic2 code.
[1] P. Giannozzi et al, J. Phys. Condens. Matter 21 (2009) 395502 <http://www.quantum-espresso.org>
[3] A. Otero-de-la-Roza and V. Luaña, "Equations of state and thermodynamics of solids using empirical corrections in the quasiharmonic approximation", Phys. Rev. B 84 (2011) 184103.
[3] A. Otero-de-la-Roza, D. Abbasi-Perez, and V. Luaña, "GIBBS2 A new version of the quasi-harmonic model code. II. Thermal models, features and implementation". Comput. Phys. Commun. 182 (2011) 2232.
[4] A. Otero-de-la-Roza, A. Martín Pendás, and V. Luaña, "CRITIC2" (To be published)
Isabelle Fourréa, Elsa Alvareza a,b, Patrick Chaquina
aLaboratoire de Chimie Théorique, UMR 7616 CNRS UPMC Sorbonne Universités, 4 place Jussieu F-75005 Paris, France.
bPresent address: Institut Lavoisier, UMR 8180 CNRS Université de Versailles St-Quentin-en-Yvelines, 45 avenue des Etats-Unis, 78035 Versailles, France.
email: isabelle.fourre@upmc.fr
Carbon-helium bonding: Substitution of H by He+ in some small organic molecules (CH4, C2H6, C2H4, C2H2, H2CO, C6H6)
The name of “noble gases” is due to their difficulty to undergo chemical bonds. Nevertheless, strong and relatively strong helium bonds were predicted in cationic compounds (see 1,2 for a review). This contribution deals with the possibility of forming such bonds in small typical organic molecules, by substitution of one or several H atoms by He+. A structural and energetical study (based on calculations up to CCSD(T)/aug-ccpVTZ level) of this unusual bonding, as well as a topological characterization of the resulting cations is undertaken. The most stable species are generally obtained by substituting about half of hydrogen atoms. The C-He bonds exhibit equilibrium distances ranging from 1.327 (C2H2He22+) to 1.129 Å in He2CO2+. Their dissociation energies range from ~5 kcal.mol-1 (He2CO2+, Z-C2H2He22+) to 25 kcal.mol-1 (C2HHe33+), but remain most frequently around 10 kcal.mol-1. However most of He+ substituted hydrocarbons are metastable with respect to C-C cleavage, except derivatives of ethane.
AIM and ELF topological descriptors classify the C-He bond as a weak charge-shift interaction rather than a polar covalent one, in agreement with a recent publication of Rzepa [3] He2CO2+ is the only investigated compound which presents a C-He bonding ELF basin, indicating a non-negligible covalent contribution to the bond. Other modifications in the electronic structure, such as the breaking of the triple CC bond in ethyne derivatives (Figure 1) or the loss of aromaticity in C6H3He33+ are also nicely revealed by the ELF topology.
[1] D. Bellert, W. H. Breckenridge Chem. Rev., 2002, 102, 1595-1622.
[2] F. Grandinetti Int. J. of Mass Spect., 2004, 237, 243-267.
[3] S. H. Rzepa Nature Chem., 2010, 2, 390-393.
J. El Houaria a, C. Barberot a, JC. Boisson b, H. Khartabil a, S. Gérard a, G. Monard c and E. Hénon a
a ICMR, UMR CNRS 7312, University of Reims Champagne-Ardenne, Moulin de la Housse, Reims, 51687, France
b CReSTIC, University of Reims Champagne-Ardenne, Moulin de la Housse, Reims, 51687, France
c SRSMC, UMR CNRS 7565, University of Lorraine, Boulevard des Aiguillettes, Vandoeuvre les Nancy, 54506, France
Non Covalent Interactions in Molecular Docking : a quantum approach using the AlgoGen 2.3 software
Molecular docking of ligand to biological macromolecule receptors frequently relies on the use of standard force field methods and high-quality ligand force fields are needed to evaluate protein-ligand interaction energies. Linear scaling methodology using quantum semi-empirical Hamiltonians can now be envisaged to perform more energetically precise docking and to avoid the problem of defining force field. AlgoGen is a program that greatly benefits from these algorithmic advances [1]. It is able to perform rigid molecular docking in order to ultimately pose the ligand in the receptor site by combining the Divide and Conquer linear scaling methodology using a genetic algorithm (GA). A new version of this program is presented, interfaced with MOZYME.
Firstly, we tested the AlgoGen-MOZYME/PM6-DH+ engine by exploring the potential energy surface of 22 relative small complexes (containing up to 30 atoms). We selected the S22 set designed by Hobza et al.[4] to test new computational tools for biologically oriented applications. These complexes have been chosen so that to be representative of components of biomacromolecules that possess characteristic stabilizing contributions. It involves multiple possible non-covalent interactions: hydrogen bonding, dispersion, π- stacking (offset parallel, T-shape). Additionally, a biological application on the docking of the zardaverine inhibitor in the active site of PDE4D leads to a pose in good agreement with known crystallographic structure. Complementary post-docking treatments using NCI [2,3] and charge decomposition analysis turn out to give insight into chemical weak interactions at protein-ligand interface.
Last, preliminary results are presented concerning the NCI method [2,3] used to estimate interaction energies. This could be an interesting alternative to the quantum evaluation function in molecular docking since it is only based on promolecular densities calculations, which are less expensive than full wave function calculations and which can fully take advantage of grid-based techniques.
[1] E Thiriot et al., THEOCHEM 898:31, 2009
[2] ER Johnson et al., JACS 132: 6498, 2010
[3] J Contreras-García et al., JCTC 7: 625, 2011
[4] P Jurecka et al. PCCP, 8(17):1985, 2006