Supporting Information
|
|
- Dwain Lucas
- 5 years ago
- Views:
Transcription
1 1 Supporting Information Hydrothermal Decomposition of Amino Acids and Origins of Prebiotic Meteoritic Organic Compounds Fabio Pietrucci 1, José C. Aponte 2,3,*, Richard Starr 2,4, Andrea Pérez-Villa 1, Jamie E. Elsila 2, Jason P. Dworkin 2, A. Marco Saitta 1,* Sorbonne Université, Muséum National d Histoire Naturelle, UMR CNRS 7590, IRD, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, F-75005, Paris, France 2 The Goddard Center for Astrobiology and Solar System Exploration Division, Code 691, NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, Maryland 20771, United States of America 3 Department of Chemistry, The Catholic University of America, 620 Michigan Ave. NE, Washington, DC 20064, United States of America 4 Physics Department, The Catholic University of America, 620 Michigan Ave. NE, Washington, DC 20064, United States of America * To whom correspondence may be addressed. jose.c.aponte@nasa.gov (tel.: ) or marco.saitta@sorbonne-universite.fr (tel.: )
2 Born-Oppenheimer molecular dynamics (BOMD) Ab initio calculations were based on the Density Functional Theory with the exchange-correlation Perdew-Burke-Ernzerhof functional [1], including Grimme s dispersion corrections [2,3] and the Martins- Troullier pseudopotentials [4] for C, N, O, and H atoms, as implemented in the code CPMD 4.1 [5]. BOMD simulations were performed in the canonical ensemble (NVT) at a temperature of 373 K, with a time step of 0.48 fs, with electronic wavefunction plane-wave cutoff of 70 Ry and convergence threshold of 10-5, in a periodically repeated cubic box of 13.2x13.1x13.8 Å for glycine (247 atoms) and 13.4x13.3x13.8 Å for isovaline (256 atoms), including 1 amino acid and 79 water molecules (density 1.03 g/cm 3 ). The temperature was controlled with a massive Nose-Hoover thermostat [6,7] with a chain length of 4 and with a frequency of 3000 cm Enhanced sampling simulations BOMD simulations were coupled with enhanced sampling techniques [8], according to the following protocol: In the first place, we performed metadynamics simulations [9, 10], employing the path collective variables s and z as the set of reaction coordinates [11]. s indicates the progress along a putative pathway composed by a discrete sequence of atomic configurations, while z represents the distance from the putative pathway. The path collective variables are defined by the following functional form: where D(R, R k ) corresponds to to the distance between the current atomic configuration and the k-th reference configuration: in this work we employed only two references: the reactants and the putative products, thus avoiding to bias the results with hypotheses on the reaction mechanism. The parameter λ is chosen such that λd(r k, R k+1 ) 2.3, conveniently localizing the free energy basins of the reference states around s 1.1 and z 1.9, respectively, and leading to smooth pathways and landscapes. A much larger λ would produce very irregular and discontinuous pathways, while a much smaller one would hamper the resolution of different free energy minima. We adopted a topological metric suited to chemical reactions in solution, introduced in Ref. [12], that tracks changes in chemical bond networks passing from reactants, through intermediates / transition states, until the products: 48
3 Atomic structures are therefore compared based on the difference between the coordination numbers (C IS (t), defined smoothly as in the formula above) of a specific atom I with all atoms of species S. Parameters r 0 SS, n, and m are related to the typical bond distances between two atomic species. We adopted n = 8, m = 14, and r 0 SS = 1.8 and 1.5 Å for X-X and X-H pairs respectively, with X = C, N, O. Reference coordination patterns for reactants and products were obtained as averages from equilibrium BOMD simulations of 10 ps. The patterns are reported in the following tables: Glycine deamination: {C} {O} {N} {H} Ref. 1 C OO C α N Ref. 2 C OO C α N Glycine decarboxylation: {C} {O} {N} {H} Ref. 1 C OO C α N Ref. 2 C OO C α N Isovaline deamination: {C} {O} {N} {H} Ref. 1 C OO C α N Ref. 2 C OO C α N
4 Isovaline decarboxylation: {C} {O} {N} {H} Ref. 1 C OO C α N Ref. 2 C OO C α N Ref. 3 C OO C α N Ref. 4 C OO C α N Ref. 5 C OO C α N Angelic acid decarboxylation: {C} {O} {N} {H} Ref. 1 C OO C α C =Cα N Ref. 2 C OO C α C =Cα N In a first stage, we performed metadynamics simulations [10] for the different reaction channels: a history-dependent bias potential is added to the potential energy of the system, defined as a sum of Gaussian functions of the collective variables. The procedure allows to quickly overcome barriers and explore the relevant free energy landscape. The Gaussian widths for the path CVs were set at σ s = and σ z =0.2, with a height of 3.1 kcal/mol and a deposition interval of 24 fs. The typical duration of each simulation was between 15 and 30 ps. In a second stage, to identify the transition state of each reaction and to get more insight about its mechanism, we carried out a commitor analysis [13]. For each reaction, we extracted a set of up to 20 atomic configurations close to the saddle point in the free energy landscape from reactive metadynamics (and umbrella sampling, see below) trajectories, and we performed 10 independent unbiased BOMD trajectories of length between 0.24 and 0.72 ps, starting from each configuration, verifying whether they fell into reactants or products basins. If the 10 trajectories split between the two basins, additional trajectories were generated until a total of 100. A transition state configuration is
5 identified as committed to both basins with a probability of 50±10%. Note that isovaline deamination (Fig. 3a in the main text) leads to a three-way transition state: we explicitly verified that the corresponding configurations relax to each of the three basins with similar probability. Overall, committor analysis relaxation trajectories are free from bias potentials, hence represent realistic reactive trajectories of the system and were analyzed to understand the reaction mechanisms. In a third stage, a series of umbrella sampling [14] simulations were performed, systematically restraining the s coordinate at different locations along the reaction pathway obtained by commitor analysis (or metadynamics) using a harmonic bias potential in the form 0.5 k (s-s 0 ) 2. For each reaction we adopted between 25 and 60 equally-spaced windows, with force constant set to k = k B T /(Δs/ 2.5 ) 2 with Δs the spacing between windows. We performed for every window a trajectory of duration between 5 and 10 ps. The weighted histogram analysis method (WHAM)[15] finally yielded the free energy profile of the reaction, with a statistical error bar of about ±1 kcal/mol. The error is estimated (in a conservative way) by comparing the profiles obtained using data in different intervals: [1, 10] ps, [1, 5.5] ps, and [5.5, 10] ps. For comparison, the autocorrelation time of the collective variables is less than 0.1 ps. The umbrella sampling simulation of isovaline decarboxylation to sec-butylamine (Fig. 3c in the main text) employed 5 reference configurations instead of only 2 in the definition of the path CV, to improve the resolution of the free energy landscape. The 3 additional intermediate reference configurations were obtained from preliminary metadynamics and umbrella sampling simulations based on 2 reference configurations only. For isovaline deamination an additional series of umbrella sampling simulations were performed along the z coordinate, obtaining through the WHAM the two-dimensional free energy landscape shown in Fig. 3a in the main text. All enhanced sampling simulations employed a modified version of the plugin Plumed 1.3 [16]. A guide to the use of the coordination pattern-based path variables employed in this work, as well as source code and example input files for Plumed 1.3 and Plumed 2.x can be freely downloaded from
6 Table S1: Summary of simulated amino acid decomposition reactions, with their corresponding free energy barriers and reaction free energy differences (values in kcal/mol) Reactant Product(s) Mechanism ΔF ǂ ΔF Glycine (1) Glycolic acid (2) + Ammonia Deamination Glycine (1) Glycolate (3) + Ammonium Deamination Glycine (1) Methylammonium (4) Decarboxylation Carbon dioxide + Hydroxide Isovaline (5) 2,2-Ethyl-methyl-α-acetolactone (6) + Ammonia Deamination Isovaline (5) cis-2-methyl-2-butenoate (7) Deamination cis-2-methyl-2- butenoate (7) + Ammonium Butene (8) + Carbon dioxide + Ammonia Decarboxylation Isovaline (5) Sec-butylamine (9) Decarboxylation
7 Figure S1. Schematic summary of glycine decomposition mechanisms. Note that the arrow convention, representing the displacement of electron pairs, provides a convenient schematization of the mechanisms, however for a precise description of the latter (that are typically concerted and may involve proton transfer through long solvent chains) we refer to the main text, as well as to Fig. 4 and 5 and to the Supplementary Movies
8 Figure S2. Schematic summary of isovaline decomposition mechanisms. The same remarks apply as for Fig. S
9 Movie S1. Deamination of (1), formation of (2) and ammonia. Movie S2. Decarboxylation of (1), formation of (4), ammonia and hydroxide. Movie S3. Deamination of (5), formation of (6) and ammonia. Movie S4. Deamination of (5), formation of (7). Movie S5. Decarboxylation of (5), formation of (9). Movie S6. Decarboxylation of (7), formation of (8), ammonia and carbon dioxide.
10 References [1] J. P. Perdew, K. Burke, and M. Ernzerhof. Generalized gradient approximation made simple. Phys. Rev. Lett., 77: , Oct [2] S. Grimme. Semiempirical gga-type density functional constructed with a long-range dispersion correction. J. Comput. Chem., 27(15): , [3] V. Barone, M. Casarin, D. Forrer, M. Pavone, M. Sambi, and A. Vittadini. Role and effective treatment of dispersive forces in materials: Polyethylene and graphite crystals as test cases. J. Comput. Chem., 30(6): , [4] N. Troullier and J. L. Martins. Efficient pseudopotentials for plane-wave calculations. Phys. Rev. B, 43: , Jan [5] CPMD Copyright IBM Corp. and by Max Planck Institute Stuttgart [6] S. Nosé. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys., 81(1): , [7] W. G. Hoover. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A, 31: , Mar [8] F. Pietrucci. Strategies for the exploration of free energy landscapes: unity in diversity and challenges ahead. Rev. Phys., 2:32 45, [9] A. Laio and M. Parrinello. Escaping free-energy minima. Proc. Natl. Acad. Sci. U.S.A., 99(20): , [10] A. Laio and F. L. Gervasio. Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science. Rep. Prog. Phys., 71(12):126601, [11] D. Branduardi, F. L. Gervasio, and M. Parrinello. From a to b in free energy space. J. Chem. Phys.,126(5):054103, [12] F. Pietrucci and A. M. Saitta. Formamide reaction network in gas phase and solution via a unified theoretical approach: Toward a reconciliation of different prebiotic scenarios. Proc. Natl. Acad. Sci. U.S.A., 112(49): , [13] P. G. Bolhuis, D. Chandler, C. Dellago, and P. L. Geissler. Transition path sampling: Throwing ropes over rough mountain passes, in the dark. Annu. Rev. Phys. Chem., 53(1): , [14] G.M. Torrie and J.P. Valleau. Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling. J. Comput. Phys., 23(2): , 1977.
11 [15] S. Kumar, J. M. Rosenberg, D. Bouzida, R. H. Swendsen, and P. A. Kollman. The weighted histogram analysis method for free-energy calculations on biomolecules. i. the method. J. Comput. Chem., 13(8): , [16] M. Bonomi, D. Branduardi, G. Bussi, C. Camilloni, D. Provasi, P. Raiteri, D. Donadio, F. Marinelli, F. Pietrucci, R. A. Broglia, and M. Parrinello. Plumed: A portable plugin for free-energy calculations with molecular dynamics. Comput. Phys. Commun., 180(10): , 2009.
Supporting 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 informationSupporting 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 informationA Tunable, Strain-Controlled Nanoporous MoS 2 Filter for Water Desalination
Supporting Information A Tunable, Strain-Controlled Nanoporous MoS 2 Filter for Water Desalination Weifeng Li 1, Yanmei Yang 1, Jeffrey K. Weber 2, Gang Zhang* 3, Ruhong Zhou* 1,2,4 1. School for Radiological
More informationSupporting Information
Supporting Information The effect of strong acid functional groups on electrode rise potential in capacitive mixing by double layer expansion Marta C. Hatzell, a Muralikrishna Raju, b Valerie J. Watson,
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 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 informationComputing free energies with PLUMED 2.0
Computing free energies with PLUMED 2.0 Davide Branduardi Formerly at MPI Biophysics, Frankfurt a.m. TyOutline of the talk Relevance of free energy computation Free-energy from histograms: issues Tackling
More informationFrom fundamental physics to the origins of life: ab initio Miller experiments
From fundamental physics to the origins of life: ab initio Miller experiments A. Marco Saitta Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie (IMPMC) Université Pierre et Marie Curie
More informationAb initio molecular dynamics
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
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 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 informationCharacterization of the free-energy landscapes of proteins by NMR-guided metadynamics
Characterization of the free-energy landscapes of proteins by NMR-guided metadynamics B Results and Discussion BIOPHYSICS ND COMPUTTIONL BIOLOGY SI Text Inset B χ χ Inset SI Text Inset C Folding of GB3
More informationa micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics FHI-aims Developers and Users Meeting Organized by Fritz Haber Institute, CECAM-mm1p,
More informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
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 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 informationa micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED Tutorial for FHI-aims users PLUMED: portable plugin for free-energy calculations with molecular dynamics Density functional theory and beyond: Computational materials science for real materials
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
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 informationa micro PLUMED2.0 Tutorial for ESPResSo users PLUMED2.0: portable plugin for free-energy calculations with molecular dynamics
a micro PLUMED2.0 Tutorial for ESPResSo users PLUMED2.0: portable plugin for free-energy calculations with molecular dynamics ESPResSo summer school 2013 Organized by CECAM, SimTech, ICP and University
More informationFlexMD Manual. ADF Program System Release 2013
FlexMD Manual ADF Program System Release 2013 Scientific Computing & Modelling NV Vrije Universiteit, Theoretical Chemistry De Boelelaan 1083; 1081 HV Amsterdam; The Netherlands WWW: www.scm.com E-mail:
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 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 informationQuick reference guide on PLUMED with Quantum ESPRESSO
Quick reference guide on PLUMED with Quantum ESPRESSO Changru Ma SISSA, Trieste March 30, 2011 Contents 1 Introduction 2 1.1 Overview................................... 2 1.2 Collective variables.............................
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 informationBacterial Outer Membrane Porins as Electrostatic Nanosieves: Exploring Transport Rules of Small Polar Molecules
Bacterial Outer Membrane Porins as Electrostatic Nanosieves: Exploring Transport Rules of Small Polar Molecules Harsha Bajaj, Silvia Acosta Gutiérrez, Igor Bodrenko, Giuliano Malloci, Mariano Andrea Scorciapino,
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 informationUnraveling the degradation of artificial amide bonds in Nylon oligomer hydrolase: From induced-fit to acylation processes
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 Supporting Information for Unraveling the degradation of artificial amide bonds
More informationMETAGUI A VMD EXTENSION TO ANALYZE AND VISUALIZE METADYNAMICS SIMULATIONS
METAGUI A VMD EXTENSION TO ANALYZE AND VISUALIZE METADYNAMICS SIMULATIONS Alessandro Laio SISSA & DEMOCRITOS, Trieste Coworkers: Xevi Biarnes Fabio Pietrucci Fabrizio Marinelli Metadynamics (Laio A. and
More informationMaking the best of a bad situation: a multiscale approach to free energy calculation
Making the best of a bad situation: arxiv:1901.04455v2 [physics.comp-ph] 15 Jan 2019 a multiscale approach to free energy calculation Michele Invernizzi,, and Michele Parrinello, Department of Physics,
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 informationMetadynamics with adaptive Gaussians
Metadynamics with adaptive Gaussians Davide Branduardi Theoretical Molecular Biophysics Group, Max Planck Institute for Biophysics, Max-von-Laue strasse 5, 60438, Frankfurt am Main, Germany Giovanni Bussi
More informationDetermining Energy Barriers and Selectivities of a Multi-Pathway System With Infrequent Metadynamics
Determining Energy Barriers and Selectivities of a Multi-Pathway System With Infrequent Metadynamics Christopher D. Fu 1, Luiz F.L. Oliveira 1 and Jim Pfaendtner 1,a) 1 Department of Chemical Engineering,
More informationA self-activated mechanism for nucleic acid polymerization catalyzed by DNA/RNA polymerases
A self-activated mechanism for nucleic acid polymerization catalyzed by DNA/RNA polymerases Vito Genna 1,2,, Pietro Vidossich 2, Emiliano Ippoliti 2, Paolo Carloni 2* and Marco De Vivo 1,2* 1 Laboratory
More informationSupporting Information
Supporting Information Structure and Dynamics of Uranyl(VI) and Plutonyl(VI) Cations in Ionic Liquid/Water Mixtures via Molecular Dynamics Simulations Katie A. Maerzke, George S. Goff, Wolfgang H. Runde,
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 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 informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationApplication of Computational Modelling to Protein Folding and Aggregation Studies
Application of Computational Modelling to Protein Folding and Aggregation Studies NEVENA TODOROVA, ANDREW HUNG AND IRENE YAROVSKY School of Applied Sciences RMIT University GPO Box 2476V Melbourne Victoria
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 informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
More informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
More informationSupplementary Materials for Oxygen-induced self-assembly of quaterphenyl molecule on metal surfaces
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2014 Supplementary Materials for Oxygen-induced self-assembly of quaterphenyl molecule on metal surfaces
More informationarxiv: v1 [q-bio.bm] 9 Jul 2007
Exploring the Protein G Helix Free Energy Surface by Solute Tempering Metadynamics Carlo Camilloni, 1, 2 Davide Provasi, 1, 2 Guido Tiana, 1, 2, and Ricardo A. Broglia 1,2,3 1 Department of Physics, University
More informationSupplementary Information. Supplementary Figure 1 Synthetic routes to the organic linker H 2 ATBDC.
Supplementary Information Supplementary Figure 1 Synthetic routes to the organic linker H 2 ATBDC. S1 Supplementary Figure 2 1 H NMR (D 2 O, 500MHz) spectrum of H 2 ATBDC. S2 Supplementary Figure 3 13
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 informationSupporting Information
Supporting Information Interface-Induced Affinity Sieving in Nanoporous Graphenes for Liquid-Phase Mixtures Yanan Hou, Zhijun Xu, Xiaoning Yang * State Key Laboratory of Material-Orientated Chemical Engineering,
More informationCrossing the barriers - simulations of activated processes
Crossing the barriers - simulations of activated processes Mgr. Ján Hreha for 6 th Student Colloquium and School on Mathematical Physics Faculty of Mathematics, Physics and Informatics Comenius University
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 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 informationDynamics of Nucleation in the Ising Model
J. Phys. Chem. B 2004, 108, 19681-19686 19681 Dynamics of Nucleation in the Ising Model Albert C. Pan and David Chandler* Department of Chemistry, UniVersity of California, Berkeley, California 94720-1460
More informationMetadynamics. Day 2, Lecture 3 James Dama
Metadynamics Day 2, Lecture 3 James Dama Metadynamics The bare bones of metadynamics Bias away from previously visited configuraaons In a reduced space of collecave variables At a sequenaally decreasing
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 informationCPMD Tutorial Atosim/RFCT 2009/10
These exercices were inspired by the CPMD Tutorial of Axel Kohlmeyer http://www.theochem.ruhruni-bochum.de/ axel.kohlmeyer/cpmd-tutor/index.html and by other tutorials. Here is a summary of what we will
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More informationSampling the free energy surfaces of collective variables
Sampling the free energy surfaces of collective variables Jérôme Hénin Enhanced Sampling and Free-Energy Calculations Urbana, 12 September 2018 Please interrupt! struct bioinform phys chem theoretical
More informationHydrogenation of Penta-Graphene Leads to Unexpected Large. Improvement in Thermal Conductivity
Supplementary information for Hydrogenation of Penta-Graphene Leads to Unexpected Large Improvement in Thermal Conductivity Xufei Wu, a Vikas Varshney, b,c Jonghoon Lee, b,c Teng Zhang, a Jennifer L. Wohlwend,
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 informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2016 Electronic Supplementary Information Active Site Structure and Absorption Spectrum of Channelrhodopsin-2
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 2 Feb 2005
Topological defects and bulk melting of hexagonal ice Davide Donadio, Paolo Raiteri, and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus,
More informationSupporting Information
Supporting Information Indirect Four-Electron Oxygen Reduction Reaction on Carbon Materials Catalysts in Acidic Solutions Guo-Liang Chai* 1, Mauro Boero 2, Zhufeng Hou 3, Kiyoyuki Terakura 3,4 and Wendan
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Electronic Supplementary Information Rational modifications on champion porphyrin
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 informationarxiv: v3 [physics.comp-ph] 30 Apr 2009
PLUMED: a portable plugin for free-energy calculations with molecular dynamics arxiv:0902.0874v3 [physics.comp-ph] 30 Apr 2009 Massimiliano Bonomi, 1 Davide Branduardi, 1 Giovanni Bussi, 2 Carlo Camilloni,
More informationMolecular Dynamics and Accelerated Molecular Dynamics
Molecular Dynamics and Accelerated Molecular Dynamics Arthur F. Voter Theoretical Division National Laboratory Lecture 3 Tutorial Lecture Series Institute for Pure and Applied Mathematics (IPAM) UCLA September
More informationOn the Reaction Rate Calculation Methods using Path Sampling
On the Reaction Rate Calculation Methods using Path Sampling Paris, France France 28 November, 2006 Titus S. van Erp Centre for Surface Chemistry and Catalysis, KULeuven Leuven, Belgium Introduction Straightforward
More informationarxiv: v1 [cond-mat.mtrl-sci] 29 Jan 2015
Interatomic potentials for ionic systems with density functional accuracy based on charge densities obtained by a neural network S. Alireza Ghasemi 1,, Albert Hofstetter 2, Santanu Saha 2 and Stefan Goedecker
More informationDescription: Supplementary Figures, Supplementary Methods, and Supplementary References
File Name: Supplementary Information Description: Supplementary Figures, Supplementary Methods, and Supplementary References File Name: Supplementary Movie 1 Description: Footage of time trace of seeds
More informationSupporting Information for. The Origin of Coupled Chloride and Proton Transport in a Cl /H + Antiporter
Supporting Information for The Origin of Coupled Chloride and Proton Transport in a Cl /H + Antiporter Sangyun Lee, Heather B. Mayes, Jessica M. J. Swanson, and Gregory A. Voth * Department of Chemistry,
More informationA self-learning algorithm for biased molecular dynamics
A self-learning algorithm for biased molecular dynamics Tribello, G. A., Ceriotti, M., & Parrinello, M. (2010). A self-learning algorithm for biased molecular dynamics. Proceedings of the National Academy
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 informationSupporting Information
Supporting Information Solid-State 17 O NMR Reveals Hydrogen-Bonding Energetics: Not All Low-Barrier Hydrogen Bonds Are Strong Jiasheng Lu, Ivan Hung, Andreas Brinkmann, Zhehong Gan, Xianqi Kong, and Gang
More informationTesting Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers
computation Article Testing Convergence of Different Free-Energy Methods in a Simple Analytical System with Hidden Barriers S. Alexis Paz 1, ID and Cameron F. Abrams 3, * ID 1 Departamento de Química Teórica
More informationMolecular dynamics (MD) and the Monte Carlo simulation
Escaping free-energy minima Alessandro Laio and Michele Parrinello* Centro Svizzero di Calcolo Scientifico, Via Cantonale, CH-6928 Manno, Switzerland; and Department of Chemistry, Eidgenössische Technische
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 informationMethods for predicting phase transition: VCNEB and TPS approaches
Methods for predicting phase transition: VCNEB and TPS approaches Guangrui Qian Aug 7, 2013 Department of Geosciences, Stony Brook University Materials at High Pressure New physics and new chemistry New
More informationarxiv: v1 [physics.comp-ph] 12 Dec 2014
Accurate multiple time step in biased molecular simulations Marco Jacopo Ferrarotti, Sandro Bottaro, Andrea Pérez-Villa, and Giovanni Bussi Scuola Internazionale Superiore di Studi Avanzati (SISSA), via
More informationSelectivity in the initial C-H bond cleavage of n-butane on PdO(101)
Supporting Information for Selectivity in the initial C-H bond cleavage of n-butane on PdO(101) Can Hakanoglu (a), Feng Zhang (a), Abbin Antony (a), Aravind Asthagiri (b) and Jason F. Weaver (a) * (a)
More informationHow wet should be the reaction coordinate for ligand unbinding?
THE JOURNAL OF CHEMICAL PHYSICS 145, 054113 (2016) How wet should be the reaction coordinate for ligand unbinding? Pratyush Tiwary a) and B. J. Berne b) Department of Chemistry, Columbia University, New
More informationA self-learning algorithm for biased molecular dynamics
A self-learning algorithm for biased molecular dynamics Gareth A. Tribello, Michele Ceriotti and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETHZ Zurich USI-Campus,
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationMolecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments)
Molecular Dynamics Simulations of Fusion Materials: Challenges and Opportunities (Recent Developments) Fei Gao gaofeium@umich.edu Limitations of MD Time scales Length scales (PBC help a lot) Accuracy of
More informationIntroduction to Numerical Simulations and High Performance Computing: From Materials Science to Biochemistry
Introduction to Numerical Simulations and High Performance Computing: From Materials Science to Biochemistry Mauro Boero Institut de Physique et Chimie des Matériaux de Strasbourg University of Strasbourg
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 informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationSupporting Information
Supporting Information ph-responsive self-assembly of polysaccharide through a rugged energy landscape Brian H. Morrow, Gregory F. Payne, and Jana Shen Department of Pharmaceutical Sciences, School of
More informationOn the efficient sampling of pathways in the transition path ensemble
On the efficient sampling of pathways in the transition path ensemble Thijs J. H. Vlugt{ and Berend Smit* Department of Chemical Engineering, University of Amsterdam, Nieuwe Achtergracht 166, 1018 WV Amsterdam,
More informationSupporting Information for. Structural and Chemical Dynamics of Pyridinic Nitrogen. Defects in Graphene
Supporting Information for Structural and Chemical Dynamics of Pyridinic Nitrogen Defects in Graphene Yung-Chang Lin, 1* Po-Yuan Teng, 2 Chao-Hui Yeh, 2 Masanori Koshino, 1 Po-Wen Chiu, 2 Kazu Suenaga
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
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 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 informationDetermination of the Electric Dipole Moment of a Molecule from. Density Functional Theory Calculations
Determination of the Electric Dipole Moment of a Molecule from Density Functional Theory Calculations Byeong June Min Department of Physics, Daegu University, Kyungsan 712-714, Korea Density functional
More informationSupplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization
Supplementary Information for Atomistic Simulation of Spinodal Phase Separation Preceding Polymer Crystallization Richard H. Gee * Naida Lacevic and Laurence E. Fried University of California Lawrence
More informationKinetic Pathways of Ion Pair Dissociation in Water
3706 J. Phys. Chem. B 1999, 103, 3706-3710 Kinetic Pathways of Ion Pair Dissociation in Water Phillip L. Geissler, Christoph Dellago, and David Chandler* Department of Chemistry, UniVersity of California,
More informationEncarna Pucheta-Martinez+, Nicola D Amelio+, Moreno Lelli, Jorge L. Martinez-Torrecuadrada, Marius Sudol, Giorgio Saladino* and Francesco L.
Supplementary Information for: Changes in the folding landscape of the WW domain provide a molecular mechanism for an inherited genetic syndrome Encarna Pucheta-Martinez+, Nicola D Amelio+, Moreno Lelli,
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 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 informationDetermination of Kamlet-Taft parameters for selected solvate ionic liquids.
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Determination of Kamlet-Taft parameters for selected solvate ionic liquids. Daniel
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
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 information