Molecular Mechanics / ReaxFF
|
|
- Noah Henderson
- 6 years ago
- Views:
Transcription
1 Molecular Dynamics simulations Lecture 09: Molecular Mechanics / ReaxFF Dr. Olli Pakarinen University of Helsinki Fall 2012 Lecture notes based on notes by Dr. Jani Kotakoski, 2010 CONTENTS Molecular mechanics (MM) Force fields for molecular interactions Force field approach in MD GROMACS ReaxFF
2 Molecular Mechanics Molecular mechanics (MM) is a method to describe complex molecular systems Structurally there are two differences between MD and MM: The constituent particles can be molecules or parts of molecules rather than single atoms Interactions are typically described via force fields instead of potentials For example, methyl ( CH 3 ) and methylene ( CH 2 ) groups are often modeled as single entities Also, large proteins are sometimes modeled as a set of attached beads Bonds are presented as springs no bond breakings occur However, the prototypical application is simple energy minimization, not time evolution different optimization methods However, many variations exist The example below used Folding@Home for about a million CPU days of simulation time [ [Christopher D. Snow, Houbi Nguyen, Vijay S. Pande and Martin Gruebele, Nature 420, 102 (2002)]
3 Force Fields Molecular interactions become easily complicated because bond energies depend on the overall structure instead of only the local neighborhood U = U bond + U angle + U torsion + U oop + U cross + U nonbond (1) Ubond is the energy related to bond length Uangle is associated with changes in the bond angles Utorsion term including variations coming from the rotation of two parts of the molecule with respect to each other Uoop out-of-plane -term, which affects the planarity of the molecule Unonbond ionic or the van der Waals terms which are not associated with covalent bonding MM2: [J. Am. Chem. Soc. 99, 8127 (1977)], AMBER: [J. Comp. Chem. 24, 1999 (2003)], CHARMm: [J. Comp. Chem. 4, 187 (2004)] From original site removed Some terms included in a force field parametrization
4 BOND TERM The term describing covalent bonding depends on the bond length variation from an ideal value b 0 It is thus a simple pair potential, which can be of Morse or Lennard-Jones type At least until recently these potentials have been too expensive for molecular interactions In fact, the simple harmonic form is often used (AMBER, CHARMm): U bond = 1 2 k b(b b 0 ) 2 (2) bonds Since bond breaking is not allowed, the harmonic approximation is more or less ok for most cases in solids this corresponds to the harmonic regime (Hooke s law) However, it fails for large inter-atomic distances Another often used form is the cubic approximation (MM2) U bond = bonds 1 2 k b(b b 0 ){1 2(b b 0 )} (3) Morse Harmonic Cubic Morse potential is the most realistic of these However, as mentioned, it is sometimes too expensive Bond energy Bond length b 0 Cubic is asymmetric, and in that sense better than harmonic Unfortunately it becomes unphysical at long b
5 ANGLE TERM Also the angle dependency is typically (AMBER, CHARMm) modeled with a harmonic function U angle = bonds 1 2 k θ(θ θ 0 ) 2 (4) A better description is offered by a more expensive sixth-order function (MM2) U angle = k θ (θ θ 0 ) 2 { (θ θ 9 ) 4 } (5) Bond angle bending force constants (k θ ) are proportionately smaller than the bond distortion force constants (k b ) because much lower energies are required for distorting angles TORSIONAL INTERACTIONS Torsional angle rotations are the most important intramolecular terms in a force field τ They differ from bond stretching and bending interactions in two important ways: Internal rotation barriers are low, and thus changes in dihedral angles can be large Torsional potential U torsion is periodic through a 360 rotation
6 The U torsion depends heavily on the atoms forming the molecule Therefore, the functional form must be chosen so that it can model a wide variety of different potentials Commonly, a Fourier series is used for this purpose Different implementations include: AMBER: U torsion = 1 2 k τ{1 + cos(nτ δ)} CHARMm: U torsion = k τ k τ cos(nτ) MM2: U torsion = 1 2 {k τ 1 (1 + cos τ) + k τ2 (1 cos 2τ) + k τ3 (1 + cos 3τ) These functions include different numbers of minima The last one is a sum of functions with 1, 2 and 3 minima OUT-OF-PLANE BENDING INTERACTIONS Out-of-plane interactions are closely related to the torsional terms For four atoms (A,B,C,D) which for three angles (A-B-C, A-B-D, C-B-D), there can be an energy cost for moving one of the atoms out of the plane spanned by the three other atoms This can be accounted for in one of two main ways: Harmonic potential Improper torsion term
7 NON-BONDING INTERACTION Non-bonded interactions act both between atoms in the same molecule and those in other molecules These are often divided to Electrostatic interactions Van der Waals interactions The electrostatic part arises from charge differences within a molecule This can be modeled by distributing point charges in the molecule Of course, the simple approach is to implement an interaction of the form (CHARMm & AMBER) U elect = 1 4πε 0 q 1 q 2 r ij (6) Also a dipole-dipole interaction can be implemented to act between polar groups A functional form is (MM2): U e = µ iµ j ε 0 rij 3 {cos χ 3 cos α i cos α j } (7) The µ i,j are dipole moments, χ is the angle between them. The α i,j are the angles between the moments The van der Waals interaction is understood to include all non-bonding interactions
8 Most typically, this is included via a LJ potential (something similar to this is used in CHARMm and AMBER) [ (σij ) 12 ( ) ] 6 σij U LJ = 4ε ij r ij r ij (8) A better description of the repulsive part would be an exponential function This leads to the Hill/Buckingham potential (MM2) U Hill vdw = A ij exp( B ij r ij ) c ij /r 6 ij (9) As mentioned during the pair potential lecture, this leads to a maximum at finite r ij and thus allows for a Buckingham explosion at close distances Also the Morse potential is occasionally used Comparison of the behavior of LJ and Buckingham potentials Note that this is a schematic, albeit with correct functional forms Buckingham U LJ r
9 ADDITIONAL TERMS The above-listed terms are common for all force fields However, there are typically many additional terms They are often used as coupling terms between the different listed terms, e.g., to join the the bond length and bond angle Also, it s possible to include terms describing electronegativity or conjugative effects Another missing feature is hydrogen bonding, which has been implemented in different ways CHOOSING A FORCE FIELD Obviously, the force field has to be chosen based on an analysis on the accuracy vs. efficiency of the method Some have been implemented for simulating bulk phases (e.g., AMBER) These force fields have generally a simple form with harmonic terms Other force fields, such as the MM-group FFs (MM4: [J. Comput. Chem. 14, 642 (1996)], [J. Comput. Chem. 14, 669 (1996)], [J. Comput. Chem. 14, 695 (1996)], [J. Comput. Chem. 14, 747 (1996)]), are designed to determine structures and vibrational frequencies
10 In this case, more accurate high order polynomials are used for bond stretches and bends and a Hill/Buckingham potential is used for the van der Waals interaction The different force fields are classified as class 1, 2 or 3 force fields, where the class 1 FFs are of the simplest form Another choise to be made is whether all the atoms are to be explicitly described (all-atom force field) The other option is to disregard the unimportant atoms, usually hydrogen (united-atom force field) The neglected atoms are accounted for by increasing the size of the atoms they are bonded to (by increasing the van der Waals radius) This obviously helps greatly in reducing the required number of computations Also entire groups of atoms can be replaced by a single site However, one needs to remember that neglecting the hydrogen atoms will yield a poorer description of electrostatic effects (e.g., charge separation) A hybrid method has also been proposed [J. Mol. Struct. (Theochem) 464, 39 (1999)], in which an all-atom description is used for hydrogen atoms bonded to phenyl rings and a united-atom approach for hydrogen atoms in aliphatic groups
11 PARAMETRIZATION OF A FORCE FIELD For a two-atomic system (here C-H), the following parameters need to be parametrized: 2 parameters for the two different bond types (H-H neglected) 2 for each three different angle types (C-C-C, C-C-H, H-C-H) 3 for each dihedral type (C-C-C-C, H-C-C-C, H-C-C-H) 2 van der Waals parameters for the three different atom combinations (C-C, C-H, H-H) up to 5 different charges So, the simplest possible two-atomic system of any interest includes at least 25 parameters For example the MM2 force field (for 30 atoms) has 3722 parameters ON CURRENT FORCE FIELDS The source: [ is from 2003 The force fields mentioned above are of an intuitive simple form More accurate description have been developed more recently The new ones are fitted largely to experimental data (i.e., heats of formation and vibrational frequencies) A complicated functional form is required for the accurate description Methods based on the MM2 aim at making accurate prediction of molecular structures and properties Two more recent ones are: MM3: [J. Am. Chem. Soc. 111, 8551 (1989)], [J. Am. Chem. Soc. 111, 8566 (1989)], [J. Am. Chem. Soc. page 8576 (1989)] MM4: [J. Comput. Chem. 14, 642 (1996)], [J. Comput. Chem. 14, 669 (1996)], [J. Comput. Chem. 14, 695 (1996)], [J. Comput. Chem. 14, 747 (1996)]
12 The other group of force fields aims at modelling large molecules (polymers/proteins) Typical force fields in this category are AMBER [J. Am. Chem. Soc. 106 (1984)], [J. Am. Chem. Soc. 117, 5179 (1995)] CHARMm [J. Comput. Chem. 4, 1234 (1983)] OPLS [J. Am. Chem. Soc. 118, (1996)] These have a simpler functional form, typically only containing the terms covered during this lecture Often they also utilize the united-atom approach For references to more advanced methods, see the mentioned online material GROMACS GROMACS is a GPL d MD code which is able to use various force fields (Groningen Machine for Chemical Simulations) It is available in many Linux distributions as a standard package Website: Manual: This software is primarily designed for biochemical molecules like proteins, lipids and nucleic acids that have many complicated bonded interactions, with solvent or implicit solvent GROMACS at least claims to be extremely fast at calculating the nonbonded interactions which typically dominate the simulation time The code is parallelized with MPI and can thus be run using several PC computers or several cores of a single PC
13 The actual manual is 370 pages long ([ Manual gives a very good introduction to the method (MD/MM) It even highlights the approximations done in the simulations: The simulations are classical, Electrons are in the ground state, Force fields are approximate, The force field is pair-additive, Long-range interactions are cut off, Boundary conditions are unnatural This manual is instructive to any MD user, even if not using GROMACS The biggest problem in using GROMACS is that no bond breaking is allowed all force fields are non-reactive An example from YouTube: [ ReaxFF ReaxFF is an approach to combine force fields with a bond order approach The name comes from Reactive Force Field This method has been developed by van Duin, Goddard and others to allow chemical reactions within a force field approach ReaxFF aims at being as general as possible, and has been parametrized and tested for hydrocarbon reactions, transition-metal-catalyzed nanotube formation and high-energy materials ReaxFF manual is available at Adri van Duin s website at (updated in 2003) ReaxFF is trying to provide a transferable potential which would be applicable to a wide range of chemical environments
14 To ensure transferability, the following guidelines have been adopted No discontinuities in energy or forces, even during reactions (Same condition as for any other inter-atomic potential!) Each element is described by just one force field atom type No pre-definition of reactive sites (although it is possible to drive reactions with restraints) ReaxFF software has following features implemented NVT and NVE dynamics (Berendsen method) Steepest descent and conjugate gradient minimization methods Optimization of simulation cell parameters Simulations using a crystal unit cells keeping track of the system between periodic images of atoms EEM charge derivation method allowing calculation of geometry dependent charge distributions From the ReaxFF manual, the following parameter sets are available (website update date: July 13, 2003): Name System Quality CH Hydrocarbons *** CH All-carbon (fullerenes, nanotubes) ** RDX Nitramines *** CONSH CONSH-systems, proteins ** SiO Si/SiO 2 clusters and condensed phase *** AlO Al/AlO *** SiN SiN clusters and condensed phase * ZrO Zr and ZrO 2 condensed phase * MoO Mo and Mo-oxides ** Pt Pt bulk metal; Pt-C, Pt-H and Pt-O systems * Quality: * means that parameterization has only been performed for against a fairly limited training set, ** indicates that the parameters have been tested against a reasonably good training set but that further modifications and improvements are expected, *** indicates that these parameters have reached the application-stage.
15 Originally, the method was described in a publication in 2001: [J. Phys. Chem. A 2001, 105, ] From a lecture by the main author found online; Current status of ReaxFF (2008): Code is (apparently) used by 70 research groups around the world List of implementations:
16 The ReaxFF approach is attractive, but due to the huge number of involved parameters, it gets difficult to make sure that the fit works in all circumstances An example from YouTube (reactions): [ According to (now removed) YouTube videos, at least someone has had problems with the Si fit: badly unstable at wrong lattice constant. Also, due to a large number of different terms, ReaxFF is much more demanding than the bond order potentials Another example from YouTube: [ EXAMPLE OF A REAXFF FIT From Development of a ReaxFF description for gold by Tommi Järvi and co-workers [Eur. Phys. J. B 66, (2008)] Gold is not an obvious target for a ReaxFF potential The main motivation of the fit was to allow atomistic simulations of a thiol-gold system A thiol is a organosulfur compound that contains a carbon-bonded sulfhydryl (-C-SH or R-SH) group which are used in many areas of nanotechnology, typically on gold surfaces Of the possible interactions, only three terms are included U system = U bond + U over + U vdw (10) These are supposed to be enough to describe a metal U over is implemented to correct over-coordinations, other terms are self-explanatory
17 APPLICATION FOR THIOL-GOLD DESCRIPTION From Binding of deposited gold clusters to thiol self-assembled monolayers on Au(111) surfaces by Leila Costelle and co-workers [Appl. Phys. Lett. 98(4) (2011)] For this simplified ReaxFF potential, we get the following equations with fitted parameters
18 Properties of different structures as given by the potential Particular goals for the potential were reasonable descriptions of surface and dimer properties For the surface energetics, the potential gives γ (100) : 1.21 J/m 2 γ (111) : 1.07 J/m 2 Experiments give γ: 1.54 J/m 2 Adatom/vacancy energetics turn out to be reasonable (LDA/PBE DFT values in parentheses): E adat(100) : 0.59 (0.30/0.14) ev, E vac(100) : 0.59 (0.40/0.31) ev, E adat(111),fcc : (0.95/0.61) ev, E adat(100),hcp : (1.00/0.64) ev, E adat(100) : 0.87 (0.88/0.61) ev
19 Of all the potentials the ReaxFF parametrization was compared to, it gave clearly the best description of the dimer dissociation energetics This is easy to understand since the other potentials were not fitted to this Taking into account to goal of the potential, the fit seems to perform well SUMMARY Molecular mechanics is a wide term partially overlapping with molecular dynamics biggest differences are that the former is typically more concerned with molecules than individual atoms, and focuses on finding minimum energy configurations or other properties by several types of optimization methods, not necessarily following the time evolution Description of molecules often requires implementing constraints For example, some atoms may be described as a single large particle (united-atom approach) Terms included in force fields include: bond term, angle term, torsion term, out-of-plane term and non-bonding interactions Additionally, cross-terms are often implemented GROMACS is a GPL d software developed for force field MD ReaxFF builds upon the force field approach, but implements continuous functions to allow bond breaking It has a large set of parameters which allows for better fits but may yield problems in cases not included in the fitting procedure
Structural Bioinformatics (C3210) Molecular Mechanics
Structural Bioinformatics (C3210) Molecular Mechanics How to Calculate Energies Calculation of molecular energies is of key importance in protein folding, molecular modelling etc. There are two main computational
More informationMolecular Mechanics. Yohann Moreau. November 26, 2015
Molecular Mechanics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Mechanics, Label RFCT 2015 November 26, 2015 1 / 29 Introduction A so-called Force-Field
More informationT6.2 Molecular Mechanics
T6.2 Molecular Mechanics We have seen that Benson group additivities are capable of giving heats of formation of molecules with accuracies comparable to those of the best ab initio procedures. However,
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
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 informationReactive Empirical Force Fields
Reactive Empirical Force Fields Jason Quenneville jasonq@lanl.gov X-1: Solid Mechanics, EOS and Materials Properties Applied Physics Division Los Alamos National Laboratory Timothy C. Germann, Los Alamos
More informationReactive potentials and applications
1.021, 3.021, 10.333, 22.00 Introduction to Modeling and Simulation Spring 2011 Part I Continuum and particle methods Reactive potentials and applications Lecture 8 Markus J. Buehler Laboratory for Atomistic
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 informationAll-atom Molecular Mechanics. Trent E. Balius AMS 535 / CHE /27/2010
All-atom Molecular Mechanics Trent E. Balius AMS 535 / CHE 535 09/27/2010 Outline Molecular models Molecular mechanics Force Fields Potential energy function functional form parameters and parameterization
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
More informationThe Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii
The Potential Energy Surface (PES) And the Basic Force Field Chem 4021/8021 Video II.iii Fundamental Points About Which to Be Thinking It s clear the PES is useful, so how can I construct it for an arbitrary
More informationMolecular mechanics. classical description of molecules. Marcus Elstner and Tomáš Kubař. April 29, 2016
classical description of molecules April 29, 2016 Chemical bond Conceptual and chemical basis quantum effect solution of the SR numerically expensive (only small molecules can be treated) approximations
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 informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More information3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D.
3rd Advanced in silico Drug Design KFC/ADD Molecular mechanics intro Karel Berka, Ph.D. Martin Lepšík, Ph.D. Pavel Polishchuk, Ph.D. Thierry Langer, Ph.D. Jana Vrbková, Ph.D. UP Olomouc, 23.1.-26.1. 2018
More informationSolutions to Assignment #4 Getting Started with HyperChem
Solutions to Assignment #4 Getting Started with HyperChem 1. This first exercise is meant to familiarize you with the different methods for visualizing molecules available in HyperChem. (a) Create a molecule
More informationMolecular Mechanics, Dynamics & Docking
Molecular Mechanics, Dynamics & Docking Lawrence Hunter, Ph.D. Director, Computational Bioscience Program University of Colorado School of Medicine Larry.Hunter@uchsc.edu http://compbio.uchsc.edu/hunter
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
More informationBiomolecular modeling I
2016, December 6 Biomolecular structure Structural elements of life Biomolecules proteins, nucleic acids, lipids, carbohydrates... Biomolecular structure Biomolecules biomolecular complexes aggregates...
More informationA Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane
A Molecular Dynamics Simulation of a Homogeneous Organic-Inorganic Hybrid Silica Membrane Supplementary Information: Simulation Procedure and Physical Property Analysis Simulation Procedure The molecular
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationThis semester. Books
Models mostly proteins from detailed to more abstract models Some simulation methods This semester Books None necessary for my group and Prof Rarey Molecular Modelling: Principles and Applications Leach,
More informationForce fields in computer simulation of soft nanomaterials
Lecture 2 Part B Force fields in computer simulation of soft nanomaterials Recommended reading: Leach Chapter 4 1 Force Field Methods Force field methods are simulation methods that use classical force
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationMolecular Mechanics. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology. January 2001
Molecular Mechanics C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology January 2001 Introduction Molecular Mechanics uses classical type models to predict the energy
More informationAtomic and molecular interaction forces in biology
Atomic and molecular interaction forces in biology 1 Outline Types of interactions relevant to biology Van der Waals interactions H-bond interactions Some properties of water Hydrophobic effect 2 Types
More informationPair Potentials and Force Calculations
Molecular Dynamics simulations Lecture 04: Pair Potentials and Force Calculations Dr. Olli Pakarinen University of Helsinki Fall 01 Lecture notes based on notes by Dr. Jani Kotakoski, 010 CONSTRUCTING
More informationSubject of the Lecture:
Subject of the Lecture: Conceptual basis for the development of force fields. Implementation/validation Water - a worked example Extensions - combining molecular mechanics and quantum mechanics (QM/MM)
More informationCharge equilibration
Charge equilibration Taylor expansion of energy of atom A @E E A (Q) =E A0 + Q A + 1 @Q A 0 2 Q2 A @ 2 E @Q 2 A 0 +... The corresponding energy of cation/anion and neutral atom E A (+1) = E A0 + @E @Q
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationNanotube AFM Probe Resolution
Influence of Elastic Deformation on Single-Wall Carbon Nanotube AFM Probe Resolution Ian R. Shapiro, Santiago D. Solares, Maria J. Esplandiu, Lawrence A. Wade, William A. Goddard,* and C. Patrick Collier*
More informationForce Fields in Molecular Mechanics
Force Fields in Molecular Mechanics Rajarshi Guha (9915607) and Rajesh Sardar (9915610) March 21, 2001 1 Introduction With the advent of computers chemists have realized the utility of carrying out simulations
More information4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro. Karel Berka, Ph.D.
4 th Advanced in silico Drug Design KFC/ADD Molecular Modelling Intro Karel Berka, Ph.D. UP Olomouc, 21.1.-25.1. 2019 Motto A theory is something nobody believes, except the person who made it An experiment
More informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 09 Quantum Mechanics/Molecular Mechanics (QM/MM) Techniques Fundamentals of Sustainable Technology These notes created by David Keffer, University of Tennessee, Knoxville,
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationScuola di Chimica Computazionale
Societa Chimica Italiana Gruppo Interdivisionale di Chimica Computazionale Scuola di Chimica Computazionale Introduzione, per Esercizi, all Uso del Calcolatore in Chimica Organica e Biologica Modellistica
More informationElectronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani
Technische Universiteit Eindhoven University of Technology Electronic excitations in conjugated polymers Electronic from excitations QM/MM simulations in conjugated pol based from on Many-Body QM/MM simulations
More informationEnergy functions and their relationship to molecular conformation. CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror
Energy functions and their relationship to molecular conformation CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror Yesterday s Nobel Prize: single-particle cryoelectron microscopy 2 Outline Energy
More informationThe broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals.
Physical Metallurgy The broad topic of physical metallurgy provides a basis that links the structure of materials with their properties, focusing primarily on metals. Crystal Binding In our discussions
More informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
More informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationWhy Is Molecular Interaction Important in Our Life
Why Is Molecular Interaction Important in ur Life QuLiS and Graduate School of Science iroshima University http://www.nabit.hiroshima-u.ac.jp/iwatasue/indexe.htm Suehiro Iwata Sept. 29, 2007 Department
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 informationChapter 11 Molecular Mechanics
Chapter 11 Molecular Mechanics Molecular Mechanics uses an analytical, differentiable, and relatively simple potential energy function, (R), for describing the interactions between a set of atoms specified
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationEnergy functions and their relationship to molecular conformation. CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror
Energy functions and their relationship to molecular conformation CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror Outline Energy functions for proteins (or biomolecular systems more generally)
More informationIntroduction to Molecular Dynamics
Introduction to Molecular Dynamics Dr. Kasra Momeni www.knanosys.com Overview of the MD Classical Dynamics Outline Basics and Terminology Pairwise interacting objects Interatomic potentials (short-range
More informationMolecular Simulation II. Classical Mechanical Treatment
Molecular Simulation II Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationSupporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes
Supporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes Details of the molecular dynamics simulations and computational procedures
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationForce Fields for Classical Molecular Dynamics simulations of Biomolecules. Emad Tajkhorshid
Force Fields for Classical Molecular Dynamics simulations of Biomolecules Emad Tajkhorshid Theoretical and Computational Biophysics Group, Beckman Institute Departments of Biochemistry and Pharmacology,
More informationAn Informal AMBER Small Molecule Force Field :
An Informal AMBER Small Molecule Force Field : parm@frosst Credit Christopher Bayly (1992-2010) initiated, contributed and lead the efforts Daniel McKay (1997-2010) and Jean-François Truchon (2002-2010)
More informationIntroduction to ReaxFF: Reactive Molecular Dynamics
Introduction to ReaxFF: Reactive Molecular Dynamics Ole Carstensen carstensen@scm.com TCCM ADF Tutorial April 21 Amsterdam Outline ReaxFF - general aspects Molecular Dynamics Intro 200 DFT 100 ReaxFF Harmonic
More informationMOLECULAR DYNAMICS SIMULATIONS OF HEAT TRANSFER ISSUES IN CARBON NANOTUBES
The st International Symposium on Micro & Nano Technology, 4-7 March, 4, Honolulu, Hawaii, USA MOLECULAR DYNAMICS SIMULATIONS OF HEAT TRANSFER ISSUES IN CARBON NANOTUBES S. Maruyama, Y. Igarashi, Y. Taniguchi
More informationThe split-charge method:
Split charge equilibration: A charge transfer method for electrolytes and other non-metallic materials The split-charge method: Martin H. Müser Materialwissenschaften und Werkstoffwissenschaften Motivation
More informationSemi Empirical Force Fields and Their Limitations. Potential Energy Surface (PES)
Semi Empirical Force Fields and Their Limitations Ioan Kosztin Beckman Institute University of Illinois at Urbana-Champaign Potential Energy Surface (PES) Schrödinger equation: H T Ψ( r, = E Ψ( r, H =
More informationJASS Modeling and visualization of molecular dynamic processes
JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St Petersburg State University, Physics faculty, Department of Computational Physics Supervisor PhD Stepanova Margarita
More informationWater models in classical simulations
Water models in classical simulations Maria Fyta Institut für Computerphysik, Universität Stuttgart Stuttgart, Germany Water transparent, odorless, tasteless and ubiquitous really simple: two H atoms attached
More informationChemistry Review: Atoms
Chemistry Review: Atoms Atoms are made up : nucleus containing protons and neutrons orbitals containing electrons (2, 8, 8,...). Valence electrons outermost electrons Chemistry Review: Atoms All atoms
More informationLecture 1. Conformational Analysis in Acyclic Systems
Lecture 1 Conformational Analysis in Acyclic Systems Learning Outcomes: by the end of this lecture and after answering the associated problems, you will be able to: 1. use Newman and saw-horse projections
More informationIonic Bonding - Electrostatic Interactions and Polarization
Ionic Bonding - Electrostatic Interactions and Polarization Chemistry 754 Solid State Chemistry Dr. Patrick Woodward Lecture #13 Born-Haber Cycle for NaCl It is energetically unfavorable for Na metal and
More informationMolecular dynamics simulation. CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror
Molecular dynamics simulation CS/CME/BioE/Biophys/BMI 279 Oct. 5 and 10, 2017 Ron Dror 1 Outline Molecular dynamics (MD): The basic idea Equations of motion Key properties of MD simulations Sample applications
More informationFrom Atoms to Materials: Predictive Theory and Simulations
From Atoms to Materials: Predictive Theory and Simulations Week 3 Lecture 4 Potentials for metals and semiconductors Ale Strachan strachan@purdue.edu School of Materials Engineering & Birck anotechnology
More informationReactive Force Field & Molecular Dynamics Simulations (Theory & Applications)
Reactive Force Field & Molecular Dynamics Simulations (Theory & Applications) Ying Li Collaboratory for Advanced Computing & Simulations Department of Chemical Engineering & Materials Science Department
More informationPH 548 Atomistic Simulation Techniques
PH 548 Atomistic Simulation Techniques Lectures: Lab: 2-0-2-6 Tuesday 12-1 PM Thursday 12-1 PM Monday 2-5 PM P. K. Padmanabhan PH548: Two Excellent Books How to do it? 1 2 F 12 F 23 3 F i = F ij j i F
More informationDownloaded from
I.I.T.Foundation - XI Chemistry MCQ #4 Time: 45 min Student's Name: Roll No.: Full Marks: 90 Chemical Bonding I. MCQ - Choose Appropriate Alternative 1. The energy required to break a chemical bond to
More informationLecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions. Let's get together.
Lecture C2 Microscopic to Macroscopic, Part 2: Intermolecular Interactions Let's get together. Most gases are NOT ideal except at very low pressures: Z=1 for ideal gases Intermolecular interactions come
More informationUnit 6: Molecular Geometry
Unit 6: Molecular Geometry Molecular Geometry [6-5] the polarity of each bond, along with the geometry of the molecule determines Molecular Polarity. To predict the geometries of more complicated molecules,
More informationChapter 3. Crystal Binding
Chapter 3. Crystal Binding Energy of a crystal and crystal binding Cohesive energy of Molecular crystals Ionic crystals Metallic crystals Elasticity What causes matter to exist in three different forms?
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 Figure 1. Schematic of rapid thermal annealing process: (a) indicates schematics and SEM cross-section of the initial layer-by-layer
Supplementary Figure 1. Schematic of rapid thermal annealing process: (a) indicates schematics and SEM cross-section of the initial layer-by-layer film configuration, (b) demonstrates schematic and cross-section
More informationWhat factors affect whether something is a solid, liquid or gas? What actually happens (breaks) when you melt various types of solids?
States of Mattter What factors affect whether something is a solid, liquid or gas? What actually happens (breaks) when you melt various types of solids? What external factors affect whether something is
More informationModeling and visualization of molecular dynamic processes
JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St-Petersburg State University Physics faculty Department of Computational Physics Supervisor PhD Stepanova Margarita
More informationLecture: P1_Wk1_L1 IntraMolecular Interactions. Ron Reifenberger Birck Nanotechnology Center Purdue University 2012
Lecture: IntraMolecular Interactions Distinguish between IntraMolecular (within a molecule) and InterMolecular (between molecules) Ron Reifenberger Birck Nanotechnology Center Purdue University 2012 1
More informationHomology modeling. Dinesh Gupta ICGEB, New Delhi 1/27/2010 5:59 PM
Homology modeling Dinesh Gupta ICGEB, New Delhi Protein structure prediction Methods: Homology (comparative) modelling Threading Ab-initio Protein Homology modeling Homology modeling is an extrapolation
More informationGeneral Physical Chemistry II
General Physical Chemistry II Lecture 13 Aleksey Kocherzhenko October 16, 2014" Last time " The Hückel method" Ø Used to study π systems of conjugated molecules" Ø π orbitals are treated separately from
More informationAdvanced Molecular Molecular Dynamics
Advanced Molecular Molecular Dynamics Technical details May 12, 2014 Integration of harmonic oscillator r m period = 2 k k and the temperature T determine the sampling of x (here T is related with v 0
More informationUnderstanding Chemistry in 3D
Contents 1 Understanding Chemistry in 3D 2 Organic Chemistry 2.1 So how do cells keep from falling apart in a water environment? 2.2 How are macromolecules assembled? 3 Where do we find macromolecules?
More informationClose agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations
Close agreement between the orientation dependence of hydrogen bonds observed in protein structures and quantum mechanical calculations Alexandre V. Morozov, Tanja Kortemme, Kiril Tsemekhman, David Baker
More informationLecture 2: Bonding in solids
Lecture 2: Bonding in solids Electronegativity Van Arkel-Ketalaar Triangles Atomic and ionic radii Band theory of solids Molecules vs. solids Band structures Analysis of chemical bonds in Reciprocal space
More informationwith the larger dimerization energy also exhibits the larger structural changes.
A7. Looking at the image and table provided below, it is apparent that the monomer and dimer are structurally almost identical. Although angular and dihedral data were not included, these data are also
More informationChemistry: The Central Science
Chemistry: The Central Science Fourteenth Edition Chapter 8 Basic Concepts of Chemical Bonding Chemical Bonds Three basic types of bonds Ionic Electrostatic attraction between ions Covalent Sharing of
More informationChemical bonding in solids from ab-initio Calculations
Chemical bonding in solids from ab-initio Calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology, University
More informationChapter 8. Basic Concepts of Chemical Bonding
Chapter 8 Basic Concepts of Chemical Bonding Chemical Bonds An attractive force that holds two atoms together in a more complex unit Three basic types of bonds Ionic Electrons are transferred from one
More informationBiophysics II. Hydrophobic Bio-molecules. Key points to be covered. Molecular Interactions in Bio-molecular Structures - van der Waals Interaction
Biophysics II Key points to be covered By A/Prof. Xiang Yang Liu Biophysics & Micro/nanostructures Lab Department of Physics, NUS 1. van der Waals Interaction 2. Hydrogen bond 3. Hydrophilic vs hydrophobic
More informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More informationCS 2, HCN, BeF 2 Trigonal planar. Cl 120 BF 3, AlCl 3, SO 3, NO 3-, CO NCl 3,PF 3,ClO 3,H 3 O + ...
Shape of molecules Name No bonding pairs No lone pairs Diagram Bond angle Examples linear 2 0 l Be l 180 2, S 2, N, Be 2 Trigonal planar 3 0 l l 120 B 3, All 3, S 3, N 3-, 2-3 B Tetrahedral 4 0 109.5 Sil
More informationTOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS
TOPOLOGIES AND FORCE FIELD PARAMETERS FOR NITROXIDE SPIN LABELS 1. Introduction The force field used in these studies is the CHARMM19 (Chemistry at HARvard Molecular Mechanics) extended atom force field
More informationMolecular Simulation II
Molecular Simulation II Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationParameterization of a reactive force field using a Monte Carlo algorithm
Parameterization of a reactive force field using a Monte Carlo algorithm Eldhose Iype (e.iype@tue.nl) November 19, 2015 Where innovation starts Thermochemical energy storage 2/1 MgSO 4.xH 2 O+Q MgSO 4
More information3. An Introduction to Molecular Mechanics
3. An Introduction to Molecular Mechanics Introduction When you use Chem3D to draw molecules, the program assigns bond lengths and bond angles based on experimental data. The program does not contain real
More informationElectronic structure and transport in silicon nanostructures with non-ideal bonding environments
Purdue University Purdue e-pubs Other Nanotechnology Publications Birck Nanotechnology Center 9-15-2008 Electronic structure and transport in silicon nanostructures with non-ideal bonding environments
More information1. Reactions can be followed by measuring changes in concentration, mass and volume of reactants and products.
Higher Chemistry - Traffic Lights Unit 1 CHEMICAL CHANGES AND STRUCTURE I know: Controlling the rate Collision theory and relative rates 1. Reactions can be followed by measuring changes in concentration,
More information3. An Introduction to Molecular Mechanics
3. An Introduction to Molecular Mechanics Introduction When you use Chem3D to draw molecules, the program assigns bond lengths and bond angles based on experimental data. The program does not contain real
More informationSolutions and Intermolecular Forces
Solutions and Intermolecular Forces REVIEW Chemical Bonds Three basic types of bonds: Ionic Electrostatic attraction between ions Covalent Sharing of electrons Metallic Metal atoms bonded to several other
More information