Electron Density of Delocalized Bonds As a Universal Tool for Assessing Global and Local Effects of Chemical Resonance

Total Page:16

File Type:pdf, Size:1020Kb

Electron Density of Delocalized Bonds As a Universal Tool for Assessing Global and Local Effects of Chemical Resonance Summary of Professional Accomplishments Electron density of delocalized bonds as a universal tool for assessing global and local effects of chemical resonance Dr. Dariusz Wojciech Szczepanik Kraków 2021 1. Name Name and Surname: Dariusz Wojciech Szczepanik Degree: Doctor of Philosophy in Chemistry Current employment: Department of Theoretical Chemistry, Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland. ORCID: 0000-0002-2013-0617 ResearcherID: E-2787-2014 ScopusID: 36835515900 Personal website: http://www.eddb.pl/aboutme Presentation of the achievements: http://www.eddb.pl/hab 2. Diplomas, degrees conferred in specific areas of science or arts, including the name of the institution which conferred the degree, year of degree conferment, title of the PhD dissertation 2008 Master of Science Department of Computational Methods in Chemistry, Faculty of Chemistry UJ Thesis: „Entropic indices of the chemical bonds from information theory” Supervisor: dr. hab. Janusz Mrozek 2013 Doctor of Philosophy in Chemistry Department of Theoretical Chemistry, Faculty of Chemistry UJ Thesis: „Probabilistic models of the chemical bond in function spaces” (de- fended with honors) Supervisor: dr. hab. Janusz Mrozek 3. Information on employment in research institutes or faculties/departments or school of arts 2015 – 2020 Department of Theoretical Chemistry, Faculty of Chemistry, Jagiellonian Uni- versity. Position: Technician 2018 – 2020 Institute of Computational Chemistry and Catalysis, University of Girona, Posi- tion: EU-researcher (MSCA-IF, postdoc). 2020 – Department of Theoretical Chemistry, Faculty of Chemistry, Jagiellonian Uni- versity. Position: Adjunct 4. Description of the achievements, set out in art. 219 para 1 point 2 of the Act The basis for the scientific achievement entitled “Electron density of delocalized bonds as a universal tool for assessing global and local effects of chemical resonance” is a series of thematically linked publications composed of 9 scientific articles (H1-H9) published in 2014-2019 in peer-reviewed jour- nals from the JCR list, and 1 monographic chapter (H10) of 2021. The subject of scientific articles 2 selected for the habilitation cycle (H1-H9) concerns development of the theoretical basis of my orig- inal computational method that enables one to “extract” from the total electron density of a molecule (obtained from quantum-chemical calculations) the so-called electron density of delocalized bonds (EDDB), which allows visualization and quantification of different effects of electron delocalization in chemical species, regardless of their size, topology and the electronic state. Publication (H10), in turn, can be regarded as a review article. All publications included in the habilitation cycle are based on my own original research ideas, and in all of these works I am the first and the corresponding author. I am the sole author of H3, H4, and H8, while in the case of H9 and H10 the second co-author is Prof. Miquel Solà, at whose invitation these works were prepared. The high percentage estimate of my contribution to the remaining articles is guided by the criterion of the amount of time spent plan- ning research study, writing software and scripts automating calculations (partially carried out by PhD students), elaboration of the results, preparation and submission of the manuscript, improving them in the review process, etc. In addition, I am the sole author of the original model of the "migrat- ing π-cycles" introduced in H7 (and being a generalization of the well-known in the literature model of “the migrating Clar’s π-sextet”). 4.1. Motivation1 It is commonly accepted that electron delocalization in aromatic rings is linked with unusual thermo- dynamic stability by means of the π-electron bookkeeping rules, like “4n+2” and “4n” (depending on the system topology and multiplicity), known from the age-old chemistry textbooks. Although these qualitative criteria of aromaticity (and antiaromaticity) adequately relate topology, symmetry, and degeneracy of molecular orbitals (MO) in the [n]annulene-like systems predominated by covalent resonance forms at their singlet or the lowest-lying triplet states, their use in a more general context regarding topologically diversified poly- and macrocyclic species (e.g. expanded porphyrins predom- inated by ionic forms), non-Kekulé molecules (e.g. radicals), etc., is not well-founded. Over the last decades an overwhelming number of quantitative ‘measures’ of aromaticity has been proposed in the literature, based on energetic, structural, magnetic, and electronic properties of molecules, thus providing a far more accurate account of aromatic stabilization than the electron-bookkeeping crite- ria. The most commonly used quantitative criteria of aromaticity within each of these groups are: 1) the aromatic stabilization energy (ASE), which is an energetic measure of π-aromaticity that emanates from the theory of valence bonds, and it can be evaluated by means of thermodynamic data for iso- desmic, (hyper)homodesmotic or isomerization reactions (also, a multitude of schemes can be found in the literature that allows one to efficiently estimate the aromatic stabilization energies for a specific class of aromatic species); 2) the harmonic-oscillator model of aromaticity (HOMA), which is a π- aromaticity index based on structural properties of molecules (being a normalized measure of devia- tions of bond lengths in aromatic molecule from the corresponding optimum bond lengths in an ide- alized non-aromatic molecule as a reference) – HOMA is close to 0 for nonaromatic species, ap- proaches 1 for highly aromatic ones, while for potentially antiaromatic rings it usually assumes neg- 1 Formal definitions, descriptions and scientific arguments adduced in this section comes from review articles in the spe- cial issue Chem. Soc. Rev 44 (2015), original publications H5, H7, and H10, as well as the following papers: A. Stanger, Chem. Comm. 2009 (2009) 1939; R. Hoffmann, Am. Sci. 103 (2015) 18; M. Solà, Front. Chem. 5 (2017) 22. 3 ative values; 3) the nucleus-independent chemical shift (NICS), which quantifies the effective mag- netic shielding at the centroid (or above) of the aromatic ring in external magnetic field – the more negative (positive) value of NICS, the more aromatic (antiaromatic) is the molecular ring in question; 4) the multicenter index (MCI), which is a non-reference index of aromaticity that can be calculated from both the ab initio molecular wave function as well as the electron density; MCI has been shown to be superior to other aromaticity descriptors as the only one that passes a set of rigorous tests for aromaticity quantifiers designed by Prof. M. Solà (University of Girona, Spain). Unfortunately, each of the above-mentioned aromaticity ‘measures’ has shortcomings and lim- itations that sometimes may lead to wrong predictions. The ASE seems to outwardly be the most adequate measure of global aromaticity since it can be evaluated by means of thermodynamical data and it emerges for the direct relationship between structural consequences of electron delocalization and the stability. However, designing of isodesmic and homodesmotic reaction scenarios is very dif- ficult in practice and opens the door to a lot of arbitrariness. In contrast, an unquestionable advantage of HOMA is its computational and interpretative simplicity – it allows one to straightforwardly clas- sify any molecular ring as aromatic, non-aromatic or (potentially) antiaromatic. Unfortunately, the principal problem with HOMA is the necessity of parametrization of bond lengths for an idealized reference molecule, which obviously cannot be chosen unambiguously. Consequently, the practical use of HOMA is limited to aromatic and heteroaromatic systems since the parameters for chemical bonds with metal atoms are not available. Furthermore, parametrization of HOMA should be per- formed using exactly the same quantum-chemical method as used in calculations of equilibrium ge- ometries of the molecule under study, since routinely computed HOMA with the experimentally de- termined parameters is bound to suffer from large unsystematic errors and strong sensitivity to the choice of the basis set and the exchange-correlation functional in the density functional theory (DFT) based calculations. Aromaticity descriptors based on magnetic properties of molecules, especially NICS(0), NICS(1), and its axial component, NICS(1)zz, dissected NICS, etc., are one the most pref- erable measures of local aromaticity due to their relation to experiment; diatropic (aromatic) and paratropic (antiaromatic) ring currents indirectly manifest itself in the NMR spectra. However, the magnetic-based measures of aromaticity have also come under bitter criticism due to their complexity (NICS relies on the condensation of potentially complicated patterns of induced currents to a single number), and methodological shortcomings. Finally, the main disadvantage connected with the cal- culation of MCI is its numerical instability and computational cost which prevent using MCI to ana- lyze systems containing more than 12-14 atoms; additionally, MCI suffers from the method depend- ence – in particular, the exchange-correlation functional selection at the density functional theory (DFT) level. To summarize, most of the aromaticity descriptors suffer from serious issues such as the arbi- trariness of choice of a reference system, lack of parametrization, ring-size extensivity issue, limited applicability,
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]
  • 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.
    [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]