Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
|
|
- Lucinda Price
- 5 years ago
- Views:
Transcription
1 Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014
2 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
3 What is MD? We call molecular dynamics (MD) a computer simulation technique where the time evolution of a set of interacting particles is followed by integrating their equations of motion. In MD particles obey classical mechanics (Newton s law): ~F i = m i ~a i In contrast with the Monte Carlo methods, MD is a deterministic method: given an initial set of positions and velocities, the time evolution of the system of interacting particles is in principle completely determined. Through a MD simulations we can compute equilibrium and transport properties of a system.
4 What is MD? The result of a MD simulation is a trajectory in a 6N-dimensional phase space. Often we are not particularly interested in the trajectory itself, but in using this trajectory as a way to sample some statistical distribution function (statistical ensemble), like in Monte Carlo simulations. The probability density function ρ in the various ensembles is given by: Microcanonical Ensemble (NVE): Canonical Ensemble (NVT): ( ) / (H( ) E) ( ) / e H( ) K B T Isothermal-Isobaric Ensemble (NPT): ( ) / e H( )+pv k B T Grand Canonical Ensemble (µvt): ( ) / e H( ) µ K B T
5 What is MD? In statistical physics, physical quantities are given by averages over configurations distributed according to specific statistical ensembles (ensemble averages). Under the ergodic hypothesis, a sufficiently long MD trajectory provides such a set for the ensemble in which it was run. Ensemble averages can hence be evaluated as time averages over an MD trajectory: Ensemble average: hai = Z A( ) ( )d Time average: Ergodicity: A T = 1 M MX n=1 lim A T!hAi t!1 A(t n )
6 Force Fields Forces F i are given in (classical) MD by empirical force fields. Examples: 2-body potentials, Lennard-Jones
7 Force Fields Forces F i are given in (classical) MD by empirical force fields. Examples: 3-body potentials, Stillinger-Weber
8 Force Fields Forces F i are given in (classical) MD by empirical force fields. Examples: Force fields for proteins (Amber, Charmm, ) (Image: Adv Drug Del Rev)
9 Integration of the E.O.M. Taylor expansion of r(t): Summing these 2 eqns:
10 Integration of the E.O.M. Good things about the Verlet algorithm: Time reversibility Symplectic algorithm Liouville s theorem (Hamiltonian dynamics): da=da Stability of the algorithm
11 MD run 1) Initialization: Assign positions and velocities to particles, usually followed by structural optimization 2) Equilibration: Physical quantities need to become stationary before we start taking averages 3) Measurement: Simulate a sufficiently long trajectory and perform time averages of the quantities of interest Instantaneous temperature: he kin i = 1 2 3NX i=1 m i v 2 i = 1 2 (3N)k BT In NVT and NPT ensembles we want to impose a certain temperature T w
12 Thermostats Velocity rescaling = p T W /T (t) Berendsen (1984) dt (t) dt = 1 (T 0 T (t)) 2 = q + t T0 T (t) 1 τ : τ=δt: velocity rescaling microcanonical ensemble Andersen (1980) Select one particle at random every τ and reassign its velocity drawing from the Maxwell-Botzmann distribution. Correct NVT sampling for τ No constant of motion. Discontinuous trajectories.
13 Thermostats Nose -Hoover S. Nosé, J. Chem. Phys. 81, 511 (1984). W. G. Hoover, Phys. Rev. A 31, 1695 (1985). Canonical sampling through velocity rescaling (CSVR) Bussi et al. J. Chem. Phys. 126, (2007) Both have a constant of motion and sample the NVT ensemble Barostats (NPT) Parrinello-Rahman barostat Parrinello and Rahman, J. Appl. Phys. 52 (1981) 7182
14 Why ab initio MD? Free of parameterization and only fundamental constants required. Bond forming and breaking processes can, in principle, be treated accurately. Transferable, unlike many classical MD methods.
15 Quantum Dynamics r electrons, R ions electron-electron repulsion ion-ion repulsion electron-ion attraction Total wavefunction (both electrons and nuclei treated as QM particles) Full quantum dynamics is extremely complicated and not needed in most cases.
16 Decoupling electrons and ions First solve for clamped nuclei Instantaneous representation for full solution Integrate out r dependence
17 Decoupling electrons and ions Non-adiabatic coupling terms: l leads to the Born-Oppenheimer approximation
18 Born-Oppenheimer Dynamics Kinetic energy Effective potential energy In the BO approximation, the dynamics of ions does not change the state of the electronic sub-system. Thus, The dependence of the electronic wave-function on the nuclear positions is now suppressed. Ek(R) defines an energy eigenvalue for every configuration R - Potential Energy Surface (of the k th electronic eigenstate) Classical ions would move on the ground-state (k=0) surface.
19 Classical dynamics for ions Let s define a classical Lagrangean for B-O dynamics Equations of motion are obtained by Thus we have The potential energy at every instant is obtained by solving the electronic problem variationally.
20 DFT Energy and forces where the first term is just the DFT energy computed from the groundstate density n(r) by methods discussed earlier in this school. Hellmann-Feynman forces : External potential on electrons
21 Thermodynamic averages
22 BO molecular dynamics At every time step we need to perform a minimization to reach self consistency. This could be time-consuming (especially around 1985 where iterative diagonalization schemes were not employed in first-principles calculations). The accuracy of the simulation critically depends on the accuracy of the SCF minimization (Energy drifts). The performance depends on the algorithms used to extrapolate the wavefunctions from the previous steps.
23 BO-MD: energy conservation The SCF threshold strongly influences the energy conservation (drift) The extrapolation algorithm strongly influences the number of iterations required to reach the desired SCF convergence Herbert and Head-Gordon, PCCP 7, 3269 (2005)
24 BO-MD: energy conservation Similarly, the error on the forces critically depends on the accuracy of the SCF calculation Herbert and Head-Gordon, PCCP 7, 3269 (2005)
25 Choosing parameters for MD A larger time step means smaller overall computational cost. The maximum time step is set by the fastest ionic motion: dt ~ dtmax where dtmax=1/ωmax=period of fastest phonon. The criterion for SCF convergence at each step should be very strict to conserve energy, avoid systematic drifts and ensure accurate forces.
26 Historically, this is the first firstprinciples MD reported. Today (July 2014) total no. of citations exceeds 6500! Trieste, June 2007
27 The CP Method The dynamics of electrons and nuclei happen of very different time scales: electrons move very fast compared to the slow nuclei. Car & Parrinello: exploit this adiabatic time scale separation Idea : introduce a fictitious dynamics for the electrons which keeps the electrons close to the ground state, by treating the orbital degrees of freedom as classical variables.
28 The CP Method CP Lagrangean : Euler-Lagrange equations : Equations of motion :
29 The CP Method Electronic degrees of freedom are the expansion coefficients of KS orbitals into PW. Forces on electrons are determined by the KS Hamiltonian calculated from current values of ϕk and RI. Forces acting on nuclei have the Hellman-Feynman form : but they slightly differ from true (i.e. BO) forces. The dynamics of the orbitals is fictitious, i.e. not real dynamics of electrons.
30 The CP Method Ions evolve in time at a certain (instantaneous) physical temperature Electrons evolve in time at a certain (instantaneous) fictitious temperature We have to ensure that during the dynamics the electrons stay cold while following adiabatically the motion of the nuclei. This is possible if the two sub-systems (electrons and nuclei) are decoupled, i.e. if their power spectra do not overlap. This therefore amounts to impose and preserve a metastability condition for a sufficiently long time, since equipartition would ensure that each degree of freedom acquires ½ kbt. The fictitious kinetic energy is a measure of the departure from the BO surface.
31 The CP Method Econs (constant of motion): accuracy of the integration of the e.o.m. Ephys and Te: degree of adiabatic separation.
32 Adiabatic separation in CP (i.e. why does it work?) Power spectrum of the dynamics of the electronic degrees of freedom for a 8-Si atom system sampled at the Γ point. Power spectrum obtained from the Fourier transform of the velocity-velocity correlation function:
33 Adiabatic separation in CP Frequency spectrum of the orbital degrees of freedom Lowest electronic frequency: To ensure adiabatic separation => reduce the fictitious mass as much as possible Highest electronic frequency: Maximum time step to integrate the fastest frequency.
34 Adiabatic separation in CP The fictitious mass µ and time step δt must be chosen in such a way that: µ is big enough to enable the use of a reasonable time step, and small enough to guarantee: o adiabaticity, i.e. no energy transfer from nuclei to electrons, which always remain close to the ground state (BO surface) o correctness of the nuclear trajectory (small error in forces) Typical values: ~ electron masses δt should be the largest value that yields a stable dynamics (no drifts, no loss of orthonormality)
35 Errors on forces The error of the CP forces compared to the BO ones show two distinct time scales: High-frequency small-amplitude oscillations. Low-frequency ( = nuclear) large amplitude fluctuations, due to the drag of the nuclei exerted on the electrons, due to the fact that they have a finite mass µ. Note the the finite mass of the electrons dampens artificially the nuclear motion. This effect can be controlled by reducing the fictitious mass.
36 Energy drift : CP vs BO The use of the Verlet algorithm for propagating the e.o.m. ensures reversibility, which can be shown to lead to no energy drift. BO dynamics, on the other hand, can suffer from systematic energy drifts, which however can be controlled by the degree of convergence of the SCF procedure.
37 Metals with CP? If the gap reduces to zero the approach breaks down, since it is not possible to decouple the electronic and ionic energy scales. Possible remedies: Thermostat the electrons to prevent them from heating up Pasquarello et al. PRL 69, 1982 (1992). Ensemble DFT: the function to minimize is a free energy Marzari et al. PRL 82, 3296 (1999). This can be cast in an extended Lagrangian formalism: Stengel et al. PRB (2000). Use BO molecular dynamics.
38 CP in Quantum-ESPRESSO CP molecular dynamics can be performed with the cp.x code. Norm-conserving as well as ultrasoft pseudopotentials can be used. Gamma-point only (wave functions are real). NVT (Nose-Hoover chains) as well as NPT (Parrinello-Rahman) ensembles. Simulations in finite electric fields. The code is interfaced with the PLUMED plug-in to perform free energy calculations (metadynamics, umbrella sampling,...)
39 References Most of the equations have been taken from D. Marx and J. Hutter, Ab initio molecular dynamics, Cambridge University Press, 2009 Part of the material has been adapted from P. Giannozzi s tutorial on QE The original CP paper is R. Car and M. Parrinello, Phys. Rev. Lett. 55, 2471 (1985).
Energy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationAb initio molecular dynamics : BO
School on First Principles Simulations, JNCASR, 2010 Ab initio molecular dynamics : BO Vardha Srinivasan IISER Bhopal Why ab initio MD Free of parametrization and only fundamental constants required. Bond
More 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 informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationAb initio molecular dynamics
Ab initio molecular dynamics Molecular dynamics Why? allows realistic simulation of equilibrium and transport properties in Nature ensemble averages can be used for statistical mechanics time evolution
More informationAb initio Molecular Dynamics Born Oppenheimer and beyond
Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator
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 informationAb initio molecular dynamics: Basic Theory and Advanced Met
Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site
More informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
More 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 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 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 informationGear methods I + 1/18
Gear methods I + 1/18 Predictor-corrector type: knowledge of history is used to predict an approximate solution, which is made more accurate in the following step we do not want (otherwise good) methods
More informationClassical Molecular Dynamics
Classical Molecular Dynamics Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~mijp1 Overview of lecture n Motivation n Types of
More informationAb initio Molecular Dynamics
Ab initio Molecular Dynamics Jürg Hutter Department of Chemistry, University of Zurich Overview Equations of motion and Integrators General Lagrangian for AIMD Car Parrinello MD and Born Oppenheimer MD
More informationAb Initio Molecular Dynamic
Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master / RFCT Ab Initio Molecular Dynamic p.1/67 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD)
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 informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More informationJavier Junquera. Statistical mechanics
Javier Junquera Statistical mechanics From the microscopic to the macroscopic level: the realm of statistical mechanics Computer simulations Thermodynamic state Generates information at the microscopic
More informationAb Initio Molecular Dynamic
Ab Initio Molecular Dynamic Paul Fleurat-Lessard AtoSim Master Ab Initio Molecular Dynamic p.1/67 Summary I. Introduction II. Classical Molecular Dynamic III. Ab Initio Molecular Dynamic (AIMD) 1 - Which
More informationLangevin Dynamics in Constant Pressure Extended Systems
Langevin Dynamics in Constant Pressure Extended Systems D. Quigley and M.I.J. Probert CMMP 2004 1 Talk Outline Molecular dynamics and ensembles. Existing methods for sampling at NPT. Langevin dynamics
More informationTime reversible Born Oppenheimer molecular dynamics
Time reversible Born Oppenheimer molecular dynamics Jianfeng Lu Mathematics Department Department of Physics Duke University jianfeng@math.duke.edu KI-Net Conference, CSCAMM, University of Maryland, May
More informationComputational Chemistry - MD Simulations
Computational Chemistry - MD Simulations P. Ojeda-May pedro.ojeda-may@umu.se Department of Chemistry/HPC2N, Umeå University, 901 87, Sweden. May 2, 2017 Table of contents 1 Basics on MD simulations Accelerated
More informationCar-Parrinello Molecular Dynamics
Car-Parrinello Molecular Dynamics Eric J. Bylaska HPCC Group Moral: A man dreams of a miracle and wakes up with loaves of bread Erich Maria Remarque Molecular Dynamics Loop (1) Compute Forces on atoms,
More informationUnderstanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry
JNCASR August 20, 21 2009 Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles Srikanth Sastry Jawaharlal Nehru Centre for Advanced Scientific Research, Bangalore
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationModeling of Nanostructures and Materials Jacek A. Majewski Summer Semester 2013 Lecture 6 Lecture Modeling of Nanostructures Molecular Dynamics
Chair of Condensed Matter Physics nstitute of Theoretical Physics Faculty of Physics, Universityof Warsaw Summer Semester 013 Lecture Modeling of Nanostructures and Materials Jacek A. Majewski E-mail:
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 14, 2012 All electrons vs pseudopotentials Classes of Basis-set Condensed phase: Bloch s th and PBC Hamann-Schlüter-Chiang pseudopotentials
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationSimulation of molecular systems by molecular dynamics
Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction
More informationConstant Pressure Langevin Dynamics: Theory and Application to the Study of Phase Behaviour in Core-Softened Systems.
Constant Pressure Langevin Dynamics: Theory and Application to the Study of Phase Behaviour in Core-Softened Systems. David Quigley A thesis submitted for the degree of Doctor of Philosophy University
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationModeling & Simulation of Glass Structure
Modeling & Simulation of Glass Structure VCG Lecture 21 John Kieffer Department of Materials Science and Engineering University of Michigan 1 Overview Historical perspective Simulation methodologies Theoretical
More informationSet the initial conditions r i. Update neighborlist. r i. Get new forces F i
Set the initial conditions r i t 0, v i t 0 Update neighborlist Get new forces F i r i Solve the equations of motion numerically over time step t : r i t n r i t n + v i t n v i t n + Perform T, P scaling
More informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Length Scales (c) ICAMS: http://www.icams.de/cms/upload/01_home/01_research_at_icams/length_scales_1024x780.png Goals for today: Background
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 informationA Brief Introduction to Ab Initio Molecular Dynamics
A Brief Introduction to Ab Initio Molecular Dynamics Matt Probert Condensed Matter Dynamics Group, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Overview of Talk In this talk
More informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationThe micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction
ISBN 978-1-84626-081-0 Proceedings of the 2010 International Conference on Application of Mathematics and Physics Volume 1: Advances on Space Weather, Meteorology and Applied Physics Nanjing, P. R. China,
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationOptimization Methods via Simulation
Optimization Methods via Simulation Optimization problems are very important in science, engineering, industry,. Examples: Traveling salesman problem Circuit-board design Car-Parrinello ab initio MD Protein
More informationCitation for published version (APA): Hess, B. (2002). Stochastic concepts in molecular simulation Groningen: s.n.
University of Groningen Stochastic concepts in molecular simulation Hess, Berk IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF) if you wish to cite from it. Please check
More informationModeling Materials. Continuum, Atomistic and Multiscale Techniques. gg CAMBRIDGE ^0 TADMOR ELLAD B. HHHHM. University of Minnesota, USA
HHHHM Modeling Materials Continuum, Atomistic and Multiscale Techniques ELLAD B. TADMOR University of Minnesota, USA RONALD E. MILLER Carleton University, Canada gg CAMBRIDGE ^0 UNIVERSITY PRESS Preface
More informationIntroduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012
Introduction to Computer Simulations of Soft Matter Methodologies and Applications Boulder July, 19-20, 2012 K. Kremer Max Planck Institute for Polymer Research, Mainz Overview Simulations, general considerations
More informationCHE3935. Lecture 4 Quantum Mechanical Simulation Methods Continued
CHE3935 Lecture 4 Quantum Mechanical Simulation Methods Continued 1 OUTLINE Review Introduction to CPMD MD and ensembles The functionals of density functional theory Return to ab initio methods Binding
More 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 informationMolecular Dynamics Simulations
Molecular Dynamics Simulations Dr. Kasra Momeni www.knanosys.com Outline Long-range Interactions Ewald Sum Fast Multipole Method Spherically Truncated Coulombic Potential Speeding up Calculations SPaSM
More informationFast Eigenvalue Solutions
Fast Eigenvalue Solutions Techniques! Steepest Descent/Conjugate Gradient! Davidson/Lanczos! Carr-Parrinello PDF Files will be available Where HF/DFT calculations spend time Guess ρ Form H Diagonalize
More informationAdvanced sampling. fluids of strongly orientation-dependent interactions (e.g., dipoles, hydrogen bonds)
Advanced sampling ChE210D Today's lecture: methods for facilitating equilibration and sampling in complex, frustrated, or slow-evolving systems Difficult-to-simulate systems Practically speaking, one is
More informationEnabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package
Enabling constant pressure hybrid Monte Carlo simulations using the GROMACS molecular simulation package Mario Fernández Pendás MSBMS Group Supervised by Bruno Escribano and Elena Akhmatskaya BCAM 18 October
More informationIV. Classical Molecular Dynamics
IV. Classical Molecular Dynamics Basic Assumptions: 1. Born-Oppenheimer Approximation 2. Classical mechanical nuclear motion Unavoidable Additional Approximations: 1. Approximate potential energy surface
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationComputer simulation methods (1) Dr. Vania Calandrini
Computer simulation methods (1) Dr. Vania Calandrini Why computational methods To understand and predict the properties of complex systems (many degrees of freedom): liquids, solids, adsorption of molecules
More information4. The Green Kubo Relations
4. The Green Kubo Relations 4.1 The Langevin Equation In 1828 the botanist Robert Brown observed the motion of pollen grains suspended in a fluid. Although the system was allowed to come to equilibrium,
More informationarxiv: v1 [cond-mat.stat-mech] 28 Mar 2008
Canonical sampling through velocity-rescaling Giovanni Bussi, Davide Donadio, and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationAnalysis of MD Results Using Statistical Mechanics Methods. Molecular Modeling
Analysis of MD Results Using Statistical Mechanics Methods Ioan Kosztin eckman Institute University of Illinois at Urbana-Champaign Molecular Modeling. Model building. Molecular Dynamics Simulation 3.
More informationFirst-principles Molecular Dynamics Simulations
First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul
More informationA (short) practical introduction to kinetic theory and thermodynamic properties of gases through molecular dynamics
A (short) practical introduction to kinetic theory and thermodynamic properties of gases through molecular dynamics Miguel A. Caro mcaroba@gmail.com March 28, 2018 Contents 1 Preface 3 2 Review of thermodynamics
More informationTHERMOSTATS AND THERMAL TRANSPORT IN SOLIDS
Hands-On Tutorial Workshop, July 19 th 2011 THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS Christian Carbogno FRITZ-HABER-INSTITUT DER MAX-PLANCK-GESELLSCHAFT, BERLIN - GERMANY Hands-On Tutorial Workshop,
More informationQUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES
QUANTUM AND THERMAL MOTION IN MOLECULES FROM FIRST-PRINCIPLES 1 Tapio T. Rantala, Department of Physics, Tampere University of Technology http://www.tut.fi/semiphys CONTENTS MOTIVATION PATH INTEGRAL APPROACH
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 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 informationLECTURE 6 : BASICS FOR MOLECULAR SIMULATIONS - Historical perspective - Skimming over Statistical Mechanics - General idea of Molecular Dynamics -
LECTURE 6 : BASICS FOR MOLECULAR SIMULATIONS - Historical perspective - Skimming over Statistical Mechanics - General idea of Molecular Dynamics - Force calculations, structure of MD, equations of motion
More information23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation M I
23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation 1. Now we will write down the Hamiltonian for a molecular system comprising N nuclei and n electrons.
More informationIntroduction to atomic scale simulations
Powder Technology course Autumn semester 2017 Introduction to atomic scale simulations Peter M Derlet Condensed Matter Theory Paul Scherrer Institut peter.derlet@psi.ch Outline of lecture The big picture
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 informationNext generation extended Lagrangian first principles molecular dynamics
Next generation extended Lagrangian first principles molecular dynamics Anders M. N. Niklasson Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 (Dated: June 1, 017) LA-UR-17-3007
More information7 To solve numerically the equation of motion, we use the velocity Verlet or leap frog algorithm. _ V i n = F i n m i (F.5) For time step, we approxim
69 Appendix F Molecular Dynamics F. Introduction In this chapter, we deal with the theories and techniques used in molecular dynamics simulation. The fundamental dynamics equations of any system is the
More informationIntro to ab initio methods
Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio
More informationMethods of Computer Simulation. Molecular Dynamics and Monte Carlo
Molecular Dynamics Time is of the essence in biological processes therefore how do we understand time-dependent processes at the molecular level? How do we do this experimentally? How do we do this computationally?
More informationBasics of Statistical Mechanics
Basics of Statistical Mechanics Review of ensembles Microcanonical, canonical, Maxwell-Boltzmann Constant pressure, temperature, volume, Thermodynamic limit Ergodicity (see online notes also) Reading assignment:
More informationQuantum Molecular Dynamics Basics
Quantum Molecular Dynamics Basics Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationHeat Transport in Glass-Forming Liquids
Heat Transport in Glass-Forming Liquids by VINAY VAIBHAV The Institute of Mathematical Sciences CIT Campus, Taramani, Chennai 600 113, India. A thesis submitted in partial fulfillment of requirements for
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 informationSupporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure
Supporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure Bladimir Ramos-Alvarado, Satish Kumar, and G. P. Peterson The George
More informationMolecular Dynamics at Constant Pressure: Allowing the System to Control Volume Fluctuations via a Shell Particle
Entropy 2013, 15, 3941-3969; doi:10.3390/e15093941 Review OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Molecular Dynamics at Constant Pressure: Allowing the System to Control Volume
More information2. Thermodynamics. Introduction. Understanding Molecular Simulation
2. Thermodynamics Introduction Molecular Simulations Molecular dynamics: solve equations of motion r 1 r 2 r n Monte Carlo: importance sampling r 1 r 2 r n How do we know our simulation is correct? Molecular
More informationNanoscale simulation lectures Statistical Mechanics
Nanoscale simulation lectures 2008 Lectures: Thursdays 4 to 6 PM Course contents: - Thermodynamics and statistical mechanics - Structure and scattering - Mean-field approaches - Inhomogeneous systems -
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 informationThermodynamics of Solids: Harmonic and Quasi-harmonic Approximations
Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations, USA, July 9-14, 2017 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it 2017 Outline -
More informationRandom Walks A&T and F&S 3.1.2
Random Walks A&T 110-123 and F&S 3.1.2 As we explained last time, it is very difficult to sample directly a general probability distribution. - If we sample from another distribution, the overlap will
More informationStatistical methods in atomistic computer simulations
Statistical methods in atomistic computer simulations Prof. Michele Ceriotti, michele.ceriotti@epfl.ch This course gives an overview of simulation techniques that are useful for the computational modeling
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationA Brief Introduction to Statistical Mechanics
A Brief Introduction to Statistical Mechanics E. J. Maginn, J. K. Shah Department of Chemical and Biomolecular Engineering University of Notre Dame Notre Dame, IN 46556 USA Monte Carlo Workshop Universidade
More informationIntroduction to model potential Molecular Dynamics A3hourcourseatICTP
Introduction to model potential Molecular Dynamics A3hourcourseatICTP Alessandro Mattoni 1 1 Istituto Officina dei Materiali CNR-IOM Unità di Cagliari SLACS ICTP School on numerical methods for energy,
More informationCPMD Tutorial Atosim/RFCT 2009/10
These exercices were inspired by the CPMD Tutorial of Axel Kohlmeyer http://www.theochem.ruhruni-bochum.de/ axel.kohlmeyer/cpmd-tutor/index.html and by other tutorials. Here is a summary of what we will
More informationTime-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
More informationComputational Physics
Computational Physics Molecular Dynamics Simulations E. Carlon, M. Laleman and S. Nomidis Academic year 015/016 Contents 1 Introduction 3 Integration schemes 4.1 On the symplectic nature of the Velocity
More informationElectronic Structure Methods for molecular excited states
Electronic Structure Methods for molecular excited states Summer School on Time dependent Density-Functional Theory: Prospects and Applications Benasque, Spain September 2008 Ivano Tavernelli EPFL CH-1015
More informationMolecular Dynamics. Park City June 2005 Tully
Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationPrinciples of Equilibrium Statistical Mechanics
Debashish Chowdhury, Dietrich Stauffer Principles of Equilibrium Statistical Mechanics WILEY-VCH Weinheim New York Chichester Brisbane Singapore Toronto Table of Contents Part I: THERMOSTATICS 1 1 BASIC
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 informationBrief Review of Statistical Mechanics
Brief Review of Statistical Mechanics Introduction Statistical mechanics: a branch of physics which studies macroscopic systems from a microscopic or molecular point of view (McQuarrie,1976) Also see (Hill,1986;
More informationMolecular dynamics: Car-Parrinello method
Formulations Potential energy Initialization Verlet algorithm MD: Steps MD: Thermo and barostats CP: Car-Parrinello Atomic units MD and CP textbo Molecular dynamics: Car-Parrinello method Víctor Luaña
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