For Peer Review Journal: Journal of Computational Chemistry

Total Page:16

File Type:pdf, Size:1020Kb

For Peer Review Journal: Journal of Computational Chemistry Journal of Computational Chemistry A Toolkit to Assist ONIOM Calculations For Peer Review Journal: Journal of Computational Chemistry Manuscript ID: Draft Wiley - Manuscript type: Software News and Updates Date Submitted by the Author: Complete List of Authors: Tao, Peng; Wayne State University, Department of Chemistry Schlegel, H. Bernhard; Wayne State University, Department of Chemistry Key Words: QM/MM, ONIOM, charge fitting, efficiency John Wiley & Sons, Inc. Page 1 of 20 Journal of Computational Chemistry 1 2 3 4 A Toolkit to Assist ONIOM Calculations 5 6 7 8 9 10 11 Peng Tao, and H. Bernhard Schlegel * 12 13 14 Department of Chemistry, Wayne State University, 5101 Cass Ave, Detroit, Michigan 48202 15 16 17 [email protected] 18 For Peer Review 19 RECEIVED DATE 20 21 22 Correspondence to: H. Bernhard Schlegel, Department of Chemistry, Wayne State University, Detroit, 23 24 Michigan 48202. Email: [email protected] 25 26 27 This article contains supplementary material available in the online version of this article. 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 1 John Wiley & Sons, Inc. Journal of Computational Chemistry Page 2 of 20 Abstract: A general procedure for QM/MM studies on biochemical systems is outlined, and a 1 2 3 collection of PERL scripts to facilitate ONIOM type QM/MM calculations is described. This toolkit is 4 5 designed to assist in the different stages of an ONIOM QM/MM study of biomolecules, including input 6 7 file preparation and checking, job monitoring, production calculations, and results analysis. An iterative 8 9 10 procedure for refitting the partial charges of QM region atoms is described and yields a more accurate 11 12 treatment of the electrostatic interaction between QM and MM regions during QM/MM calculations. 13 14 15 The toolkit fully supports this partial charge refitting procedure. By using this toolkit for file 16 17 conversions, structure manipulation, input sanity checks, parameter lookup, charge refitting, tracking 18 For Peer Review 19 optimizations and analyzing results, QM/MM studies of large size biochemical systems can be much 20 21 22 more convenient and practical. 23 24 25 26 Keywords: QM/MM, ONIOM, charge fitting, efficiency 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 2 John Wiley & Sons, Inc. Page 3 of 20 Journal of Computational Chemistry 1 Introduction 2 3 Computational chemistry and biology are becoming routine methods to study many biological 4 5 processes. In computational studies of biochemical reactions involving macromolecules, one of the most 6 7 powerful tools available is combined quantum mechanics and molecular mechanics (QM/MM) 8 1-5 9 methods. In QM/MM methods, the large system of interest is divided into two (or more) parts. The 10 larger part, which is not directly involved with the chemical reaction, is treated at a low cost level of 11 12 theory, such as molecular mechanics (MM). The smaller part where the chemical reaction occurs is 13 14 treated with a more accurate level of theory, such as quantum mechanics (QM), with a higher 15 16 computational cost. Using this scheme, the entire system, e.g. an enzyme with substrate and solvent 17 molecules, can be treated computationally. The reaction of interest in the active site of the enzyme can 18 For Peer Review 19 be studied at an appropriate level of theory with an affordable computational cost. 20 21 Although different cases may vary, the general procedure for a typical QM/MM study of a 22 23 biochemical system can be described as follows. First, reliable experimental structures, usually crystal 24 25 structures, are identified and preprocessed. The preprocessing includes adding hydrogen atoms, 26 determining the most probable protonation states of titratable amino acids, add counter ions and solvent 27 28 molecules, geometry optimization with MM and equilibration by molecular dynamics (MD) using MM. 29 30 Then, representative geometries are selected from the MD simulations and are subject to optimization 31 32 again before the QM/MM study is initiated. Third, QM/MM geometry optimizations are set up from 33 selected MD geometries and conducted for the initial structures, usually reactant complexes in the 34 35 reaction sequence under study. Fourth, transition states (TS), intermediate structures and products of 36 37 interest are generated from reactant structures and optimized. Constructing and optimizing 38 39 representative structures along the reaction path is the most time consuming part of the study and needs 40 41 to be carried out very carefully. Fifth, the structures identified in along reaction path are subject to 42 production calculations at a higher level of theory than used for the geometry optimizations. This 43 44 general process is illustrated in Figure 1. 45 46 Based on the mechanism of coupling between the QM and MM region, two types of QM/MM 47 6 48 methods are generally available, additive and subtractive. ONIOM is one popular and robust 49 subtractive QM/MM method. 7,8 The total ONIOM energy for a two layer system is 8 50 51 E ONIOM = E real , MM + E model , QM − E model ,MM 52 53 In a two-layer ONIOM QM/MM calculation, the real system contains all the atoms (including both QM 54 real,MM 55 and MM region) and is calculated at the MM level ( E ). The model system is the QM region treated 56 model,QM 57 at the QM level ( E ). To obtain the total ONIOM energy, the model system also needs to be 58 treated at MM level (Emodel,MM) and be subtracted from real system MM energy. 59 60 3 John Wiley & Sons, Inc. Journal of Computational Chemistry Page 4 of 20 Currently, the QM/MM methods implemented in the Gaussian package 9 are two- and three-layer 1 2 ONIOM models for calculating the energy, gradient and Hessian for geometry optimization and 3 4 vibrational frequency analysis. ONIOM calculations using Gaussian have been applied in numerous 5 10-18 6 studies in computational biology, chemistry and material science. In a subtractive QM/MM method 7 such as ONIOM, the MM force field is needed for entire system, including QM region. The MM force 8 9 field includes parameters for bond and non-bond interactions and partial charges. Due to its subtractive 10 11 nature, most bond and short-range non-bonded terms in the QM region cancel in the ONIOM QM/MM 12 13 method. 14 15 Two schemes are generally applied to treat the electrostatic interaction between the QM and MM 16 regions. 19 In the mechanical embedding (ME) scheme, this interaction is treated at MM level. In the 17 18 electrostatic embedding (EE)For scheme, Peer the partial charges Review of MM atoms are included in the Hamiltonian 19 20 of the QM region. In the ME scheme, the electrostatic interaction between the QM and MM regions do 21 22 not cancel in the QM/MM energy. The partial charges of QM atoms play an important role in this 23 scheme, therefore they need to be treated carefully in QM/MM calculations. 24 25 Since the size of biochemical systems studied by QM/MM can be rather large, some of the structure 26 27 manipulation and processing of these systems can be cumbersome. During a QM/MM calculation, it is a 28 29 good practice to monitor running jobs in terms of both energies and geometries. If problems in running 30 31 jobs can be spotted at an early stage, much time and resources can be saved. However, since the sizes of 32 structures and files in QM/MM studies can be rather large, it is not always convenient to check running 33 34 QM/MM jobs on the fly, especially geometries. Although there are many graphical interfaces available 35 36 to facilitate QM/MM calculations, the efficiency and convenience of QM/MM studies can be 37 38 substantially increased when appropriate auxiliary tools are available for job preparation, monitoring, 39 and analysis. In this paper, we describe a Toolkit to Assist ONIOM calculations (TAO). The toolkit is 40 41 written in PERL using an object oriented design and was developed during our previous QM/MM 42 18 43 studies. With these tools, QM/MM studies using the ONIOM method in Gaussian become much easier 44 45 and more efficient, especially for large enzymatic systems. To illustrate the use of the toolkit and to 46 guide a user through an ONIOM QM/MM study of an enzymatic system, we also provide a tutorial in 47 48 the Supplementary Information. 49 50 In this paper, the general procedure of a QM/MM study is discussed. The tools available for the 51 52 various stages of this study are listed in Table I. 53 54 55 56 57 58 59 60 4 John Wiley & Sons, Inc. Page 5 of 20 Journal of Computational Chemistry ONIOM input preparation 1 2 Since many large scale QM/MM studies are focused on biological systems, such as proteins, this 3 4 paper is directed mainly toward calculations of these systems. However, this toolkit can also be applied 5 6 to any other system suitable for QM/MM studies. The Protein Data Bank (PDB) is a common resource 7 for initial structures of biomolecules. The PDB file format is one of the most commonly used formats 8 9 for biomolecule structure storage.
Recommended publications
  • The Norbornene Mystery Revealed
    Supplementary Material (ESI) for Chemical Communications This journal is (c) The Royal Society of Chemistry 2010 Electronic Supporting Information The Norbornene Mystery Revealed Stephan N. Steinmann, Pierre Vogel, Yirong Mo, and Clémence Corminboeuf* To be published in Chemical Communication. The supplementary Data were prepared on June 22, 2010 and contains 14 pages. Contents: Method Section Figure S1: Molecules used for the comparisons between known experimental 1 Jcc-coupling constants (in Hertz) and the computed values. 1 Table S1 JCC-coupling comparison between experiment and computation. 1 Table S2: JCC-coupling constants for all compounds considered. Cartesian coordinates: B3LYP/6-311+G(d,p) (BLW and canonical) optimized structures for all compounds considered in the article. Contents: Computational Details 1 Table S1: JCC-coupling constants for all compounds considered. 1 Figure S1: Molecules for JCC-coupling comparison between experiment and computation. 1 Table S2: JCC-coupling comparison between experiment and computation. B3LYP/6-311+G(d,p) (BLW and canonical) optimized structures for all compounds considered in the article. Computational Details Standard geometries were optimized at the B3LYP1, 2/6-311+G** level using Gaussian 09.3 The BLW- constrained B3LYP/6-311+G** geometries (‘loc’) were optimized using a modified version of GAMESS-US (release 2008)4interfaced with the BLW-module.5, 6 Both the density and the J-coupling constants have been computed at the PBE7/IGLO-III level in Dalton 2.0.8 The BLW-eigenvalues and eigenvectors were optimized at the given DFT level with the SCF module of Dalton, SIRIUS, that has been interfaced with the BLW- module.
    [Show full text]
  • Computational Chemistry for Chemistry Educators Shawn C
    Volume 1, Issue 1 Journal Of Computational Science Education Computational Chemistry for Chemistry Educators Shawn C. Sendlinger Clyde R. Metz North Carolina Central University College of Charleston Department of Chemistry Department of Chemistry and Biochemistry 1801 Fayetteville Street, 66 George Street, Durham, NC 27707 Charleston, SC 29424 919-530-6297 843-953-8097 [email protected] [email protected] ABSTRACT 1. INTRODUCTION In this paper we describe an ongoing project where the goal is to The majority of today’s students are technologically savvy and are develop competence and confidence among chemistry faculty so often more comfortable using computers than the faculty who they are able to utilize computational chemistry as an effective teach them. In order to harness the student’s interest in teaching tool. Advances in hardware and software have made technology and begin to use it as an educational tool, most faculty research-grade tools readily available to the academic community. members require some level of instruction and training. Because Training is required so that faculty can take full advantage of this chemistry research increasingly utilizes computation as an technology, begin to transform the educational landscape, and important tool, our approach to chemistry education should reflect attract more students to the study of science. this. The ability of computer technology to visualize and manipulate objects on an atomic scale can be a powerful tool to increase both student interest in chemistry as well as their level of Categories and Subject Descriptors understanding. Computational Chemistry for Chemistry Educators (CCCE) is a project that seeks to provide faculty the J.2 [Physical Sciences and Engineering]: Chemistry necessary knowledge, experience, and technology access so that they can begin to design and incorporate computational approaches in the courses they teach.
    [Show full text]
  • Car-Parrinello Molecular Dynamics
    CPMD Car-Parrinello Molecular Dynamics An ab initio Electronic Structure and Molecular Dynamics Program The CPMD consortium WWW: http://www.cpmd.org/ Mailing list: [email protected] E-mail: [email protected] January 29, 2019 Send comments and bug reports to [email protected] This manual is for CPMD version 4.3.0 CPMD 4.3.0 January 29, 2019 Contents I Overview3 1 About this manual3 2 Citation 3 3 Important constants and conversion factors3 4 Recommendations for further reading4 5 History 5 5.1 CPMD Version 1.....................................5 5.2 CPMD Version 2.....................................5 5.2.1 Version 2.0....................................5 5.2.2 Version 2.5....................................5 5.3 CPMD Version 3.....................................5 5.3.1 Version 3.0....................................5 5.3.2 Version 3.1....................................5 5.3.3 Version 3.2....................................5 5.3.4 Version 3.3....................................6 5.3.5 Version 3.4....................................6 5.3.6 Version 3.5....................................6 5.3.7 Version 3.6....................................6 5.3.8 Version 3.7....................................6 5.3.9 Version 3.8....................................6 5.3.10 Version 3.9....................................6 5.3.11 Version 3.10....................................7 5.3.12 Version 3.11....................................7 5.3.13 Version 3.12....................................7 5.3.14 Version 3.13....................................7 5.3.15 Version 3.14....................................7 5.3.16 Version 3.15....................................8 5.3.17 Version 3.17....................................8 5.4 CPMD Version 4.....................................8 5.4.1 Version 4.0....................................8 5.4.2 Version 4.1....................................8 5.4.3 Version 4.3....................................9 6 Installation 10 7 Running CPMD 11 8 Files 12 II Reference Manual 15 9 Input File Reference 15 9.1 Basic rules........................................
    [Show full text]
  • Ab Initio Search for Novel Bxhy Building Blocks with Potential for Hydrogen Storage
    Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 12-2010 Ab Initio Search for Novel BxHy Building Blocks with Potential for Hydrogen Storage Jared K. Olson Utah State University Follow this and additional works at: https://digitalcommons.usu.edu/etd Part of the Physical Chemistry Commons Recommended Citation Olson, Jared K., "Ab Initio Search for Novel BxHy Building Blocks with Potential for Hydrogen Storage" (2010). All Graduate Theses and Dissertations. 844. https://digitalcommons.usu.edu/etd/844 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]. AB INITIO SEARCH FOR NOVEL BXHY BUILDING BLOCKS WITH POTENTIAL FOR HYDROGEN STORAGE by Jared K. Olson A dissertation submitted in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY in Chemistry (Physical Chemistry) Approved: Dr. Alexander I. Boldyrev Dr. Steve Scheiner Major Professor Committee Member Dr. David Farrelly Dr. Stephen Bialkowski Committee Member Committee Member Dr. T.C. Shen Byron Burnham Committee Member Dean of Graduate Studies UTAH STATE UNIVERSITY Logan, Utah 2010 ii Copyright © Jared K. Olson 2010 All Rights Reserved iii ABSTRACT Ab Initio Search for Novel BxHy Building Blocks with Potential for Hydrogen Storage by Jared K. Olson, Doctor of Philosophy Utah State University, 2010 Major Professor: Dr. Alexander I. Boldyrev Department: Chemistry and Biochemistry On-board hydrogen storage presents a challenging barrier to the use of hydrogen as an energy source because the performance of current storage materials falls short of platform requirements.
    [Show full text]
  • Surface Modification of Mgni by Perylene
    Materials Transactions, Vol. 43, No. 11 (2002) pp. 2711 to 2716 Special Issue on Protium New Function in Materials c 2002 The Japan Institute of Metals Surface Modification of MgNi by Perylene Tiejun Ma, Yuji Hatano, Takayuki Abe and Kuniaki Watanabe Hydrogen Isotope Research Center, Toyama University, Toyama 930-8555, Japan Amorphous MgNi was modified by ball milling with perylene and its effects on the charge/discharge capacity were examined by using a conventional two-electrode cell. It was found that both the ball milling time and perylene/MgNi ratio had great influence on the discharge capacity and cycle life of MgNi. Three types of effects were identified, depending on ball-milling conditions. One of them was the increase in the discharge capacity at the first cycle, the second type was the deceleration of the degradation of the discharge capacity with charge/discharge cycle, and the last type was the reduction in the charge/discharge capacity. Chemical states of modified surfaces were analyzed by Auger electron spectroscopy (AES) as well as ab-inito calculation. Both AES and ab-initio calculations indicated that carbon atoms can form bonding with both magnesium and nickel atoms, but bonding with magnesium atoms is most preferable. The change in the charge/discharge capacity is attributed to such kind of reactions, and the three distinct effects are ascribed to the presence of different MgNi-perylene composites, formed during the ball milling on the surface, resulting in the retardation or acceleration of Mg(OH)2 formation on the electrode. (Received May 20, 2002; Accepted July 16, 2002) Keywords: magnesium, nickel, aromatic compound, surface modification, battery, electrode 1.
    [Show full text]
  • Application of Hartree-Fock Method for Modeling of Bioactive Molecules Using SAR and QSPR
    Computational Molecular Bioscience, 2014, 4, 1-24 Published Online March 2014 in SciRes. http://www.scirp.org/journal/cmb http://dx.doi.org/10.4236/cmb.2014.41001 Application of Hartree-Fock Method for Modeling of Bioactive Molecules Using SAR and QSPR Cleydson B. R. Santos1,2,3*, Cleison C. Lobato1, Francinaldo S. Braga1, Sílvia S. S. Morais4, Cesar F. Santos1, Caio P. Fernandes1,3, Davi S. B. Brasil1,5, Lorane I. S. Hage-Melim1, Williams J. C. Macêdo1,2, José C. T. Carvalho1,2,3 1Laboratory of Modeling and Computational Chemistry, Federal University of Amapá, Macapá, Brazil 2Postgraduate Program in Biotechnology and Biodiversity-Network BIONORTE, Campus Universitário Marco Zero, Macapá, Brazil 3Laboratory of Drug Research, School of Pharmaceutical Sciences, Federal University of Amapá, Macapá, Brazil 4Laboratory of Statistical and Computational Modeling, University of the State of Amapá, Macapá, Brazil 5Institute of Technology, Federal University of Pará, Belém, Brazil Email: *[email protected] Received 12 January 2014; revised 12 February 2014; accepted 20 February 2014 Copyright © 2014 by authors and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/ Abstract The central importance of quantum chemistry is to obtain solutions of the Schrödinger equation for the accurate determination of the properties of atomic and molecular systems that occurred from the calculation of wave functions accurate for many diatomic and polyatomic molecules, us- ing Self Consistent Field method (SCF). The application of quantum chemical methods in the study and planning of bioactive compounds has become a common practice nowadays.
    [Show full text]
  • Bringing Open Source to Drug Discovery
    Bringing Open Source to Drug Discovery Chris Swain Cambridge MedChem Consulting Standing on the shoulders of giants • There are a huge number of people involved in writing open source software • It is impossible to acknowledge them all individually • The slide deck will be available for download and includes 25 slides of details and download links – Copy on my website www.cambridgemedchemconsulting.com Why us Open Source software? • Allows access to source code – You can customise the code to suit your needs – If developer ceases trading the code can continue to be developed – Outside scrutiny improves stability and security What Resources are available • Toolkits • Databases • Web Services • Workflows • Applications • Scripts Toolkits • OpenBabel (htttp://openbabel.org) is a chemical toolbox – Ready-to-use programs, and complete programmer's toolkit – Read, write and convert over 110 chemical file formats – Filter and search molecular files using SMARTS and other methods, KNIME add-on – Supports molecular modeling, cheminformatics, bioinformatics – Organic chemistry, inorganic chemistry, solid-state materials, nuclear chemistry – Written in C++ but accessible from Python, Ruby, Perl, Shell scripts… Toolkits • OpenBabel • R • CDK • OpenCL • RDkit • SciPy • Indigo • NumPy • ChemmineR • Pandas • Helium • Flot • FROWNS • GNU Octave • Perlmol • OpenMPI Toolkits • RDKit (http://www.rdkit.org) – A collection of cheminformatics and machine-learning software written in C++ and Python. – Knime nodes – The core algorithms and data structures are written in C ++. Wrappers are provided to use the toolkit from either Python or Java. – Additionally, the RDKit distribution includes a PostgreSQL-based cartridge that allows molecules to be stored in relational database and retrieved via substructure and similarity searches.
    [Show full text]
  • A DFT/ECP-Small Basis Set Modelling of Cisplatin: Molecular Structure and Vibrational Spectrum
    Computational Molecular Bioscience, 2012, 2, 35-44 http://dx.doi.org/10.4236/cmb.2012.22004 Published Online June 2012 (http://www.SciRP.org/journal/cmb) A DFT/ECP-Small Basis Set Modelling of Cisplatin: Molecular Structure and Vibrational Spectrum Nicolay I. Dodoff Acad. R. Tsanev Institute of Molecular Biology, Bulgarian Academy of Sciences, Sofia, Bulgaria Email: [email protected] Received March 25, 2012; revised April 25, 2012; accepted May 2, 2012 ABSTRACT A DFT conformational and vibrational analysis of a single molecule of cisplatin (cis-[Pt(NH3)2Cl2]) was performed by means of PW91 functional and LANL08 ECP basis set for the Pt atom. 3-21G and 3-21G* Basis sets were used for the remaining atoms. All the initially chosen conformations were found to converge to the global minimum conformation of C2v symmetry with H atoms lying in the coordination plane and pointing to the Cl atoms. The computational results were compared with the newest experimental structural data and with the vibrational spectroscopic data for cisplatin, obtained by other workers. The chosen level of theory was found to describe satisfactory the molecular structure (r. m. s. of the relative deviations ≤ 6%) and the harmonic vibrational frequencies (r. m. s. of the relative deviations ≤ 5%) of cisplatin. Keywords: Cisplatin; DFT Calculations; Basis Set; Effective Core Potential; Molecular Structure; Vibrational Analysis 1. Introduction and effectiveness of different computational schemes in predicting molecular geometry and vibrational frequen- Cisplatin (cis-diamminedichloroplatinum(II)) is the first cies [22-25]. inorganic compound introduced in clinical use for the Although the vibrational spectra of cisplatin have been treatment of cancer [1-4].
    [Show full text]
  • 1 Introduction – the Basic Computational Environment This Chapter Provides You with Some Basic Introduction to Your 'Computational Work Desk'
    The Basic Computational Environment 1 1 Introduction – The Basic Computational Environment This chapter provides you with some basic introduction to your 'Computational Work desk'. It is structured into three parts. The first one consists of some quick tutorial how to login into your machine. The second part deals with basic steps in a UNIX-like environment and introduces the 'SHELL', the UNIX command line. In the third part two of the most used editors are presented in order to enable the user to manipulate text files at will. 1.1 Login Into Your Computer Before you can login it is essential that you haven been provided with a username and a password. Usually these are provided by your chemistry department or your system administrator. The first thing you will see is the 'login-screen', asking you to type in your username. Figure 1: a typical login screen The screen might look somewhat different, but it should be similar to the one shown above. Once you typed in your username and pressed the <RETURN> key, the system will ask you for your password, which you are expected to type, again followed by <RETURN>. (The pressing of the <RETURN> key after typing in some commands or other things is assumed from now on.) Having succeeded in providing the correct password you are provided with your graphical desktop. All your actions will take place in this environment. The next thing to do is to open a window containing a command- line, which is usally called a 'shell'. This is nothing else than a window where you can type in your commands.
    [Show full text]
  • Visualization of Polymer Processing at the Continuum Level Jeremy Hicks Clemson University, [email protected]
    Clemson University TigerPrints All Theses Theses 12-2006 Visualization of Polymer Processing at the Continuum Level Jeremy Hicks Clemson University, [email protected] Follow this and additional works at: https://tigerprints.clemson.edu/all_theses Part of the Computer Sciences Commons Recommended Citation Hicks, Jeremy, "Visualization of Polymer Processing at the Continuum Level" (2006). All Theses. 15. https://tigerprints.clemson.edu/all_theses/15 This Thesis is brought to you for free and open access by the Theses at TigerPrints. It has been accepted for inclusion in All Theses by an authorized administrator of TigerPrints. For more information, please contact [email protected]. VISUALIZATION OF POLYMER PROCESSING AT THE CONTINUUM LEVEL A Thesis Presented to the Graduate School of Clemson University In Partial Fulfillment of the Requirements for the Degree Master of Fine Arts Digital Production Arts by Jeremy L. Hicks December 2006 Accepted by: Dr. Timothy Davis, Committee Chair Dr. John Kundert-Gibbs Dr. Christopher Cox ABSTRACT Computer animation, coupled with scientific experimentation and modeling, allows scientists to produce detailed visualizations that potentially enable more comprehensive perception of physical phenomena and ultimately, new discoveries. With the use of Maya, an animation and modeling program that incorporates the natural laws of physics to control the behavior of virtual objects in computer animation, data from the modeling of physical processes such as polymer fibers and films can be explored in the visual realm. Currently, few attempts have been made at the continuum level to represent polymer properties via computer animation using advanced graphics. As a result, scientists may be unable to recognize patterns and trends in a specific polymer quickly and efficiently, and thus, lose time and money commonly required for further experiments.
    [Show full text]
  • An Interdisciplinary Study on Monosubstituted Benzene Involving Computation, Statistics and Chemistry
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Elsevier - Publisher Connector PARA SABER, EXPERIMENTAR Y SIMULAR Educ. quím., 25(4), 418–424, 2014. © Universidad Nacional Autónoma de México, ISSN 0187-893-X Publicado en línea el 8 de agosto de 2014, ISSNE 1870-8404 An interdisciplinary study on monosubstituted benzene involving computation, statistics and chemistry João E. V. Ferreira,¹ Clauber H. S. da Costa,¹ Ricardo Morais de Miranda,¹ Antonio F. Figueiredo¹ and Yoanna M. A. Ginarte² ABSTRACT Computation and statistics have given a great contribution to chemical research and education be- cause both sciences make possible to deal with many different types of information about atoms and molecules. This paper describes the use of computation and statistics to study 20 different groups in monosubstituted benzene. Firstly, compounds had their geometries optimized and molecular de- scriptors were computed. Then, two exploratory methods, principal component analysis and hierar- chical cluster analysis, were employed to classify these groups according to their activating and deac- tivating effects on rate of electrophilic aromatic substitution reactions. The pedagogical objective is to show an interdisciplinary application in organic chemistry and motivate students and teachers to apply this strategy in the classroom. KEYWORDS: graduate education, interdisciplinary, computer-based learning, statistics, chemo- metrics Resumen (Un estudio interdisciplinario sobre el benceno monosustituido empleando cómputo, estadística y química) El cómputo y la estadística han dado una gran contribución a la investigación y la educación quími- cas, porque ambas han hecho posible el manejo de diversos tipos de información sobre átomos y moléculas.
    [Show full text]
  • Nwchem: Past, Present, and Future
    1 NWChem: Past, Present, and Future NWChem: Past, Present, and Future a) E. Apra,1 E. J. Bylaska,1 W. A. de Jong,2 N. Govind,1 K. Kowalski,1, T. P. Straatsma,3 M. Valiev,1 H. J. 4 5 6 5 7 8 9 J. van Dam, Y. Alexeev, J. Anchell, V. Anisimov, F. W. Aquino, R. Atta-Fynn, J. Autschbach, N. P. Bauman,1 J. C. Becca,10 D. E. Bernholdt,11 K. Bhaskaran-Nair,12 S. Bogatko,13 P. Borowski,14 J. Boschen,15 J. Brabec,16 A. Bruner,17 E. Cauet,18 Y. Chen,19 G. N. Chuev,20 C. J. Cramer,21 J. Daily,1 22 23 9 24 11 25 M. J. O. Deegan, T. H. Dunning Jr., M. Dupuis, K. G. Dyall, G. I. Fann, S. A. Fischer, A. Fonari,26 H. Fruchtl,27 L. Gagliardi,21 J. Garza,28 N. Gawande,1 S. Ghosh,29 K. Glaesemann,1 A. W. Gotz,30 6 31 32 33 34 2 10 J. Hammond, V. Helms, E. D. Hermes, K. Hirao, S. Hirata, M. Jacquelin, L. Jensen, B. G. Johnson,35 H. Jonsson,36 R. A. Kendall,11 M. Klemm,6 R. Kobayashi,37 V. Konkov,38 S. Krishnamoorthy,1 M. 19 39 40 41 42 43 44 Krishnan, Z. Lin, R. D. Lins, R. J. Little eld, A. J. Logsdail, K. Lopata, W. Ma, A. V. Marenich,45 J. Martin del Campo,46 D. Mejia-Rodriguez,47 J. E. Moore,6 J. M. Mullin,48 T. Nakajima,49 D. R. Nascimento,1 J.
    [Show full text]