DALTON2011 Program Manual

Total Page:16

File Type:pdf, Size:1020Kb

DALTON2011 Program Manual DALTON2011 Program Manual C. Angeli, K. L. Bak, V. Bakken, O. Christiansen, R. Cimiraglia, S. Coriani, P. Dahle, E. K. Dalskov, T. Enevoldsen, B. Fernandez, L. Ferrighi, L. Frediani, C. H¨attig, K. Hald, A. Halkier, H. Heiberg, T. Helgaker, H. Hettema, B. Jansik, H. J. Aa. Jensen, D. Jonsson, P. Jørgensen, S. Kirpekar, W. Klopper, S. Knecht, R. Kobayashi, J. Kongsted, H. Koch, A. Ligabue, O. B. Lutnæs, K. V. Mikkelsen, C. B. Nielsen, P. Norman, J. Olsen, A. Osted, M. J. Packer, T. B. Pedersen, Z. Rinkevicius, E. Rudberg, T. A. Ruden, K. Ruud, P. Sa lek,C. C. M. Samson, A. Sanchez de Meras, T. Saue, S. P. A. Sauer, B. Schimmelpfennig, A. H. Steindal, K. O. Sylvester-Hvid, P. R. Taylor, O. Vahtras, D. J. Wilson, and H. Agren,˚ Contents Preface ix 1 Introduction 1 1.1 General description of the manual . 2 1.2 Acknowledgments . 3 2 New features in the Dalton releases 4 2.1 New features in Dalton2011 . 4 2.2 New features in Dalton 2.0 (2005) . 7 2.3 New features in Dalton 1.2 . 10 I DALTON Installation Guide 13 3 Installation 14 3.1 Hardware/software supported . 14 3.2 Source files . 14 3.3 Installing the program using the Makefile . 15 3.4 Running the dalton2011 test suite . 18 4 Maintenance 20 4.1 Memory requirements . 20 4.1.1 Redimensioning dalton2011 ...................... 20 4.2 New versions, patches . 21 4.3 Reporting bugs and user support . 22 II DALTON User’s Guide 23 5 Getting started with dalton2011 24 5.1 The DALTON.INP file . 24 i CONTENTS ii 5.1.1 A CASSCF geometry optimization . 24 5.1.2 A RASSCF calculation of NMR parameters . 25 5.1.3 A parallel cubic response calculation . 26 5.2 General structure of the DALTON.INP file . 27 5.3 The molecule input file . 29 5.4 The first calculation with dalton2011 ..................... 31 6 Getting the wave function you want 35 6.1 Necessary input to SIRIUS . 36 6.2 An input example for SIRIUS . 36 6.3 Hints on the structure of the **WAVE FUNCTIONS input . 39 6.4 How to restart a wave function calculation . 41 6.5 Transfer of molecular orbitals between different computers . 42 6.6 Wave function input examples . 42 7 Potential energy surfaces 51 7.1 Locating stationary points . 52 7.1.1 Equilibrium geometries . 52 7.1.2 Transition states using the image method . 57 7.1.3 Transition states using first-order methods . 59 7.1.4 Transition states following a gradient extremal . 60 7.1.5 Level-shifted mode-following . 62 7.2 Trajectories and Dynamics . 63 7.2.1 Intrinsic reaction coordinates . 63 7.2.2 Doing a dynamical walk . 64 7.2.3 Calculating relative translational energy release . 67 7.3 Geometry optimization using non-variational wave functions . 67 8 Molecular vibrations 69 8.1 Vibrational frequencies . 69 8.2 Infrared (IR) intensities . 70 8.3 Dipole-gradient based population analysis . 71 8.4 Raman intensities . 72 8.5 Vibrational g factor . 74 9 Electric properties 77 9.1 Dipole moment . 77 9.2 Quadrupole moment . 77 9.3 Nuclear quadrupole coupling constants . 78 CONTENTS iii 9.4 Static and frequency dependent polarizabilities . 79 10 Calculation of magnetic properties 81 10.1 Magnetizabilities . 82 10.2 Nuclear shielding constants . 84 10.3 Rotational g tensor . 85 10.4 Nuclear spin–rotation constants . 86 10.5 Indirect nuclear spin–spin coupling constants . 87 10.6 Hyperfine Coupling Tensors . 89 10.7 Electronic g-tensors . 91 10.8 Zero field splitting . 91 10.9 CTOCD-DZ calculations . 92 10.9.1 General considerations . 92 10.9.2 Input description . 93 11 Calculation of optical and Raman properties 96 11.1 Electronic excitation energies and oscillator strengths . 96 11.2 Vibrational Circular Dichroism calculations . 98 11.3 Electronic circular dichroism (ECD) . 100 11.4 Optical Rotation . 103 11.5 Vibrational Raman Optical Activity (VROA) . 105 12 Getting the property you want 109 12.1 General considerations . 109 12.2 Input description . 110 12.2.1 Linear response . 110 12.2.2 Quadratic response . 113 12.2.3 Cubic response . 115 13 Direct and parallel calculations 117 13.1 Direct methods . 117 13.2 Parallel methods . 118 14 Finite field calculations 119 14.1 General considerations . 119 14.2 Input description . 120 15 Solvent calculations 122 15.1 General considerations . 122 15.2 Input description . 123 CONTENTS iv 15.2.1 Geometry optimization . 125 15.2.2 Non-equilibrium solvation . 125 16 Vibrational corrections 128 16.1 Effective geometries . 128 16.2 Vibrational averaged properties . 130 16.3 Vibrationally averaged spin–spin coupling constants . 132 17 Relativistic Effects 134 18 SOPPA, SOPPA(CC2), SOPPA(CCSD) and RPA(D) 136 18.1 General considerations . 136 18.2 Input description molecular orbital based SOPPA . 138 18.3 Input description atomic orbital based SOPPA module . 141 19 NEVPT2 calculations 145 19.1 General considerations . 145 19.2 Input description . 146 20 Examples of generalized active space CI calculations 147 20.1 Energy calculation with a GAS-type active space decomposition I . 147 20.2 Energy calculation with a GAS-type active space decomposition II . 149 20.3 Energy calculation with a RAS-type active space decomposition . 150 21 Examples of coupled cluster calculations 152 21.1 Multiple model energy calculations . 152 21.2 First-order property calculation . 153 21.3 Static and frequency-dependent dipole polarizabilities and corresponding dis- persion coefficients . 153 21.4 Static and frequency-dependent dipole hyperpolarizabilities and correspond- ing dispersion coefficients . 154 21.5 Excitation energies and oscillator strengths . 155 21.6 Gradient calculation, geometry optimization . 156 21.7 R12 methods . 157 22 Examples of Cholesky decomposition-based calculations 158 22.1 Hartree-Fock energy and polarizability . 158 22.2 KT3 magnetic properties using London orbitals . 159 22.3 MP2 energy . 159 22.4 Restart of MP2 energy . 160 CONTENTS v 22.5 CC2 magnetic properties using the CTOCD/DZ method . 162 22.6 CCSD(T) energy calculation using decomposed energy denominators . 163 III DALTON Reference Manual 164 23 General input module 165 23.1 General input to DALTON : **DALTON ..................... 165 23.1.1 Geometry optimization module 1: *OPTIMIZE ............. 168 23.1.2 Parallel calculations : *PARALLEL .................... 178 23.1.3 PCM environment model: *PCM ..................... 178 23.1.4 QM/MM environment model: *QM3 .................. 178 23.1.5 Geometry optimization module 2: *WALK ................ 179 23.1.6 Molecule geometry and basis sets, *MOLBAS .............. 184 23.2 Numerical differentiation : **NMDDRV ...................... 185 23.2.1 Vibrational averaging of molecular properties: *PROPAV ....... 188 23.2.2 Vibrational analysis: *VIBANA ...................... 189 23.3 Decomposition of two-electron integrals : **CHOLES .............. 190 24 Integral evaluation, hermit 192 24.1 General . 192 24.2 **INTEGRALS directives . 193 24.2.1 General: **INTEGRALS .......................... 193 24.2.2 One-electron integrals: *ONEINT ..................... 210 24.2.3 Two-electron integrals using twoint: *TWOINT ............ 211 24.2.4 Two-electron integrals using eri: *ER2INT ............... 212 24.2.5 Integral sorting: *SORINT ........................ 214 24.2.6 Construction of the supermatrix file: *SUPINT ............. 215 25 molecule input style 216 25.1 General molecule input ............................ 217 25.2 Cartesian geometry input . 220 25.3 Z-matrix input . 223 25.4 Using basis set libraries . 224 25.5 Auxiliary basis sets . 227 25.6 The basis sets supplied with dalton2011 ................... 228 26 Molecular wave functions, sirius 236 26.1 General notes for the sirius input reference manual . 236.
Recommended publications
  • Density Functional Theory
    Density Functional Approach Francesco Sottile Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) 22 October 2010 Density Functional Theory 1. Any observable of a quantum system can be obtained from the density of the system alone. < O >= O[n] Hohenberg, P. and W. Kohn, 1964, Phys. Rev. 136, B864 Density Functional Theory 1. Any observable of a quantum system can be obtained from the density of the system alone. < O >= O[n] 2. The density of an interacting-particles system can be calculated as the density of an auxiliary system of non-interacting particles. Hohenberg, P. and W. Kohn, 1964, Phys. Rev. 136, B864 Kohn, W. and L. Sham, 1965, Phys. Rev. 140, A1133 Density Functional ... Why ? Basic ideas of DFT Importance of the density Example: atom of Nitrogen (7 electron) 1. Any observable of a quantum Ψ(r1; ::; r7) 21 coordinates system can be obtained from 10 entries/coordinate ) 1021 entries the density of the system alone. 8 bytes/entry ) 8 · 1021 bytes 4:7 × 109 bytes/DVD ) 2 × 1012 DVDs 2. The density of an interacting-particles system can be calculated as the density of an auxiliary system of non-interacting particles. Density Functional ... Why ? Density Functional ... Why ? Density Functional ... Why ? Basic ideas of DFT Importance of the density Example: atom of Oxygen (8 electron) 1. Any (ground-state) observable Ψ(r1; ::; r8) 24 coordinates of a quantum system can be 24 obtained from the density of the 10 entries/coordinate ) 10 entries 8 bytes/entry ) 8 · 1024 bytes system alone. 5 · 109 bytes/DVD ) 1015 DVDs 2.
    [Show full text]
  • Electronic Supporting Information a New Class of Soluble and Stable Transition-Metal-Substituted Polyoxoniobates: [Cr2(OH)4Nb10o
    Electronic Supplementary Material (ESI) for Dalton Transactions This journal is © The Royal Society of Chemistry 2012 Electronic Supporting Information A New Class of Soluble and Stable Transition-Metal-Substituted Polyoxoniobates: [Cr2(OH)4Nb10O30]8-. Jung-Ho Son, C. André Ohlin and William H. Casey* A) Computational Results The two most interesting aspects of this new polyoxoanion is the inclusion of a more redox active metal and the optical activity in the visible region. For these reasons we in particular wanted to explore whether it was possible to readily and routinely predict the location of the absorbance bands of this type of compound. Also, owing to the difficulty in isolating the oxidised and reduced forms of the central polyanion dichromate species, we employed computations at the density functional theory (DFT) level to get a better idea of the potential characteristics of these forms. Structures were optimised for anti-parallel and parallel electronic spin alignments, and for both high and low spin cases where relevant. Intermediate spin configurations were 8- II III also tested for [Cr2(OH)4Nb10O30] . As expected, the high spin form of the Cr Cr complex was favoured (Table S1, below). DFT overestimates the bond lengths in general, but gives a good picture of the general distortion of the octahedral chromium unit, and in agreement with the crystal structure predicts that the shortest bond is the Cr-µ2-O bond trans from the central µ6- oxygen in the Lindqvist unit, and that the (Nb-)µ2-O-Cr-µ4-O(-Nb,Nb,Cr) angle is smaller by ca 5˚ from the ideal 180˚.
    [Show full text]
  • Computer-Assisted Catalyst Development Via Automated Modelling of Conformationally Complex Molecules
    www.nature.com/scientificreports OPEN Computer‑assisted catalyst development via automated modelling of conformationally complex molecules: application to diphosphinoamine ligands Sibo Lin1*, Jenna C. Fromer2, Yagnaseni Ghosh1, Brian Hanna1, Mohamed Elanany3 & Wei Xu4 Simulation of conformationally complicated molecules requires multiple levels of theory to obtain accurate thermodynamics, requiring signifcant researcher time to implement. We automate this workfow using all open‑source code (XTBDFT) and apply it toward a practical challenge: diphosphinoamine (PNP) ligands used for ethylene tetramerization catalysis may isomerize (with deleterious efects) to iminobisphosphines (PPNs), and a computational method to evaluate PNP ligand candidates would save signifcant experimental efort. We use XTBDFT to calculate the thermodynamic stability of a wide range of conformationally complex PNP ligands against isomeriation to PPN (ΔGPPN), and establish a strong correlation between ΔGPPN and catalyst performance. Finally, we apply our method to screen novel PNP candidates, saving signifcant time by ruling out candidates with non‑trivial synthetic routes and poor expected catalytic performance. Quantum mechanical methods with high energy accuracy, such as density functional theory (DFT), can opti- mize molecular input structures to a nearby local minimum, but calculating accurate reaction thermodynamics requires fnding global minimum energy structures1,2. For simple molecules, expert intuition can identify a few minima to focus study on, but an alternative approach must be considered for more complex molecules or to eventually fulfl the dream of autonomous catalyst design 3,4: the potential energy surface must be frst surveyed with a computationally efcient method; then minima from this survey must be refned using slower, more accurate methods; fnally, for molecules possessing low-frequency vibrational modes, those modes need to be treated appropriately to obtain accurate thermodynamic energies 5–7.
    [Show full text]
  • Camcasp 5.9 Alston J. Misquitta† and Anthony J
    CamCASP 5.9 Alston J. Misquittay and Anthony J. Stoneyy yDepartment of Physics and Astronomy, Queen Mary, University of London, 327 Mile End Road, London E1 4NS yy University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW September 1, 2015 Abstract CamCASP is a suite of programs designed to calculate molecular properties (multipoles and frequency- dependent polarizabilities) in single-site and distributed form, and interaction energies between pairs of molecules, and thence to construct atom–atom potentials. The CamCASP distribution also includes the programs Pfit, Casimir,Gdma 2.2, Cluster, and Process. Copyright c 2007–2014 Alston J. Misquitta and Anthony J. Stone Contents 1 Introduction 1 1.1 Authors . .1 1.2 Citations . .1 2 What’s new? 2 3 Outline of the capabilities of CamCASP and other programs 5 3.1 CamCASP limits . .7 4 Installation 7 4.1 Building CamCASP from source . .9 5 Using CamCASP 10 5.1 Workflows . 10 5.2 High-level scripts . 10 5.3 The runcamcasp.py script . 11 5.4 Low-level scripts . 13 6 Data conventions 13 7 CLUSTER: Detailed specification 14 7.1 Prologue . 15 7.2 Molecule definitions . 15 7.3 Geometry manipulations and other transformations . 16 7.4 Job specification . 18 7.5 Energy . 23 7.6 Crystal . 25 7.7 ORIENT ............................................... 26 7.8 Finally, . 29 8 Examples 29 8.1 A SAPT(DFT) calculation . 29 8.2 An example properties calculation . 30 8.3 Dispersion coefficients . 34 8.4 Using CLUSTER to obtain the dimer geometry . 36 9 CamCASP program specification 39 9.1 Global data .
    [Show full text]
  • Dalton2018.0 – Lsdalton Program Manual
    Dalton2018.0 – lsDalton Program Manual V. Bakken, A. Barnes, R. Bast, P. Baudin, S. Coriani, R. Di Remigio, K. Dundas, P. Ettenhuber, J. J. Eriksen, L. Frediani, D. Friese, B. Gao, T. Helgaker, S. Høst, I.-M. Høyvik, R. Izsák, B. Jansík, F. Jensen, P. Jørgensen, J. Kauczor, T. Kjærgaard, A. Krapp, K. Kristensen, C. Kumar, P. Merlot, C. Nygaard, J. Olsen, E. Rebolini, S. Reine, M. Ringholm, K. Ruud, V. Rybkin, P. Sałek, A. Teale, E. Tellgren, A. Thorvaldsen, L. Thøgersen, M. Watson, L. Wirz, and M. Ziolkowski Contents Preface iv I lsDalton Installation Guide 1 1 Installation of lsDalton 2 1.1 Installation instructions ............................. 2 1.2 Hardware/software supported .......................... 2 1.3 Source files .................................... 2 1.4 Installing the program using the Makefile ................... 3 1.5 The Makefile.config ................................ 5 1.5.1 Intel compiler ............................... 5 1.5.2 Gfortran/gcc compiler .......................... 9 1.5.3 Using 64-bit integers ........................... 12 1.5.4 Using Compressed-Sparse Row matrices ................ 13 II lsDalton User’s Guide 14 2 Getting started with lsDalton 15 2.1 The Molecule input file ............................ 16 2.1.1 BASIS ................................... 17 2.1.2 ATOMBASIS ............................... 18 2.1.3 User defined basis sets .......................... 19 2.2 The LSDALTON.INP file ............................ 21 3 Interfacing to lsDalton 24 3.1 The Grand Canonical Basis ........................... 24 3.2 Basisset and Ordering .............................. 24 3.3 Normalization ................................... 25 i CONTENTS ii 3.4 Matrix File Format ................................ 26 3.5 The lsDalton library bundle .......................... 27 III lsDalton Reference Manual 29 4 List of lsDalton keywords 30 4.1 **GENERAL ................................... 31 4.1.1 *TENSOR ...............................
    [Show full text]
  • Quantum Chemistry on a Grid
    Quantum Chemistry on a Grid Lucas Visscher Vrije Universiteit Amsterdam 1 Outline Why Quantum Chemistry needs to be gridified Computational Demands Other requirements What may be expected ? Program descriptions Dirac • Organization • Distributed development and testing Dalton Amsterdam Density Functional (ADF) An ADF-Dalton-Dirac grid driver Design criteria Current status Quantum Chemistry on the Grid, 08-11-2005 L. Visscher, Vrije Universiteit Amsterdam 2 Computational (Quantum) Chemistry QC is traditionally run on supercomputers Large fraction of the total supercomputer time annually allocated Characteristics Compute Time scales cubically or higher with number of atoms treated. Typical : a few hours to several days on a supercomputer Memory Usage scales quadratically or higher with number of atoms. Typical : A few hundred MB till tens of GB Data Storage variable. Typical: tens of Mb till hundreds of GB. Communication • Often serial (one-processor) runs in production work • Parallelization needs sufficient band width. Typical lower limit 100 Mb/s. Quantum Chemistry on the Grid, 08-11-2005 L. Visscher, Vrije Universiteit Amsterdam 3 Changing paradigms in quantum chemistry Computable molecules are often already “large” enough Many degrees of freedom Search for global minimum energy is less important Statistical sampling of different conformations is desired Let the nuclei move, follow molecules that evolve in time 2 steps : (Quantum) dynamics on PES 1 step : Classical dynamics with forces computed via QM Changing characteristics General: more similar calculations Compute time : Increases : more calculations Memory usage : Stays constant : typical size remains the same Data storage : Decreases : intermediate results are less important Communication: Decreases : only start-up data and final result Quantum Chemistry on the Grid, 08-11-2005 L.
    [Show full text]
  • Computational Studies of the Unusual Water Adduct [Cp2time (OH2)]: the Roles of the Solvent and the Counterion
    Dalton Transactions View Article Online PAPER View Journal | View Issue Computational studies of the unusual water + adduct [Cp2TiMe(OH2)] : the roles of the Cite this: Dalton Trans., 2014, 43, 11195 solvent and the counterion† Jörg Saßmannshausen + The recently reported cationic titanocene complex [Cp2TiMe(OH2)] was subjected to detailed computational studies using density functional theory (DFT). The calculated NMR spectra revealed the importance of including the anion and the solvent (CD2Cl2) in order to calculate spectra which were in good agreement Received 29th January 2014, with the experimental data. Specifically, two organic solvent molecules were required to coordinate to Accepted 19th March 2014 the two hydrogens of the bound OH2 in order to achieve such agreement. Further elaboration of the role DOI: 10.1039/c4dt00310a of the solvent led to Bader’s QTAIM and natural bond order calculations. The zirconocene complex + www.rsc.org/dalton [Cp2ZrMe(OH2)] was simulated for comparison. Creative Commons Attribution-NonCommercial 3.0 Unported Licence. Introduction Group 4 metallocenes have been the workhorse for a number of reactions for some decades now. In particular the cationic + compounds [Cp2MR] (M = Ti, Zr, Hf; R = Me, CH2Ph) are believed to be the active catalysts for a number of polymeriz- – ation reactions such as Kaminsky type α-olefin,1 10 11–15 16–18 This article is licensed under a carbocationic or ring-opening lactide polymerization. For these highly electrophilic cationic compounds, water is usually considered a poison as it leads to catalyst decompo- sition. In this respect it was quite surprising that Baird Open Access Article. Published on 21 March 2014.
    [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]
  • Thomas–Fermi–Dirac–Von Weizsäcker Models in Finite Systems Garnet Kin-Lic Chan, Aron J
    View metadata, citation and similar papers at core.ac.uk brought to you by CORE provided by Caltech Authors - Main Thomas–Fermi–Dirac–von Weizsäcker models in finite systems Garnet Kin-Lic Chan, Aron J. Cohen, and Nicholas C. Handy Citation: 114, (2001); doi: 10.1063/1.1321308 View online: http://dx.doi.org/10.1063/1.1321308 View Table of Contents: http://aip.scitation.org/toc/jcp/114/2 Published by the American Institute of Physics JOURNAL OF CHEMICAL PHYSICS VOLUME 114, NUMBER 2 8 JANUARY 2001 Thomas–Fermi–Dirac–von Weizsa¨cker models in finite systems Garnet Kin-Lic Chan,a) Aron J. Cohen, and Nicholas C. Handy Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, United Kingdom ͑Received 7 June 2000; accepted 8 September 2000͒ To gain an understanding of the variational behavior of kinetic energy functionals, we perform a numerical study of the Thomas–Fermi–Dirac–von Weizsa¨cker theory in finite systems. A general purpose Gaussian-based code is constructed to perform energy and geometry optimizations on polyatomic systems to high accuracy. We carry out benchmark studies on atomic and diatomic systems. Our results indicate that the Thomas–Fermi–Dirac–von Weizsa¨cker theory can give an approximate description of matter, with atomic energies, binding energies, and bond lengths of the correct order of magnitude, though not to the accuracy required of a qualitative chemical theory. We discuss the implications for the development of new kinetic functionals. © 2001 American Institute of Physics. ͓DOI: 10.1063/1.1321308͔ I. INTRODUCTION
    [Show full text]
  • 1 Installation Guide for Columbus Version 7.1 with Openmolcas Support (Colmol-2.1.0)
    1 Installation Guide for Columbus Version 7.1 with OpenMolcas Support (colmol-2.1.0) 1.1 Overview Columbus and Molcas cooperate by exchanging well-defined, transferable quantities, specif- ically AO integrals, AO density matrices and MO coefficients in the native Molcas format by using library functions from the Molcas library. Additionally, the RunFile is processed, which is used by Molcas to store additional information often required to be passed in between dif- ferent Molcas modules. Molcas provides some mechanism to incorporate external programs and making them accessible through the standard Molcas input file. Since it is necessary to link a portion of the Molcas library, both Columbus and Molcas must be compiled and linked in a compatible way, that is to say (i) the same compilers and (ii) compatible compiler options. Hence, both codes cooperate on a binary level and do not rely on some file conversion utilities. Please note, that the Molcas library is frequently restructured, so that files produced with a previous version may or may not be compatible with the current version of Molcas (cf. page 4). Additionally, there are now two molcas driver utilities: molcas.exe in Molcas-8.2 developer version and pymolcas in the OpenMolcas version, which may not behave exactly the same way. 1.2 New Features Primarily, rewritten versions of the SCF, MCSCF, AO-MO transformation and CI-gradient code are introduced. With exception of the CI gradient code, all codes are now parallel and make extensive use of shared memory and much more efficient intermediate data storage. The same driver (columbus.exe) now also supports an old modified version of the dalton code (affecting primarily the integral and density matrix handling and storage format in the old HERMIT and ABACUS code portions).
    [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]