Enhanced Free Energy Based Structure Prediction in Materials Science
|
|
- Myra Barton
- 6 years ago
- Views:
Transcription
1 Enhanced Free Energy Based Structure Prediction in Materials Science Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY NYU-ECNU Center for Computational Chemistry at NYU Shanghai , China 纽约大学 - 华东师范大学计算化学中心在上海纽约大学 NEW YORK UNIVERSITY MRSEC
2 Crystallization and the Pharmaceutical Industry: Materials Science? Pharmaceutical Materials Science: An Active New Frontier in Materials Research J. Elliott and B. Hancock Guest Editors MRS Bulletin 31, Issue 11 (2006) Pharmaceutical Industry $500 billion global market $250 billion US market Majority of pharmaceuticals are small molecule crystals One of the continuing scandals in the physical sciences is that it remains impossible to predict the structure of even the simplest crystalline solids from their chemical composition. -- John Maddox, Nature (1998)
3 Targets of the sixth blind structure prediction test
4 Accuracy of intermolecular interactions Standard force fields qiq j aij b ij U( R) U 12 6 bonded ( R) ij rij rij rij Kohn-Sham density functional theory (plus dispersion): 1 1 n( r) n( r) U d d E n E n U 2 ( R) min i i xc[ ] ext[, ] NN ( ) { } 2 i 2 r r R R rr' Empirical potentials fit to SAPT(DFT) calculations ( n) ( l) l ijr q ij (1) iq j C ( n) ij U( R) 1 Aij rij e f1( ij, rij ) fn( ij, rij ) Ubonded ( ) n R ij l rij n6,8,10 rij m n r r fn(, r) 1e m! Machine-learning ( mathematical ) potentials 2 n( r) ( r) i i i j ij m0
5 Structure prediction protocol for Target XXII Optimize monomer geometry using DFT [PBE0 + D3, 6-311G**(d,p) Generate Unit Cells (UPACK) Run MD at exptl. P,T 25,000 generated structures 4,042 unique after clustering 107 within 10 kj/mol of lowest energy Lowest Energy Structures C2/c Pbca Pbca P P2 1 /n P2 1 /n
6 Summary of 50 predicted structures based on MD averaged energies
7 I: P2 1 /n II: Pbca, E = 3.46 kj/mol III: P2 1 /c, E = 5.37 kj/mol IV: P2 1 /n, E = 5.99 kj/mol V: P2 1 /c, E = 6.31 kj/mol
8 Prediction of the experimental structure Final ranking based on POLY1 SAPT(DFT) potential At 300 K: RMSD 20 from experimental structure = 0.18 Å. Rank = 4 on our list At 150 K: RMSD 20 from experimental structure = 0.15 Å. Rank = 4 on our list Spc grp a b c β Expt. P2 1 /n Calc. P2 1 /n Attempted predictions = 21, Times predicted = 12, Times ranked 5 or lower = 7
9 Polymorphism refers to the ability of a compound to form multiple crystal structures. Ritonavir: HIV protease inhibitor NIH: $3,500,000 investment Abbott Labs: $200,000,000 investment Originally dispensed in 1996 as ordinary capsules, no refrigeration required Converted to lower energy (hitherto unknown) polymorph (form I to form II) on the shelf Form II: Poor solubility, lower bioavailability Required recall (1998) and reformulation as gel caps (2002). Pharmaceutical polymorphism HIV-1 protease
10 Pharmaceutical Polymorphism C 13 H 22 N 4 O 3 S Ranitidine hydrochloride (Zantac) Acetylsalicylic acid Bond et al. Amer. Pharma. Rev. (2007)
11 Rough energy landscapes in proteins and crystals Dunitz and Scheraga PNAS 101, (2004)
12 Sampling via collective variables in molecular simulation Peptides: φ ψ Crystals: Cell shape: Packing: Q 4 l 2 l f rb Ylm b b (2l1) N ml b1 N b ( ) (, ) h a b c a b c a b c x x x y y y z z z Orientation:
13 Adiabatic free energy (temperature-accelerated molecular) dynamics (AFED/TAMD) L. Rosso, P. Minary, Z. Zhu, MET J. Chem. Phys. 116, 4389 (2002) Margliano and Vanden-Eijnden, Chem. Phys. Lett. 426, 168 (2006); J. B. Abrams and MET, J. Phys. Chem. B 112, (2008) For prediction of crystal structures: Yu and MET Phys. Rev. Lett. 107, (2011); Yu et al. J. Chem. Phys. 150, (2014) Suppose n collective variables characterize a free energy landscape of interest Canonical probability distribution and free energy surface: q q (,..., ) q( ) 1,..., n r1 rn r n H (, ) ( 1,..., n) pr d p d r ( ( r) ) 1 P s s e q s A( s,..., s ) k T ln P( s,..., s ) 1 n B 1 n Write δ-functions as product of Gaussians: 1/2 2 { } q ( r) s lim exp q ( r) s 2 2 Dynamically sample Gaussian centers s 1,,s n by introducing a kinetic energy for them And extending the phase space, which gives the following Hamiltonian: n 2 p n s H ( p, r, s, ps) H( p, r) q( r) s 2m
14 Adiabatic free energy (temperature-accelerated molecular) dynamics (AFED/TAMD) Introduce high temperature T T for extended variables and high masses m m s i Adiabatically decoupled equations of motion: H q ( ) miri s q r heat bath( T) r r i m s s q r T s ( ) heat bath( ) Under adiabatic conditions, we generate a distribution { } ({ }) lim kts ln Padb ( s1,..., sn, Ts, T) { } { } i P ( s,..., s, T, T) ({ }) adb 1 / lim P ( s,..., s, T, T) P( s,..., s, T) ({ }) adb 1 n s 1 n T lim kts ln P( s1,..., sn, T) A( s1,..., sn, T) Ts TT s n s
15 Other developments: 1. Combining with metadynamics used as a bias: M. Chen. M. Cuendet, MET J. Chem. Phys. 137, (2012) 2. Optimal free energy reconstruction protocols: M. Cuendet and MET J. Chem. Theor. Comput. 10, 2975 (2014) Alanine Decamer (gas phase) Force field: CHARMM22 20 CVs: All (φ,ψ) pairs CV Temperature: T s = 900 K Physical Temp: T = 300 K CV mass: m (ψ,φ) = 600m H Harmonic coupling: kcal/mol/rad 2
16
17 Constant pressure version with cell matrix as CVs GROMOS FF, T 100 K, T s 32,000 40,000 K N 216 (333), P 2 GPa, Total run time = 500 ps, 5 ns
18 X-ray powder diffraction patterns Yonetani et al. (2001). II(98) II(01) 98 01
19 Free energies of benzene polymorphs Six-dimensional free energy surface analyzed using clustering based on a Gaussian mixture model. (see, also, J. Yang et al. Science 345, 640 (2014)) P2 1 /c Cmca P C2/c P2 1 /c II(98) II(01) MURI/ARO
20 Polymorphs of resorcinol Resorcinol (C 6 H 6 O 2 ) is a simple benzene derivative. At ambient conditions, forms two stable forms (α and β) depending on the relative orientations of the OH bonds, both having space group Pna2 11. High pressure phase (δ phase) exists and arises from an α δ at pressures above 3 GPa 2. Structure unknown. Another high pressure phase (ε phase) can be reached from β phase above 5.6 GPa. This phase is also thought to be disordered based on Raman scattering 2. Another ambient polymorph (γ phase) in the R b conformation recently found 3. Its exact structure is unknown 4. Recent new experimental activity in study of resorcinol. 5 1 Robertson and Ubbelohde Proc. Royal Soc. London (1938) 2 Rao et al. Phys. Rev. B (2002). 3 Kahr et al. Cryst. Growth Design (2011). 4 Zhu et al. Angew. Chem. Intl. Ed. (submitted) 5 J. Phys. Chem. B (2015); J. Mol. Struct. (2015).
21 polymorph
22 P21
23 transition at 3 GPa
24 Preliminary free energy estimates for resorcinol Supercell length A form δ form a/a b/a c/a a/a b/a c/a Cell length Cell angle Space GP Pna2 1 Pna2 1
25 Energetic Co-crystals: DADP.TxTNB DADP TCTNB Polarizable system without pre-determined force field parameters = opportunity for ab initio enhanced sampling! Experimental Structures α β ( chloro- ) ( iodo- ) ( bromo- ) 1 3 Exp. Structures (Maztger Group): J. Am. Chem. Soc. 2015, 137, Angew. Chem. Int. Ed. 2013, 52,
26 Energetic Co-crystals: DADP.TxTNB α ( chloro- ) ( bromo- ) β ( iodo- ) For optimized unit cells within PBE+D3, DZVP, 450 Ry in CP2K Cl α < Cl β Br α > Br β ΔE 1.6 kj/mol ΔE 4.9 kj/mol I α >> I β ΔE 15.8 kj/mol Observed DADP.TBTNB co-crystal is a kinetic product aggregation promotes co-crystal growth? Can we obtain the actual free energy difference between the two forms for each co-crystal? Implementation of crystal-afed in CP2K CP2K 2.5: PBE+D3/DZVP-MOLOPT-(SR)-GTH
27 Ab initio molecular dynamics M R I I I min E[{ }, R] { } Start with nuclei Add electrons Compute i, n, i n F Add electrons Propagate nuclei a short time Δt with F e.g. Verlet: 2 t RI ( t) 2 RI (0) RI ( t) FI (0) M I
28 11.0 picosecond of ab initio crystal-afed with CP2K for TCTNB/DADP co-crystal
29 Experimental structure: α New Polymorph Alternate polymorph: β a b c a b c a b c g/cm g/cm g/cm
30 Crystal structures of Xenon at 2700 K and 25 GPa T. Q. Yu, P. Y. Chen, A. Samanta, E. Vanden-Eijnden, MET J. Chem. Phys. (2014) Experimental evidence mixed fcc-hcp structures in pressure-induced Martensitic transitions from fcc to hcp as intermediates: Jephcoat et al. Phys. Rev. Lett. 59, 2670 (1987) Cynn et al. Phys. Rev. Lett. 86, 4552 (2001) Using Q 4 and Q 6 only as collective variables T s = 1.5 x 10 5 K m s = 10 8 amu Using Q 4, Q 6, and the cell matrix h d as collective variables T s = T h = 10 5 K, m s = 10 8 amu MURI/ARO
31
32 The mechanism of equilibrium melting of a solid A. Samanta, T.-Q. Yu, W. E, MET Science (November, 2014) Classical nucleation theory: 4 G r r r ( ) 4 s, < 0 r* r 2 s Gr ( ) s 2
33 Crystal d-afed simulations of the melting of Cu A. Samanta, T.-Q. Yu, W. E, MET Science (November, 2014) Collective variables: h d, Q 4, Q 6, T = T m, T s =T h = 10 7 K Potential model: Embedded Atom Model (EAM) of Mishin et al. Phys. Rev. B 63, (2001).
34 Free energy surfaces One-dimensional profiles constructed using string method [Maragliano et al. JCP (2006)]
35
36 From a minimum, we ask for the nearest accessible saddle points. These can be found using gentlest ascent dynamics (GAD) on the FES. Activation: Navigating on a high-dimensional surface A. Samanta, M. Chen, T. Q. Yu, W. E, MET J. Chem. Phys. 140, (2014). M. Chen, T. Q. Yu, MET PNAS 112, 3235 (2015) Combine with relaxation to next minimum: s F Termed STochastic Activation-Relaxation Technique or START
37 Graph representation of free energy surfaces Landmarks determined using density-based clustering methods. Minima and saddles are assigned as vertices of different colors. An edge connects a saddle vertex to a minimum vertex if a relaxation initiated at the saddle ends up at the corresponding minimum. Example: Alanine tripeptide
38 Full network diagram for met-enkephalin using 10 Ramachandran angles as CVs: 1081 minima, 1431 saddles Graph nodes sorted based on Sketch Map analysis 2 F( Rij ) f ( rij ) i j 2 M. Ceriotti, G. A. Tribello, M. Parrinello Proc. Natl. Acad. Sci. 108, (2011)
39 If we want to compute free energies, use the minima obtained in a START run to indicate where a set of bins should be placed and use AFED sampling to obtain the populations in those bins. Can handle explicit solvent! Can be adapted for crystal structure prediction! [M. Chen, T. Q. Yu, MET (in prep)]
40
41 M. Chen, T. Q. Yu, MET (in prep) Spc grp a b c β Expt. P2 1 /c Calc. P2 1 /c Experimental structure first found in 1 ns using orientational CVs. RMSD 4 of carbons = 0.16 Å, full RMSD 4 = 0.44 Å using a simple, universal force field
42 Conclusions and Perspectives 1. Protocol for crystal structure prediction including specialized force field generation [SAPT(DFT)] combined with molecular dynamics refinement. 2. Enhanced sampling approaches further allow for polymorph prediction and ranking based on free energy, which is the proper figure of merit. 3. Surface navigation combined with enhanced sampling looks to be a promising approach going forward, however, it is important to have a robust set of internal collective variables. 4. Considerable challenges remain including: i. Highly flexible molecules. ii. High-temperature and/or high-pressure polymorphs. iii. Stacking faults, screw axes, mixed phases,. iv. Isotopic polymorphism (requires, e.g., path integrals). 5. From Nature press release based on CCDC 6 th blind test, It is fair to claim that, to a large extent, this blind test has shown that the problem of organic structure prediction has been solved. M. Neumann. This claim might be exaggerated, and the field of crystal structure prediction far from dead.
43 Acknowledgments External Krzysztof Szalewicz (U. Del.) Michael Metz (U. Del.) Andrew Rohl (Curtin U.) Bart Kahr (NYU) Qiang Zhu (Stony Brook) Postdocs Elia Schneider Linas Vilciauskas Leslie Vogt Students Ming Chen Tang-Qing Yu Manav Kumar Hongxing Song
Free 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 informationEnhanced sampling via molecular dynamics II: Unbiased approaches
Enhanced sampling via molecular dynamics II: Unbiased approaches Mark E. Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Sciences New York University, 100 Washington Square East, NY
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 11 I. PHYSICAL AND CHEMICAL RELEVANCE OF FREE ENERGY In this section, we will consider some examples showing the significance of free energies.
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 informationPolymorphism. 6 June 2012 Lecture 4 RUB
Polymorphism 1 Polymorphism It is the phenomenon in which the same chemical substance exhibits different internal crystal packing arrangements. Polymorphism is an exclusively solid state phenomenon. Polymorphs
More informationAb initio phonon calculations in mixed systems
Ab initio phonon calculations in mixed systems Andrei Postnikov apostnik@uos.de Outline: Experiment vs. ab initio theory Ways of theory: linear response and frozen phonon approaches Applications: Be x
More informationStructure of Cement Phases from ab initio Modeling Crystalline C-S-HC
Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model
More 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 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationMachine learning the Born-Oppenheimer potential energy surface: from molecules to materials. Gábor Csányi Engineering Laboratory
Machine learning the Born-Oppenheimer potential energy surface: from molecules to materials Gábor Csányi Engineering Laboratory Interatomic potentials for molecular dynamics Transferability biomolecular
More informationWhy polymorphism? An Evaluation using Experimental Charge Densities Analysis
Why polymorphism? An Evaluation using Experimental Charge Densities Analysis T. N. Guru Row Solid State and Structural Chemistry Unit Indian Institute of Science Bangalore 560012 INDIA Email: ssctng@sscu.iisc.ernet.in
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 informationOptimisation: From theory to better prediction
Optimisation: From theory to better prediction Claire Adjiman Claire Adjiman 23 May 2013 Christodoulos A. Floudas Memorial Symposium 6 May 2017 Why look for the global solution? Deterministic global optimisation
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 informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More 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 informationSupporting Information. Ze-Min Zhang, Lu-Yi Pan, Wei-Quan Lin, Ji-Dong Leng, Fu-Sheng Guo, Yan-Cong Chen, Jun-Liang Liu, and Ming-Liang Tong*
Supporting Information Wheel-shaped nanoscale 3d-4f {Co II 16Ln III 24} clusters (Ln = Dy and Gd) Ze-Min Zhang, Lu-Yi Pan, Wei-Quan Lin, Ji-Dong Leng, Fu-Sheng Guo, Yan-Cong Chen, Jun-Liang Liu, and Ming-Liang
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 informationCrystal Structure Prediction using CRYSTALG program
Crystal Structure Prediction using CRYSTALG program Yelena Arnautova Baker Laboratory of Chemistry and Chemical Biology, Cornell University Problem of crystal structure prediction: - theoretical importance
More informationHeadache because of Aspirin! Roland Boese
Headache because of Aspirin! Roland Boese Headache because of Aspirin! About 80% of all APIs exhibit polymorphism! What is so special about Aspirin so that it seems to be monomorphic in spite of millions
More informationSpeed of light c = m/s. x n e a x d x = 1. 2 n+1 a n π a. He Li Ne Na Ar K Ni 58.
Physical Chemistry II Test Name: KEY CHEM 464 Spring 18 Chapters 7-11 Average = 1. / 16 6 questions worth a total of 16 points Planck's constant h = 6.63 1-34 J s Speed of light c = 3. 1 8 m/s ħ = h π
More informationSupplementary Online Materials: Formation of Stoichiometric CsF n Compounds
1 2 3 4 5 6 7 8 9 1 11 12 13 Supplementary Online Materials: Formation of Stoichiometric CsF n Compounds Qiang Zhu, 1, a) Artem R. Oganov, 1, 2, 3 and Qingfeng Zeng 4 1) Department of Geosciences, Stony
More informationFrom Dynamics to Thermodynamics using Molecular Simulation
From Dynamics to Thermodynamics using Molecular Simulation David van der Spoel Computational Chemistry Physical models to describe molecules Software to evaluate models and do predictions - GROMACS Model
More informationAn EAM potential for the dynamical simulation of Ni-Al alloys
J. At. Mol. Sci. doi: 10.4208/jams.022310.031210a Vol. 1, No. 3, pp. 253-261 August 2010 An EAM potential for the dynamical simulation of Ni-Al alloys Jian-Hua Zhang, Shun-Qing Wu, Yu-Hua Wen, and Zi-Zhong
More informationHow Accurately do Current Force Fields Predict. Experimental Peptide Conformations? An Adiabatic Free Energy Dynamics Study:
How Accurately do Current Force Fields Predict Experimental Peptide Conformations? An Adiabatic Dynamics Study: Supplementary Information (SI) Alexandar T. Tzanov, Michel A. Cuendet,, and Mark E. Tuckerman,,
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 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 informationPhysical Chemistry I CHEM 4641 Final Exam 13 questions, 30 points
Physical Chemistry I CHEM 4641 Final Exam 13 questions, 30 points Name: KEY Gas constant: R = 8.314 J mol -1 K -1 = 0.008314 kj mol -1 K -1. Boltzmann constant k = 1.381 10-23 J/K = 0.6950 cm -1 /K h =
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 informationCasos de Estudio usando Materials Studio
Casos de Estudio usando Materials Studio Javier Ramos Biophysics of Macromolecular Systems group (BIOPHYM) Departamento de Física Macromolecular Instituto de Estructura de la Materia CSIC j.ramos@iem.cfmac.csic.es
More informationControlling fluctuations
Controlling fluctuations Michele Parrinello Department of Chemistry and Applied Biosciences ETH Zurich and ICS, Università della Svizzera Italiana, Lugano, Switzerland Today s menu Introduction Fluctuations
More informationMany-body dispersion interactions between semiconducting wires
Many-body dispersion interactions between semiconducting wires Alston J. Misquitta TCM Group, Cavendish Laboratory University of Cambridge 26 th July, 2010 Collaborators James Spencer (Imperial College,
More informationSoot - Developing anisotropic potentials from first principles for PAH molecules. Tim Totton, Alston Misquitta and Markus Kraft 12/11/2009
Soot - Developing anisotropic potentials from first principles for PAH molecules. Tim Totton, Alston Misquitta and 12/11/2009 HRTEM images of soot Some evidence for different soot structures based on different
More informationMolecular dynamics simulations of EXAFS in germanium
Cent. Eur. J. Phys. 93 2011 710-715 DOI: 10.2478/s11534-010-0074-0 Central European Journal of Physics Molecular dynamics simulations of EXAFS in germanium Research Article Janis Timoshenko Alexei Kuzmin
More informationOne hydrogen bond doesn t make a separation or does it? Resolution of amines by diacetoneketogulonic acid
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2015 Electronic Supplementary Information One hydrogen bond doesn t make a separation or does it? Resolution
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
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 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 informationMolecular Dynamics Simulations of Glass Formation and Crystallization in Binary Liquid Metals
Citation & Copyright (to be inserted by the publisher ) Molecular Dynamics Simulations of Glass Formation and Crystallization in Binary Liquid Metals Hyon-Jee Lee 1,2, Tahir Cagin 2, William A. Goddard
More informationarxiv:physics/ v1 [physics.atm-clus] 21 Jun 2004
Conformational properties of neutral and charged alanine and glycine chains arxiv:physics/0406093v1 [physics.atm-clus] 21 Jun 2004 Alexander V. Yakubovitch, Ilia A. Solov yov, Andrey V. Solov yov, and
More informationTheory and Applications of Residual Dipolar Couplings in Biomolecular NMR
Theory and Applications of Residual Dipolar Couplings in Biomolecular NMR Residual Dipolar Couplings (RDC s) Relatively new technique ~ 1996 Nico Tjandra, Ad Bax- NIH, Jim Prestegard, UGA Combination of
More informationStructure and interactions in benzamide molecular crystals
Structure and interactions in benzamide molecular crystals Philipp Ectors 1, Dominique Ectors 2, Dirk Zahn 1 1) Lehrstuhl für Theoretische Chemie/Computer-Chemie-Centrum Friedrich-Alexander- Universität
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 informationMaking Electronic Structure Theory work
Making Electronic Structure Theory work Ralf Gehrke Fritz-Haber-Institut der Max-Planck-Gesellschaft Berlin, 23 th June 2009 Hands-on Tutorial on Ab Initio Molecular Simulations: Toward a First-Principles
More informationAnalysis of the Lennard-Jones-38 stochastic network
Analysis of the Lennard-Jones-38 stochastic network Maria Cameron Joint work with E. Vanden-Eijnden Lennard-Jones clusters Pair potential: V(r) = 4e(r -12 - r -6 ) 1. Wales, D. J., Energy landscapes: calculating
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 information1. What is the difference between intermolecular forces and intramolecular bonds? Variations in the Boiling Point of Noble Gases
NAME: DATE: Chemical Bonding Forces Assignment 1. What is the difference between intermolecular forces and intramolecular bonds? 2. Use your data booklet to fill in the following chart Variations in the
More informationMolecular Dynamics Simulation Study of the Ionic Mobility of OH Using the OSS2 Model
1154 Bull. Korean Chem. Soc. 2006, Vol. 27, No. 8 Song Hi Lee Molecular Dynamics Simulation Study of the Ionic Mobility of OH Using the OSS2 Model Song Hi Lee Department of Chemistry, Kyungsung University,
More informationRepresenting High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks
Energy Landscapes, Chemnitz 200 Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks Jörg Behler Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum D-44780 Bochum,
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
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 informationSimulations of Self-Assembly of Polypeptide-Based Copolymers
East China University of Science and Technology Theory, Algorithms and Applications of Dissipative Particle Dynamics Simulations of Self-Assembly of Polypeptide-Based Copolymers Jiaping LIN ( 林嘉平 ) East
More informationCh. 11: Liquids and Intermolecular Forces
Ch. 11: Liquids and Intermolecular Forces Learning goals and key skills: Identify the intermolecular attractive interactions (dispersion, dipole-dipole, hydrogen bonding, ion-dipole) that exist between
More informationAtom in Molecules a Quantum Theory (AIM)
Atom in Molecules a Quantum Theory (AIM) Richard F. W. Bader* The primary purpose in postulating the existence of atoms in molecules is a consequence of the observation that atoms or groupings of atoms
More informationDimer Dissociation of a Photoreceptor Protein from QM/MM and MD Simulations
Dimer Dissociation of a Photoreceptor Protein from QM/MM and MD Simulations IMA University of Minnesota Minneapolis, MN, July 20, 2015 Haisheng Ren Advisor: Prof. Jiali Gao Department of Chemistry, University
More informationHKUST-Shenzhen Research Institute, No. 9 Yuexing 1st RD, South Area, Hi-tech Park, Nanshan, Shenzhen , China
Electronic Supporting Information Aggregation-Induced Emission and Photocyclization of Poly(hexaphenyl-1,3-butadiene)s Synthesized from 1+2 Polycoupling of Internal Alkynes and Arylboronic Acids Yajing
More informationScreening for cocrystals of succinic acid and 4-aminobenzoic acid. Supplementary Information
Screening for cocrystals of succinic acid and 4-aminobenzoic acid Nizar Issa, Sarah A. Barnett, Sharmarke Mohamed, Doris E. Braun, Royston C. B. Copley, Derek A. Tocher, Sarah L Price* Supplementary Information
More informationChemistry 101 Chapter 9 CHEMICAL BONDING. Chemical bonds are strong attractive force that exists between the atoms of a substance
CHEMICAL BONDING Chemical bonds are strong attractive force that exists between the atoms of a substance Chemical Bonds are commonly classified into 3 types: 1. IONIC BONDING Ionic bonds usually form between
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
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 informationSupplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water
Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical
More informationAb initio Berechungen für Datenbanken
J Ab initio Berechungen für Datenbanken Jörg Neugebauer University of Paderborn Lehrstuhl Computational Materials Science Computational Materials Science Group CMS Group Scaling Problem in Modeling length
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationON PREDICTING THE CRYSTAL STRUCTURE OF ENERGETIC MATERIALS FROM QUANTUM MECHANICS
ON PREDICTING THE CRYSTAL STRUCTURE OF ENERGETIC MATERIALS FROM QUANTUM MECHANICS Betsy M. Rice* U. S. Army Research Laboratory, AMSRD-ARL-WM-BD Aberdeen Proving Ground, MD 21005-5069 Rafal Podeszwa University
More informationVERSION A White. perm. version. name. Chem 1C - Spring exam 2
Chem 1C - Spring 2010 - exam 2 VERSION A White ON YOUR SCANTRON: BUBBLE IN YOUR 7 DIDGIT PERM (Leaving the last three digits blank) BUBBLE IN THE VERSION: A WRITE YOUR NAME AT THE END, HAND IN ONLY THE
More informationChemistry 112 Spring 2007 Prof. Metz Practice Exam 1 Solutions
Chemistry 112 Spring 2007 Prof. Metz Practice Exam 1 Solutions 1. The intermolecular attractive forces would be greatest in which of the following molecules: CH 4, CH 2 Cl 2 and CO 2. (A) CH 4 (B) CH 2
More informationCoupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation
Coupling the Level-Set Method with Variational Implicit Solvent Modeling of Molecular Solvation Bo Li Math Dept & CTBP, UCSD Li-Tien Cheng (Math, UCSD) Zhongming Wang (Math & Biochem, UCSD) Yang Xie (MAE,
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 informationChapter 6 Cyclic urea - a new central unit in bent-core compounds
82 Chapter 6 Cyclic urea - a new central unit in bent-core compounds A new class of five-ring bent-core molecules with a cyclic urea group as a central unit was synthesized [94]. A significant difference
More 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 informationTopics in the November 2009 Exam Paper for CHEM1101
November 2009 Topics in the November 2009 Exam Paper for CHEM1101 Click on the links for resources on each topic. 2009-N-2: 2009-N-3: 2009-N-4: 2009-N-5: 2009-N-6: 2009-N-7: 2009-N-8: 2009-N-9: 2009-N-10:
More informationApplications of the CSD to Structure Determination from Powder Data
Applications of the CSD to Structure Determination from Powder Data Prof. Alastair J. Florence Strathclyde Institute of Pharmacy and Biomedical Sciences, Glasgow, UK ACS, Salt Lake City, March 2009 Overview
More informationVERSION B Yellow. perm. version. name. Chem 1C - Spring exam 2
Chem 1C - Spring 2010 - exam 2 VERSION B Yellow ON YOUR SCANTRON: BUBBLE IN YOUR 7 DIDGIT PERM (Leaving the last three digits blank) BUBBLE IN THE VERSION: B WRITE YOUR NAME AT THE END, HAND IN ONLY THE
More informationQuantum Monte Carlo simulation of spin-polarized tritium
Higher-order actions and their applications in many-body, few-body, classical problems Quantum Monte Carlo simulation of spin-polarized tritium I. Bešlić, L. Vranješ Markić, University of Split, Croatia
More informationCharacterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi
Characterizing Structural Transitions of Membrane Transport Proteins at Atomic Detail Mahmoud Moradi NCSA Blue Waters Symposium for Petascale Science and Beyond Sunriver, Oregon May 11, 2015 Outline Introduction
More informationAspects of nonautonomous molecular dynamics
Aspects of nonautonomous molecular dynamics IMA, University of Minnesota, Minneapolis January 28, 2007 Michel Cuendet Swiss Institute of Bioinformatics, Lausanne, Switzerland Introduction to the Jarzynski
More informationCrystal and molecular structure of N-(p-nitrobenzylidene)- 3-chloro-4-fluoroaniline
PRAMANA cfl Indian Academy of Sciences Vol. 55, No. 3 journal of September 2000 physics pp. 441 446 Crystal and molecular structure of N-(p-nitrobenzylidene)- 3-chloro-4-fluoroaniline K V ARJUNA GOWDA,
More information1. Increasing the pressure above a liquid will cause the boiling point of the liquid to:
JASPERSE CHEM 210 PRACTICE TEST 1 VERSION 2 Ch. 11 Liquids, Solids, and Materials Ch. 10 Gases Ch. 15 The Chemistry of Solutes and Solutions Ch. 13 Chemical Kinetics 1 Constants and/or Formulas Formulas
More informationCurrent approaches to predicting molecular organic crystal structures
Crystallography Reviews ISSN: 0889-311X (Print) 1476-3508 (Online) Journal homepage: http://www.tandfonline.com/loi/gcry20 Current approaches to predicting molecular organic crystal structures Graeme M.
More informationChem 112 Dr. Kevin Moore
Chem 112 Dr. Kevin Moore Gas Liquid Solid Polar Covalent Bond Partial Separation of Charge Electronegativity: H 2.1 Cl 3.0 H Cl δ + δ - Dipole Moment measure of the net polarity in a molecule Q Q magnitude
More informationGround-state Structure and Dynamics
Ground-state Structure and Dynamics Jonathan Yates jonathan.yates@materials.ox.ac.uk Materials Modelling Laboratory, Oxford Materials For a given set of atomic positions the ions will experience a force
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 informationApplications: Molecular crystals Graphite MgO(001)/CO MIL-53(Al) 2
Bartolomeo Civalleri Voice: Loredana Valenzano B3LYP augmented with an empirical dispersion term (B3LYP-D*) as applied to solids Università di Torino Dipartimento di Chimica IFM & NIS Torino - MSSC2009-10/09/09
More informationComputing free energy: Replica exchange
Computing free energy: Replica exchange 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 Macroscopic i a t e regime
More informationOverview & Applications. T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 04 June, 2015
Overview & Applications T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 4 June, 215 Simulations still take time Bakan et al. Bioinformatics 211. Coarse-grained Elastic
More informationAdaptive algorithm for saddle point problem for Phase Field model
Adaptive algorithm for saddle point problem for Phase Field model Jian Zhang Supercomputing Center, CNIC,CAS Collaborators: Qiang Du(PSU), Jingyan Zhang(PSU), Xiaoqiang Wang(FSU), Jiangwei Zhao(SCCAS),
More informationWatching crystals grow and transform
Watching crystals grow and transform Elias Vlieg Radboud University, Institute for Molecules and Materials crystals transforming crystals solution-mediated solid-solid transition Institute for Molecules
More informationLiquids, Solids and Phase Changes
Chapter 10 Liquids, Solids and Phase Changes Chapter 10 1 KMT of Liquids and Solids Gas molecules have little or no interactions. Molecules in the Liquid or solid state have significant interactions. Liquids
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 informationJoint ICTP-IAEA Workshop on Physics of Radiation Effect and its Simulation for Non-Metallic Condensed Matter.
2359-16 Joint ICTP-IAEA Workshop on Physics of Radiation Effect and its Simulation for Non-Metallic Condensed Matter 13-24 August 2012 Introduction to atomistic long time scale methods Roger Smith Loughborough
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 informationImproving the accuracy of lattice energy calculations in crystal structure prediction using experimental data
Improving the accuracy of lattice energy calculations in crystal structure prediction using experimental data Christina-Anna Gatsiou A thesis submitted for the Doctor of Philosophy degree of Imperial College
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 informationEnergy Profiles and Chemical Reactions
Energy rofiles and Chemical eactions + B C D E egensburg er-la Norrby Modeling Kinetics Molecular Modeling Stationary points Energies Barriers eaction rates DFT, ab initio, MM Ground & Transition states
More informationMultiple Choice. Multiple Choice
1. At what temperature in degree Celcius is the value in degree Fahrenheit twice of that in degree Celcius? A) 160 o C B) -24.6 o C C) 6.4 o C D) 22.2 o C E) 32 o C 2. The correct name for NaOCl is, A)
More information2. As gas P increases and/or T is lowered, intermolecular forces become significant, and deviations from ideal gas laws occur (van der Waal equation).
A. Introduction. (Section 11.1) CHAPTER 11: STATES OF MATTER, LIQUIDS AND SOLIDS 1. Gases are easily treated mathematically because molecules behave independently. 2. As gas P increases and/or T is lowered,
More informationTemperature and pressure dependence of the Raman frequency shifts in anthracene
Indian Journal of Pure & Applied Physics Vol. 54, August 2016, pp. 489-494 Temperature and pressure dependence of the Raman frequency shifts in anthracene H Özdemir & H Yurtseven* Department of Physics,
More information7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
More informationIAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.
IAP 2006: From nano to macro: Introduction to atomistic modeling techniques and application in a case study of modeling fracture of copper (1.978 PDF) http://web.mit.edu/mbuehler/www/teaching/iap2006/intro.htm
More information