Overview to Molecular Modeling Computational Chemistry

Total Page:16

File Type:pdf, Size:1020Kb

Overview to Molecular Modeling Computational Chemistry Overview to Molecular Modeling Computational Chemistry E of a molecular structure Geometry optimization Related properties vibrational frequencies nmr e ) density Energy method / Energy basis set // Geometry method / Geometry basis set Computational Chemistry Molecular Mechanics Atoms obey laws of classical physics No e) structure MM2, MM3, MM+, others Useful Large (bio) molecules Small molecules N O energy value Computational Chemistry Molecular Mechanics 3 E = E i Large number of parameters C 2 H 6 C-C, 6 @ C-H 6 @ C - C - H 9 @ H - C - C - H C 6 H 6 6 @ C - H, 6 @ C -/= C (not C - C or C = C) 6 @ C - C - H, 24 torsion Parameters determined empirically Computational Chemistry Molecular Quantum Mechanics Mechanics Electronic structure based on , ø = E ø , is known exactly ø is unknown except for simple systems (H-like atoms, SHO, RR, particles in boxes, etc.) Computational Chemistry Problems Molecular Quantum Mechanics Mechanics Overlap Integral Exchange Integral Exchange Functional (HF theory) Correlation Functional Computational Chemistry Molecular Quantum Mechanics Mechanics Semiempirical Methods Ignore part of , Hückel molecular orbital theory MOPAC theory ZINDO theory Computational Chemistry Molecular Quantum Hückel Theory Mechanics Mechanics Semiempirical Methods HMO Hückel molecular orbital theory Applied to conjugated hydrocarbons Assumes ALL overlap integrals are zero EHT Extended Hückel theory Applied to any molecule type Useful for “quick and dirty” calculations and starting point for more advanced calculations Computational Chemistry Molecular Quantum MOPAC Mechanics Mechanics Molecular Orbital Package Semiempirical Methods CNDO Complete Neglect of Differential Overlap INDO Intermediate Neglect of ... NDDO Neglect of Diatomic ... MINDO Modified INDO MINDO/3 MNDO Modified Neglect of ... AM1 Austin Model 1 PM3 Parameterized Model Series 3 AM1/d and MNDO-d (MOPAC 2000, d e- ’s) Useful for ground state energy and geometry Computational Chemistry Molecular Quantum ZINDO Mechanics Mechanics Zerner’s INDO Semiempirical Methods ZINDO/1, ZINDO/3, ZINDO-d, etc Useful for Transition states Energies Spectroscopy Transition elements Not useful for optimizations Computational Chemistry Molecular Quantum Mechanics Mechanics Semiempirical ab initio Methods Methods Use complete , Estimate ø $ Variation Principle (E trial E experimental) Computational Chemistry Molecular Quantum Mechanics Mechanics Level of Theory Semiempirical ab initio Methods Methods HF-SCF Hartree-Fock Self-Consistent Field B3LYP Density Function Theory (DFT) Becke Exchange with Lee-Yang-Parr Correlation MP2/MP4 Second/Fourth Order Møller-Plesset perturbation theory QCISD(T) Quadratic configuration interaction Computational Chemistry Molecular Quantum Mechanics Mechanics Trial Wave Functions Semiempirical ab initio (Basis Sets) Methods Methods Computational Chemistry Molecular Quantum Mechanics Mechanics Electron Spin Semiempirical ab initio Methods Methods Open Shell (unrestricted) Odd number of electrons Excited states 2 or more unpaired electrons Bond dissociation processes Closed Shell (restricted) Even number of electrons--all paired Comparison of ab initio Methods (p 94) Comparison of Models (F/F p 96) Comparison of Commercial Software Basis Sets Basis Set Criteria Capable of describing actual wave function well enough to give chemically useful results Can be used to evaluate I ’s accurately and “cheaply” Basis Functions Hydrogenlike Orbitals Hydrogenlike Orbitals ( n - l - 1) nodes Basis Functions Slater-type Orbitals (STO’s) Basis Functions Gaussian-type Orbitals (GTO’s) GTO’s Advantages Complete Favorable math properties Disadvantages Not mutually orthogonal Poor representation of electron probability near and far away from nucleus (overcome using large number of GTO’s Use of STO’s One or more STO on each nucleus Accuracy of calculation increases as Orbital exponents chosen well Number of STO’s used increases Number of STO’s used Minimal Basis Set Use STO for occupied AO’s Examples H 1s C 1s 2s 2px 2 p y 2 p z Number of STO’s used Split (Double Zeta æ ) Basis Set Linear combination of two similar orbitals with different orbital exponents (different sizes) ö 2 p = a ö 2p,inner + b ö 2p,outer If a > b charge cloud contracted around nucleus If b > a diffuse cloud Number of STO’s used Split (Double Zeta æ ) Basis Set Examples: H 1s, 1sN N N N N C 1s, 2s, 2s , 2px , 2py , 2pz , 2px , 2py , 2pz Triple Zeta basis sets are also used Number of STO’s used Diffuse Basis Set Extra s and p wave functions included that are significantly larger than usual ones Useful for Distant electrons Molecules with lone pairs Anions Species with significant negative charge Excited states Species with low ionization potentials Describing acidities Number of STO’s used Polarized Basis Set Linear combination of different types of orbitals Examples H 1s and 2p C 1s, 2s, 2p and 3d Shifts charge in/out of bonding regions Number of STO’s used Other attempts Place STO’s in center of bonds instead of on only nuclei Problems with increasing number of STO’s used Number of I ’s increases as N4 where N is the number of basis functions As minimization occurs, orbital exponents change thus defining a new basis set to rebegin the calculation Use of STO’s/GTO’s Wrong shape of GTO’s accounted for by Choosing several á ’s to get set of “primitive” gaussians for compact and diffuse Linear combination of primitives (usually 1-7) to get STO Optimize “Freeze” as “contracted” gaussian function Use minimal, split/double zeta, polarization, diffuse sets Use of STO’s / GTO’s Jargon: minimal basis set STO-NG where N is the number of primitive gaussians STO-3G 3 primitve gaussians per basis set not the simplest minimal basis set popular Use of STO’s / GTO’s Jargon: split basis set K-LMG where K is the number of sp type inner shell primitive gaussians L is the number of inner valence s and p primitive gaussians M is the number of outer valence s and p primitive gaussians Use of STO’s / GTO’s Jargon: split basis set 3-21G 3 primitives for inner shell 2 sizes of basis functions for each valence orbital 6-311G 6 primitives for inner shell 3 sizes of basis functions for each valence orbital Use of STO’s / GTO’s Jargon: polarization * d-type orbital added to atoms with Z > 2 ** d-type orbital added to atoms with Z > 2 and p-type orbital added to H and He d’s added: STO-NG are 5 regular 3d’s L-KMG are 6 3d’s dx x , d y y , d z z , d x y , d y z , d x z (formed by linear combination of 5 regular 3d’s and 3s) Use of STO’s / GTO’s Jargon: polarization 6-31G* or 6-31G(d) 6-31G with d added for Z > 2 (FF choice) 6-31G** or 6-31G(d,p) 6-31G with d added for Z > 2 and p added to H 6-31G(2d) 6-31G with 2d functions added for Z > 2 Use of STO’s / GTO’s Jargon: diffuse + diffuse function included for Z > 2 ++ diffuse function included for Z > 2 and for H 6-31+G(d) 6-31G(d) with diffuse function added for Z > 2 6-31++G(d) 6-31+G(d) with diffuse function added for H Some Recommended Standard Basis Sets (F/F p 102) ~DZVP ~TZVP Common Basis Sets.
Recommended publications
  • Semiempirical Implementations of Molecular
    Semiempirical Implementations of Molecular Orbital Theory How one can make Hartree-Fock theory less computationally intensive without much sacrificing its accuracy? The most demanding step – calculation of two-electron (four-index) integrals (J and K integrals) appearing in the Fock matrix elements (N4 where N is the number of basis functions). One way to save time – to estimate their value accurately in an a priori fashion and thus to avoid numerical integration. Coulomb integrals measure the repulsion between electrons in regions of space defined by the basis functions. When the basis functions in the integral for one electron are very far from the basis functions for the other, the value of that integral will approach zero. In a large molecule, one might be able to avoid the calculation of a very large number of integrals simply by assuming them to be zero. HF theory is intrinsically inaccurate as it does not include correlation energy. Therefore, modifications of the theory introduced in order to simplify its formalism may actually improve it, provided the new approximations somehow introduce an accounting for correlation energy. Most typically, such approximations involve the adoption of a parametric form for some aspect of the calculation where the parameters involved are chosen so as best reproduce experimental data – ‘semiempirical’. Another motivation for introducing semiempirical approximation into HF theory was to facilitate the computation of derivatives (gradients, Hessians) so that geometries could be more efficiently optimized. Extended Hückel Theory Before considering semiempirical methods we revisit Hückel theory: H11 − ES11 H12 − ES12 L H1N − ES1N H21 − ES21 H22 − ES22 L H2N − ES2N = 0 M M O M H − ES H − ES H − ES N1 N1 N2 N 2 L NN NN The dimension of the secular determinant depends on the choice of the basis set.
    [Show full text]
  • Semiempirical Quantum-Chemical Methods Max-Planck-Institut Für
    Max-Planck-Institut für Kohlenforschung This is the peer reviewed version of the following article: WIREs Comput. Mol. Sci. 4, 145-157 (2014), which has been published in final form at https://doi.org/10.1002/wcms.1161. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Use of Self- Archived Versions. Semiempirical quantum-chemical methods Walter Thiel Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1, 45470 Mülheim, Germany [email protected] Abstract The semiempirical methods of quantum chemistry are reviewed, with emphasis on established NDDO-based methods (MNDO, AM1, PM3) and on the more recent orthogonalization-corrected methods (OM1, OM2, OM3). After a brief historical overview, the methodology is presented in non- technical terms, covering the underlying concepts, parameterization strategies, and computational aspects, as well as linear scaling and hybrid approaches. The application section addresses selected recent benchmarks and surveys ground-state and excited-state studies, including recent OM2- based excited-state dynamics investigations. Introduction Quantum mechanics provides the conceptual framework for understanding chemistry and the theoretical foundation for computational methods that model the electronic structure of chemical compounds. There are three types of such approaches: Quantum-chemical ab initio methods provide a convergent path to the exact solution of the Schrödinger equation and can therefore give “the right answer for the right reason”, but they are costly and thus restricted to relatively small molecules (at least in the case of the highly accurate correlated approaches). Density functional theory (DFT) has become the workhorse of computational chemistry because of its favourable price/performance ratio, allowing for fairly accurate calculations on medium-size molecules, but there is no systematic path of improvement in spite of the first-principles character of DFT.
    [Show full text]
  • Universit`A Di Parma
    UNIVERSITA` DI PARMA DEPARTMENT OF CHEMISTRY, LIFE SCIENCES AND ENVIRONMENTAL SUSTAINABILITY Doctoral Programme in Chemical Sciences XXIX Cycle DYES AND NANOPARTICLES FOR 2PA APPLICATIONS: MODELS AND COMPUTATIONS PhD Student: Somananda Sanyal Supervisors: Prof. Anna Painelli Prof. Swapan K. Pati Coordinator: Prof. Roberto Corradini 2014-2017 Dedicated to My Family List of Abbreviations BFC BF2 complex of Curcumin CT Charge Transfer CNT Carbon Nanotubes D/A Electron Donor / Acceptor 0 DANS p,p -dimethylamino-nitrostilbene DFT Density Functional Theory DMRG Density Matrix Renormalization Group Theory eV electron-volt unit ESM Essential State Model ESP Electrostatic Potential FMO Frontier Molecular Orbital GM G¨oppert-Mayer units (1GM ≡ 10−50 cm 4 s photon−1) HOMO Highest Occupied Molecular Orbital HRS Hyper Rayleigh Scattering LUMO Lowest Unoccupied Molecular Orbital MO Molecular Orbital OPA One Photon Absorption PCM Polarizable Continuum Model PPP Pariser-Parr-Pople TPA Two Photon Absorption ZINDO Zerner's Intermediate Neglect of Differential Overlap i Acknowledgements It's time to say Thank You to many people who have either actively or silently been supporting me in successfully completing the PhD thesis! These three years of my PhD tenure in Parma, Italy has been an eye opener for me, as I slowly learnt to adapt as a girl in a new city. I have been extremely fortunate to meet some people whose influence has been profound and worth cherishing, and I want to grab this opportunity to Thank them All! My parents have been my biggest support system throughout my life and from them I have learnt the first lessons of morality and humanity.
    [Show full text]
  • Biochem Press
    Internet Electronic Journal of Molecular Design 2003, 2, 86–95 ISSN 1538–6414 BioChem Press http://www.biochempress.com Internet Electronic Journal of Molecular Design February 2003, Volume 2, Number 2, Pages 86–95 Editor: Ovidiu Ivanciuc Special issue dedicated to Professor Haruo Hosoya on the occasion of the 65th birthday Part 6 Guest Editor: Jun–ichi Aihara Partial Structures of C60 Responsible for Its Lowest Electronic Absorption Band: Corannulene or Triphenylene? Tai–ichi Shibuya,1 Susumu Narita,1 Yasushi Nomura,1 and Tetsuo Morikawa2 1 Department of Chemistry, Faculty of Textile Science and Technology, Shinshu University, Ueda, Nagano–ken, 386–8567, Japan 2 Department of Chemistry, Joetsu University of Education, Joetsu, Niigata–ken, 943–8512, Japan Received: October 21, 2002; Revised: November 22, 2002; Accepted: December 15, 2002; Published: February 28, 2003 Citation of the article: T. Shibuya, S. Narita, Y. Nomura, and T. Morikawa, Partial Structures of C60 Responsible for Its Lowest Electronic Absorption Band: Corannulene or Triphenylene?, Internet Electron. J. Mol. Des. 2003, 2, 86–95, http://www.biochempress.com. Copyright © 2003 BioChem Press T. Shibuya, S. Narita, Y. Nomura, and T. Morikawa Internet Electronic Journal of Molecular Design 2003, 2, 86–95 Internet Electronic Journal BioChem Press of Molecular Design http://www.biochempress.com Partial Structures of C60 Responsible for Its Lowest Electronic Absorption Band: Corannulene or Triphenylene?# Tai–ichi Shibuya,1,* Susumu Narita,1 Yasushi Nomura,1 and Tetsuo Morikawa2 1 Department of Chemistry, Faculty of Textile Science and Technology, Shinshu University, Ueda, Nagano–ken, 386–8567, Japan 2 Department of Chemistry, Joetsu University of Education, Joetsu, Niigata–ken, 943–8512, Japan Received: October 21, 2002; Revised: November 22, 2002; Accepted: December 15, 2002; Published: February 28, 2003 Internet Electron.
    [Show full text]
  • Α-Acetylpyridine Ketone'un Ve Argerol'lü Bileşiğinin
    α-ACETYLPYRIDINE KETONE’UN VE ARGEROL'LÜ BİLEŞİĞİNİN TİTREŞİMLERİNİN DENEYSEL VE TEORİK OLARAK İNCELENMESİ ABDULLAH ATILGAN YÜKSEK LİSANS TEZİ FİZİK GAZİ ÜNİVERSİTESİ FEN BİLİMLERİ ENSTİTÜSÜ ŞUBAT 2013 ANKARA Abdullah ATILGAN tarafından hazırlanan “ α - ACETYLPYRIDINE KETONE’UN VE ARGEROL'LÜ BİLEŞİĞİNİN TİTREŞİMLERİNİN DENEYSEL VE TEORİK OLARAK İNCELENMESİ” adlı bu tezin Yüksek Lisans tezi olarak uygun olduğunu onaylarım. Prof. Dr. Şenay YURDAKUL ………………………………….. Tez Danışmanı, Fizik Anabilim Dalı Bu çalışma, jürimiz tarafından oy birliği ile Fizik Anabilim Dalında Yüksek Lisans tezi olarak kabul edilmiştir. Prof. Dr. Atike Semra BİLGİÇ ………………………………….. Kimya Anabilim Dalı, Ankara Üniversitesi Prof. Dr. Şenay YURDAKUL ………………………………….. Fizik Anabilim Dalı, Gazi Üniversitesi Doç. Dr. Semran SAĞLAM ………………………………….. Fizik Anabilim Dalı, Gazi Üniversitesi Tez SavunmaTarihi: 04/02/2013 Bu tez ile G.Ü. Fen Bilimleri Enstitüsü Yönetim Kurulu Yüksek Lisans derecesini onamıştır. Prof.Dr. Şeref SAĞIROĞLU ………………………………….. Fen Bilimleri Enstitüsü Müdürü TEZ BİLDİRİMİ Tez içindeki bütün bilgilerin etik davranış ve akademik kurallar çerçevesinde elde edilerek sunulduğunu, ayrıca tez yazım kurallarına uygun olarak hazırlanan bu çalışmada bana ait olmayan her türlü kaynağa eksiksiz atıf yapıldığını bildiririm. ABDULLAH ATILGAN iv α-ACETYLPYRİDİNE KETONE’ÜN VE ARGEROL'LÜ BİLEŞİĞİNİN TİTREŞİMLERİNİN DENEYSEL VE TEORİK OLARAK İNCELENMESİ (Yüksek Lisans Tezi) Abdullah ATILGAN GAZİ ÜNİVERSİTESİ FEN BİLİMLERİ ENSTİTÜSÜ Şubat 2013 ÖZET Bu çalışmada C7H7NO (α-acetylpyridin ketone) genel formülü ile verilen ligandın ve argerollü bileşiğinin IR spektrumu kaydedilip bu bileşiklerin titreşim frekans ve kipleri saptandı. Ligandın ve bileşiğin titreşim frekanslarına işaretleme yapıldı. Gaussian 09 bilgisayar programı yardımıyla teorik frekans değerleri, geometrik parametreler ve NBO yük dağılımı elde edildi ve bu değerler deneysel sonuçlar ile karşılaştırıldı. Bilim Kodu : 202.1.008 Anahtar Kelimeler : α-acetylpyridine ketone, argerol bileşiği, IR, Ra, NBO Sayfa Adedi : 65 Tez Yöneticisi : Prof.
    [Show full text]
  • Computational Chemistry: a Practical Guide for Applying Techniques to Real-World Problems
    Computational Chemistry: A Practical Guide for Applying Techniques to Real-World Problems. David C. Young Copyright ( 2001 John Wiley & Sons, Inc. ISBNs: 0-471-33368-9 (Hardback); 0-471-22065-5 (Electronic) COMPUTATIONAL CHEMISTRY COMPUTATIONAL CHEMISTRY A Practical Guide for Applying Techniques to Real-World Problems David C. Young Cytoclonal Pharmaceutics Inc. A JOHN WILEY & SONS, INC., PUBLICATION New York . Chichester . Weinheim . Brisbane . Singapore . Toronto Designations used by companies to distinguish their products are often claimed as trademarks. In all instances where John Wiley & Sons, Inc., is aware of a claim, the product names appear in initial capital or all capital letters. Readers, however, should contact the appropriate companies for more complete information regarding trademarks and registration. Copyright ( 2001 by John Wiley & Sons, Inc. All rights reserved. No part of this publication may be reproduced, stored in a retrieval system or transmitted in any form or by any means, electronic or mechanical, including uploading, downloading, printing, decompiling, recording or otherwise, except as permitted under Sections 107 or 108 of the 1976 United States Copyright Act, without the prior written permission of the Publisher. Requests to the Publisher for permission should be addressed to the Permissions Department, John Wiley & Sons, Inc., 605 Third Avenue, New York, NY 10158-0012, (212) 850-6011, fax (212) 850-6008, E-Mail: PERMREQ @ WILEY.COM. This publication is designed to provide accurate and authoritative information in regard to the subject matter covered. It is sold with the understanding that the publisher is not engaged in rendering professional services. If professional advice or other expert assistance is required, the services of a competent professional person should be sought.
    [Show full text]
  • UFRRJ Catarina De Nigris Del Cistia
    UFRRJ INSTITUTO DE CIÊNCIAS EXATAS PROGRAMA DE PÓS-GRADUAÇÃO EM QUÍMICA TESE MODELAGEM MOLECULAR APLICADA AO ESTUDO DE REAÇÕES DE INIBIÇÃO ENZIMÁTICA COM APLICAÇÃO POTENCIAL NO CONTROLE DE Leishmania amazonensis Catarina De Nigris Del Cistia 2010 Livros Grátis http://www.livrosgratis.com.br Milhares de livros grátis para download. UNIVERSIDADE FEDERAL RURAL DO RIO DE JANEIRO INSTITUTO DE CIÊNCIAS EXATAS PROGRAMA DE PÓS-GRADUAÇÃO EM QUÍMICA MODELAGEM MOLECULAR APLICADA AO ESTUDO DE REAÇÕES DE INIBIÇÃO ENZIMÁTICA COM APLICAÇÃO POTENCIAL NO CONTROLE DE Leishmania amazonensis CATARINA DE NIGRIS DEL CISTIA Sob a orientação do Professor Carlos Maurício Rabello de Sant’Anna e Co-orientação da Professora Aurea Echevarria Tese submetida como requisito parcial para obtenção do grau de Doutor em Ciências, no Programa de Pós- Graduação em Química, Área de Concentração em Química Orgânica Seropédica, RJ Junho de 2010 UFRRJ / Biblioteca Central / Divisão de Processamentos Técnicos 547.21 D331m Del Cistia, Catarina de Nigris, 1980-. T Modelagem molecular aplicada ao estudo de reações de inibição enzimática com aplicação potencial no controle de Leishmania amazonensis / Catarina de Nigris Del Cistia – 2010. 193 f.: il. Orientador: Carlos Maurício Rabello de Sant’Anna. Tese (doutorado) – Universidade Federal Rural do Rio de Janeiro, Programa de Pós- Graduação em Química. Bibliografia: f. 144-159. 1. Síntese orgânica – Teses. 2. Leishmania - Teses. 3. Leishmaniose - Teses. 4. Inibidores enzimáticos – Teses. I. Sant’Anna, Carlos Maurício Rabello de, 1965-. II. Universidade Federal Rural do Rio de Janeiro. Programa de Pós-Graduação em Química. III. Título. “Os loucos abrem os caminhos que mais tarde os sábios percorrem.” James Joyce AGRADECIMENTOS .
    [Show full text]
  • Faster. Smarter. Just Better
    AMPAC Faster. Smarter. Just better. Feature AMPAC 9 MOPAC 2007 Graphical User Interface (GUI) Molecule Building / Viewing ▪ Visual display of properties, surfaces, MOs ▪ Animation of reaction coordinates and vibrations ▪ Gaussian03 – shares common interface ▪ Graphic images for publication ▪ Submit and manage jobs from the GUI ▪ Read pdb files and accurately complete hydrogens ▪ Semiempirical Methods AM1, MNDO, MINDO3, PM3, MNDO/d, RM1, PM6 ▪ ▪ SAM1 (transition metals Fe and Cu), MNDOC ▪ Geometry Optimization and SCF Convergence Automatic Heuristic SCF Convergence ∆ ▪ RHF/UHF ▪ ▪ TRUSTE and TRUSTG geometry optimization ∆ ▪ ▪ Eigenvector Following (EF) ▪ CHN and LTRD transition state location methods ∆ ▪ PATH, IRC for potential surfaces and reaction pathways ▪ Reaction pathway definition ▪ GRID 2D reaction pathway investigation ▪ LFORCE for rapid characterization of TSs and minima ▪ Sparse matrix method for large molecules * ▪ Additional Methods Configuration Interaction (CI) ∆ ▪ Selected State Optimized CI ∆ ▪ Analytic CI Gradients and higher spin multiplicities (20) ∆ ▪ Simulated Annealing for Multiple Minima Searches ▪ AMSOL Method for Solvated Molecules ▪ COSMO Solvation Method ▪ ▪ Tomasi Solvation Method ▪ Property Calculations Thermodynamic Properties ▪ ▪ Unpaired Electron Spin Density ▪ ▪ Population Analysis: Coulson / Mulliken / ESP ▪ ▪ Non-linear Optical Properties ▪ Polymers and solid states ▪ Limited molecular dynamics ▪ Compatibility and Support Fully Compatible with CODESSA™ QSAR Program ▪ Multiple File formats for read/write (mol, mol2, G03, pdf, CIF) ▪ Manual: new, fully updated and indexed in hypertext format ▪ Updating and addition of new features regularly ▪ Customer Support - knowledgeable and available ▪ Generous site-licensing for academics ▪ * Under active development. Common feature enhanced and improved in AMPAC. ∆ AMPAC much faster and more reliable .
    [Show full text]
  • 11. MOLECULAR MODELING Chemists Frequently
    11. MOLECULAR MODELING Chemists frequently use models to better understand chemical phenomena. In this laboratory models are used to analyze chemical bonding and molecular structure. We will start with the Lewis electron dot (Localized Electron) model1 and use it to predict the electronic structure of some simple molecules. We will then use ball-and-stick models to construct three-dimensional models of those molecules. The three-dimensional models will be refined using Gillespie’s Valence Shell Electron Pair Repulsion (VSEPR) model2. Finally, we will use Spartan3, a computational chemistry program, to calculate the molecular geometry. The results from Spartan will then be compared with the VSEPR results. Computational chemistry. In theory, all of chemistry can be understood as solutions of Schrödinger’s wave equation. In practice, solution of the Schrödinger equation is so complicated that chemically accurate calculations for molecules involving more than three electrons and three nuclei, e.g. H3, are beyond the reach of current computer hardware/algorithms. To model the properties of larger (and more interesting) molecules, chemists must use some form of approximate computational method. The simplest, and thus fastest, computational method is called molecular mechanics. In molecular mechanics the bonds of a molecule are treated as springs and the force constants of the springs are determined from experimental measurements. The overall molecular energy is the sum of the energies of stretching, bending and torsion (the molecular force field) for all of the bonds in the molecule. Although molecular mechanics ignores electronic properties, it gives quite good structural results for large molecules such as proteins. Drug companies performed much of the pioneering work on molecular mechanics.
    [Show full text]
  • Fundamentals of Quantum Mechanics for Chemistry
    Fundamentals of Quantum Mechanics for Chemistry G´eraldMONARD Equipe de Chimie et Biochimie Th´eoriques UMR 7565 CNRS - Universit´eHenri Poincar´e Facult´edes Sciences - B.P. 239 54506 Vandœuvre-les-Nancy Cedex - FRANCE http://www.monard.info/ Outline . 1. Fundamentals of Quantum Mechanics for Chemistry 3 Hartree-Fock methods 3 Density Functional Theory 3 The QM scaling problem 3 Semiempirical methods 3 Molecular Mechanics 2. Fundamentals of QM/MM methods 3 Partionning 3 QM/MM interactions 3 Cutting covalent bonds 3 ONIOM 3 Some available software 3. Selected QM/MM applications 3 Solvent effects 3 Spectroscopy 3 Biochemistry . Outline 4. Fundamentals of Linear Scaling methods 3 QM Bottlenecks 3 General ideas and solutions 3 Some available software 5. Focus on some Linear Scaling methods 3 CG-DMS 3 Mozyme 3 Divide & Conquer 6. Selected Linear Scaling applications 3 Energy Decomposition; Charge Transfer & Polarization 3 Born-Oppenheimer Molecular Dynamics 7. Parallelization of QM/MM and Linear Scaling methods Fundamentals of Quantum Mechanics for Chemistry (1) Some approximations And there was the Schr¨odinger equation. H0Ψ0 = E0Ψ0 where: H0 is an Hamiltonian operator that describes a molecular system Ψ0 is a wavefunction (solution of the Schr¨odingerequation) that describe a state of the system E0 the energy associated to Ψ0 + 1 equation + 2 unknowns (given H0) = +¥ solutions! ­ From now on: 3 ground state 3 closed shell 3 non-relativistic Fundamentals of Quantum Mechanics for Chemistry (2) Born-Oppenheimer approximation ­ nuclei are fixed point charges ­ only electrons are represented by a wavefunction Ψ HΨ = Eelec Ψ (1) H = Te + VeN + Vee 1 −Z 1 K (2) = − ∑∆i + ∑∑ + ∑∑ 2 i i K riK i i>j rij | {z } | {z } | {z } kinetic energy e−-nuclei inter.
    [Show full text]
  • Semiemprical Methods Chapter 5 Semiemprical Methods
    Semiemprical Methods Chapter 5 Semiemprical Methods Because of the difficulties in applying ab initio methods to medium and large molecules, many semiempirical methods have been developed to treat such molecules. The earliest semiempirical methods treated only the π electrons of conjugated molecules. In the π-electron approximation, the nπ π electrons are treated separately by incorporating the effects of the σ electrons and the nuclei into some sort of effective π- electron Hamiltonian Ĥπ (recall the similar valence-electron approximation where Vi is the potential energy of the ith π electron in the field produced by the nuclei and the σ electrons. The core is everything except the π electrons. The most celebrated semiempirical π-electron theory is the Hückel molecular-orbital (HMO) method, developed in the 1930s to treat planar conjugated hydrocarbons. Here the π-electron Hamiltonian is approximated by the simpler form where Ĥeff(i) incorporates the effects of the π-electron repulsions in an average way. In fact the Hückel method does not specify any explicit form for Ĥeff(i). Since the Hückel π-electron Hamiltonian is the sum of one-electron Hamiltonians, a separation of variables is possible. We have The wave function takes no account of spin or the antisymmetry requirement. To do so, we must put each electron in a spin-orbital ui=ii. The wave function π is then written as a Slater determinant of spin-orbitals. Since Ĥeff(i) is not specified, there is no point in trying to solve above eq. directly. Instead, the variational method is used. The next assumption in the HMO method is to approximate the π MOs as LCAOs.
    [Show full text]
  • Winmostar™ User Manual Release 10.7.0
    Winmostar™ User Manual Release 10.7.0 X-Ability Co., Ltd. Sep 30, 2021 Contents 1 Introduction 2 2 Installation Guide 26 3 Main Window 30 4 Basic Operation Flow 33 5 Structure Building 36 6 Main Menu and Subwindows 43 7 Remote job 183 8 Add-On 191 9 Integration with other software 200 10 Other topics 202 11 Known problems 208 12 Frequently asked questions · Troubleshooting 212 Bibliography 242 i Winmostar™ User Manual, Release 10.7.0 This manual describes the operation method of each function of Winmostar (TM). The latest version of this document is available from Official site. If you are using Winmostar (TM) for the first time, please refer to Quick Manual. If there is an uncertain point or it does not move as expected, please confirm Frequently asked questions · Troubleshooting which is updated from time to time. For specific operational procedures for each purpose, such as chemical reaction analysis and calculation of specific physical properties, see various tutorials. Contents 1 CHAPTER 1 Introduction Winmostar (TM) provides a graphical user interface that can efficiently manipulate quantum chemical calculations, first principles calculations, and molecular dynamics calculations. From the creation of the initial structure, from the calculation execution to the result analysis, you can carry out the one operation required for the simulation on Winmostar (TM). For molecular modeling function it has been confirmed to operate up to 100,000 atoms. The function of MD calculation has been confirmed in a larger system. 1.1 About quotation When announcing data created using Winmostar (TM) in academic presentations, articles, etc., please describe the Winmostar (TM) main body as follows, for example.
    [Show full text]