Density Functional Theory for Transition Metals and Transition Metal Chemistry

Total Page:16

File Type:pdf, Size:1020Kb

Density Functional Theory for Transition Metals and Transition Metal Chemistry REVIEW ARTICLE www.rsc.org/pccp | Physical Chemistry Chemical Physics Density functional theory for transition metals and transition metal chemistry Christopher J. Cramer* and Donald G. Truhlar* Received 8th April 2009, Accepted 20th August 2009 First published as an Advance Article on the web 21st October 2009 DOI: 10.1039/b907148b We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids. 1. Introduction chemical systems, in part because its cost scales more favorably with system size than does the cost of correlated WFT, and yet Density functional theory (DFT) describes the electronic states it competes well in accuracy except for very small systems. This of atoms, molecules, and materials in terms of the three- is true even in organic chemistry, but the advantages of DFT dimensional electronic density of the system, which is a great are still greater for metals, especially transition metals. The simplification over wave function theory (WFT), which involves reason for this added advantage is static electron correlation. a3N-dimensional antisymmetric wave function for a system It is now well appreciated that quantitatively accurate 1 with N electrons. Although DFT is sometimes considered the electronic structure calculations must include electron correlation. B ‘‘new kid on the block’’, it is now 45 years old in its modern It is convenient to recognize two types of electron correlation, 2 formulation (more than half as old as quantum mechanics the first called dynamical electron correlation and the second 3,4 itself), and it has roots that are almost as ancient as the called static correlation, near-degeneracy correlation, or non- Schro¨ dinger equation. In practical work, DFT is almost always dynamical correlation. Dynamical correlation is a short-range 5 applied in the form introduced by Kohn and Sham, including effect by which electrons avoid one another to reduce electron 6,7 its spin-polarized extension. The basic quantity in DFT is the repulsion. It is a very general effect for all finite systems many-electron spin density, r. The spin-polarized Kohn–Sham containing two or more electrons. Accounting for dynamical formalism involves a determinant formed from a set of N correlation by a configuration interaction wave function is fictitious single-particle spin-orbitals corresponding to a non- very slowly convergent and requires a very large number of interacting system of electrons with the same spin densities, ra configurations. Other correlation effects, which are very system and rb, as the real system, where r is the sum of ra and rb,the specific and can be either medium ranged or long ranged, can spin density ra is the 3-dimensional electron density of all be accounted for to a large extent by mixing a small number spin-up electrons, and rb is the same for spin-down electrons. (sometimes two, sometimes more) of configurations that are In the original Kohn–Sham formalism (applicable to closed- ‘‘nearly’’ degenerate.8,9 Such correlation effects are called shell molecules and nonmagnetic solids), ra equals rb.(The static or near-degeneracy correlation and systems exhibiting original Kohn–Sham formalism may also be labeled spin- significant static correlation effects are often called multi- restricted Kohn–Sham or restricted Kohn–Sham, and the reference systems;10 likewise, WFT methods based on spin-polarized version may be called spin-unrestricted or multi-configurational zero-order states are often called multi- unrestricted.) We note that ‘‘spin density’’ is a generic term for reference methods. Due to partially filled d subshells, and the density associated with the subset of electrons characterized nearly degenerate (n +1)s and nd subshells, systems containing by the same definite value of Sz, i.e.,eithera or b (thus one might transition metals often have a plethora of low-lying nearly say ‘‘spin-up density and spin-down density’’ rather than ‘‘spin degenerate states, and near-degeneracy correlation effects on densities’’). In a many-electron system comprised of both spin-up the ground-state structure, electron distribution, and energy of and spin-down electrons, the term spin-density is also sometimes transition metal systems can be very large. Even though static used to refer to the position-dependent difference between the up correlation can often be accounted for to a zero-order anddownspindensities,ortothevector analog of this quantity. approximation by a small number of configurations, it is often To avoid confusion, in this article we will refer to the difference very difficult to include correlation effects in a well balanced density as the spin polarization density. way in WFT calculations on multi-reference systems.11,12 Density functional theory (DFT) has now become the DFT, however, remains simple for such systems and is often preferred method for electronic structure theory for complex surprisingly accurate. This empirical fact adds to the advantage of computational efficiency in making DFT a preferred method for transition metal chemistry. Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis, MN 55455-0431, USA. The subject of transition-metal DFT is too large for any E-mail: [email protected], [email protected] single review to be complete. We selected recent papers that This journal is c the Owner Societies 2009 Phys.Chem.Chem.Phys., 2009, 11, 10757–10816 | 10757 best illustrate the promise of DFT in a number of very active includes only dynamic correlation, and DFT exchange areas of transition-metal chemistry, and we also include some includes not only exchange but also some static correlation, older references that help to put them in perspective. We although the latter is present in an unspecified and largely exclude biological applications, which deserve and uncontrolled way.18,19 In subsection 2.2 we introduce receive13–17 their own reviews. some exchange and correlation density functionals, and in Section 2 reviews the theory, section 3 contains some subsections 2.4 and 2.5 we describe a greater number of them. comments on methodology, and section 4 reviews validation Transition metal chemistry often involves open-shell systems studies. Section 5 reviews recent applications and section 6 and excited states. Kohn–Sham theory is not general enough to contains concluding remarks. treat all open-shell systems or excited states.20–22 Although generalized Kohn–Sham theories have been advanced to 2. Overview of DFT and functionals overcome this limitation, they are not in widespread use. Thus one encounters difficulties not just due to inaccurate density 2.1 Fundamental background functionals but due to Kohn–Sham theory itself. There is a rich theoretical literature on these questions, but it contains more Kohn–Sham spin-orbitals, c where s is a or b and j denotes js than one point of view, and much practical work in DFT the other quantum numbers, are obtained by a self-consistent involves using approximate functionals in a rough-and-ready field (SCF) calculation and are formally functions of the exact way (see, e.g., section 5.1). Despite such issues, Kohn–Sham density of the system. Then theory, even with approximate functionals, is often the most Xocc accurate available approach for practical work on a given r ¼ jc j2 ð1Þ s js system, especially for complex systems. j One type of generalized Kohn–Sham theory that is in where s is the spin component (a or b), the spin-orbitals are widespread use is often called hybrid DFT, and it involves normalized, and the sum is over occupied orbitals of a given combining Hartree–Fock exchange, which is orbital-dependent, spin component. The electronic energy of the system is with explicit functions of local spin densities and their gradients. approximated as a sum of four terms, Tn, ene, eee, and exc. (Note that one can also treat Hartree–Fock exchange within Tn is the kinetic energy of a system of noninteracting electrons the ungeneralized Kohn–Sham framework by using the 23,24 with the same spin densities as the real system; ene is the ‘‘optimized effective potential’’ (OEP) method, but we interaction of the electron distribution with the nuclear frame- shall not consider this formalism in detail within the present work; eee is the classical Coulomb energy of the spin densities review.) Further generalizing hybrid functionals to include interacting with each other and with themselves; and exc, called dependences on the local Laplacians of the spin densities or the exchange–correlation energy, is everything else (everything on the local spin kinetic energies computed from the spin- except Tn, ene, and eee). Therefore exc includes the interaction orbitals yields what are called hybrid meta functionals, and correction to Tn, the correction to Vee for the fact that real these are the most powerful functionals now available. electrons do not interact with themselves, the exchange energy Although one should distinguish between Kohn–Sham and (due to the indistinguishability of electrons exchanging their generalized Kohn–Sham methods, and making fundamental space and spin variables), and the correlation energy (due to progress at extending DFT to arbitrarily complex systems the fact that the many-electron spin densities are not certainly will require careful attention to the fundamentals 25 26 uncorrelated products of spin-orbital densities). exc is written of the theory, especially for magnetic properties, most as a functional, called the spin-density functional, of the spin practical algorithms in popular use ignore such issues and densities. Since the Kohn–Sham spin-orbitals are functions of the address the theory from a computational perspective with spin-densities, exc can depend explicitly on the spin densities and approximate methods and procedures being accepted or also implicitly on them by depending on the spin-orbitals.
Recommended publications
  • Benchmark Study of the Performance of Density
    DOI: 10.1002/open.201300012 Benchmark Study of the Performance of Density Functional Theory for Bond Activations with (Ni,Pd)-Based Transition-Metal Catalysts Marc Steinmetz and Stefan Grimme*[a] The performance of 23 density functionals, including one LDA, due to partial breakdown of the perturbative treatment for the four GGAs, three meta-GGAs, three hybrid GGAs, eight hybrid correlation energy in cases with difficult electronic structures meta-GGAs, and ten double-hybrid functionals, was investigat- (partial multi-reference character). Only double hybrids either ed for the computation of activation energies of various cova- with very low amounts of perturbative correlation (e.g., PBE0- lent main-group single bonds by four catalysts: Pd, PdClÀ , DH) or that use the opposite-spin correlation component only PdCl2, and Ni (all in the singlet state). A reactant complex, the (e.g., PWPB95) seem to be more robust. We also investigated barrier, and reaction energy were considered, leading to 164 the effect of the D3 dispersion correction. While the barriers energy data points for statistical analysis. Extended Gaussian are not affected by this correction, significant and mostly posi- AO basis sets were used in all calculations. The best functional tive results were observed for reaction energies. Furthermore, for the complete benchmark set relative to estimated CCSD(T)/ six very recently proposed double-hybrid functionals were ana- CBS reference data is PBE0-D3, with an MAD value of 1.1 kcal lyzed regarding the influence of the amount of Fock exchange molÀ1 followed by PW6B95-D3, the double hybrid PWPB95-D3, as well as the type of perturbative correlation treatment.
    [Show full text]
  • An Ab Initio Materials Simulation Code
    CP2K: An ab initio materials simulation code Lianheng Tong Physics on Boat Tutorial , Helsinki, Finland 2015-06-08 Faculty of Natural & Mathematical Sciences Department of Physics Brief Overview • Overview of CP2K - General information about the package - QUICKSTEP: DFT engine • Practical example of using CP2K to generate simulated STM images - Terminal states in AGNR segments with 1H and 2H termination groups What is CP2K? Swiss army knife of molecular simulation www.cp2k.org • Geometry and cell optimisation Energy and • Molecular dynamics (NVE, Force Engine NVT, NPT, Langevin) • STM Images • Sampling energy surfaces (metadynamics) • DFT (LDA, GGA, vdW, • Finding transition states Hybrid) (Nudged Elastic Band) • Quantum Chemistry (MP2, • Path integral molecular RPA) dynamics • Semi-Empirical (DFTB) • Monte Carlo • Classical Force Fields (FIST) • And many more… • Combinations (QM/MM) Development • Freely available, open source, GNU Public License • www.cp2k.org • FORTRAN 95, > 1,000,000 lines of code, very active development (daily commits) • Currently being developed and maintained by community of developers: • Switzerland: Paul Scherrer Institute Switzerland (PSI), Swiss Federal Institute of Technology in Zurich (ETHZ), Universität Zürich (UZH) • USA: IBM Research, Lawrence Livermore National Laboratory (LLNL), Pacific Northwest National Laboratory (PNL) • UK: Edinburgh Parallel Computing Centre (EPCC), King’s College London (KCL), University College London (UCL) • Germany: Ruhr-University Bochum • Others: We welcome contributions from
    [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]
  • Exploring the Limit of Accuracy of the Global Hybrid Meta Density Functional for Main-Group Thermochemistry, Kinetics, and Noncovalent Interactions
    J. Chem. Theory Comput. XXXX, xxx, 000 A Exploring the Limit of Accuracy of the Global Hybrid Meta Density Functional for Main-Group Thermochemistry, Kinetics, and Noncovalent Interactions Yan Zhao and Donald G. Truhlar* Department of Chemistry and Supercomputing Institute, UniVersity of Minnesota, 207 Pleasant Street S.E., Minneapolis, Minnesota 55455-0431 Received June 25, 2008 Abstract: The hybrid meta density functionals M05-2X and M06-2X have been shown to provide broad accuracy for main group chemistry. In the present article we make the functional form more flexible and improve the self-interaction term in the correlation functional to improve its self-consistent-field convergence. We also explore the constraint of enforcing the exact forms of the exchange and correlation functionals through second order (SO) in the reduced density gradient. This yields two new functionals called M08-HX and M08-SO, with different exact constraints. The new functionals are optimized against 267 diverse main-group energetic data consisting of atomization energies, ionization potentials, electron affinities, proton affinities, dissociation energies, isomerization energies, barrier heights, noncovalent complexation energies, and atomic energies. Then the M08-HX, M08-SO, M05-2X, and M06-2X functionals and the popular B3LYP functional are tested against 250 data that were not part of the original training data for any of the functionals, in particular 164 main-group energetic data in 7 databases, 39 bond lengths, 38 vibrational frequencies, and 9 multiplicity-changing electronic transition energies. These tests include a variety of new challenges for complex systems, including large-molecule atomization energies, organic isomerization energies, interaction energies in uracil trimers, and bond distances in crowded molecules (in particular, cyclophanes).
    [Show full text]
  • Performance of DFT Procedures for the Calculation of Proton-Exchange
    Article pubs.acs.org/JCTC Performance of Density Functional Theory Procedures for the Calculation of Proton-Exchange Barriers: Unusual Behavior of M06- Type Functionals Bun Chan,*,†,‡ Andrew T. B. Gilbert,¶ Peter M. W. Gill,*,¶ and Leo Radom*,†,‡ † School of Chemistry, University of Sydney, Sydney, NSW 2006, Australia ‡ ¶ ARC Centre of Excellence for Free Radical Chemistry and Biotechnology and Research School of Chemistry, Australian National University, Canberra, ACT 0200, Australia *S Supporting Information ABSTRACT: We have examined the performance of a variety of density functional theory procedures for the calculation of complexation energies and proton-exchange barriers, with a focus on the Minnesota-class of functionals that are generally highly robust and generally show good accuracy. A curious observation is that M05-type and M06-type methods show an atypical decrease in calculated barriers with increasing proportion of Hartree−Fock exchange. To obtain a clearer picture of the performance of the underlying components of M05-type and M06-type functionals, we have investigated the combination of MPW-type and PBE-type exchange and B95-type and PBE-type correlation procedures. We find that, for the extensive E3 test set, the general performance of the various hybrid-DFT procedures improves in the following order: PBE1-B95 → PBE1-PBE → MPW1-PBE → PW6-B95. As M05-type and M06-type procedures are related to PBE1-B95, it would be of interest to formulate and examine the general performance of an alternative Minnesota DFT method related to PW6-B95. 1. INTRODUCTION protocol.13,14 In a recent study,15 however, it was shown that, Density functional theory (DFT) procedures1 are undeniably for complexation energies of water, ammonia and hydrogen the workhorse for many computational chemists.
    [Show full text]
  • Kepler Gpus and NVIDIA's Life and Material Science
    LIFE AND MATERIAL SCIENCES Mark Berger; [email protected] Founded 1993 Invented GPU 1999 – Computer Graphics Visual Computing, Supercomputing, Cloud & Mobile Computing NVIDIA - Core Technologies and Brands GPU Mobile Cloud ® ® GeForce Tegra GRID Quadro® , Tesla® Accelerated Computing Multi-core plus Many-cores GPU Accelerator CPU Optimized for Many Optimized for Parallel Tasks Serial Tasks 3-10X+ Comp Thruput 7X Memory Bandwidth 5x Energy Efficiency How GPU Acceleration Works Application Code Compute-Intensive Functions Rest of Sequential 5% of Code CPU Code GPU CPU + GPUs : Two Year Heart Beat 32 Volta Stacked DRAM 16 Maxwell Unified Virtual Memory 8 Kepler Dynamic Parallelism 4 Fermi 2 FP64 DP GFLOPS GFLOPS per DP Watt 1 Tesla 0.5 CUDA 2008 2010 2012 2014 Kepler Features Make GPU Coding Easier Hyper-Q Dynamic Parallelism Speedup Legacy MPI Apps Less Back-Forth, Simpler Code FERMI 1 Work Queue CPU Fermi GPU CPU Kepler GPU KEPLER 32 Concurrent Work Queues Developer Momentum Continues to Grow 100M 430M CUDA –Capable GPUs CUDA-Capable GPUs 150K 1.6M CUDA Downloads CUDA Downloads 1 50 Supercomputer Supercomputers 60 640 University Courses University Courses 4,000 37,000 Academic Papers Academic Papers 2008 2013 Explosive Growth of GPU Accelerated Apps # of Apps Top Scientific Apps 200 61% Increase Molecular AMBER LAMMPS CHARMM NAMD Dynamics GROMACS DL_POLY 150 Quantum QMCPACK Gaussian 40% Increase Quantum Espresso NWChem Chemistry GAMESS-US VASP CAM-SE 100 Climate & COSMO NIM GEOS-5 Weather WRF Chroma GTS 50 Physics Denovo ENZO GTC MILC ANSYS Mechanical ANSYS Fluent 0 CAE MSC Nastran OpenFOAM 2010 2011 2012 SIMULIA Abaqus LS-DYNA Accelerated, In Development NVIDIA GPU Life Science Focus Molecular Dynamics: All codes are available AMBER, CHARMM, DESMOND, DL_POLY, GROMACS, LAMMPS, NAMD Great multi-GPU performance GPU codes: ACEMD, HOOMD-Blue Focus: scaling to large numbers of GPUs Quantum Chemistry: key codes ported or optimizing Active GPU acceleration projects: VASP, NWChem, Gaussian, GAMESS, ABINIT, Quantum Espresso, BigDFT, CP2K, GPAW, etc.
    [Show full text]
  • Dr. David Danovich
    Dr. David Danovich ! 14 January 2021 (h-index: 34; 124 publications) Personal Information Name: Contact information Danovich David Private: Nissan Harpaz st., 4/9, Jerusalem, Date of Birth: 9371643 Israel June 29, 1959 Mobile: +972-54-4768669 E-mail: [email protected] Place of Birth: Irkutsk, USSR Work: Institute of Chemistry, The Hebrew Citizenship: University, Edmond Safra Campus, Israeli Givat Ram, 9190401 Jerusalem, Israel Marital status: Phone: +972-2-6586934 Married +2 FAX: +972-2-6584033 E-mail: [email protected] Current position 1992- present: Senior computational chemist and scientific programmer in the Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem, Israel Scientific skills Very experienced scientific researcher with demonstrated ability to implement ab initio molecular electronic structure and valence bond calculations for obtaining insights into physical and chemical properties of molecules and materials. Extensive expertise in implementation of Molecular Orbital Theory based on ab initio methods such as Hartree-Fock theory (HF), Density Functional Theory (DFT), Moller-Plesset Perturbation theory (MP2-MP5), Couple Cluster theories [CCST(T)], Complete Active Space (CAS, CASPT2) theories, Multireference Configuration Interaction (MRCI) theory and modern ab initio valence bond theory (VBT) methods such as Valence Bond Self Consistent Field theory (VBSCF), Breathing Orbital Valence Bond theory (BOVB) and its variants (like D-BOVB, S-BOVB and SD-BOVB), Valence Bond Configuration Interaction theory (VBCI), Valence Bond Perturbation theory (VBPT). Language skills Advanced Conversational English, Russian Hebrew Academic background 1. Undergraduate studies, 1978-1982 Master of Science in Physics, June 1982 Irkutsk State University, USSR Dr. David Danovich ! 2. Graduate studies, 1984-1989 Ph.
    [Show full text]
  • An Ion Mobility-Aided, Quantum Chemical Benchmark Analysis of H+GPGG Conformers Daniel Beckett, Tarick J
    Article Cite This: J. Chem. Theory Comput. 2018, 14, 5406−5418 pubs.acs.org/JCTC Electronic Energies Are Not Enough: An Ion Mobility-Aided, Quantum Chemical Benchmark Analysis of H+GPGG Conformers Daniel Beckett, Tarick J. El-Baba, David E. Clemmer, and Krishnan Raghavachari* Department of Chemistry, Indiana University, Bloomington, Indiana 47405, United States *S Supporting Information ABSTRACT: Peptides and proteins exist not as single, lowest energy structures but as ensembles of states separated by small barriers. In order to study these species we must be able to correctly identify their gas phase conformational distributions, and ion mobility spectrometry (IMS) has arisen as an experimental method for assessing the gas phase energetics of flexible peptides. Here, we present a thorough exploration and benchmarking of the low energy conformers of the small, hairpin peptide H+GPGG with the aid of ion mobility ). spectrometry against a wide swath of density functionals (35 dispersion-corrected and uncorrected functionals represented by 21 unique exchange−correlation functionals) and wave function theory methods (15 total levels of theory). The three experimentally resolved IMS peaks were found to correspond to three distinct pairs of conformers, each pair composed of species differing only by the chair or boat configuration of the proline. Two of the H+GPGG conformer pairs possess a cis configuration about the Pro−Gly1 peptide bond while the other adopts a trans configuration. While the experimental spectrum reports a higher intensity for the o legitimately share published articles. cis-1 conformer than the trans conformer, 13 WFT levels of theory, including a complete basis set CCSD(T) extrapolation, obtain trans to be favored in terms of the electronic energy.
    [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]
  • Jaguar 5.5 User Manual Copyright © 2003 Schrödinger, L.L.C
    Jaguar 5.5 User Manual Copyright © 2003 Schrödinger, L.L.C. All rights reserved. Schrödinger, FirstDiscovery, Glide, Impact, Jaguar, Liaison, LigPrep, Maestro, Prime, QSite, and QikProp are trademarks of Schrödinger, L.L.C. MacroModel is a registered trademark of Schrödinger, L.L.C. To the maximum extent permitted by applicable law, this publication is provided “as is” without warranty of any kind. This publication may contain trademarks of other companies. October 2003 Contents Chapter 1: Introduction.......................................................................................1 1.1 Conventions Used in This Manual.......................................................................2 1.2 Citing Jaguar in Publications ...............................................................................3 Chapter 2: The Maestro Graphical User Interface...........................................5 2.1 Starting Maestro...................................................................................................5 2.2 The Maestro Main Window .................................................................................7 2.3 Maestro Projects ..................................................................................................7 2.4 Building a Structure.............................................................................................9 2.5 Atom Selection ..................................................................................................10 2.6 Toolbar Controls ................................................................................................11
    [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]
  • Trends in Atomistic Simulation Software Usage [1.3]
    A LiveCoMS Perpetual Review Trends in atomistic simulation software usage [1.3] Leopold Talirz1,2,3*, Luca M. Ghiringhelli4, Berend Smit1,3 1Laboratory of Molecular Simulation (LSMO), Institut des Sciences et Ingenierie Chimiques, Valais, École Polytechnique Fédérale de Lausanne, CH-1951 Sion, Switzerland; 2Theory and Simulation of Materials (THEOS), Faculté des Sciences et Techniques de l’Ingénieur, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 3National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland; 4The NOMAD Laboratory at the Fritz Haber Institute of the Max Planck Society and Humboldt University, Berlin, Germany This LiveCoMS document is Abstract Driven by the unprecedented computational power available to scientific research, the maintained online on GitHub at https: use of computers in solid-state physics, chemistry and materials science has been on a continuous //github.com/ltalirz/ rise. This review focuses on the software used for the simulation of matter at the atomic scale. We livecoms-atomistic-software; provide a comprehensive overview of major codes in the field, and analyze how citations to these to provide feedback, suggestions, or help codes in the academic literature have evolved since 2010. An interactive version of the underlying improve it, please visit the data set is available at https://atomistic.software. GitHub repository and participate via the issue tracker. This version dated August *For correspondence: 30, 2021 [email protected] (LT) 1 Introduction Gaussian [2], were already released in the 1970s, followed Scientists today have unprecedented access to computa- by force-field codes, such as GROMOS [3], and periodic tional power.
    [Show full text]