On the Covalent Character of Rare Gas Bonding Interactions: a New Kind of Weak Interaction” by Wenli Zou, Davood Nori-Shargh, and James E

Total Page:16

File Type:pdf, Size:1020Kb

On the Covalent Character of Rare Gas Bonding Interactions: a New Kind of Weak Interaction” by Wenli Zou, Davood Nori-Shargh, and James E Supporting Information to the article entitled “On the covalent character of rare gas bonding interactions: a new kind of weak interaction” by Wenli Zou, Davood Nori-Shargh, and James E. Boggs Table S1. Equilibrium geometries (bond length in Å and bond angle in degree) of some rare gas containing systems. All of them have C ∞v symmetry except HeB(CH) 3 with C 3v . No. Mol. State Structure parameters Ref. 1 FHeO - 1Σ+ He-O=1.100, He-F=1.621 61 2 HHeF 1Σ+ He-H=0.818, He-F=1.418 63 3 HArF 1Σ+ Ar-H=1.329, Ar-F=1.967 6 4 HeCuF 1Σ+ He-Cu=1.659, Cu-F=1.732 26 5 HeAgF 1Σ+ He-Ag=2.143, Ag-F=1.967 26 6 HeAuF 1Σ+ He-Au=1.841, Au-F=1.904 26 7 NeCuF 1Σ+ Ne-Cu=2.172, Cu-F=1.738 26 8 NeAgF 1Σ+ Ne-Ag=2.700, Ag-F=1.974 26 9 NeAuF 1Σ+ Ne-Au=2.444, Au-F=1.917 26 10 ArCuF a) 1Σ+ Ar-Cu=2.219, Cu-F=1.753 12 11 ArAgF a) 1Σ+ Ar-Ag=2.558, Ag-F=1.986 13 12 ArAuF a) 1Σ+ Ar-Au=2.391, Au-F=1.918 14 13 KrCuF a) 1Σ+ Kr-Cu=2.316, Cu-F=1.745 11 14 KrAgF a) 1Σ+ Kr-Ag=2.594, Ag-F=1.984 b) 10 15 KrAuF a) 1Σ+ Kr-Au=2.460, Au-F=1.918 10 16 XeCuF a) 1Σ+ Xe-Cu=2.430, Cu-F=1.745 9 17 XeAgF a) 1Σ+ Xe-Ag=2.666, Ag-F=1.983 8 18 XeAuF a) 1Σ+ Xe-Au=2.543, Au-F=1.918 7 19.a HePtF 2Σ+ He-Pt=1.828, Pt-F=1.881 31 19.b 2Π He-Pt=1.860, Pt-F=1.862 19.c 2∆ He-Pt=1.798, Pt-F=1.900 20 HePtXe 1Σ+ He-Pt=1.818, Xe-Pt=2.509 33 1 21 HeB(CH) 3 A1 He-B=1.355, B-C=1.493, C-H=1.066, 32 C-B-He=143.7, H-C-B=161.4 22 HeBeO 1Σ+ He-Be=1.522, Be-O=1.330 25 23 NeBeO 1Σ+ Ne-Be=1.798, Be-O=1.339 37 24 ArBeO 1Σ+ Ar-Be=2.076, Be-O=1.341 37 25 KrBeO 1Σ+ Kr-Be=2.211, Be-O=1.343 37 26 XeBeO 1Σ+ Xe-Be=2.385, Be-O=1.344 37 a) Experimental geometries. b) R 0 is used here. 1.957 Å in ref 10 may be a typo. Natural charge and WBI. It has been reported that natural charge and WBI may have problems in the case of some metal compounds. S1-S4 For example, the WBI value of the Be-O bond is only 0.6, being much smaller than the Mayer bond order (MBO) of 1.8 and the empirical bond order of 2.0. This problem also exists in MgO, CaO, ZnO, and so on. For rare gas bonds, however, NBO analysis provides reasonable results, and the WBI values are qualitatively in line with the MBO ones. So NBO S1 analysis can be used safely in this research. References (S1) Maseras, F.; Morokuma, K. Chem. Phys. Lett. 1992 , 195, 500-504. (S2) Clark, A. E.; Sonnenberg, J. L.; Hay, P. J.; Martin, R. L. J. Chem. Phys. 2004 , 121, 2563-2570. (S3) Landis, C. R.; Weinhold, F. J. Comput. Chem. 2007 , 28, 198-203. (S4) Lu, T.; Chen, F.-W. Acta Phys.-Chim. Sin. 2012 , 28, 1-18. S2 Programs used in this research Molden2AIM It converts the data format from Molden to AIM's WFN. The latter format can be used in AIMPAC, AIM2000, AIMALL, PAMoC, TopMoD, Xaim, MORPHY98, Multiwfn, ProMolDen, and so on. Molden file can be generated using the following quantum chemistry softwares, but some of them do not work with Molden2AIM. 1. MOLPRO 2. MOLCAS (spherical basis functions only) 3. deMon2k 4. Q-Chem (spherical S, P, D, and F or Cartesian S, P, and D basis functions) 5. CFour (S, P, D, and F basis functions; insert a command line [Program] Cfour into the Molden file) 6. Turbomole (insert a command line [Program] Turbomole into the Molden file) 7. ORCA (insert a command line [Program] Orca into the Molden file) 8. Columbus (thanks to Dr. Marat Talipov for testing) 9. Priroda (thanks to Dr. Evgeniy Pankratyev for testing) 10.StoBe (not tested) 11.TeraChem (not tested) 12.ACES II and III (not supported) 13.DALTON (not supported) 14.Jaguar (not supported) Limitations: 1. In general, only S, P, D, F, and G Gaussian basis functions are supported by Molden2AIM. AIMALL can also use H functions, but the Molden format cannot. 2. ECP or MCP is not supported because there is no ECP data block in the Molden format. See AIMALL for the analysis using ECP. 3. The total energy and virial ratio printed on the last line of the WFN file do not make sense. If they are used in your AIM analysis, please modify them manually. NBO2Molden It saves NBO's plot files (*.31 ~ *.40) into Molden format, which can be plotted using Molden, GabEdit, MacMolPlt, and so on. MolBO It reads basis sets, overlap matrix, Fork matrix, density matrix, ... from MOLPRO's output file (use GPRINT and MATROP commands to obtain these data; see the examples), and prepares NBO's *.47 file. It also calculates Mayer bond orders (MBO). Electronic density can be computed using the following theoretical methods: S3 1. SCF: RHF, ROHF, UHF, RDFT, RODFT, UDFT 2. Density Fitting (nosym): DF-HF, DF-KS 3. Post-HF: MP2, MP3, CCSD, QCISD, QCISD(T), EOM-CCSD, Full-CI 4. Local Post-HF (nosym): LMP2 5. MC/MR: MCSCF, CASVB, MRCI (including SRCI), CASPT2 (including SRMP2) Supports: 1. MOLPRO 2006 - 2010. Lower MOLPRO versions were not tested. 2. NBO 3.0 and 5.x. Version 4.x has not been tested, but there should be no problem. Limitations: 1. Symmetry equivalent atoms are not allowed. It is necessary to decrease the symmetry until all the atoms are symmetry unique. For example, use C 2v instead of D 2h for CO 2, C s instead of C 2v for H 2O, and C 1 instead of C s for CF 4. 2. For the all-electron relativistic calculations (DKHn), do not print kinetic energy and potential energy matrices because they are non-relativistic. Please read the manual for other questions. S4 .
Recommended publications
  • Aqueous Pka Prediction for Tautomerizable Compounds Using Equilibrium Bond Lengths
    ARTICLE https://doi.org/10.1038/s42004-020-0264-7 OPEN Aqueous pKa prediction for tautomerizable compounds using equilibrium bond lengths Beth A. Caine1,2, Maddalena Bronzato3, Torquil Fraser3, Nathan Kidley3, Christophe Dardonville 4 & ✉ Paul L.A. Popelier 1,2 1234567890():,; The accurate prediction of aqueous pKa values for tautomerizable compounds is a formidable task, even for the most established in silico tools. Empirical approaches often fall short due to a lack of pre-existing knowledge of dominant tautomeric forms. In a rigorous first-principles approach, calculations for low-energy tautomers must be performed in protonated and deprotonated forms, often both in gas and solvent phases, thus representing a significant computational task. Here we report an alternative approach, predicting pKa values for her- bicide/therapeutic derivatives of 1,3-cyclohexanedione and 1,3-cyclopentanedione to within just 0.24 units. A model, using a single ab initio bond length from one protonation state, is as accurate as other more complex regression approaches using more input features, and outperforms the program Marvin. Our approach can be used for other tautomerizable spe- cies, to predict trends across congeneric series and to correct experimental pKa values. 1 Department of Chemistry, University of Manchester, Manchester, UK. 2 Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester, UK. 3 Syngenta AG, Jealott’s Hill, Warfield, Bracknell RG42 6E7, UK. 4 Instituto de Química Médica, IQM–CSIC, C/Juan de la Cierva 3, Madrid 28006, Spain. ✉ email: [email protected] COMMUNICATIONS CHEMISTRY | (2020) 3:21 | https://doi.org/10.1038/s42004-020-0264-7 | www.nature.com/commschem 1 ARTICLE COMMUNICATIONS CHEMISTRY | https://doi.org/10.1038/s42004-020-0264-7 pproximately 21% of the compounds that make up (http://vcclab.org), Marvin (http://www.chemaxon.com) and pharmaceutical databases are said to exist in two or more Pallas (www.compudrug.com)) on 248 compounds of the Gold A 1 12 tautomeric forms .
    [Show full text]
  • DFT and QTAIM Study of the Tetra-Tert ... -.:. Michael Pittelkow
    Journal of Molecular Structure 1026 (2012) 127–132 Contents lists available at SciVerse ScienceDirect Journal of Molecular Structure journal homepage: www.elsevier.com/locate/molstruc DFT and QTAIM study of the tetra-tert-butyltetraoxa[8]circulene regioisomers structure ⇑ Gleb V. Baryshnikov a, Boris F. Minaev a, , Valentina A. Minaeva a,b, Alina T. Baryshnikova a, Michael Pittelkow c a Bohdan Khmelnytsky National University, 18031 Cherkasy, Ukraine b Theoretical Chemistry, School of Biotechnology, Royal Institute of Technology, SE-10691 Stockholm, Sweden c Department of Chemistry, University of Copenhagen, Universitetsparken 5, DK-2100 Copenhagen Ø, Denmark highlights " Tetra-tert-butyltetraoxa[8]circulene regioisomers were studied by DFT method. " Electronic density distribution was calculated by the QTAIM method. " The presence of stabilizing non-valence bonds is detected by X-ray experiment. " The HÁÁÁH contacts are dynamically unstable due to high ellipticity. " The energy of the HÁÁÁH and CHÁÁÁO contacts was estimated by the Espinosa equation. article info abstract Article history: The recently synthesized tetra-tert-butyltetraoxa[8]circulene regioisomers characterized by unusual Received 6 March 2012 solution-state aggregation behavior are calculated at the density functional theory (DFT) level with the Received in revised form 24 May 2012 quantum theory of atoms in molecules (QTAIMs) approach to the electron density distribution analysis. Accepted 24 May 2012 The presence of stabilizing intramolecular hydrogen bonds and hydrogen–hydrogen interactions in the Available online 31 May 2012 studied molecules is predicted and the energies of these interactions are estimated with QTAIM. Occur- rence of the CHÁÁÁO bonds is detected by the single-crystal X-ray analysis for two regioisomers, obtained Keywords: in high purity.
    [Show full text]
  • Formation of Β-Cyclodextrin Complexes in an Anhydrous Environment
    JMolModel (2016) 22:207 DOI 10.1007/s00894-016-3061-6 ORIGINAL PAPER Formation of β-cyclodextrin complexes in an anhydrous environment Hocine Sifaoui1 & Ali Modarressi2 & Pierre Magri2 & Anna Stachowicz-Kuśnierz3 & Jacek Korchowiec3 & Marek Rogalski2 Received: 16 December 2015 /Accepted: 3 July 2016 # The Author(s) 2016. This article is published with open access at Springerlink.com Abstract The formation of inclusion complexes of β- Keywords β-Cyclodextrin . Differential scanning cyclodextrin was studied at the melting temperature of guest calorimetry . Inclusion complexes . Polyaromatics . Organic compounds by differential scanning calorimetry. The com- acids . Molecular modeling plexes of long-chain n-alkanes, polyaromatics, and organic acids were investigated by calorimetry and IR spectroscopy. The complexation ratio of β-cyclodextrin was compared with results obtained in an aqueous environment. The stability and Introduction structure of inclusion complexes with various stoichiometries were estimated by quantum chemistry and molecular dynam- Many lipophilic drug molecules display low bioavailability ics calculations. Comparison of experimental and theoretical that hinder their efficacy [1]. Cyclodextrins, which form stable results confirmed the possible formation of multiple inclusion complexes with numerous molecules, are used to improve complexes with guest molecules capable of forming hydrogen aqueous solubility and bioavailability. β-Cyclodextrin (β- bonds. This finding gives new insight into the mechanism of CD) is composed of seven glucopyranose units forming a formation of host–guest complexes and shows that hydropho- cyclic, cone-shaped cavity with a hydrophilic outer surface bic interactions play a secondary role in this case. and a relatively hydrophobic inner surface [2]. This molecular structure of β-CD allows the formation of complexes with a large variety of organic molecules of different shape and po- larity [3–8].
    [Show full text]
  • Predikce Hodnot Pka Na Základě EEM Atomových Nábojů
    MASARYKOVA UNIVERZITA PRˇ I´RODOVEˇ DECKA´ FAKULTA NA´ RODNI´ CENTRUM PRO VY´ZKUM BIOMOLEKUL Diplomova´pra´ce BRNO 2013 STANISLAV GEIDL MASARYKOVA UNIVERZITA PRˇ I´RODOVEˇ DECKA´ FAKULTA NA´ RODNI´ CENTRUM PRO VY´ZKUM BIOMOLEKUL Predikce hodnot pKa na za´kladeˇ EEM atomovy´ch na´boju˚ Diplomova´pra´ce Stanislav Geidl Vedoucı´pra´ce: prof. RNDr. Jaroslav Kocˇa, DrSc. Konzultant: RNDr. Radka Svobodova´Varˇekova´, Ph.D. Brno 2013 Bibliograficky´za´znam Autor: Bc. Stanislav Geidl Prˇı´rodoveˇdecka´fakulta, Masarykova univerzita Na´rodnı´centrum pro vy´zkum biomolekul Na´zev pra´ce: Predikce hodnot pKa na za´kladeˇEEM atomovy´ch na´boju˚ Studijnı´program: Biochemie Studijnı´obor: Chemoinformatika a bioinformatika Vedoucı´pra´ce: prof. RNDr. Jaroslav Kocˇa, DrSc. Akademicky´rok: 2012/2013 Pocˇet stran: xi + 78 Klı´cˇova´slova: disociacˇnı´konstanta, pKa, QSPR, vı´cerozmeˇrna´linea´rnı´ regrese, atomove´ na´boje, kvantova´ mechanika, EEM, fenoly, karboxylove´kyseliny Bibliographic Entry Author: Bc. Stanislav Geidl Faculty of Science, Masaryk University National Centre for Biomolecular Research Title of Thesis: Predicting pKa values from EEM atomic charges Degree Programme: Biochemistry Field of Study: Chemoinformatics and Bioinformatics Supervisor: prof. RNDr. Jaroslav Kocˇa, DrSc. Academic Year: 2012/2013 Number of Pages: xi + 78 Keywords: dissociation constant, pKa, QSPR, multilinear regres- sion, partial atomic charge, quantum mechanics, EEM, phenols, carboxilic acids Abstrakt Disociacˇnı´ konstanta pKa je velmi du˚lezˇitou vlastnostı´ molekuly, a proto je vy´voj spolehlivy´ch a rychly´ch metod pro predikci pKa jednou z klı´cˇovy´ch oblastı´ vy´zkumu. Tato pra´ce zjisˇt’uje, jestli lze predikovat pKa pomocı´QSPR modelu˚ vyuzˇı´vajı´cı´ch empiricke´atomove´na´boje, konkre´tneˇna´boje vypocˇı´tane´metodou EEM (Electronegativity Equalization Method).
    [Show full text]
  • Multiwfn - a Multifunctional Wavefunction Analyzer
    Multiwfn - A Multifunctional Wavefunction Analyzer - Software Manual with abundant tutorials and examples in Chapter 4 Version 3.3.8 2015-Dec-1 http://Multiwfn.codeplex.com Tian Lu [email protected] Beijing Kein Research Center for Natural Sciences (www.keinsci.com) !!!!!!!!!! ALL USERS MUST READ !!!!!!!!!! I know, most people, including myself, are unwilling to read lengthy manual. Since Multiwfn is a heuristic and very user-friendly program, it is absolutely unnecessary to read through the whole manual before using it. However, you should never skip reading following content! 1. The BEST way to get started quickly is directly reading Chapter 1 and the tutorials in Chapter 4. After that if you want to learn more about Multiwfn, then read Chapters 2 and 3. Note that the tutorials and examples given in Chapter 4 only cover the most frequently used functions of Multiwfn. Some application skills are described in Chapter 5, which may be useful for you. 2. Different functions of Multiwfn require different type of input file, please read Section 2.5 for explanation. 3. If you do not know how to copy the output of Multiwfn from command-line window to a plain text file, consult Section 5.4. If you do not know how to enlarge screen buffer size of command-line window of Windows system, consult Section 5.5. 4. If the error “No executable for file l1.exe” appears in screen when Multiwfn is invoking Gaussian, you should set up Gaussian environment variable first. For Windows version, you can refer Appendix 1. (Note: Most functions in Multiwfn DO NOT require Gaussian installed on your local machine) 5.
    [Show full text]
  • Solving the Problem of Aqueous Pka Prediction for Tautomerizable Compounds Using Equilibrium Bond Lengths
    Solving the Problem of Aqueous pKa Prediction for Tautomerizable Compounds Using Equilibrium Bond Lengths Beth A. Cainea,b, Maddalena Bronzatoc, Torquil Fraserc, Nathan Kidleyc, Christophe Dardonvilled and Paul L. A. Popeliera,b aSchool of Chemistry, University of Manchester, Great Britain, bManchester Institute of Biotechnology (MIB), 131 Princess Street, Great Britain, c Syngenta AG, Jealott’s Hill, Warfield, Bracknell, RG42 6E7, Great Britain d Instituto de Química Médica, IQM–CSIC, C/ Juan de la Cierva 3, 28006 Madrid, Spain The accurate prediction of aqueous pKa values for tautomerizable compounds is a formidable task, even for the most established in silico tools. Empirical approaches often fall short due to a lack of pre-existing knowledge of dominant tautomeric forms. In a rigorous first-principles approach, calculations for low- energy tautomers must be performed, in protonated and deprotonated forms, both in gas and solvent phase, thus representing a significant computational task. Here we report an alternative approach, predicting pKa values for herbicide/therapeutic derivatives of 1,3-cyclohexanedione and 1,3- cyclopentanedione to within just 0.24 units. A model, with as input feature a single ab initio bond length from one protonation state, is as accurate as other, more complex machine learning approaches (SVR, RFR, GPR, PLS) using more input features, and outperforms the program Marvin. Our approach can be used for other tautomerizable species, to predict trends across congeneric series and to correct experimental pKa values. 1 Approximately 21% of the compounds that make up pharmaceutical databases are said to exist in two or more tautomeric forms1. Tautomerism is a form of structural isomerism that is characterized by a species having two or more structural representations, between which interconversion can be achieved by “proton hopping” from one atom to another.
    [Show full text]
  • New Tools for Chemical Bonding Analysis
    New tools for chemical bonding analysis Eduard Matito February 19, 2019 Contents 1 Theoretical Background 4 1.1 TheElectronDensity......................... 5 1.2 ThePairDensity ........................... 7 1.3 ElectronCorrelation . 11 2 The Atom in the Molecule 13 2.1 TheOriginoftheAtom ....................... 13 2.2 ModernTheoryoftheAtom. 15 2.3 Thedefinitionofanatominamolecule . 15 2.4 Hilbertspacepartition . 16 2.5 Realspacepartition . .. .. .. .. .. .. .. .. .. .. 17 2.5.1 The quantum theory of atoms in molecules (QTAIM) . 18 2.5.2 Otherreal-spacepartitions . 26 2.5.3 References........................... 28 2.6 Populationanalysis. 28 2.6.1 Mullikenpopulationanalysis . 28 2.6.2 Populationfromrealspacepartitions. 29 2.6.3 The electron-sharing indices (bond orders) . 30 2.6.4 MulticenterIndices. 32 3 Aromaticity 38 3.1 GlobalandLocalAromaticity. 40 3.2 Geometricalindices. 40 3.3 Magneticindices ........................... 41 3.4 Electronicindices. .. .. .. .. .. .. .. .. .. .. .. 42 3.4.1 The Aromatic Fluctuation Index: FLU . 43 3.4.2 A Multicenter based index: Iring .............. 43 3.4.3 TheMulticenterIndex: MCI . 44 3.4.4 Thepara-DelocalizationIndex: PDI . 44 3.4.5 TheBond-OrderAlternation . 44 3.4.6 EXERCISE8(Aromaticity). 45 4 Oxidation state 51 4.1 Definition ............................... 51 4.2 Computational calculation of Oxidation States . ... 52 4.2.1 BondValenceSum . .. .. .. .. .. .. .. .. .. 52 4.2.2 AtomPopulationAnalysis. 53 4.2.3 SpinDensities......................... 53 1 4.2.4 Orbital-localization Methods: Effective Oxidation State . 56 5 The Electron Localization Function (ELF) 59 5.1 The topological analysis of the electron density . .... 60 6 The local Spin 63 7 Appendices 64 7.1 AppendixI:Howtoperformcalculations. 64 7.2 AppendixII:ManualofESI-3D . 65 7.2.1 Howtocitetheprogram. 65 7.2.2 Compulsorykeywords . 66 7.3 Optionalkeywords .........................
    [Show full text]
  • WHAT INFLUENCE WOULD a CLOUD BASED SEMANTIC LABORATORY NOTEBOOK HAVE on the DIGITISATION and MANAGEMENT of SCIENTIFIC RESEARCH? by Samantha Kanza
    UNIVERSITY OF SOUTHAMPTON Faculty of Physical Sciences and Engineering School of Electronics and Computer Science What Influence would a Cloud Based Semantic Laboratory Notebook have on the Digitisation and Management of Scientific Research? by Samantha Kanza Thesis for the degree of Doctor of Philosophy 25th April 2018 UNIVERSITY OF SOUTHAMPTON ABSTRACT FACULTY OF PHYSICAL SCIENCES AND ENGINEERING SCHOOL OF ELECTRONICS AND COMPUTER SCIENCE Doctor of Philosophy WHAT INFLUENCE WOULD A CLOUD BASED SEMANTIC LABORATORY NOTEBOOK HAVE ON THE DIGITISATION AND MANAGEMENT OF SCIENTIFIC RESEARCH? by Samantha Kanza Electronic laboratory notebooks (ELNs) have been studied by the chemistry research community over the last two decades as a step towards a paper-free laboratory; sim- ilar work has also taken place in other laboratory science domains. However, despite the many available ELN platforms, their uptake in both the academic and commercial worlds remains limited. This thesis describes an investigation into the current ELN landscape, and its relationship with the requirements of laboratory scientists. Market and literature research was conducted around available ELN offerings to characterise their commonly incorporated features. Previous studies of laboratory scientists examined note-taking and record-keeping behaviours in laboratory environments; to complement and extend this, a series of user studies were conducted as part of this thesis, drawing upon the techniques of user-centred design, ethnography, and collaboration with domain experts. These user studies, combined with the characterisation of existing ELN features, in- formed the requirements and design of a proposed ELN environment which aims to bridge the gap between scientists' current practice using paper lab notebooks, and the necessity of publishing their results electronically, at any stage of the experiment life cycle.
    [Show full text]
  • Computational Chemistry of Non-Covalent Interaction and Its Application in Chemical Catalysis
    Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 5-2016 Computational Chemistry of Non-Covalent Interaction and its Application in Chemical Catalysis Vincent de Paul Nzuwah Nziko Utah State University Follow this and additional works at: https://digitalcommons.usu.edu/etd Part of the Chemistry Commons Recommended Citation Nziko, Vincent de Paul Nzuwah, "Computational Chemistry of Non-Covalent Interaction and its Application in Chemical Catalysis" (2016). All Graduate Theses and Dissertations. 5248. https://digitalcommons.usu.edu/etd/5248 This Dissertation is brought to you for free and open access by the Graduate Studies at DigitalCommons@USU. It has been accepted for inclusion in All Graduate Theses and Dissertations by an authorized administrator of DigitalCommons@USU. For more information, please contact [email protected]. COMPUTATIONAL CHEMISTRY OF NON-COVALENT INTERACTION AND ITS APPLICATION IN CHEMICAL CATALYSIS by Vincent de Paul Nzuwah Nziko A dissertation submitted in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY in Chemistry Approved: Steve Scheiner Alexander I. Boldyrev Major Professor Committee Member Alvan Hengge T.C. Shen Committee Member Committee Member Cheng-Wei Tom Chang Mark R. McLellan Committee Member Vice President for Research and Dean of the School of Graduate Studies UTAH STATE UNIVERSITY Logan, Utah 2016 ii Copyright © Vincent de Paul Nzuwah Nziko 2016 All Rights Reserved iii ABSTRACT Computational Chemistry of Non-Covalent Interaction and its Application in Chemical Catalysis by Vincent de Paul Nzuwah Nziko, Doctor of Philosophy Utah State University, 2016 Major Professor: Dr. Steve Scheiner Department: Chemistry and Biochemistry Unconventional non-covalent interactions such as halogen, chalcogen, and tetrel bonds are gaining interest in several domains including but not limited to drug design, as well as novel catalyst design.
    [Show full text]
  • An Exploration of the Ozone Dimer Potential Energy Surface
    An Exploration of the Ozone Dimer Potential Energy Surface Luis Miguel Azofra, † Ibon Alkorta † and Steve Scheiner ‡, * †Instituto de Química Médica, CSIC, Juan de la Cierva, 3, E-28006, Madrid, Spain ‡Department of Chemistry and Biochemistry, Utah State University, Logan, UT 84322-0300, USA *Author to whom correspondence should be addressed. Fax: (+1) 435-797-3390 E-mail: [email protected] ABSTRACT The (O 3)2 dimer potential energy surface is thoroughly explored at the ab initio CCSD(T) computational level. Five minima are characterized with binding energies between 0.35 and 2.24 kcal/mol. The most stable may be characterized as slipped parallel, with the two O 3 monomers situated in parallel planes. Partitioning of the interaction energy points to dispersion and exchange as the prime contributors to the stability, with varying contributions from electrostatic energy, which is repulsive in one case. Atoms in Molecules (AIM) analysis of the wavefunction presents specific O···O bonding interactions, whose number is related to the overall stability of each dimer. All internal vibrational frequencies are shifted to the red by dimerization, particularly the antisymmetric stretching mode whose shift is as high as 111 cm –1. In addition to the five minima, 11 other higher-order stationary points are identified. KEYWORDS: chalcogen bonds; weak interactions; dispersion; frequency shifts; AIM 1 INTRODUCTION The paradigmatic molecule of ozone (O 3) was the first allotrope of a chemical element to be recognized. Its composition was determined in 1865 by Soret,[1] and subsequently, in 1867, confirmed by Schönbein.[2] O3 is found in several atmospheric layers, especially in the stratosphere, where it is most concentrated.
    [Show full text]
  • Rhorix: an Interface Between Quantum Chemical Topology and the 3D Graphics Program Blender
    Lawrence Berkeley National Laboratory Recent Work Title Rhorix: An interface between quantum chemical topology and the 3D graphics program blender. Permalink https://escholarship.org/uc/item/18n454gq Journal Journal of computational chemistry, 38(29) ISSN 0192-8651 Authors Mills, Matthew JL Sale, Kenneth L Simmons, Blake A et al. Publication Date 2017-11-01 DOI 10.1002/jcc.25054 Peer reviewed eScholarship.org Powered by the California Digital Library University of California SOFTWARE NEWS AND UPDATES WWW.C-CHEM.ORG Rhorix: An Interface between Quantum Chemical Topology and the 3D Graphics Program Blender Matthew J. L. Mills ,[a,b] Kenneth L. Sale,[a,b] Blake A. Simmons,[a,c] and Paul L. A. Popelier [d] Chemical research is assisted by the creation of visual repre- topological objects and their 3D representations. Possible sentations that map concepts (such as atoms and bonds) to mappings are discussed and a canonical example is suggested, 3D objects. These concepts are rooted in chemical theory that which has been implemented as a Python “Add-On” named predates routine solution of the Schrodinger€ equation for sys- Rhorix for the state-of-the-art 3D modeling program Blender. tems of interesting size. The method of Quantum Chemical This allows chemists to use modern drawing tools and artists Topology (QCT) provides an alternative, parameter-free means to access QCT data in a familiar context. A number of exam- to understand chemical phenomena directly from quantum ples are discussed. VC 2017 The Authors. Journal of Computa- mechanical principles. Representation of the topological ele- tional Chemistry Published by Wiley Periodicals, Inc.
    [Show full text]
  • Laboratory for Molecular Simulation Subscription Rates for TAMU Researchers FY20
    Laboratory for Molecular Simulation Subscription Rates for TAMU Researchers FY20 Subscription Number of Users with Access to LMS Consulting Cost Type Software and Hardware Hours Included per year 1 PI & 1 researcher 5 $500 2 PI & 2 researchers 10 $1000 3 PI & 3 researchers 15 $1500 4 PI & 4+ researchers 20 $2000 Consulting Fee for non-subscribers and additional hours for subscribers: $65/hour In addition to the consulting hours, an annual subscription to the LMS gives the PI and designated researcher(s) access to the licensed software and hardware outlined below. Licensed Software Company Software BIOVIA Materials Studio Visualizer and the following modules: Conformers, Forcite Plus Parallel, Gaussian Interface, QSAR+, Reflex, VAMP, MS Pipeline Pilot Collection, Adsorption Locator, Amorphous Cell, Blends, Compass, GULP, Mesocite, Mesodyn, Sorption, Synthia, CASTEP, DFTB+, DMOL3, NMR CASTEP, ONETEP, QMERA BIOVIA Discovery Studio Visualizer and the following modules: Analysis, Biopolymer, Catalyst Conformation, Catalyst Score, CHARMm, DMOL3 Molecular, MMFF (Merk Molecular Force Field), Protein Refine, QUANTUMm (QM/MM - DMOL3/CHARMm), Catalyst DB Build, Catalyst DB Search, Catalyst Hypothesis, Catalyst SBP, Catalyst Shape, CFF (Consistent Force-Field), De Novo Evolution, De Novo Ligand Builder, Flexible Docking, LibDock, LigandFit, LigandScore, Ludi, MCSS (Multiple Copy Simultaneous Search), MODELER, Protein Docking: ZDOCK and RDOCK, Protein Families, Protein Health, Sequence Analysis, X-ray. Schrödinger Schrödinger suite of software: Maestro, CombiGlide, Glide, Liaison, Strike, QikProp, Canvas, LigPrep, BioLuminate GUI, Prime, Qsite, MacroModel, ConfGen, Jaguar, pKa Predictor, Epik, SiteMap, and PIPER. CCG MOE: Molecular Operating Environment, One fully integrated drug discovery software package, including structure-based design, fragment-based design, pharmacophore discovery, medicinal and biologics applications, protein and antibody modeling, molecular mechanics/dynamics, cheminformatics and QSAR.
    [Show full text]