Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
|
|
- Kerry Hines
- 6 years ago
- Views:
Transcription
1 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 techniques. Without giving detailed mathematical descriptions, briefly outline the advantages and disadvantages of these two techniques. Outline Answer: The steepest descents method of energy minimization is a simple and effective method of improving the geometry of a molecular model (i.e. finding a lower energy structure). It moves straight downhill on the energy surface, taking steps in the direction of the negative of the gradient. This involves calculating the forces (gradients) on the atoms - the first derivative of the energy). The gradient is a vector (g), with one component for each coordinate. The step, s, is in the direction g. The steepest descents step is determined mostly by the largest forces. It is therefore good for relieving strain in a built model. Only first derivatives are required - this is an advantage for large molecules. Steepest descents provides fast minimization (big steps) when far from minimum, but small steps (slow minimization) when close to minimum. It has slow convergence (i.e. it will find the minimum energy structure only very slowly). It may also oscillate around a minimum and never actually reach the minimum. Because it only moves downhill on the energy surface, it is no use for TS searches. The Newton-Raphson method uses in addition second derivatives of the energy and can be a better minimization method. It is efficient for small molecules, and converges quickly, in contrast to steepest descents minimization. However, it makes the approximation that the energy surface has a quadratic form, which is often poor approximation for molecular potential energy surfaces, particularly far from minimum The calculation, inversion and storage of the Hessian matrix (i.e. the matrix of second derivatives) is computationally difficult for large molecules. It is possible to locate stationary points other than minima on potential energy surfaces with the Newton- Raphson method, and this can be useful for transition state optimization. The Newton-Raphson method works well if the initial geometry is close to a stationary point, but may not work well if it is not. Also, technically the method may be difficult to use for some types of problem. Technical difficulties arise because of the need to calculate the Hessian, which is a potentially very large matrix (no. of variables squared). The calculation of the second derivatives can be demanding (requiring a lot of computer time), and their storage can be difficult (a lot of computer memory is required). The Hessian must be inverted (diagonalized), which is a technical problem for more than 1000 variables ( matrix), i.e. approx. 300 atoms. It is possible to avoid calculating the Hessian matrix at every step by either using a guessed matrix (e.g. the identity matrix) or by using a Hessian update scheme (constructs an approximate Hessian based on e.g. the gradient and coordinates at current and new points. Calculating the Hessian at a lower level of theory can also reduce the demands of the calculation. When dealing with large numbers of atoms, though, the Newton-Raphson method is not suitable because of the large computational requirements. A first-order method has to be used instead (first-order methods use first derivatives of the energy; second-order methods use first and second derivatives). 2) Question: Give reasons why calculation and analysis of the harmonic vibrational frequencies of molecules may be useful. Outline the process involved in calculation of harmonic vibrational frequencies.
2 There are several reasons for calculating vibrational frequencies of molecules, including: To identify the nature of a stationary point found by geometry optimization (e.g. a minimum or a saddlepoint). At a minimum, all the 3N-6 eigenvalues of the Hessian are positive, so all the calculated frequencies will be positive. To calculate the infrared and Raman spectra of a molecule To calculate thermodynamic properties of a molecule Vibrational frequencies are usually calculated by finding the normal modes of the molecule. 3) Question: Give an expression for a molecular mechanics potential energy function of the type used for standard macromolecular simulations. Include intra- and intermolecular terms. State clearly what each energy term represents, and its physical origin. How, typically, would the parameters of such a potential energy function be found? Discuss the range of applicability of a potential energy function of this type, its strengths and its limitations. Give brief examples of applications. What improvements could be made to a potential energy function of this type? Answer: The answer should discuss the harmonic function used to model bond stretching, the limitations of the harmonic potential (e.g. no dissociation; poor behaviour away from minimum); the harmonic valence angle bending; the cosine series for torsion angles, mentioning in each case the meaning of the symbols (i.e. force constants, barrier heights, phase factors), and the need for improper torsion terms. The final terms should be identified as the non-bonded interactions in the potential function, van der Waals and electrostatic terms. The answer should discuss which non-bonded interactions would be calculated (i.e. 1-4 interactions etc). In discussing the van der Waals terms, the repulsive and attractive interactions should be identified, and the cause (dispersion interactions and exchange repulsion) should be outlined. Atomic partial charges should be described and the ability of these models to reproduce hydrogen bond structures. Limitations of molecular mechanics force fields such as incomplete and potentially incorrect description of electrostatic interactions, lack of electronic polarizability (fixed atomic charges), inability to model reactions (incorrect functional form), dependence on parameterization (e.g. limited dataset), should be discussed. Possible improvements include: inclusion of atomic polarizability, different forms for bonding interactions; and combined quantum mechanics/molecular mechanics method (outline briefly) the latter allow modelling of reactions. 4) Question i. Write down the Lennard-Jones potential energy function typically used to model van der Waals interactions in molecular modelling of proteins, in terms of the collision diameter,, well-depth,, and interatomic separation, r.
3 ii. State briefly the meaning of and. iii. Sketch the form of this energy equation, labelling the axes as energy and interatomic separation. Indicate and on the sketch. Outline answer: i. V r ii. σ 4ε r 12 σ r This function contains two parameters: the well depth, (the depth of the energy minimum for the interaction), and the collision diameter, (the separation at which the interaction energy is zero). 6 iii. Sketch the curve: 5) What types of points on a potential energy surface are particularly relevant in understanding a chemical reaction? Transition state structures are saddlepoints on the potential energy surface; stable structures are minima on the potential energy surface. 6) Give reasons why calculation and analysis of the harmonic vibrational frequencies of molecules may be useful. Outline the process involved in calculation of harmonic vibrational frequencies. Identify the type of stationary point found by geometry optimization (e.g. is it a minimum or a saddlepoint?) IR and Raman spectra of a molecule Calculate zero-point energy important for some reactions Thermodynamic properties (e.g. vibrational energy contribution to energy barrier) Contribution of quantum tunnelling to a reaction barrier
4 7) Describe how you could confirm that a structure is a transition state structure. Also outline how you could show that the transition state really is the right one for the reaction of interest. Outline answer: Transition state structures are saddlepoints They can be transition states for chemical changes, or for conformational changes For a transition state structure (saddlepoint), one of the vibrational eigenvalues will be negative - it is a maximum in one direction In all other directions, it is a minimum The frequency is given by the square root of the eigenvalue, so formally there will be one imaginary frequency If there is one, and only one, imaginary frequency, the structure really is a transition state The atomic displacements for the imaginary frequency are very important - this normal mode corresponds to motion along the reaction coordinate For a TS structure, the atomic displacements for the imaginary frequency show which atoms are involved are involved in motion along the reaction coordinate A more detailed check is to calculate the steepest descent path from the TS: it should lead (in one direction) to the reactants and (and in the other direction) to the products. This path (calculated in mass-weighted coordinates from the TS) is called the intrinsic reaction coordinate 8) Briefly describe the following methods for calculating molecular energies and geometries (outline their advantages and disadvantages): (i) ab initio molecular orbital methods (ii) semiempirical molecular orbital methods (iii) molecular mechanics methods (iv) combined quantum mechanics/molecular mechanics methods Outline answer: Describe each type of method briefly. The answer should mention e.g. that semiempirical methods are more approximate than ab initio methods; both are quantum mechanical electronic structure methods; and that both can be applied to model chemical reactions unlike molecular mechanics. QM/MM methods combine QM and approaches and thus allow modelling of reactions in large systems such as enzymes. 9) Discuss the applicability and accuracy of the following methods to calculation of barriers for (a) conformational changes and (b) chemical reactions: molecular mechanics, Hartree-Fock ab initio molecular orbital calculations, semiempirical QM calculations; MP2 ab initio calculations. (e.g. for each method, discuss whether it is likely to give calculated rate constants close to experiment for conformational changes, and (separately) for chemical reactions. Outline answer: Molecular mechanics methods can calculate conformational properties, including energy barriers, quite well if suitably parameterized. They are not useful for modelling chemical reactions: for example, the use of a harmonic term for bonds means that they cannot break. Hartree-Fock ab initio molecular orbital methods can deal well with conformational barriers. This is because correlation energy does not change much in a conformational change because chemical bonding does not change. In a chemical reaction, on the other hand, the correlation energy changes significantly, so Hartree-Fock calculations are not very useful for chemical reactions: they tend to give barriers which are much too high.
5 Semiempirical molecular orbital methods can be applied to model conformational or chemical changes, but are more approximate than ab initio methods and so are prone to large errors. MP2 ab initio methods improve on Hartree-Fock by including electron correlation effects. Rate constants (for chemical or conformational changes) can be calculated using the Eyring equation. MP2 calculations would be expected to give rate constants close to experimental values. The other electronic structure (QM) methods would be likely to give results in poorer agreement with experiment. 10) Outline the essential features the TIP3P model used for simulations of liquid water. Discuss its strengths and weaknesses. Outline the features of more sophisticated water models for molecular simulations. Outline answer: TIP3P is a standard model of water for atomistic simulations. It is an invariant point charge model. Each atom has a fixed charge (q(o) = 0.834, (H)=0.417)). There are Lennard-Jones parameters for the oxygen only. Fixed geometry (based on experimental structure, O-H distance = nm, angle H-O- Hw = º. Dipole moment of TIP3P (2.35D) is higher than the experimental result for a water molecule in the gas phase (1.85D). Closer to effective dipole of liquid water (2.6D). Interaction energy for TIP3P water dimer is 6.5 kcal/mol, higher than experiment ( 5.4 ± 0.7 kcal/mol). Liquid water is polarized, each molecule by its neighbours. TIP3P was developed by fitting to experimental properties of liquid water. It is a poor model for small numbers of water molecules, but reasonable for modelling liquid water. Better models can include more point charges to give a better description of the electrostatic interactions of water (e.g. TIP4P, TIP5P). More sophisticated still are models that include electronic polarization explicitly, i.e. polarizable water models. 10) Classical molecular dynamics simulations are an important technique in molecular modelling (e.g. of proteins). Describe the essential features of such molecular dynamics simulations. Outline answer: Classical MD determines atomic positions through solving Newton's equations of motion. Outline integration algorithms for MD, typical timesteps, use of SHAKE algorithm for increasing timestep. A molecular dynamics simulation requires the definition of a potential function, or a description of the terms by which the particles in the simulation will interact. This is usually referred to as a force field. Potentials may be defined at many levels of physical accuracy; those most commonly used in chemistry are based on molecular mechanics and embody a classical treatment of particle-particle interactions that can reproduce structural and conformational changes but usually cannot treat chemical reactions. The reduction from a fully quantum description to a classical potential entails two main approximations. The first one is the Born-Oppenheimer approximation, which states that the dynamics of electrons is so fast that they can be considered to react instantaneously to the motion of their nuclei. As a consequence, they may be treated separately. The second one treats the nuclei, which are much heavier than electrons, as point particles that follow classical Newtonian dynamics. In classical molecular dynamics the effect of the electrons is approximated as a single potential energy surface, usually representing the ground state. When treatment of electronic properties is required, potentials based on quantum mechanics are used; some techniques attempt to create hybrid classical/quantum ( QM/MM ) potentials where the bulk of the system is treated classically but a small region is treated as a quantum system, e.g. for modelling a chemical transformation.
6 11) Hartree-Fock ab initio calculations give an energy of Hartree for ethane in its staggered conformation, and Hartree for the eclipsed form. From these results, calculate the barrier to rotation in ethane in kj/mol, given that 1 H is equivalent to kj/mol. Comment on the fact that the barrier obtained by a Hartree-Fock calculation is found to be close to the experimental result. Outline answer: Energy barrier = H = 13.1 kj/mol. A Hartree-Fock calculation gives a good answer in this case because the electron correlation energy does not change significantly during rotation around the C-C bond. Hartree-Fock calculations do not include the effects of electron correlation, and so overestimate electron-electron repulsion. Therefore the absolute energy calculated by a Hartree-Fock method is too high. However, for processes of conformational change, the correlation energy remains roughly constant, so HF calculations can give a useful prediction of relative energies. This is not the case for processes involving breaking or making of bonds.
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 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 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 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 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 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 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 informationComputational Chemistry. An Introduction to Molecular Dynamic Simulations
Computational Chemistry An Introduction to Molecular Dynamic Simulations Computational chemistry simulates chemical structures and reactions numerically, based in full or in part on the fundamental laws
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 informationComputational Modeling of Protein-Ligand Interactions
Computational Modeling of Protein-Ligand Interactions Steven R. Gwaltney Department of Chemistry Mississippi State University Mississippi State, MS 39762 Auguste Comte, 1830 Every attempt to refer chemical
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 informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
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 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 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 informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
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 informationFigure 1: Transition State, Saddle Point, Reaction Pathway
Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces
More informationIntroduction to Computational Chemistry Computational (chemistry education) and/or. (Computational chemistry) education
Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools to help increase student understanding of material
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 informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
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 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 informationIntroduction to Geometry Optimization. Computational Chemistry lab 2009
Introduction to Geometry Optimization Computational Chemistry lab 9 Determination of the molecule configuration H H Diatomic molecule determine the interatomic distance H O H Triatomic molecule determine
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 informationTECHNIQUES TO LOCATE A TRANSITION STATE
32 III. TECHNIQUES TO LOCATE A TRANSITION STATE In addition to the location of minima, the calculation of transition states and reaction pathways is an interesting task for Quantum Chemistry (QC). The
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 informationChem 1140; Molecular Modeling
P. Wipf 1 Chem 1140 $E = -! (q a + q b )" ab S ab + ab!k
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 information( ) 2 75( ) 3
Chemistry 380.37 Dr. Jean M. Standard Homework Problem Set 3 Solutions 1. The part of a particular MM3-like force field that describes stretching energy for an O-H single bond is given by the following
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
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 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 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 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 informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationUsing Molecular Dynamics to Compute Properties CHEM 430
Using Molecular Dynamics to Compute Properties CHEM 43 Heat Capacity and Energy Fluctuations Running an MD Simulation Equilibration Phase Before data-collection and results can be analyzed the system
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 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 information( ) + ( ) + ( ) = 0.00
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 32, April 14, 2006 (Some material in this lecture has been adapted from Cramer, C.
More informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More information16.1 Molecular Vibrations
16.1 Molecular Vibrations molecular degrees of freedom are used to predict the number of vibrational modes vibrations occur as coordinated movement among many nuclei the harmonic oscillator approximation
More informationIntroduction to Computational Chemistry for Experimental Chemists... (Part 2/2)
12 th PhD seminar, Garching, October 31 st 2008 Introduction to Computational Chemistry for Experimental Chemists... (Part 2/2) Dr. Markus Drees, TU München Universität Regensburg Universität Augsburg
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 informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationCHEMISTRY 4021/8021 MIDTERM EXAM 1 SPRING 2014
CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy
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 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 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 informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
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 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 informationNH 3 inversion: Potential energy surfaces and transition states CH342L March 28, 2016
N 3 inversion: Potential energy surfaces and transition states C342L March 28, 2016 Last week, we used the IR spectrum of ammonia to determine the splitting of energy levels due to inversion of the umbrella
More informationDispersion relation for transverse waves in a linear chain of particles
Dispersion relation for transverse waves in a linear chain of particles V. I. Repchenkov* It is difficult to overestimate the importance that have for the development of science the simplest physical and
More informationSame idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)
XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low
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 informationBiomolecules are dynamic no single structure is a perfect model
Molecular Dynamics Simulations of Biomolecules References: A. R. Leach Molecular Modeling Principles and Applications Prentice Hall, 2001. M. P. Allen and D. J. Tildesley "Computer Simulation of Liquids",
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 informationMulti-scale modelling of transition metal enzymes
Degree Project C in Chemistry, 15 c Multi-scale modelling of transition metal enzymes Nathalie Proos Vedin 18th June 2015 Supervisor: Marcus Lundberg Uppsala University Department of Chemistry Ångström
More informationP. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION
P. W. Atkins and R. S. Friedman Molecular Quantum Mechanics THIRD EDITION Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1997 Introduction and orientation 1 Black-body radiation 1 Heat capacities 2 The
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationFree Energy Simulation Methods
Free Energy Simulation Methods Free energy simulation methods Many methods have been developed to compute (relative) free energies on the basis of statistical mechanics Free energy perturbation Thermodynamic
More informationComputational Material Science Part II
Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Aim of Part II Get familiar with the computational methodologies often used and properties often predicted in
More informationR BC. reaction coordinate or reaction progress R. 5) 8pts) (a) Which of the following molecules would give an infrared spectrum? HCl O 2 H 2 O CO 2
Physical Chemistry Spring 2006, Prof. Shattuck Final Name Part Ia. Answer 4 (four) of the first 5 (five) questions. If you answer more than 4, cross out the one you wish not to be graded. 1) 8pts) Of absorption
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationInteratomic Potentials. The electronic-structure problem
Interatomic Potentials Before we can start a simulation, we need the model! Interaction between atoms and molecules is determined by quantum mechanics: Schrödinger Equation + Born-Oppenheimer approximation
More informationGeometry Optimisation
Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background
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 informationDielectric polarization of 2-pyrrolidinone molecules in benzene solution - a quantum-chemical study
Dielectric polarization of 2-pyrrolidinone molecules in benzene solution - a quantum-chemical study L. Gorb* ), J. Jadżyn $) and K. W. Wojciechowski #) Institute of Molecular Physics, Polish Academy of
More informationOutline for Today. Monday, Nov. 26. Chapter 11: Intermolecular Forces and Liquids. Intermolecular Foces. Comparing States of Matter
Outline for Today Monday, Nov. 26 Chapter 11: Intermolecular Forces and Liquids Intermolecular Foces Comparing States of Matter Properties of Liquids Phase Changes!1 Exam 3 Reminders When: Wednesday at
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
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 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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationHomework Problem Set 1 Solutions
Chemistry 380.37 Dr. Jean M. Standard omework Problem Set 1 Solutions 1. A student investigates a bond between atoms A and B in a molecule using a software package for molecular mechanics. The student
More informationCHEM 344 Molecular Modeling
CHEM 344 Molecular Modeling The Use of Computational Chemistry to Support Experimental Organic Chemistry Part 1: Molecular Orbital Theory, Hybridization, & Formal Charge * all calculation data obtained
More informationMolecular Symmetry 10/25/2018
Symmetry helps us understand molecular structure, some chemical properties, and characteristics of physical properties (spectroscopy). Predict IR spectra or Interpret UV-Vis spectra Predict optical activity
More information74 these states cannot be reliably obtained from experiments. In addition, the barriers between the local minima can also not be obtained reliably fro
73 Chapter 5 Development of Adiabatic Force Field for Polyvinyl Chloride (PVC) and Chlorinated PVC (CPVC) 5.1 Introduction Chlorinated polyvinyl chloride has become an important specialty polymer due to
More informationINPUT DESCRIPTION FOR SQM version 2.0
INPUT DESCRIPTION FOR SQM version 2.0 INTRODUCTION SQM is an add-on module for the PQS program which scales force constants to produce a Scaled Quantum Mechanical (SQM) Force Field. This can correct for
More informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer. Lecture 10, February 10, / 4
Chem 350/450 Physical Chemistry II (Quantum Mechanics 3 Credits Spring Semester 006 Christopher J. Cramer Lecture 10, February 10, 006 Solved Homework We are asked to find and for the first two
More informationIntermolecular Forces
Intermolecular Forces Molecular Compounds The simplest molecule is H 2 : Increased electron density draws nuclei together The pair of shared electrons constitutes a covalent bond. Intermolecular Forces
More informationTools for QM studies of large systems
Tools for QM studies of large systems Automated, hessian-free saddle point search & characterization QM/MM implementation for zeolites Shaama Mallikarjun Sharada Advisors: Prof. Alexis T Bell, Prof. Martin
More informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationNew Perspective on structure and bonding in water using XAS and XRS
New Perspective on structure and bonding in water using XAS and XRS Anders Nilsson Stanford Synchrotron Radiation Laboratory (SSRL) and Stockholm University, Sweden R. Ludwig Angew. Chem. 40, 1808 (2001)
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 informationTHE VIBRATIONAL SPECTRA OF A POLYATOMIC MOLECULE (Revised 3/27/2006)
THE VIBRATIONAL SPECTRA OF A POLYATOMIC MOLECULE (Revised 3/27/2006) 1) INTRODUCTION The vibrational motion of a molecule is quantized and the resulting energy level spacings give rise to transitions in
More informationIdeal Gas Laws Empirical Gas Laws The Mole Equations of State Dalton's Law The Mole Fraction Extensive and Intensive Variables Graham's Law of
Ideal Gas Laws Empirical Gas Laws The Mole Equations of State Dalton's Law The Mole Fraction Extensive and Intensive Variables Graham's Law of Effusion The Maxwell-Boltzmann Distribution A Digression on
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 informationChem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester Christopher J. Cramer. Lecture 30, April 10, 2006
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 20056 Christopher J. Cramer Lecture 30, April 10, 2006 Solved Homework The guess MO occupied coefficients were Occupied
More informationChimica Farmaceutica
Chimica Farmaceutica Drug Targets Why should chemicals, some of which have remarkably simple structures, have such an important effect «in such a complicated and large structure as a human being? The answer
More informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationSection 8.1 The Covalent Bond
Section 8.1 The Covalent Bond Apply the octet rule to atoms that form covalent bonds. Describe the formation of single, double, and triple covalent bonds. Contrast sigma and pi bonds. Relate the strength
More informationIntroduction to Molecular Vibrations and Infrared Spectroscopy
hemistry 362 Spring 2017 Dr. Jean M. Standard February 15, 2017 Introduction to Molecular Vibrations and Infrared Spectroscopy Vibrational Modes For a molecule with N atoms, the number of vibrational modes
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 informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 3 The Born-Oppenheimer approximation C.-K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the Born-Oppenheimer
More informationComputer Simulation of Shock Waves in Condensed Matter. Matthew R. Farrow 2 November 2007
Computer Simulation of Shock Waves in Condensed Matter Matthew R. Farrow 2 November 2007 Outline of talk Shock wave theory Results Conclusion Computer simulation of shock waves Shock Wave Theory Shock
More information