DOCSLIB.ORG

Accurate Molecular Polarizabilities with Coupled Cluster Theory And

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Accurate Molecular Polarizabilities with Coupled Cluster Theory And Accurate molecular polarizabilities with coupled cluster theory and machine learning David M. Wilkinsa, Andrea Grisafia, Yang Yangb, Ka Un Laob, Robert A. DiStasio Jr.b,1, and Michele Ceriottia,1 aLaboratory of Computational Science and Modeling, Institut des Materiaux,´ Ecole´ Polytechnique Fed´ erale´ de Lausanne, 1015 Lausanne, Switzerland; and bDepartment of Chemistry and Chemical Biology, Cornell University, Ithaca, NY 14853 Edited by Michael L. Klein, Institute of Computational Molecular Science, Temple University, Philadelphia, PA, and approved December 26, 2018 (received for review September 23, 2018) The molecular dipole polarizability describes the tendency of a junction with a sufficiently large (diffuse) one-particle basis set molecule to change its dipole moment in response to an applied (20–23). However, such a prediction is accompanied by a sub- electric field. This quantity governs key intra- and intermolecu- stantially larger computational cost (scaling with the sixth power lar interactions, such as induction and dispersion; plays a vital of the system size), which can become quite prohibitive even role in determining the spectroscopic signatures of molecules; and when treating molecules with as few as 10 − 15 atoms. is an essential ingredient in polarizable force fields. Compared In the last few years, machine learning (ML) has gained trac- with other ground-state properties, an accurate prediction of tion as an alternative approach to the prediction of molecular the molecular polarizability is considerably more difficult, as this properties, substituting or complementing electronic structure response quantity is quite sensitive to the underlying electronic methods (24–26). In particular, it has been shown that accu- structure description. In this work, we present highly accurate racy on par with (or even better than) DFT can be achieved quantum mechanical calculations of the static dipole polarizabil- in the prediction of many molecular properties (27, 28) and ity tensors of 7,211 small organic molecules computed using that DFT (29) or coupled cluster (30) accuracy can be reached linear response coupled cluster singles and doubles theory (LR- more easily when using a less accurate but more computation- CCSD). Using a symmetry-adapted machine-learning approach, ally efficient electronic structure method as a stepping stone. we demonstrate that it is possible to predict the LR-CCSD molec- The polarizability, however, poses an additional challenge to ular polarizabilities of these small molecules with an error that ML. Due to its tensorial nature, the predicted α must trans- CHEMISTRY is an order of magnitude smaller than that of hybrid density form according to the symmetries of the SO(3) rotation group. functional theory (DFT) at a negligible computational cost. The For rigid molecules, this is easily achieved by learning the com- resultant model is robust and transferable, yielding molecular ponents of the tensor written in the reference frame of the polarizabilities for a diverse set of 52 larger molecules (includ- molecule (31, 32). However, to obtain a transferable model ing challenging conjugated systems, carbohydrates, small drugs, that would also be suitable for flexible molecules—as well as amino acids, nucleobases, and hydrocarbon isomers) at an accu- different compounds—this line of thought would require a cum- racy that exceeds that of hybrid DFT. The atom-centered decom- bersome and inelegant fragment decomposition. To avoid these position implicit in our machine-learning approach offers some complications, a symmetry-adapted Gaussian process regression insight into the shortcomings of DFT in the prediction of this fundamental quantity of interest. Significance dipole polarizability j machine learning j coupled cluster theory j The dipole polarizability of molecules and materials is cen- density functional theory j Gaussian process regression tral to several physical phenomena, modeling techniques, and the interpretation of many experiments. Its accurate evalua- he last decade has seen great progress in the first princi- tion from first principles requires quantum chemistry methods Tples evaluation of the structures, stabilities, and properties of that are often too demanding for routine use. The highly molecules and materials. Kohn–Sham density functional theory accurate calculations reported herein provide a much-needed (DFT) has played a pivotal role in this endeavor by providing benchmark of the accuracy of hybrid density functional theory ground-state properties with an accuracy that is sufficient for (DFT) as well as training data for a machine-learning model many useful applications at a manageable computational cost that can predict the polarizability tensor with an error that is (1–3). However, DFT is not equally accurate for every property about 50% smaller than DFT. This framework provides an accu- of interest. For instance, an accurate and reliable description of rate, inexpensive, and transferable strategy for estimating the the molecular dipole polarizability α, a tensor that describes how polarizabilities of molecules containing dozens of atoms, and the molecular dipole changes in the presence of an applied elec- therefore removes a considerable obstacle to accurate and tric field E, can be quite difficult to obtain (4). This is primarily reliable atomistic-based modeling of matter. due to the fact that α is a response property that is particularly sensitive to the quantum mechanical description of the underly- Author contributions: D.M.W., A.G., Y.Y., K.U.L., R.A.D., and M.C. designed research; D.M.W., A.G., Y.Y., and K.U.L. performed research; D.M.W., A.G., Y.Y., K.U.L., R.A.D., and ing electronic structure. As such, nontrivial electron correlation M.C. analyzed data; and D.M.W., A.G., Y.Y., K.U.L., R.A.D., and M.C. wrote the paper.y effects and basis set incompleteness error must be simultane- The authors declare no conflict of interest.y ously accounted for when determining α. For these reasons and This article is a PNAS Direct Submission.y in light of the fact that α is a fundamental quantity of interest that underlies induction and dispersion interactions (5–7) and Published under the PNAS license.y Raman and sum frequency generation spectroscopy (8–11), and Data deposition: Coordinates and dipole polarizabilities of the QM7b database, the showcase dataset, and the alkenes and acenes have been deposited in the Materials represents a key ingredient in the development of next genera- Cloud Archive (https://archive.materialscloud.org/2019.0002/v1).y tion polarizable force fields (12–16), it is important to provide 1 To whom correspondence may be addressed. Email: [email protected] or michele. benchmark values for α beyond the accuracy of DFT. In this ceriotti@epfl.ch.y regard, linear response coupled cluster singles and doubles the- This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10. ory (LR-CCSD) (17–19) has been shown to provide considerably 1073/pnas.1816132116/-/DCSupplemental.y more accurate and reliable predictions for α when used in con- www.pnas.org/cgi/doi/10.1073/pnas.1816132116 PNAS Latest Articles j 1 of 6 Downloaded by guest on September 23, 2021 (SA-GPR) scheme has recently been derived to naturally incor- LR-CCSD results, which will be referred to as DFT and cou- porate this SO(3) covariance into an ML scheme that is suitable pled cluster singles and doubles theory (CCSD), respectively, to predict tensorial quantities of arbitrary order (33). In this throughout the remainder of the manuscript. paper, we present comprehensive coupled cluster-level bench- marks for the polarizabilities of the ∼ 7, 000 small organic Improved SA-GPR. The formalism underlying the SA-GPR molecules contained in the QM7b database (34). We use scheme in general and the λ-SOAP (smooth overlap of atomic these reference calculations to assess the accuracy of different positions) descriptors on which our model is based have been hybrid DFT schemes and train an SA-GPR–based ML scheme introduced elsewhere (33) and are summarized in Materials and (AlphaML) that can accurately predict the polarizability tensor Methods. In this work, we include several substantial improve- with nominal computational cost. We then test the extrapolative ments that increase the accuracy and speed of the SA-GPR prediction capabilities of AlphaML on a showcase dataset com- model, and these are worth a separate discussion. For one, eval- posed of 52 larger molecules and demonstrate that this approach uation of the λ-SOAP representation is greatly accelerated by provides a viable alternative to state-of-the-art and computa- choosing the most significant few hundred spherical harmonic tionally prohibitive electronic structure methods for predicting components (of several tens of thousands) using farthest point molecular polarizabilities. sampling (FPS) (43). The calculation of the kernel in Eq. 1 can be carried out with essentially the same result as if all Results components were retained but with a much lower computa- Electronic Structure Calculations. The QM7b database (26, 34, 35) tional cost. A second improvement is the generalization of the includes N = 7,211 molecules containing up to seven “heavy” λ-SOAP kernels beyond the linear kernels used in ref. 33. It atoms (i.e., C, N, O, S, Cl) with varying levels of H satura- has been shown that, in many cases, taking an integer power tion. This dataset is based on a systematic enumeration of small of the scalar SOAP kernel improves the performance of the organic compounds (35) and contains a rich diversity of chemical associated ML model. This can be understood in terms of the groups, making it a challenging test of the accuracy associated order (two body, three body,: : :) of the interatomic correla- with DFT and quantum chemical methodologies. DFT-based tions that are described by different kernels (44, 45). In the molecular polarizabilities were obtained by (numerical) differ- tensorial case, one should be careful, as the linear nature of entiation of the molecular dipole moment, µ, with respect to the kernel is essential to ensure the correct covariant behav- an external electric field E using the hybrid B3LYP (36, 37) ior.
Load more

Recommended publications
  • Appendix 1 Table of Contents
    @ECHA EUROPEAN CHEMICALS AGENCY Appendix 1 Table of Contents: Supplementary Table C:Substances on Annex II of the Cosmetic Products Regulation ... 1 Supplementary Table D: Substances proposed to be restricted due their use restriction in CPR Annex IV column 9..,,.. 119 Supplementary Table E: Substances allowed in tattoo inks subject to the conditions specified in CPR Annex IV columns h and i ......r29 Annankatu 18. P.O. Box 400, FI-00121 Helsinki, Finland I Tel. +358 9 686180 | Fax +358 9 68618210 | echa.europa.eu ANNEX XV RESTRICTION REPORT - SUBSTANCES IN TATTOO INKS AND PERMANENT MAKE UP Su lemen Table C:Substances on Annex II of the Cosmetic Products Re ulationl Substance EC# cAs # Substance EC# cAs # R T T T A- s c R I T Name Name e b b b II s M 7 c I s I I I #4 5 6 D 9 2 1 2 3 a 3 3 3 N-(s- Chlorobenzoxa zol-2- 35783- vl)acetamide 57-4 1 (2- Acetoxyethyl)t rimethylammo (2- nium acetoxyethyl hydroxide )trimethyla 200- (Acetylcholine) 200- mmonium 124-9 5 1-84-3 and its salts t2a-9 51-84-3 2 Deanol Deanol 222- 3342- aceglumate 222- 3342- aceqlumate 085-5 61-8 (INN) 085-5 61-8 3 Spironolacto 200- Spironolactone 200- ne 133-6 52-O1-7 rINN) 133-6 52-0L-7 4 14-(4- Hydroxy-3- iodophenoxy)- 3,5- diiodophenylla cetic acid (Tiratricol 200- (INN)) and its 200- Tiratricol 086- 1 5r-24-7 salts 086- 1 5l-24-7 5 Methotrexat 200- Methotrexate 200- e 413-8 59-05-2 (INN) 413-8 59-05-2 6 Aminocaproic Aminocaproi 200- acid (INN) and 200- c acid 469-3 60-32-2 its salts 469-3 60-32-2 7 Cinchophen (rNN), its salts, derivatives and salts of 205- 132-60- these 205- 132-60- Cinchophen 067-r 5 derivatives 067-l 5 Thyropropic acid (INN) and its salts 5L-26-3 9 Trichloroacet 200- Trichloroacetic 200- ic acid 927-2 75-03-9 acid 927-2 76-03-9 l0 Aconitum napellus L.
    [Show full text]
  • Correlation.Pdf
    INTRODUCTION TO THE ELECTRONIC CORRELATION PROBLEM PAUL E.S. WORMER Institute of Theoretical Chemistry, University of Nijmegen, Toernooiveld, 6525 ED Nijmegen, The Netherlands Contents 1 Introduction 2 2 Rayleigh-SchrÄodinger perturbation theory 4 3 M¿ller-Plesset perturbation theory 7 4 Diagrammatic perturbation theory 9 5 Unlinked clusters 14 6 Convergence of MP perturbation theory 17 7 Second quantization 18 8 Coupled cluster Ansatz 21 9 Coupled cluster equations 26 9.1 Exact CC equations . 26 9.2 The CCD equations . 30 9.3 CC theory versus MP theory . 32 9.4 CCSD(T) . 33 A Hartree-Fock, Slater-Condon, Brillouin 34 B Exponential structure of the wavefunction 37 C Bibliography 40 1 1 Introduction These are notes for a six hour lecture series on the electronic correlation problem, given by the author at a Dutch national winterschool in 1999. The main purpose of this course was to give some theoretical background on the M¿ller-Plesset and coupled cluster methods. Both computational methods are available in many quantum chemical \black box" programs. The audi- ence consisted of graduate students, mostly with an undergraduate chemistry education and doing research in theoretical chemistry. A basic knowledge of quantum mechanics and quantum chemistry is pre- supposed. In particular a knowledge of Slater determinants, Slater-Condon rules and Hartree-Fock theory is a prerequisite of understanding the following notes. In Appendix A this theory is reviewed briefly. Because of time limitations hardly any proofs are given, the theory is sketchily outlined. No attempt is made to integrate out the spin, the theory is formulated in terms of spin-orbitals only.
    [Show full text]
  • Effective Hamiltonians Derived from Equation-Of-Motion Coupled-Cluster
    Effective Hamiltonians derived from equation-of-motion coupled-cluster wave-functions: Theory and application to the Hubbard and Heisenberg Hamiltonians Pavel Pokhilkoa and Anna I. Krylova a Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482 Effective Hamiltonians, which are commonly used for fitting experimental observ- ables, provide a coarse-grained representation of exact many-electron states obtained in quantum chemistry calculations; however, the mapping between the two is not triv- ial. In this contribution, we apply Bloch's formalism to equation-of-motion coupled- cluster (EOM-CC) wave functions to rigorously derive effective Hamiltonians in the Bloch's and des Cloizeaux's forms. We report the key equations and illustrate the theory by examples of systems with electronic states of covalent and ionic characters. We show that the Hubbard and Heisenberg Hamiltonians are extracted directly from the so-obtained effective Hamiltonians. By making quantitative connections between many-body states and simple models, the approach also facilitates the analysis of the correlated wave functions. Artifacts affecting the quality of electronic structure calculations such as spin contamination are also discussed. I. INTRODUCTION Coarse graining is commonly used in computational chemistry and physics. It is ex- ploited in a number of classic models serving as a foundation of modern solid-state physics: tight binding[1, 2], Drude{Sommerfeld's model[3{5], Hubbard's [6] and Heisenberg's[7{9] Hamiltonians. These models explain macroscopic properties of materials through effective interactions whose strengths are treated as model parameters. The values of these pa- rameters are determined either from more sophisticated theoretical models or by fitting to experimental observables.
    [Show full text]
  • Synthesis and Properties of Heterocyclic Acene Diimides
    ORGANIC LETTERS XXXX Synthesis and Properties of Heterocyclic Vol. XX, No. XX Acene Diimides 000–000 Cheng Li,†,‡ Chengyi Xiao,†,‡ Yan Li,*,† and Zhaohui Wang*,† Beijing National Laboratory for Molecular Sciences, Key Laboratory of Organic Solids, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China, and University of Chinese Academy of Sciences, Beijing 100049, China [email protected]; [email protected] Received December 27, 2012 ABSTRACT A series of heterocyclic acene diimides were synthesized effectively based on the condensation of o-phenylenediamine, 1,2-benzenedithiol, and 2-aminothiophenol with 2,3,6,7-tetrabromo-1,4,5,8-naphthalene tetracarboxylic diimide. The diimides exhibit interesting optical and electrical properties with one of them showing a hole mobility up to 0.02 cm2 VÀ1 sÀ1. In recent years, naphthalene tetracarboxylic diimides has been applied as high-performance solution-deposited (NDIs, 1, Figure 1) and their core-expanded derivatives ambipolar organic transistors.2 In contrast, the expansion have attracted a great deal of attention due to their of the π-system along the lateral position of NDIs has been interesting electro-optical properties and potential appli- demonstrated only recently due to the synthetic difficulties. cations as organic semiconductors in organic electronics.1 Recently, the synthesis of tetracene tetracarboxylic The π-skeleton of NDIs could be expanded along two diimides based on two methods of direct double ring directions: the peri position (1, 4, 5, 8) and the lateral extension of electron-deficient NDIs involving metal- position (2, 3, 6, 7). The expansion of the π-system along lacyclopentadienes and bismuth-triflate-mediated double- the peri position has been studied for many decades cyclization reaction of acid chlorides and isocyanates has because it induces impressive bathochromic shifts and been reported, which displays dramatic bathochromic shifts, smaller energy band gaps, and is a promising can- 3 † Institute of Chemistry, Chinese Academy of Sciences.
    [Show full text]
  • Chemistry of Acenes, [60]Fullerenes, Cyclacenes and Carbon Nanotubes
    University of New Hampshire University of New Hampshire Scholars' Repository Doctoral Dissertations Student Scholarship Spring 2011 Chemistry of acenes, [60]fullerenes, cyclacenes and carbon nanotubes Chandrani Pramanik University of New Hampshire, Durham Follow this and additional works at: https://scholars.unh.edu/dissertation Recommended Citation Pramanik, Chandrani, "Chemistry of acenes, [60]fullerenes, cyclacenes and carbon nanotubes" (2011). Doctoral Dissertations. 574. https://scholars.unh.edu/dissertation/574 This Dissertation is brought to you for free and open access by the Student Scholarship at University of New Hampshire Scholars' Repository. It has been accepted for inclusion in Doctoral Dissertations by an authorized administrator of University of New Hampshire Scholars' Repository. For more information, please contact [email protected]. CHEMISTRY OF ACENES, [60]FULLERENES, CYCLACENES AND CARBON NANOTUBES BY CHANDRANI PRAMANIK B.Sc., Jadavpur University, Kolkata, India, 2002 M.Sc, Indian Institute of Technology Kanpur, India, 2004 DISSERTATION Submitted to the University of New Hampshire in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Materials Science May 2011 UMI Number: 3467368 All rights reserved INFORMATION TO ALL USERS The quality of this reproduction is dependent upon the quality of the copy submitted. In the unlikely event that the author did not send a complete manuscript and there are missing pages, these will be noted. Also, if material had to be removed, a note will indicate the deletion. UMI Dissertation Publishing UMI 3467368 Copyright 2011 by ProQuest LLC. All rights reserved. This edition of the work is protected against unauthorized copying under Title 17, United States Code. ProQuest LLC 789 East Eisenhower Parkway P.O.
    [Show full text]
  • Community Soil Study
    Tonawanda Community Fund An affiliate of “The Wellness Institute of Greater Buffalo” Community Soil Study Fact Sheet In November 2012, Tonawanda Community Fund (TCF) led a preliminary study in a Tonawanda neighborhood to find out if air toxins were contaminating surface soil; and, if pollutants were present, do levels pose a health risk. TCF tested the soil for toxic chemicals on Kaufman Ave/James Ave/Sawyers Ave. streets in the Town of Tonawanda, NY. This neighborhood is located off of River Rd in the heart of Tonawanda’s industrial corridor. The study was initiated due to recent community complaints of a black gooey substance apparently coming from a close by air source and depositing on soil, vegetation and vehicles. The study gained more significance when Jackie James Creedon, TCF member and former director of the Clean Air Coalition of WNY, learned that a similar community’s surface soil in Birmingham, AL, where two foundry coke plants reside, was contaminated with a dangerous chemical called Benzo[a]pyrene (BAP). How the Soil was Sampled and Tested Ken‐Ton residents and TCF members, Charles Matteliano from TEQ Solutions, and Andrew Baumgartner, a student at University of Buffalo mentoring under Chemistry professor Dr .Joe Gardella, collected soil samples from a playground and 5 homes in the neighborhood. They also collected a background sample at Beaver Island State Park. Each soil sample, a composite of 6 sub samples, was tested at Test America (Amherst, NY) for 16 polycyclic aromatic hydrocarbon (PAH) compounds and 8 heavy metals. What Was Found Harmful Chemicals: Many chemicals were elevated above background sample.
    [Show full text]
  • Journal of Fluorine Chemistry 221 (2019) 1–7
    Journal of Fluorine Chemistry 221 (2019) 1–7 Contents lists available at ScienceDirect Journal of Fluorine Chemistry journal homepage: www.elsevier.com/locate/fluor Steric and electronic effects of3 CF conformations in acene(CF3)n derivatives T Nicholas J. DeWeerda, Eric V. Bukovskya, Karlee P. Castroa, Igor V. Kuvychkoa, ⁎ ⁎ ⁎ Alexey A. Popovb, , Steven H. Straussa, , Olga V. Boltalinaa, a Department of Chemistry, Colorado State University, Fort Collins, CO 80523, United States b Liebniz Institute for Solid State and Materials Research (IFW), Dresden 01069, Germany ARTICLE INFO ABSTRACT Keywords: The X-ray and DFT-optimized structures of the electron acceptor 2,3,6,7,9,10-anthra-cene(CF3)6 (ANTH-6-1) and Trifluoromethyl the 1/2 donor/acceptor co-crystal pyrene/(ANTH-6-1)2 are reported. These structures, along with extensive DFT PAH calculations on various conformations of ANTH-6-1 and 9,10-ANTH(CF3)2, suggest that the degree of bending of Acene the aromatic core in ANTH-6-1 and the DFT-predicted energies of ANTH-6-1 and 9,10-ANTH(CF3)2 are strongly Conformations 2 correlated with the relative eclipsed vs. staggered conformations of the CF3 groups attached to the central C(sp ) Bent pi systems carbon atoms (C9 and C10). Other literature X-ray and DFT-optimized structures of anthracene and pentacene X-ray structures 2 DFT calculations derivatives with n-C8F17 and/or CF3 substituents on the central C(sp ) atoms are analyzed and show (i) that the eclipsed vs. staggered correlation with the degree of bending of the aromatic cores may be a general phenom- enon and (ii) that molecules of this type are probably more stable when the acene core is bent than when it is planar.
    [Show full text]
  • Fullerene-Acene Chemistry
    University of New Hampshire University of New Hampshire Scholars' Repository Doctoral Dissertations Student Scholarship Spring 2007 Fullerene-acene chemistry: Part I Studies on the regioselective reduction of acenes and acene quinones; Part II Progress toward the synthesis of large acenes and their Diels-Alder chemistry with [60]fullerene Andreas John Athans University of New Hampshire, Durham Follow this and additional works at: https://scholars.unh.edu/dissertation Recommended Citation Athans, Andreas John, "Fullerene-acene chemistry: Part I Studies on the regioselective reduction of acenes and acene quinones; Part II Progress toward the synthesis of large acenes and their Diels-Alder chemistry with [60]fullerene" (2007). Doctoral Dissertations. 363. https://scholars.unh.edu/dissertation/363 This Dissertation is brought to you for free and open access by the Student Scholarship at University of New Hampshire Scholars' Repository. It has been accepted for inclusion in Doctoral Dissertations by an authorized administrator of University of New Hampshire Scholars' Repository. For more information, please contact [email protected]. FULLERENE-ACENE CHEMISTRY: PART I: STUDIES ON THE REGIOSELECTIVE REDUCTION OF ACENES AND ACENE QUINONES; PART II: PROGRESS TOWARD THE SYNTHESIS OF LARGE ACENES AND THEIR DIELS- ALDER CHEMISTRY WITH [60]FULLERENE VOLUME 1 CHAPTERS 1-5 BY ANDREAS JOHN ATHANS B.S. University of New Hampshire, 2001 DISSERTATION Submitted to the University of New Hampshire in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy m Chemistry May, 2007 Reproduced with permission of the copyright owner. Further reproduction prohibited without permission. UMI Number: 3 2 6 0 5 8 6 INFORMATION TO USERS The quality of this reproduction is dependent upon the quality of the copy submitted.
    [Show full text]
  • Møller–Plesset Perturbation Theory: from Small Molecule Methods to Methods for Thousands of Atoms Dieter Cremer∗
    Advanced Review Møller–Plesset perturbation theory: from small molecule methods to methods for thousands of atoms Dieter Cremer∗ The development of Møller–Plesset perturbation theory (MPPT) has seen four dif- ferent periods in almost 80 years. In the first 40 years (period 1), MPPT was largely ignored because the focus of quantum chemists was on variational methods. Af- ter the development of many-body perturbation theory by theoretical physicists in the 1950s and 1960s, a second 20-year long period started, during which MPn methods up to order n = 6 were developed and computer-programed. In the late 1980s and in the 1990s (period 3), shortcomings of MPPT became obvious, espe- cially the sometimes erratic or even divergent behavior of the MPn series. The physical usefulness of MPPT was questioned and it was suggested to abandon the theory. Since the 1990s (period 4), the focus of method development work has been almost exclusively on MP2. A wealth of techniques and approaches has been put forward to convert MP2 from a O(M5) computational problem into a low- order or linear-scaling task that can handle molecules with thousands of atoms. In addition, the accuracy of MP2 has been systematically improved by introducing spin scaling, dispersion corrections, orbital optimization, or explicit correlation. The coming years will see a continuously strong development in MPPT that will have an essential impact on other quantum chemical methods. C 2011 John Wiley & Sons, Ltd. WIREs Comput Mol Sci 2011 1 509–530 DOI: 10.1002/wcms.58 INTRODUCTION Hylleraas2 on the description of two-electron prob- erturbation theory has been used since the early lems with the help of perturbation theory.
    [Show full text]
  • Temperature-Induced Oligomerization of Polycyclic Aromatic Hydrocarbons
    www.nature.com/scientificreports OPEN Temperature-induced oligomerization of polycyclic aromatic hydrocarbons at ambient Received: 7 June 2017 Accepted: 10 July 2017 and high pressures Published: xx xx xxxx Artem D. Chanyshev 1,2, Konstantin D. Litasov1,2, Yoshihiro Furukawa3, Konstantin A. Kokh1,2 & Anton F. Shatskiy1,2 Temperature-induced oligomerization of polycyclic aromatic hydrocarbons (PAHs) was found at 500–773 K and ambient and high (3.5 GPa) pressures. The most intensive oligomerization at 1 bar and 3.5 GPa occurs at 740–823 K. PAH carbonization at high pressure is the fnal stage of oligomerization and occurs as a result of sequential oligomerization and polymerization of the starting material, caused by overlapping of π-orbitals, a decrease of intermolecular distances, and fnally the dehydrogenation and polycondensation of benzene rings. Being important for building blocks of life, PAHs and their oligomers can be formed in the interior of the terrestrial planets with radii less than 2270 km. High-pressure transformations of polycyclic aromatic hydrocarbons (PAHs) and benzene become extremely important due to wide applications for example in graphene- and graphene-based nanotechnology1–3, synthesis of organic superconductors4, 5, petroleum geoscience, origin of organic molecules in Universe and origin of life. In particular, PAHs were found in many space objects: meteorites6–8, cometary comae9, interstellar clouds and planetary nebulas10–12. Although the prevalent hypothesis for the formation of these PAHs is irradiation-driven polymerization of smaller hydrocarbons13, alternative explanation could be shock fragmentation of carbonaceous solid material11. PAH-bearing carbonaceous material could contribute to the delivery of extraterrestrial organic materials to the prebiotic Earth during the period of heavy bombardment of the inner Solar System from 4.5 to 3.8 Ga ago14–16.
    [Show full text]
  • Coupled Cluster Methods in Quantum Chemistry
    Title: Coupled Cluster methods in quantum chemistry Name: Thorsten Rohwedder, Reinhold Schneider Al./Addr.: Technische UniversitatBerlin Coupled Cluster methods in quantum chemistry Short description Coupled Cluster methods applied in quantum chemistry reformulate the electronic Schrodinger equation as a nonlinear equation, enabling the computation of size- consistent high-precision approximations of the ground state solution for weakly cor- related systems. Introduction The Coupled Cluster method (CC method) is one of the most successful and frequently used approaches for the computation of atomic and molecular electronic structure, i.e. for the solution of the stationary electronic Schrodinger equation, whenever high ac- curacy is required. In contrast to Hartree-Fock type methods (HF) or methods from density functional theory (see the respective entries in this work), high accuracy meth- ods have to account in particular for the quantum-mechanical phenomenon of electronic correlation. If a preliminarily calculated reference solution { usually provided by a HF calculation { already is a good approximation to the sought ground state wave func- tion, the problem is said to be weakly correlated. CC as a post-Hartree-Fock method (also see the entry by M. Lewin) then enables an ecient, accurate and size-extensive description of solutions of the electronic Schrodingerequation. In this context, the 2 size-extensivity of the CC approach is a key aspect, re ecting the correct scaling of correlation energy with respect to the number of electrons. CC methods were initially developed for the treatment of many body quantum systems in nuclear physics in the 1950s and were used for quantum chemical calculations since the 1966 initial work by Paldus and Czek,see (Czek,1991) for a historical overview.
    [Show full text]
  • The Pennsylvania State University the Eberly College of Science
    The Pennsylvania State University The Eberly College of Science Department of Chemistry SYNTHESIS OF CARBON MATERIALS VIA THE COLD COMPRESSION OF AROMATIC MOLECULES AND CARBON NANOSTRUCTURES A Dissertation in Chemistry by Thomas C. Fitzgibbons 2014 Thomas C. Fitzgibbons Submitted in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy August 2014 ii The dissertation of Thomas C. Fitzgibbons was reviewed and approved* by the following: John V. Badding Professor of Chemistry Dissertation Advisor Chair of Committee Ayusman Sen Distinguished Professor of Chemistry A. Welford Castleman, Jr Evan Pugh Professor of Chemistry and Physics William B. White Professor Emeritus of Geosciences Barbara J. Garrison Shapiro Professor of Chemistry Head of the Department of Chemistry *Signatures are on file in the Graduate School iii ABSTRACT Carbon’s ability for catenation makes it a remarkable element and allows for many interesting and surprising properties and structures. Carbon can exist in one of its two thermodynamically stable bulk crystals, graphite or diamond, one of its several nanostructures: fullerene, nanotube, or graphene, or as an amorphous material with a mixed bonding pattern. Carbon also has an ability to bond heteroatoms such as hydrogen which can increase its properties and structures even further. Pressure has been shown to be able to drastically change the bonding in and structure of carbon based materials. In this dissertation I will present how pressure can be used to synthesize new amorphous hydrogenated carbons and how a battery of analytical techniques can be used to elicit the microstructure of the carbon networks. This microstructure can then be related back to the reaction conditions and more importantly the starting small molecule.
    [Show full text]