Benchmarking and Application of Density Functional Methods In

Total Page:16

File Type:pdf, Size:1020Kb

Benchmarking and Application of Density Functional Methods In BENCHMARKING AND APPLICATION OF DENSITY FUNCTIONAL METHODS IN COMPUTATIONAL CHEMISTRY by BRIAN N. PAPAS (Under Direction the of Henry F. Schaefer III) ABSTRACT Density Functional methods were applied to systems of chemical interest. First, the effects of integration grid quadrature choice upon energy precision were documented. This was done through application of DFT theory as implemented in five standard computational chemistry programs to a subset of the G2/97 test set of molecules. Subsequently, the neutral hydrogen-loss radicals of naphthalene, anthracene, tetracene, and pentacene and their anions where characterized using five standard DFT treatments. The global and local electron affinities were computed for the twelve radicals. The results for the 1- naphthalenyl and 2-naphthalenyl radicals were compared to experiment, and it was found that B3LYP appears to be the most reliable functional for this type of system. For the larger systems the predicted site specific adiabatic electron affinities of the radicals are 1.51 eV (1-anthracenyl), 1.46 eV (2-anthracenyl), 1.68 eV (9-anthracenyl); 1.61 eV (1-tetracenyl), 1.56 eV (2-tetracenyl), 1.82 eV (12-tetracenyl); 1.93 eV (14-pentacenyl), 2.01 eV (13-pentacenyl), 1.68 eV (1-pentacenyl), and 1.63 eV (2-pentacenyl). The global minimum for each radical does not have the same hydrogen removed as the global minimum for the analogous anion. With this in mind, the global (or most preferred site) adiabatic electron affinities are 1.37 eV (naphthalenyl), 1.64 eV (anthracenyl), 1.81 eV (tetracenyl), and 1.97 eV (pentacenyl). In later work, ten (scandium through zinc) homonuclear transition metal trimers were studied using one DFT 2 functional. Electronic ground states and bond lengths for equilateral triangles were: Sc3 ( A1′ , 2.83 Å), Ti3 7 2 17 16 11 6 3 ( E ′ , 2.32), V3 ( E ″ , 2.06), Cr3 ( E ′ , 2.92), Mn3 ( A2′ , 2.73), Fe3 ( E ″ , 2.24), Co3 ( E ″ , 2.18), Ni3 ( A2″ , 2 1 2.23), Cu3 ( E ′ , 2.37), Zn3 ( A1′ , 2.93). Vibrational frequencies, several low-lying electronic states, and trends in bond lengths and atomization energies are discussed. The predicted dissociation energies ∆E -1 -1 -1 -1 (M3 → M2 + M) are 49.4 kcal mol (Sc3), 64.3 kcal mol (Ti3), 60.7 kcal mol (V3), 11.5 kcal mol (Cr3), -1 -1 -1 -1 -1 32.4 kcal mol (Mn3), 61.5 kcal mol (Fe3), 78.0 kcal mol (Co3), 86.1 kcal mol (Ni3), 26.8 kcal mol -1 (Cu3), and 4.5 kcal mol (Zn3). INDEX WORDS: Density Functional Theory, precision, G2-97, Polycyclic Aromatic Hydrocarbon, PAH, electron affinity, transition metal trimers, M3. BENCHMARKING AND APPLICATION OF DENSITY FUNCTIONAL METHODS IN COMPUTATIONAL CHEMISTRY by BRIAN N. PAPAS B.S., Michigan State University, 2002 A Dissertation Submitted to the Graduate Faculty of The University of Georgia in Partial Fulfillment of the Requirements for the Degree DOCTOR OF PHILOSOPHY ATHENS, GEORGIA 2006 © 2006 Brian N. Papas All Rights Reserved BENCHMARKING AND APPLICATION OF DENSITY FUNCTIONAL METHODS IN COMPUTATIONAL CHEMISTRY by BRIAN N. PAPAS Major Professor: Henry F. Schaefer III Committee: Geoff Smith Nigel Adams Electronic Version Approved: Maureen Grasso Dean of the Graduate School The University of Georgia May 2006 iv TABLE OF CONTENTS Page LIST OF TABLES......................................................................................................................... vi LIST OF FIGURES ...................................................................................................................... vii CHAPTER 1 INTRODUCTION AND BACKGROUND MATERIAL.............................................1 1.1 INTRODUCTION .....................................................................................................1 1.2 THEORY AND IMPLEMENTATIONS .........................................................................1 1.3 PROSPECTUS .........................................................................................................3 1.4 REFERENCES.........................................................................................................3 2 CONCERNING THE PRECISION OF STANDARD DENSITY FUNCTIONAL PROGRAMS: GAUSSIAN, MOLPRO, NWCHEM, Q-CHEM, AND GAMESS .4 2.1 ABSTRACT............................................................................................................5 2.2 INTRODUCTION.....................................................................................................5 2.3 METHODS .............................................................................................................5 2.4 RESULTS...............................................................................................................7 2.5 CONCLUSIONS ....................................................................................................10 2.6 ACKNOWLEDGEMENTS .......................................................................................10 2.7 TABLES AND FIGURES ........................................................................................10 2.7 REFERENCES.......................................................................................................21 3 THE NAPHTHALENYL, ANTHRACENYL, TETRACENYL, AND PENTACENYL RADICALS, AND THEIR ANIONS ..........................................23 v 3.1 ABSTRACT..........................................................................................................24 3.2 INTRODUCTION...................................................................................................24 3.3 METHODS ...........................................................................................................25 3.4 RESULTS AND DISCUSSION .................................................................................27 3.5 CONCLUSIONS ....................................................................................................31 3.6 ACKNOWLEDGMENTS .........................................................................................32 3.7 TABLES AND FIGURES ........................................................................................32 3.8 REFERENCES.......................................................................................................37 4 HOMONUCLEAR TRANSITION METAL TRIMERS ............................................40 4.1 ABSTRACT..........................................................................................................41 4.2 INTRODUCTION...................................................................................................41 4.3 METHODS ...........................................................................................................42 4.4 RESULTS.............................................................................................................44 4.5 CONCLUSIONS ....................................................................................................61 4.6 ACKNOWLEDGEMENTS .......................................................................................62 4.7 TABLES AND FIGURES ........................................................................................62 4.7 REFERENCES.......................................................................................................67 5 CONCLUSIONS..........................................................................................................79 vi LIST OF TABLES Page Table 2.1: G2/97 subset of neutral molecules considered in the present research.. ......................10 Table 2.2: Specification of the numerical integration grids used in this study..............................11 Table 2.3: “Exact” energies (Hartrees), computed as averages of energies from all programs with the “Large” grid quadratures...........................................................................................12 Table 3.1: Ipso C-C-C angles (in degrees) for the radicals and their anions.................................33 Table 3.2: C-C bond lengths (in Å) adjacent to the removed hydrogen atom...............................34 Table 3.3: EAs (in eV) with ZPVE corrected values in parentheses.............................................35 Table 3.4: ZPVE corrected energies (kcal/mol) relative to lowest energy species.. .....................36 Table 3.5: Global AEAs (eV), with ZPVE corrected values in parentheses.. ...............................36 Table 4.1: Low-lying electronic states of first row transition metal trimers.. ..............................62 Table 4.2: Harmonic vibrational frequencies (in cm-1) for transition metal trimers.....................63 Table 4.3: Summary of electronic ground states for the transition metal trimers in D3h symmetry.........................................................................................................................64 vii LIST OF FIGURES Page Figure 2.1: Gamess error curves. ...................................................................................................14 Figure 2.2: Molpro error curves.....................................................................................................15 Figure 2.3: NWChem error curves.................................................................................................16 Figure 2.4: QChem error curves. ...................................................................................................17 Figure 2.5: Gaussian 94 error curves.. ...........................................................................................18 Figure 2.6: BLYP results by molecule type...................................................................................19
Recommended publications
  • United States Patent Office Patented June 17, 1969
    3,450,782 United States Patent Office Patented June 17, 1969 2 1,3-dibromocyclobutane with sodium in refluxing dioxane, 3,450,782 PROCESS FOR THE PREPARATION OF K. B. Wiberg, G. M. Lampman, R. P. Ciula, D. S. Con CYCLICALKANES nor, P. Scherter, and J. Lavanesh, Tetrahedron, 21, 2749 Daniel S. Connor, Cincinnati, Ohio, assignor to The (1695), bicyclo[1.1.1 pentane by treatment of 3-(bromo Procter & Gamble Company, Cincinnati, Ohio, a cor methyl)-cyclobutyl bromide with sodium metal at a 0.5% poration of Ohio yield, with lithium amalgam in refluxing dioxane at a No Drawing. Filed Nov. 29, 1967, Ser. No. 686,738 4.2% yield and with sodium/naphtahalene at a 8% yield, Int, C. C07c I/28 K. B. Wiberg, D. S. Connor, and G. M. Lampman, Tetra U.S. C. 260-666 8 Claims hedron Letters, 531 (1964); K. B. Wiberg and D. S. Con 0 nor, J. Am. Chem. Soc., 88, 4437 (1966). On examination of the literature hereinbefore cited ABSTRACT OF THE DISCLOSURE on cyclization, it is apparent that the yields obtained were This invention concerns the preparation of cyclic al extremely low, that the conditions for reaction when it kanes from dihalogenated alkanes using lithium amalgam. did occur were quite severe, and that the starting materials 5 used were less than common. OBJECTS OF THE INVENTION SUMMARY OF THE INVENTION The use of the method of this invention to obtain cyclic The object of this invention is to prepare cyclic al compounds provides a valuable synthetic tool heretofore kanes from the halogenated straight chain starting ma 20 unknown to chemists, thus enabling the production of terials using lithium amalgam to remove the halogens valuable end products.
    [Show full text]
  • Supporting Information
    Electronic Supplementary Material (ESI) for RSC Advances. This journal is © The Royal Society of Chemistry 2020 Supporting Information How to Select Ionic Liquids as Extracting Agent Systematically? Special Case Study for Extractive Denitrification Process Shurong Gaoa,b,c,*, Jiaxin Jina,b, Masroor Abroc, Ruozhen Songc, Miao Hed, Xiaochun Chenc,* a State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Beijing, 102206, China b Research Center of Engineering Thermophysics, North China Electric Power University, Beijing, 102206, China c Beijing Key Laboratory of Membrane Science and Technology & College of Chemical Engineering, Beijing University of Chemical Technology, Beijing 100029, PR China d Office of Laboratory Safety Administration, Beijing University of Technology, Beijing 100124, China * Corresponding author, Tel./Fax: +86-10-6443-3570, E-mail: [email protected], [email protected] 1 COSMO-RS Computation COSMOtherm allows for simple and efficient processing of large numbers of compounds, i.e., a database of molecular COSMO files; e.g. the COSMObase database. COSMObase is a database of molecular COSMO files available from COSMOlogic GmbH & Co KG. Currently COSMObase consists of over 2000 compounds including a large number of industrial solvents plus a wide variety of common organic compounds. All compounds in COSMObase are indexed by their Chemical Abstracts / Registry Number (CAS/RN), by a trivial name and additionally by their sum formula and molecular weight, allowing a simple identification of the compounds. We obtained the anions and cations of different ILs and the molecular structure of typical N-compounds directly from the COSMObase database in this manuscript.
    [Show full text]
  • Automated Construction of Quantum–Classical Hybrid Models Arxiv:2102.09355V1 [Physics.Chem-Ph] 18 Feb 2021
    Automated construction of quantum{classical hybrid models Christoph Brunken and Markus Reiher∗ ETH Z¨urich, Laboratorium f¨urPhysikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Z¨urich, Switzerland February 18, 2021 Abstract We present a protocol for the fully automated construction of quantum mechanical-(QM){ classical hybrid models by extending our previously reported approach on self-parametri- zing system-focused atomistic models (SFAM) [J. Chem. Theory Comput. 2020, 16 (3), 1646{1665]. In this QM/SFAM approach, the size and composition of the QM region is evaluated in an automated manner based on first principles so that the hybrid model describes the atomic forces in the center of the QM region accurately. This entails the au- tomated construction and evaluation of differently sized QM regions with a bearable com- putational overhead that needs to be paid for automated validation procedures. Applying SFAM for the classical part of the model eliminates any dependence on pre-existing pa- rameters due to its system-focused quantum mechanically derived parametrization. Hence, QM/SFAM is capable of delivering a high fidelity and complete automation. Furthermore, since SFAM parameters are generated for the whole system, our ansatz allows for a con- venient re-definition of the QM region during a molecular exploration. For this purpose, a local re-parametrization scheme is introduced, which efficiently generates additional clas- sical parameters on the fly when new covalent bonds are formed (or broken) and moved to the classical region. arXiv:2102.09355v1 [physics.chem-ph] 18 Feb 2021 ∗Corresponding author; e-mail: [email protected] 1 1 Introduction In contrast to most protocols of computational quantum chemistry that consider isolated molecules, chemical processes can take place in a vast variety of complex environments.
    [Show full text]
  • Popular GPU-Accelerated Applications
    LIFE & MATERIALS SCIENCES GPU-ACCELERATED APPLICATIONS | CATALOG | AUG 12 LIFE & MATERIALS SCIENCES APPLICATIONS CATALOG Application Description Supported Features Expected Multi-GPU Release Status Speed Up* Support Bioinformatics BarraCUDA Sequence mapping software Alignment of short sequencing 6-10x Yes Available now reads Version 0.6.2 CUDASW++ Open source software for Smith-Waterman Parallel search of Smith- 10-50x Yes Available now protein database searches on GPUs Waterman database Version 2.0.8 CUSHAW Parallelized short read aligner Parallel, accurate long read 10x Yes Available now aligner - gapped alignments to Version 1.0.40 large genomes CATALOG GPU-BLAST Local search with fast k-tuple heuristic Protein alignment according to 3-4x Single Only Available now blastp, multi cpu threads Version 2.2.26 GPU-HMMER Parallelized local and global search with Parallel local and global search 60-100x Yes Available now profile Hidden Markov models of Hidden Markov Models Version 2.3.2 mCUDA-MEME Ultrafast scalable motif discovery algorithm Scalable motif discovery 4-10x Yes Available now based on MEME algorithm based on MEME Version 3.0.12 MUMmerGPU A high-throughput DNA sequence alignment Aligns multiple query sequences 3-10x Yes Available now LIFE & MATERIALS& LIFE SCIENCES APPLICATIONS program against reference sequence in Version 2 parallel SeqNFind A GPU Accelerated Sequence Analysis Toolset HW & SW for reference 400x Yes Available now assembly, blast, SW, HMM, de novo assembly UGENE Opensource Smith-Waterman for SSE/CUDA, Fast short
    [Show full text]
  • Starting SCF Calculations by Superposition of Atomic Densities
    Starting SCF Calculations by Superposition of Atomic Densities J. H. VAN LENTHE,1 R. ZWAANS,1 H. J. J. VAN DAM,2 M. F. GUEST2 1Theoretical Chemistry Group (Associated with the Department of Organic Chemistry and Catalysis), Debye Institute, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands 2CCLRC Daresbury Laboratory, Daresbury WA4 4AD, United Kingdom Received 5 July 2005; Accepted 20 December 2005 DOI 10.1002/jcc.20393 Published online in Wiley InterScience (www.interscience.wiley.com). Abstract: We describe the procedure to start an SCF calculation of the general type from a sum of atomic electron densities, as implemented in GAMESS-UK. Although the procedure is well known for closed-shell calculations and was already suggested when the Direct SCF procedure was proposed, the general procedure is less obvious. For instance, there is no need to converge the corresponding closed-shell Hartree–Fock calculation when dealing with an open-shell species. We describe the various choices and illustrate them with test calculations, showing that the procedure is easier, and on average better, than starting from a converged minimal basis calculation and much better than using a bare nucleus Hamiltonian. © 2006 Wiley Periodicals, Inc. J Comput Chem 27: 926–932, 2006 Key words: SCF calculations; atomic densities Introduction hrstuhl fur Theoretische Chemie, University of Kahrlsruhe, Tur- bomole; http://www.chem-bio.uni-karlsruhe.de/TheoChem/turbo- Any quantum chemical calculation requires properly defined one- mole/),12 GAMESS(US) (Gordon Research Group, GAMESS, electron orbitals. These orbitals are in general determined through http://www.msg.ameslab.gov/GAMESS/GAMESS.html, 2005),13 an iterative Hartree–Fock (HF) or Density Functional (DFT) pro- Spartan (Wavefunction Inc., SPARTAN: http://www.wavefun.
    [Show full text]
  • D:\Doc\Workshops\2005 Molecular Modeling\Notebook Pages\Software Comparison\Summary.Wpd
    CAChe BioRad Spartan GAMESS Chem3D PC Model HyperChem acd/ChemSketch GaussView/Gaussian WIN TTTT T T T T T mac T T T (T) T T linux/unix U LU LU L LU Methods molecular mechanics MM2/MM3/MM+/etc. T T T T T Amber T T T other TT T T T T semi-empirical AM1/PM3/etc. T T T T T (T) T Extended Hückel T T T T ZINDO T T T ab initio HF * * T T T * T dft T * T T T * T MP2/MP4/G1/G2/CBS-?/etc. * * T T T * T Features various molecular properties T T T T T T T T T conformer searching T T T T T crystals T T T data base T T T developer kit and scripting T T T T molecular dynamics T T T T molecular interactions T T T T movies/animations T T T T T naming software T nmr T T T T T polymers T T T T proteins and biomolecules T T T T T QSAR T T T T scientific graphical objects T T spectral and thermodynamic T T T T T T T T transition and excited state T T T T T web plugin T T Input 2D editor T T T T T 3D editor T T ** T T text conversion editor T protein/sequence editor T T T T various file formats T T T T T T T T Output various renderings T T T T ** T T T T various file formats T T T T ** T T T animation T T T T T graphs T T spreadsheet T T T * GAMESS and/or GAUSSIAN interface ** Text only.
    [Show full text]
  • Modern Quantum Chemistry with [Open]Molcas
    Modern quantum chemistry with [Open]Molcas Cite as: J. Chem. Phys. 152, 214117 (2020); https://doi.org/10.1063/5.0004835 Submitted: 17 February 2020 . Accepted: 11 May 2020 . Published Online: 05 June 2020 Francesco Aquilante , Jochen Autschbach , Alberto Baiardi , Stefano Battaglia , Veniamin A. Borin , Liviu F. Chibotaru , Irene Conti , Luca De Vico , Mickaël Delcey , Ignacio Fdez. Galván , Nicolas Ferré , Leon Freitag , Marco Garavelli , Xuejun Gong , Stefan Knecht , Ernst D. Larsson , Roland Lindh , Marcus Lundberg , Per Åke Malmqvist , Artur Nenov , Jesper Norell , Michael Odelius , Massimo Olivucci , Thomas B. Pedersen , Laura Pedraza-González , Quan M. Phung , Kristine Pierloot , Markus Reiher , Igor Schapiro , Javier Segarra-Martí , Francesco Segatta , Luis Seijo , Saumik Sen , Dumitru-Claudiu Sergentu , Christopher J. Stein , Liviu Ungur , Morgane Vacher , Alessio Valentini , and Valera Veryazov J. Chem. Phys. 152, 214117 (2020); https://doi.org/10.1063/5.0004835 152, 214117 © 2020 Author(s). The Journal ARTICLE of Chemical Physics scitation.org/journal/jcp Modern quantum chemistry with [Open]Molcas Cite as: J. Chem. Phys. 152, 214117 (2020); doi: 10.1063/5.0004835 Submitted: 17 February 2020 • Accepted: 11 May 2020 • Published Online: 5 June 2020 Francesco Aquilante,1,a) Jochen Autschbach,2,b) Alberto Baiardi,3,c) Stefano Battaglia,4,d) Veniamin A. Borin,5,e) Liviu F. Chibotaru,6,f) Irene Conti,7,g) Luca De Vico,8,h) Mickaël Delcey,9,i) Ignacio Fdez. Galván,4,j) Nicolas Ferré,10,k) Leon Freitag,3,l) Marco Garavelli,7,m) Xuejun Gong,11,n) Stefan Knecht,3,o) Ernst D. Larsson,12,p) Roland Lindh,4,q) Marcus Lundberg,9,r) Per Åke Malmqvist,12,s) Artur Nenov,7,t) Jesper Norell,13,u) Michael Odelius,13,v) Massimo Olivucci,8,14,w) Thomas B.
    [Show full text]
  • Revised Group Additivity Values for Enthalpies of Formation (At 298 K) of Carbon– Hydrogen and Carbon–Hydrogen–Oxygen Compounds
    Revised Group Additivity Values for Enthalpies of Formation (at 298 K) of Carbon– Hydrogen and Carbon–Hydrogen–Oxygen Compounds Cite as: Journal of Physical and Chemical Reference Data 25, 1411 (1996); https://doi.org/10.1063/1.555988 Submitted: 17 January 1996 . Published Online: 15 October 2009 N. Cohen ARTICLES YOU MAY BE INTERESTED IN Additivity Rules for the Estimation of Molecular Properties. Thermodynamic Properties The Journal of Chemical Physics 29, 546 (1958); https://doi.org/10.1063/1.1744539 Critical Evaluation of Thermochemical Properties of C1–C4 Species: Updated Group- Contributions to Estimate Thermochemical Properties Journal of Physical and Chemical Reference Data 44, 013101 (2015); https:// doi.org/10.1063/1.4902535 Estimation of the Thermodynamic Properties of Hydrocarbons at 298.15 K Journal of Physical and Chemical Reference Data 17, 1637 (1988); https:// doi.org/10.1063/1.555814 Journal of Physical and Chemical Reference Data 25, 1411 (1996); https://doi.org/10.1063/1.555988 25, 1411 © 1996 American Institute of Physics for the National Institute of Standards and Technology. Revised Group Additivity Values for Enthalpies of Formation (at 298 K) of Carbon-Hydrogen and Carbon-Hydrogen-Oxygen Compounds N. Cohen Thermochemical Kinetics Research, 6507 SE 31st Avenue, Portland, Oregon 97202-8627 Received January 17, 1996; revised manuscript received September 4, 1996 A program has been undertaken for the evaluation and revision of group additivity values (GAVs) necessary for predicting, by means of Benson's group additivity method, thermochemical properties of organic molecules. This review reports on the portion of that program dealing with GAVs for enthalpies of formation at 298.15 K (hereinafter abbreviated as 298 K) for carbon-hydrogen and carbon-hydrogen-oxygen compounds.
    [Show full text]
  • A Summary of ERCAP Survey of the Users of Top Chemistry Codes
    A survey of codes and algorithms used in NERSC chemical science allocations Lin-Wang Wang NERSC System Architecture Team Lawrence Berkeley National Laboratory We have analyzed the codes and their usages in the NERSC allocations in chemical science category. This is done mainly based on the ERCAP NERSC allocation data. While the MPP hours are based on actually ERCAP award for each account, the time spent on each code within an account is estimated based on the user’s estimated partition (if the user provided such estimation), or based on an equal partition among the codes within an account (if the user did not provide the partition estimation). Everything is based on 2007 allocation, before the computer time of Franklin machine is allocated. Besides the ERCAP data analysis, we have also conducted a direct user survey via email for a few most heavily used codes. We have received responses from 10 users. The user survey not only provide us with the code usage for MPP hours, more importantly, it provides us with information on how the users use their codes, e.g., on which machine, on how many processors, and how long are their simulations? We have the following observations based on our analysis. (1) There are 48 accounts under chemistry category. This is only second to the material science category. The total MPP allocation for these 48 accounts is 7.7 million hours. This is about 12% of the 66.7 MPP hours annually available for the whole NERSC facility (not accounting Franklin). The allocation is very tight. The majority of the accounts are only awarded less than half of what they requested for.
    [Show full text]
  • Porting the DFT Code CASTEP to Gpgpus
    Porting the DFT code CASTEP to GPGPUs Toni Collis [email protected] EPCC, University of Edinburgh CASTEP and GPGPUs Outline • Why are we interested in CASTEP and Density Functional Theory codes. • Brief introduction to CASTEP underlying computational problems. • The OpenACC implementation http://www.nu-fuse.com CASTEP: a DFT code • CASTEP is a commercial and academic software package • Capable of Density Functional Theory (DFT) and plane wave basis set calculations. • Calculates the structure and motions of materials by the use of electronic structure (atom positions are dictated by their electrons). • Modern CASTEP is a re-write of the original serial code, developed by Universities of York, Durham, St. Andrews, Cambridge and Rutherford Labs http://www.nu-fuse.com CASTEP: a DFT code • DFT/ab initio software packages are one of the largest users of HECToR (UK national supercomputing service, based at University of Edinburgh). • Codes such as CASTEP, VASP and CP2K. All involve solving a Hamiltonian to explain the electronic structure. • DFT codes are becoming more complex and with more functionality. http://www.nu-fuse.com HECToR • UK National HPC Service • Currently 30- cabinet Cray XE6 system – 90,112 cores • Each node has – 2×16-core AMD Opterons (2.3GHz Interlagos) – 32 GB memory • Peak of over 800 TF and 90 TB of memory http://www.nu-fuse.com HECToR usage statistics Phase 3 statistics (Nov 2011 - Apr 2013) Ab initio codes (VASP, CP2K, CASTEP, ONETEP, NWChem, Quantum Espresso, GAMESS-US, SIESTA, GAMESS-UK, MOLPRO) GS2NEMO ChemShell 2%2% SENGA2% 3% UM Others 4% 34% MITgcm 4% CASTEP 4% GROMACS 6% DL_POLY CP2K VASP 5% 8% 19% http://www.nu-fuse.com HECToR usage statistics Phase 3 statistics (Nov 2011 - Apr 2013) 35% of the Chemistry software on HECToR is using DFT methods.
    [Show full text]
  • The Molpro Quantum Chemistry Package
    The Molpro Quantum Chemistry package Hans-Joachim Werner,1, a) Peter J. Knowles,2, b) Frederick R. Manby,3, c) Joshua A. Black,1, d) Klaus Doll,1, e) Andreas Heßelmann,1, f) Daniel Kats,4, g) Andreas K¨ohn,1, h) Tatiana Korona,5, i) David A. Kreplin,1, j) Qianli Ma,1, k) Thomas F. Miller, III,6, l) Alexander Mitrushchenkov,7, m) Kirk A. Peterson,8, n) Iakov Polyak,2, o) 1, p) 2, q) Guntram Rauhut, and Marat Sibaev 1)Institut f¨ur Theoretische Chemie, Universit¨at Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany 2)School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom 3)School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom 4)Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany 5)Faculty of Chemistry, University of Warsaw, L. Pasteura 1 St., 02-093 Warsaw, Poland 6)Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States 7)MSME, Univ Gustave Eiffel, UPEC, CNRS, F-77454, Marne-la- Vall´ee, France 8)Washington State University, Department of Chemistry, Pullman, WA 99164-4630 1 Molpro is a general purpose quantum chemistry software package with a long devel- opment history. It was originally focused on accurate wavefunction calculations for small molecules, but now has many additional distinctive capabilities that include, inter alia, local correlation approximations combined with explicit correlation, highly efficient implementations of single-reference correlation methods, robust and efficient multireference methods for large molecules, projection embedding and anharmonic vibrational spectra.
    [Show full text]
  • Computational Chemistry
    Computational Chemistry as a Teaching Aid Melissa Boule Aaron Harshman Rexford Hoadley An Interactive Qualifying Project Report Submitted to the Faculty of the Worcester Polytechnic Institute In partial fulfillment of the requirements for the Bachelor of Science Degree Approved by: Professor N. Aaron Deskins 1i Abstract Molecular modeling software has transformed the capabilities of researchers. Molecular modeling software has the potential to be a helpful teaching tool, though as of yet its efficacy as a teaching technique has yet to be proven. This project focused on determining the effectiveness of a particular molecular modeling software in high school classrooms. Our team researched what topics students struggled with, surveyed current high school chemistry teachers, chose a modeling software, and then developed a lesson plan around these topics. Lastly, we implemented the lesson in two separate high school classrooms. We concluded that these programs show promise for the future, but with the current limitations of high school technology, these tools may not be as impactful. 2ii Acknowledgements This project would not have been possible without the help of many great educators. Our advisor Professor N. Aaron Deskins provided driving force to ensure the research ran smoothly and efficiently. Mrs. Laurie Hanlan and Professor Drew Brodeur provided insight and guidance early in the development of our topic. Mr. Mark Taylor set-up our WebMO server and all our other technology needs. Lastly, Mrs. Pamela Graves and Mr. Eric Van Inwegen were instrumental by allowing us time in their classrooms to implement our lesson plans. Without the support of this diverse group of people this project would never have been possible.
    [Show full text]