Ab initio molecular dynamics
|
|
- Elijah Lambert
- 5 years ago
- Views:
Transcription
1 Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field that use quantum mechanical methods to study molecules properties and reactions In most of the calculations the studied molecules are in vacuum which is seldom the case with real molecules We need computational chemistry in realistic environment The molecules also moves so often we need to simulate the molecular dynamics (MD) The main advantage of AIMD is that chemical reactions can be studied.
2 Ab Initio Molecular Dynamics combining periodic DFT-GGA and MD. the atoms are treated classically, (Born-Oppenheimer approximation) F = Ma, F = - d V/d R electrons are included and they are treated using DFT (also HF and hybrid functionals are possible) E KS (,R) forces are calculated in DFT level F = - d E KS (,R) / d R the core electrons are described with effective core potentials we need to either optimize the electrons at every time step or to use the Car-Parrinello method expensive calculations Wavefunctions do not change much between time steps the previous wf s are an excellent start for the new ones. => Fast convergence (5-10 iterations)
3 Ab Initio Molecular Dynamics CPMD smooth effective core potentials and plane wave basis set. Car- Parrinello algorithm, time step ca. 0.1 fs CP2K hard effective core potentials and gaussian basis set. Time step ca. 1 fs. Note: the wavefunctions has to be computed at every time step! accuracy of GGA is usually good Van der Waals interactions are missing. We use DFT + empirical corrections a la Grimme full arsenal of MD techniques and electronic structure analysis methods are implemented thermostats, constraints, thermodynamic integration, Wannier functions, TD-DFT
4 Computational aspects CP2K code Developed mostly in Zurich, Prof. Hutter s group Free to download (from the address above) Tutorials: 2nd CP2K Tutorial: Enabling the Power of Imagination in MD Simulations. Very complex code ( lines of code, Fortran 95) Huge amount of features Difficult to compile and difficult to learn to use Important help from CSC (compilation)
5 Computational aspects CP2K code Efficient code but in normal application not very good parallel scaling. We are interested of some simulation of 200+ waters and few MD steps. One simulation should not take more month (wall clock time, Cray XX). CP2K scaling A supercomputer is essential 100 Good scaling up to 1024 cores (10 s/step is makes few 100 ps simulations possible) Time 10 64water 128water NaCl+475w 1024water Cores
6 Ab Initio Molecular Dynamics examples: Al 2 O n H m Cl water, Al 5 O n H m Cl k water, PBE-GGA various simulations, simulation time scale ca. 100 ps CPMD or CP2K code, computations with Cray XT5/XT4
7 Ab initio molecular dynamics of acidic systems Hydrofluoric Acid HF acid is a very interesting because at low concentrations it is a weak acid (ph ca. 3) but at high concentration it become very strong acid. HF and water mix with all concentrations. Sillanpää, A., Simon, C., Klein, M.L., Laasonen, K. J. Phys. Chem. B 106, , (2002). Simon, C., Klein, M.L., ChemPhysChem, 6, , (2005). Hydrochloric Acid HCl is a typical strong acid. It s soluability limit is ca. 30 mol % We have done DCl:D 2 O simulations with concentrations of 4:28, 7:25 and 10:22 at ca. 300 K and 10:22 at 470 K and 910 K. Sillanpää, A., Laasonen, K. Phys.Chem.Chem.Phys., 6, , (2004). Heuft, J.M., Meijer, E.J., Phys.Chem.Chem.Phys., 8, , (2006).
8 Mixture of Hydrofluoric and Hydrochloric Acid Simulations: 1 HF, 3 HCl, 28 waters,15 ps, NVT, 330 K 3 HF, 4 HCl, 25 w, 12 ps + 25 ps NVT, 330 K 6 HF, 8 HCl, 18 w, 50 ps, NVT, 350 K, 20 ps, NVE, 320 K 6 HF, 8 HCl, 18 w, 40 ps annealing + 60 ps, NVT, 350 K 14 HF, 0 HCl, 18 w, 2x50 ps, NVT 350 K (in all simulations dt = fs (5 au), = 500 au) PBE GGA and Vanderbilt pseudopotentials, cut-off 35 Ry all hydrogens were replaced with deuteriums we wanted to see how HF is behaving in acid environment K. Laasonen, J. Larrucea, A. Sillanpää, J. Phys. Chem. B, 120, (2006)
9 Structure heavy atoms O-O distances are much shorter than in water (similar as in concentrated HF(aq) and HCl(aq)) O coordination is very sensitive to the acid concentration: 2.7 (1/3), 2.2 (3/4), 1.1 (6/8), 1.4 (14/0) (const. cut-off 3.0 Å) 3.2 (1/3), 2.3 (3/4), 0.9 (6/8), 1.6 (14/0) (integral to the first minimum) only 1 oxygen around water in the 6/8 simulation! (in pure water ca. 4)
10 OH pair correlation functions Good correlation with hydronium concentration (hydronium is defined by using OH cutoff of 1.25) The water molecules start loosing their meaning!! O-O coordination ca 1. H H H H Protons are getting very close to heavy atoms.
11 Free energy profiles for hydrogen. F( q) P( q) RT ln max ( P( q)) In acids the protons are very mobile, OHO barrier is around 1 kj/mol in the mixture simulation
12 Some general comments on acid AIMD simulations You do not do such simulations with empirical potentials Doable with AIMD Limitations: time scale, accuracy of GGA
13 Ab initio molecular dynamics aluminum oxide chemistry in aqueous solution Al oxides are widely used chemicals for water cleaning (coagulation) Not much are known of their formation chemistry We have a lot of new mass spectrometer data of these complexes which needs computations to resolve the molecular structures
14 Ab initio molecular dynamics aluminum oxide chemistry in aqueous solution loosely bound proton (acidic) normal proton (water) 1.48 Å Distance (Å) Distance (Å) 1.35 Å Time (ps) Time (ps)
15 Ab initio molecular dynamics aluminum oxide chemistry in aqueous solution fast proton jump correlated to Cl escape attempt O-Cl Distance (Å) O-H 1.65 Å 1.98 Å ** 1.68 Å 3.0 Å Time (ps)
16 Constrained MD simulation One can fix some geometrical parameters and compute the force to this constraint. MetaDynamics allows treatment of more complex reaction. Free energy difference is an integral of this force Tedious calculations since they need long simulation to get good averaging. The constraint can slowly grow or it can be fixed (the later turned out to be more efficient) Al - Cl H H O Al --- Cl H H O
17 Constrained MD simulation Test the hysteresis grow and reduce the constraint. The result should be the same
18 Constrained MD simulation Also the static calculations need long simulations fs
19 Aluminum oxide chemistry in aqueous solution reaction barriers very large ligand effect: Al 1 ClOHw, Al 2 Clw 2 small barriers Cl s will dissociate 14 ± 3 kj mol ± 5 kj mol J. Saukkoriipi and K. Laasonen, J. Phys. Chem. A, 112, (2008), G (kj/mol) ,2 2,4 2,6 2,8 3,0 3,2 Distance (Å)
20 Metadynamics Laio, Parrinello, PNAS 99, 12562, 2002 Biased dynamics Choose few collective variables that will describe the problem of interest Run dynamics Add small repulsive gaussians to the system in places where you have been At the end the you will get the free energy profile (from the gaussians) A good review: Laio and Gervasio, Rep. Prog. Phys (2008)
21 Metadynamics The critical point, the choice of the collective variable The Z shape potential: 1 CV going from A to B cause easily too high barrier since CV do not go backwards 2 CV better can map this problem Deep problem: real system is very high dimensional, most of them are not relevant, BUT how many are. And how many CV s are need. The MetaDynamics (in CP2K) works with 2 CV s
22 Example: Al3On reactions (Giorgio Lanzani) CV:s 1: d(al1-al2)-d(al1-al3) 2: root mean displacement of the atoms forming the Al-O ring (reference: ring like Al3On) Why such CV s Alchemy making different molecules and finding good descriptors of them (experience in Hutter s group, especially M. Iannuzzi, has been very useful) After long (few 100 s ps) we receive very interesting and complex potential Surface. To produce this data with distance constrains would have been impossible.
23 Free energy surface of the considered Al--trimeric cluster at 300 K. The geometries corresponding at the local minima (marked in red) and saddle points (blue) are reported as well.
24 Conclusions The AIMD is a very poverfull tool for small systems (less than 1000 atoms and less than 5000 electrons) Time scale is limited to 1 ns The accuracy is dermined by the DFT-GGA (+vdw) Access to structural, dynamic and electronic properties Access to energy barriers
25
Supporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface
Supporting Information for Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite Edge Water Interface Zhen Jiang, Konstantin Klyukin, and Vitaly Alexandrov,, Department
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More informationIs there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich Chemistry Chemistry is the science of atomic matter, especially its chemical reactions,
More informationSupporting Information for. Dynamics Study"
Supporting Information for "CO 2 Adsorption and Reactivity on Rutile TiO 2 (110) in Water: An Ab Initio Molecular Dynamics Study" Konstantin Klyukin and Vitaly Alexandrov,, Department of Chemical and Biomolecular
More informationSupporting Information: Surface Polarons Reducing Overpotentials in. the Oxygen Evolution Reaction
Supporting Information: Surface Polarons Reducing Overpotentials in the Oxygen Evolution Reaction Patrick Gono Julia Wiktor Francesco Ambrosio and Alfredo Pasquarello Chaire de Simulation à l Echelle Atomique
More informationChapter 2. Ab Initio Molecular Dynamics. Kari Laasonen. Abstract. 1. Background
Chapter 2 Ab Initio Molecular Dynamics Kari Laasonen Abstract In this chapter, an introduction to ab initio molecular dynamics (AIMD) has been given. Many of the basic concepts, like the Hellman Feynman
More informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
More informationMetadynamics-biased ab initio Molecular Dynamics Study of
Metadynamics-biased ab initio Molecular Dynamics Study of Heterogeneous CO 2 Reduction via Surface Frustrated Lewis Pairs Mireille Ghoussoub 1, Shwetank Yadav 1, Kulbir Kaur Ghuman 1, Geoffrey A. Ozin
More informationSupplementary Information
Pt II as proton shuttle during C H bond activation in the Shilov process Pietro Vidossich,* a Gregori Ujaque, a Agustí Lledós a a Departament de Química, Universitat Autònoma de Barcelona, 08193 Cerdanyola
More informationSupporting Information
1 Supporting Information 2 3 4 Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds 5 6 7 Fabio Pietrucci 1, José C. Aponte 2,3,*, Richard Starr 2,4, Andrea Pérez-Villa
More informationDehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways from ab initio simulations
2nd CP2K Tutorial, ETHZ 7-11 February 2011 Materials Science & Technology Nanostructures and adsorption on metallic surfaces: Dehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways
More informationCHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued
CHE3935 Lecture 4 Quantum Mechanical Simulation Methods Continued 1 OUTLINE Review Introduction to CPMD MD and ensembles The functionals of density functional theory Return to ab initio methods Binding
More informationSupplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water
Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationCurrent address: Department of Chemistry, Hong Kong Baptist University, Kowloon Tong, Hong Kong,
Hydrolysis of Cisplatin - A Metadynamics Study Supporting Information Justin Kai-Chi Lau a and Bernd Ensing* b Department of Chemistry and Applied Bioscience, ETH Zurich, USI Campus, Computational Science,
More informationSupporting Information
Prebiotic NH 3 formation: Insights from simulations András Stirling 1, Tamás Rozgonyi 2, Matthias Krack 3, Marco Bernasconi 4 1 Institute of Organic Chemistry, Research Centre for Natural Sciences of the
More informationChemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase
Supplementary Information Chemical Dynamics of the First Proton Coupled Electron Transfer of Water Oxidation on TiO 2 Anatase Jia Chen, Ye-Fei Li, Patrick Sit, and Annabella Selloni Department of Chemistry,
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationSupporting Information
alladium monophosphine d(h 3 ): does it really exist in solution? ietro Vidossich*, Gregori Ujaque,* and Agustí Lledós* Supporting Information Contents: 1. Schemes S1 and S2 2. Computational details 3.
More informationSupporting Information Synthesis of (D)-erythrose from glycolaldehyde aqueous solutions under electric field
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supporting Information Synthesis of (D)-erythrose from glycolaldehyde aqueous solutions under electric
More informationElectronic structure simulations of water solid interfaces
Electronic structure simulations of water solid interfaces Angelos Michaelides London Centre for Nanotechnology & Department of Chemistry, University College London www.chem.ucl.ac.uk/ice Main co-workers:
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationCP2K: Past, Present, Future. Jürg Hutter Department of Chemistry, University of Zurich
CP2K: Past, Present, Future Jürg Hutter Department of Chemistry, University of Zurich Outline Past History of CP2K Development of features Present Quickstep DFT code Post-HF methods (RPA, MP2) Libraries
More informationExperiment Section Fig. S1 Fig. S2
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supplementary Materials Experiment Section The STM experiments were carried out in an ultrahigh
More informationExplanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph.
Supplementary Materials Explanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph. Tao Cheng,, Lu Wang, Boris V Merinov, and William
More informationBasic introduction of NWChem software
Basic introduction of NWChem software Background NWChem is part of the Molecular Science Software Suite Designed and developed to be a highly efficient and portable Massively Parallel computational chemistry
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationarxiv: v2 [physics.comp-ph] 4 Feb 2016
Sampling Free Energy Surfaces as Slices by Combining Umbrella Sampling and Metadynamics Shalini Awasthi, Venkat Kapil and Nisanth N. Nair 1 arxiv:1507.06764v2 [physics.comp-ph] 4 Feb 2016 Department of
More informationIntroduction to CPMD 1
Introduction to CPMD 1 CPMD program ab initio electronic structure (DFT) and molecular dynamics program plane wave basis set (PBC), pseudopotentials massively parallelized, linear scaling up to thousands
More informationMODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
More informationResonances in Chemical Reactions : Theory and Experiment. Toshiyuki Takayanagi Saitama University Department of Chemistry
Resonances in Chemical Reactions : Theory and Experiment Toshiyuki Takayanagi Saitama University Department of Chemistry What is Chemical Reaction? Collision process between molecules (atoms) containing
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More informationElectronic Supporting Information Topological design of porous organic molecules
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2017 Electronic Supporting Information Topological design of porous organic molecules Valentina Santolini,
More informationStudy of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal
Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationWavelets for density functional calculations: Four families and three. applications
Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Interpolating
More informationMolecular and Electronic Dynamics Using the OpenAtom Software
Molecular and Electronic Dynamics Using the OpenAtom Software Sohrab Ismail-Beigi (Yale Applied Physics) Subhasish Mandal (Yale), Minjung Kim (Yale), Raghavendra Kanakagiri (UIUC), Kavitha Chandrasekar
More informationSupporting Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 214 Infrared Spectroscopy from Ab Initio Molecular Dynamics - the MeCN-HCl Molecular
More informationDay 1 : Introduction to AIMD and PIMD
Day 1 : Introduction to AIMD and PIMD Aug 2018 In today s exercise, we perform ab-initio molecular dynamics (AIMD) and path integral molecular dynamics (PIMD) using CP2K[?]. We will use the Zundel s cation
More informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationAccounting for Solvation in Quantum Chemistry. Comp chem spring school 2014 CSC, Finland
Accounting for Solvation in Quantum Chemistry Comp chem spring school 2014 CSC, Finland In solution or in a vacuum? Solvent description is important when: Polar solvent: electrostatic stabilization Solvent
More informationYuan Ping 1,2,3*, Robert J. Nielsen 1,2, William A. Goddard III 1,2*
Supporting Information for the Reaction Mechanism with Free Energy Barriers at Constant Potentials for the Oxygen Evolution Reaction at the IrO2 (110) Surface Yuan Ping 1,2,3*, Robert J. Nielsen 1,2, William
More informationRepresenting High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks
Energy Landscapes, Chemnitz 200 Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks Jörg Behler Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum D-44780 Bochum,
More informationAb initio molecular dynamics: Basic Theory and Advanced Met
Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationMonte Carlo strategies for first-principles simulations of elemental systems
Monte Carlo strategies for first-principles simulations of elemental systems Lev Gelb Department of Materials Science and Engineering, University of Texas at Dallas XSEDE12 Lev Gelb (UT Dallas) Monte Carlo
More informationSelf assembly of Carbon Atoms in Interstellar Space and Formation Mechanism of Naphthalene from Small Precursors: A Molecular Dynamics Study
Self assembly of Carbon Atoms in Interstellar Space and Formation Mechanism of Naphthalene from Small Precursors: A Molecular Dynamics Study Niladri Patra UIC Advisor: Dr. H. R. Sadeghpour ITAMP Large
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationRustam Z. Khaliullin University of Zürich
Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationJoint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012
2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.
More informationAb initio molecular dynamics
Ab initio molecular dynamics Molecular dynamics Why? allows realistic simulation of equilibrium and transport properties in Nature ensemble averages can be used for statistical mechanics time evolution
More informationAb initio molecular dynamics simulation of hydrogen fluoride at several thermodynamic states
Downloaded from orbit.dtu.dk on: Sep 24, 2018 Ab initio molecular dynamics simulation of hydrogen fluoride at several thermodynamic states Kreitmeir, M.; Bertagnolli, H.; Mortensen, Jens Jørgen; Parrinello,
More informationSupplementary material. From cellulose to kerogen: molecular simulation. of a geological process
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2017 Supplementary material From cellulose to kerogen: molecular simulation of a geological
More informationXYZ file format Protein Data Bank (pdb) file format Z - matrix
Chemistry block (exercise 1) In this exercise, students will be introduced how to preform simple quantum chemical calculations. Input files for Gaussian09. Output file structure. Geometry optimization,
More informationInteraction Induced Localisation of Protons Temperature and Superfluid Helium Solvation Effects Łukasz Walewski Ruhr-University Bochum
Interaction Induced Localisation of Protons Temperature and Superfluid Helium Solvation Effects Łukasz Walewski Ruhr-University Bochum VI Workshop in Electronvolt Neutron Spectroscopy, Abingdon 2014 Interaction
More informationA Posteriori Quantization of the Nuclear Motion by Backfill Calculations. Mateusz Brela Wydział Chemii Uniwersytet Jagielloński
A Posteriori Quantization of the Nuclear Motion by Backfill Calculations Mateusz Brela Wydział Chemii Uniwersytet Jagielloński Introduction The computations of the IR spectra of the hydrogen bonded systems
More informationThe PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics
The PLUMED plugin and free energy methods in electronic-structure-based molecular dynamics Davide Branduardi, Theoretical Molecular Biophysics Group Max Planck for Biophysics, Frankfurt am Main, Germany
More informationAri P Seitsonen CNRS & Université Pierre et Marie Curie, Paris
Self-organisation on noble metal surfaces Ari P Seitsonen CNRS & Université Pierre et Marie Curie, Paris Collaborations Alexandre Dmitriev, Nian Lin, Johannes Barth, Klaus Kern,... Thomas Greber, Jürg
More informationProbing the Origins of Intermolecular Vibrational and Relaxational Dynamics in Organic Solids with CP2K
Probing the Origins of Intermolecular Vibrational and Relaxational Dynamics in Organic Solids with CP2K Michael Ruggiero Department of Chemical Engineering and Biotechnology, University of Cambridge CP2K
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationModeling of Nanostructures and Materials Jacek A. Majewski Summer Semester 2013 Lecture 6 Lecture Modeling of Nanostructures Molecular Dynamics
Chair of Condensed Matter Physics nstitute of Theoretical Physics Faculty of Physics, Universityof Warsaw Summer Semester 013 Lecture Modeling of Nanostructures and Materials Jacek A. Majewski E-mail:
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationChapter 6 Cyclic urea - a new central unit in bent-core compounds
82 Chapter 6 Cyclic urea - a new central unit in bent-core compounds A new class of five-ring bent-core molecules with a cyclic urea group as a central unit was synthesized [94]. A significant difference
More informationMaterial Surfaces, Grain Boundaries and Interfaces: Structure-Property Relationship Predictions
Material Surfaces, Grain Boundaries and Interfaces: Structure-Property Relationship Predictions Susan B. Sinnott Department of Materials Science and Engineering Penn State University September 16, 2016
More informationDesign of Efficient Catalysts with Double Transition Metal. Atoms on C 2 N Layer
Supporting Information Design of Efficient Catalysts with Double Transition Metal Atoms on C 2 N Layer Xiyu Li, 1, Wenhui Zhong, 2, Peng Cui, 1 Jun Li, 1 Jun Jiang 1, * 1 Hefei National Laboratory for
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationUniversity of Zurich. Proton transfer in imidazole-based molecular crystals. Zurich Open Repository and Archive. Iannuzzi, M.
University of Zurich Zurich Open Repository and Archive Winterthurerstr. 190 CH-8057 Zurich http://www.zora.uzh.ch Year: 2006 Proton transfer in imidazole-based molecular crystals Iannuzzi, M Iannuzzi,
More informationReaction pathway and free energy barrier for defect-induced water dissociation on the 101 surface of TiO 2 -anatase
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 14 8 OCTOBER 2003 Reaction pathway and free energy barrier for defect-induced water dissociation on the 101 surface of TiO 2 -anatase Antonio Tilocca a) and
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More informationCar-Parrinello Molecular Dynamics
Car-Parrinello Molecular Dynamics Eric J. Bylaska HPCC Group Moral: A man dreams of a miracle and wakes up with loaves of bread Erich Maria Remarque Molecular Dynamics Loop (1) Compute Forces on atoms,
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationThe mutagenic action of 5-bromouracil
The mutagenic action of 5-bromouracil The structure of DNA DNA bases Cytosine (C) Thymine (T) oxygen carbon nitrogen hydrogen Guanine (G) Adenine (A) Hydrogen bonds picture from wikipedia picture from
More informationThe structure and spectroscopy of monolayers of water on MgO: An ab initio study
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 8 22 AUGUST 2000 The structure and spectroscopy of monolayers of water on MgO: An ab initio study L. Delle Site, a) A. Alavi, and R. M. Lynden-Bell b) Atomistic
More informationComputing free energy: Thermodynamic perturbation and beyond
Computing free energy: Thermodynamic perturbation and beyond Extending the scale Length (m) 1 10 3 Potential Energy Surface: {Ri} 10 6 (3N+1) dimensional 10 9 E Thermodynamics: p, T, V, N continuum ls
More informationQuantum chemical calculations and MD simulations for Be
start Quantum chemical calculations and MD simulations for Be Institute of Ion Physics University of Innsbruck Michael Probst Michael Stefan Alexander Ivan Andreas Györoek Huber Kaiser Sukuba Mauracher
More informationUnderstanding hydrogen storage in metal organic frameworks from first principles
Understanding hydrogen storage in metal organic frameworks from first principles Understanding hydrogen storage in metal organic frameworks from first principles Sohrab Ismail-Beigi (Yale Applied Physics)
More informationDouble Bond: C 2 H 4. An sp 2 hybridized C atom has one electron in each of the three sp 2 lobes
Double Bond: C 2 H 4 An sp 2 hybridized C atom has one electron in each of the three sp 2 lobes Top view of the sp 2 hybrid Side view of the sp 2 hybrid + the unhybridized p orbital 1 Double Bond: C 2
More informationFree energy sampling for electrochemical systems
Free energy sampling for electrochemical systems Mira Todorova, Anoop Kishore Vatti, Suhyun Yoo and Jörg Neugebauer Department of Computational Materials Design Düsseldorf, Germany m.todorova@mpie.de IPAM,
More informationSupporting information. Realizing Two-Dimensional Magnetic Semiconductors with. Enhanced Curie Temperature by Antiaromatic Ring Based
Supporting information Realizing Two-Dimensional Magnetic Semiconductors with Enhanced Curie Temperature by Antiaromatic Ring Based Organometallic Frameworks Xingxing Li and Jinlong Yang* Department of
More informationPost Hartree-Fock: MP2 and RPA in CP2K
Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,
More informationMolecular Interactions F14NMI. Lecture 4: worked answers to practice questions
Molecular Interactions F14NMI Lecture 4: worked answers to practice questions http://comp.chem.nottingham.ac.uk/teaching/f14nmi jonathan.hirst@nottingham.ac.uk (1) (a) Describe the Monte Carlo algorithm
More informationPhotoinduced Water Oxidation at the Aqueous. GaN Interface: Deprotonation Kinetics of. the First Proton-Coupled Electron-Transfer Step
Supporting Information Photoinduced Water Oxidation at the Aqueous Interface: Deprotonation Kinetics of the First Proton-Coupled Electron-Transfer Step Mehmed Z. Ertem,,,* eerav Kharche,,* Victor S. Batista,
More informationPerforming Metadynamics Simulations Using NAMD
Performing Metadynamics Simulations Using NAMD Author: Zhaleh Ghaemi Contents 1 Introduction 2 1.1 Goals of this tutorial.......................... 2 2 Standard metadynamics simulations 3 2.1 Free energy
More informationMultiscale modelling of D trapping in W
CMS Multiscale modelling of D trapping in W Kalle Heinola, Tommy Ahlgren and Kai Nordlund Department of Physics and Helsinki Institute of Physics University of Helsinki, Finland Contents Background Plasma-wall
More informationSchool of Chemistry, University of KwaZulu-Natal, Howard College Campus, Durban. CHEM191 Tutorial 1: Buffers
School of Chemistry, University of KwaZulu-Natal, Howard College Campus, Durban CHEM191 Tutorial 1: Buffers Preparing a Buffer 1. How many moles of NH 4 Cl must be added to 1.0 L of 0.05 M NH 3 to form
More informationBCIT Fall Chem Exam #1
BCIT Fall 2012 Chem 3615 Exam #1 Name: Attempt all questions in this exam. Read each question carefully and give a complete answer in the space provided. Part marks given for wrong answers with partially
More informationAdsorption of Iodine on Pt(111) surface. Alexandre Tkachenko Marcelo Galván Nikola Batina
Adsorption of Iodine on Pt(111) surface Alexandre Tkachenko Marcelo Galván Nikola Batina Outline Motivation Experimental results Geometry Ab initio study Conclusions Motivation Unusual structural richness
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationBenchmark of the CPMD code on CRESCO HPC Facilities for Numerical Simulation of a Magnesium Nanoparticle.
Benchmark of the CPMD code on CRESCO HPC Facilities for Numerical Simulation of a Magnesium Nanoparticle. Simone Giusepponi a), Massimo Celino b), Salvatore Podda a), Giovanni Bracco a), Silvio Migliori
More informationAb initio molecular dynamics : BO
School on First Principles Simulations, JNCASR, 2010 Ab initio molecular dynamics : BO Vardha Srinivasan IISER Bhopal Why ab initio MD Free of parametrization and only fundamental constants required. Bond
More informationMolecular Dynamics on the Angstrom Scale
Probing Interface Reactions by STM: Molecular Dynamics on the Angstrom Scale Zhisheng Li Prof. Richard Osgood Laboratory for Light-Surface Interactions, Columbia University Outline Motivation: Why do we
More informationEfficient use of CP2K
Efficient use of CP2K Fawzi Mohamed 1 1 Department of Chemisty Humboldt University of Berlin June 10, 2009 cp2k Fortran 95 opensource (GPL) general framework for different methods DFT no official release
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationYali Liu, Pengfei Zhang, Junmin Liu, Tao Wang, Qisheng Huo, Li Yang, Lei. Sun,*, Zhen-An Qiao,*, and Sheng Dai *, ASSOCIATED CONTENT
ASSOCIATED CONTENT Supporting Information Gold Cluster-CeO 2 Nanostructured Hybrid Architectures as Catalysts for Selective Oxidation of Inert Hydrocarbons Yali Liu, Pengfei Zhang, Junmin Liu, Tao Wang,
More informationExploring the anomalous behavior of metal nanocatalysts with finite temperature AIMD and x-ray spectra
Exploring the anomalous behavior of metal nanocatalysts with finite temperature AIMD and x-ray spectra F.D. Vila DOE grant DE-FG02-03ER15476 With computer support from DOE - NERSC. Importance of Theoretical
More informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
More information