Ab Ini'o Molecular Dynamics (MD) Simula?ons

Size: px
Start display at page:

Download "Ab Ini'o Molecular Dynamics (MD) Simula?ons"

Transcription

1 Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA

2 What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties of many-body systems Given initial condition (positions & momenta of N particles/ nuclei), evolve system in time following equation of motion

3 Why MD Simulations? We want to do statistical mechanics! Compute observables to connect with experiment Thermodynamics, structure, dynamics, electronic properties

4 Why MD Simulations? We want to do statistical mechanics! Compute observables to connect with experiment ho(r)i = R dro(r)e Z 3N-dimensional configurational integrals U(R) P hoi = n O ne Q Summation over all states " n Partition Function: Q = X n e " n H n = " n n

5 How to compute observables? Work under assumption of the Ergodic Hypothesis: Z hoi lim!1 0 dto(t)

6 How to compute observables? Work under assumption of the Ergodic Hypothesis: Z hoi lim!1 0 dto(t)

7 How to compute observables? Work under assumption of the Ergodic Hypothesis: Z hoi lim!1 0 dto(t) Long enough trajectories will adequately sample phase space

8 Simulating the Time Evolution Need to evaluate the forces on particles/nuclei Classical MD uses empirical potentials Numerous approaches exist for ab initio MD Ehrenfest dynamics Car-Parrinello MD (CPMD) Born-Oppenheimer MD (BOMD)

9 Born-Oppenheimer MD Work within the Born-Oppenheimer approximation Separation of nuclear & electronic DoFs Nuclei evolve in time Electronic structure problems is solution to time-independent Schrodinger equation: static problem BOMD equation of motion: m i r i (t) = r ri min 0 {h 0 H e 0i} H e 0 = E 0 0

10 Born-Oppenheimer MD Work within the Born-Oppenheimer approximation Separation of nuclear & electronic DoFs Nuclei evolve in time Electronic structure problems is solution to time-independent Schrodinger equation: static problem BOMD equation of motion: m i r i (t) = r ri min 0 {h 0 H e 0i} H e 0 = E 0 0 F

11 Integrating the equations of motion: the velocity Verlet algorithm Compute forces Update positions Update velocities Compute new Forces Update velocities Global error O( t 2 ) Need to carefully choose t

12 Exercise 1, Timesteps Run MD simulations of 8-atom liquid Si cell with t=0.5-3 fs Simulations in microcanonical (NVE) ensemble Energy is constant =0 Integration scheme should conserve energy, error in discrete integration will lead to buildup of errors

13 Results: Exercise 1, Timesteps Energy Fluctuations (and drift) increase with size of timestep Smaller t conserves E better Tradeoff between accuracy & efficiency

14 Exercise 2, Energy Convergence BOMD: Energy (and forces) depend on minimization of energy functional Sensitive to convergence criteria for self-consistent iteration Perform NVE BOMD for convergences of 1E-4 to 1E-7 ev for energy difference between iterations

15 Results: Exercise 2, Energy Convergence Energy drift increases significantly with decreasing tolerance Tradeoff between accuracy & efficiency

16 Exercise 3, Averages Compute average thermodynamic, structural, and dynamic observables Temperature in NVE ensemble Pair correlation functions Diffusion CoefFicient: Green-Kubo vs. Einstein relations

17 Results: Exercise 3a, Temperature <T> ~ 3000 K

18 Results: Exercise 3b, Structure Disordered liquid, ordered solid

19 Finite size effects in l-si

20 Exercise 3c, Dynamics MD also enables computation of dynamic properties Green-Kubo: Transport coefficients related to integral of time correlation function (TCF) Z 1 D = dt hv(0) v(t)i

21 Exercise 3c, Dynamics MD also enables computation of dynamic properties Green-Kubo: Transport coefficients related to integral of time correlation function (TCF) Z 1 D = dt hv(0) v(t)i Einstein: 6Dt = lim t!1 h r(t) r(0) i Should yield equivalent D values

22 Results: Exercise 3c, Dynamics from TCFs SCAN DSCAN ~ 0.4 Å 2 / ps DLDA ~ 1.4 Å 2 / ps LDA Green-Kubo type relations typically take longer to converge than MSD, for example

23 Results: Exercise 3d, Dynamics from MSD DSCAN ~ 0.76 Å 2 / ps SCAN DLDA ~ 1.35 Å 2 / ps Consistent with VACF results Faster dynamics in LDA LDA due to lack of good covalent bonds

24 Results: Exercise 3e, Phonon Density of States SCAN Phonon DOS is Fourier transform of velocity autocorrelation function, C(t) LDA Qualitatively similar Need better statistics Can compare with solid in phonon tutorial

25 Exercise 4, Canonical (NVT) Ensemble Perform simulations at constant temperature Use Andersen thermostat (preserves canonical distribution) Mimics coupling to heat bath (with a prescribed collision frequency) Particle momenta are reset randomly, following Maxwell- Boltzmann distribution Thermodynamics? Dynamics???

26 Results: Exercise 4, Canonical (NVT) Ensemble E is not the conserved quantity Can Fluctuate T Fluctuates about constant <T> Why isn t T constant? Small system, constant T as N tends to infinity

27 Results: Exercise 4, Canonical (NVT) Ensemble Structure is independent of ensemble

28 Results: Exercise 4, Canonical (NVT) Ensemble Dynamics are severely altered by thermostat Extreme case, others thermostats (and weak coupling) may not significantly alter dynamics Be careful!

29 SCAN description of silicon with MD

30 Background: Metallic and Covalent Bonding in Solid Silicon Si transitions from semiconducting diamond to metallic ß-Sn at P~12 GPa SCAN captures structure and energetics with quantitative accuracy Distinguishes J. Sun, et al. Nature Chem covalent from metallic bonds

31 Phase Transitions Coexistence: semiconducting solid and metallic liquid Tm ~ 1550 K, closer to Exp. (1680 K) than previous estimates using LDA: K RCR, J. Sun, M. L. Klein, to be submitted

32 Mixed Interactions: Liquid Silicon l-si is largely metallic; some covalent character 2nd peak is absent in LDA/ PBE descriptions LDA/PBE lack adequate discrimination between metallic and covalent bonding J. Sun, et al. Nature Chem. 2016

33 More phase transitions Thermodynamics Dynamics RCR, J. Sun, M. L. Klein, to be submitted

34 Changes in Electronic Structure NF DOS follows density dependence on T Transition from covalent and metallic bonding dominating RCR, J. Sun, M. L. Klein, to be submitted Covalent Heat Metallic

35 Q & A

MD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky

MD 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 information

Ab initio molecular dynamics and nuclear quantum effects

Ab 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 information

Ab initio Molecular Dynamics Born Oppenheimer and beyond

Ab 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 information

What is Classical Molecular Dynamics?

What 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 information

Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014

Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014 Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid

More information

Accelerated Quantum Molecular Dynamics

Accelerated 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 information

Energy and Forces in DFT

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 information

Density Functional Theory: from theory to Applications

Density 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 information

Introduction to molecular dynamics

Introduction 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 information

Gear methods I + 1/18

Gear 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 information

Ab initio molecular dynamics

Ab 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 information

Understanding Molecular Simulation 2009 Monte Carlo and Molecular Dynamics in different ensembles. Srikanth Sastry

Understanding 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 information

Ab initio molecular dynamics: Basic Theory and Advanced Met

Ab 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 information

Ab initio molecular dynamics simulation on temperature-dependent properties of Al Si liquid alloy

Ab initio molecular dynamics simulation on temperature-dependent properties of Al Si liquid alloy INSTITUTE OF PHYSICSPUBLISHING JOURNAL OFPHYSICS: CONDENSED MATTER J. Phys.: Condens. Matter 16 (4) 57 514 PII: S953-8984(4)7691-8 Ab initio molecular dynamics simulation on temperature-dependent properties

More information

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC

Structure 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 information

4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania

4/18/2011. Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania 1. Introduction Titus Beu University Babes-Bolyai Department of Theoretical and Computational Physics Cluj-Napoca, Romania Bibliography Computer experiments Ensemble averages and time averages Molecular

More information

Ab initio molecular dynamics

Ab 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 information

Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases

Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song

More information

Quantum (Path Integral) Molecular Dynamics

Quantum (Path Integral) Molecular Dynamics Quantum (Path Integral) 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 Background

More information

Computer simulation methods (1) Dr. Vania Calandrini

Computer 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 information

Simulation of molecular systems by molecular dynamics

Simulation 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 information

Classical Molecular Dynamics

Classical 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 information

Introduction 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 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 information

Temperature and Pressure Controls

Temperature 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 information

Ab initio Molecular Dynamics

Ab 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 information

Temperature and Pressure Controls

Temperature 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 information

Molecular Dynamics Simulations

Molecular 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 information

Set the initial conditions r i. Update neighborlist. r i. Get new forces F i

Set 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 information

A 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 information

A Brief Introduction to Ab Initio Molecular Dynamics

A 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 information

André Schleife Department of Materials Science and Engineering

André 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 information

THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS

THERMOSTATS 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 information

Optimization Methods via Simulation

Optimization 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 information

Ab Initio Molecular Dynamic

Ab 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 information

Langevin Dynamics in Constant Pressure Extended Systems

Langevin 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 information

Molecular 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 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 information

Analysis of the simulation

Analysis of the simulation Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities

More information

Electronic structure simulations of water solid interfaces

Electronic 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 information

Ab Initio Molecular Dynamic

Ab 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 information

Modeling of Nanostructures and Materials Jacek A. Majewski Summer Semester 2013 Lecture 6 Lecture Modeling of Nanostructures Molecular Dynamics

Modeling 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 information

arxiv: v1 [physics.chem-ph] 24 Apr 2018

arxiv: v1 [physics.chem-ph] 24 Apr 2018 arxiv:1804.08913v1 [physics.chem-ph] 24 Apr 2018 Fast-Forward Langevin Dynamics with Momentum Flips Mahdi Hijazi, 1 David M. Wilkins, 1, a) and Michele Ceriotti 1 Laboratory of Computational Science and

More information

Modeling & Simulation of Glass Structure

Modeling & 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 information

Generalized valence bond molecular dynamics at constant temperature

Generalized valence bond molecular dynamics at constant temperature M o l e c u l a r Ph y s i c s, 1996, Vo l. 89, No. 5, 1265± 1276 Generalized valence bond molecular dynamics at constant temperature By DOUGLAS A. GIBSON and EM ILY A. CARTER D epartment of Chemistry

More information

Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate)

Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining. Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) Supporting Information for Glass-Transition and Side-Chain Dynamics in Thin Films: Explaining Dissimilar Free Surface Effects for Polystyrene and Poly(methyl methacrylate) David D. Hsu, Wenjie Xia, Jake

More information

Introduction Statistical Thermodynamics. Monday, January 6, 14

Introduction 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 information

FUNDAMENTAL ISSUES IN NANOSCALE HEAT TRANSFER: FROM COHERENCE TO INTERFACIAL RESISTANCE IN HEAT CONDUCTION

FUNDAMENTAL ISSUES IN NANOSCALE HEAT TRANSFER: FROM COHERENCE TO INTERFACIAL RESISTANCE IN HEAT CONDUCTION HEFAT2014 10 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics 14 16 July 2014 Orlando, Florida FUNDAMENTAL ISSUES IN NANOSCALE HEAT TRANSFER: FROM COHERENCE TO INTERFACIAL

More information

arxiv: v1 [cond-mat.stat-mech] 28 Mar 2008

arxiv: 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 information

Computer 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 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

MD simulation: output

MD 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 information

Statistical Mechanics

Statistical Mechanics Franz Schwabl Statistical Mechanics Translated by William Brewer Second Edition With 202 Figures, 26 Tables, and 195 Problems 4u Springer Table of Contents 1. Basic Principles 1 1.1 Introduction 1 1.2

More information

Comparative study on methodology in molecular dynamics simulation of nucleation

Comparative study on methodology in molecular dynamics simulation of nucleation THE JOURNAL OF CHEMICAL PHYSICS 126, 224517 2007 Comparative study on methodology in molecular dynamics simulation of nucleation Jan Julin, Ismo Napari, and Hanna Vehkamäki Department of Physical Sciences,

More information

The micro-properties of [hmpy+] [Tf 2 N-] Ionic liquid: a simulation. study. 1. Introduction

The 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 information

Practical calculations using first-principles QM Convergence, convergence, convergence

Practical calculations using first-principles QM Convergence, convergence, convergence Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................

More information

Phonon calculations with SCAN

Phonon calculations with SCAN Workshop on the SCAN density functional: Fundamentals, practices, and extensions Temple university, Philadelphia May 18th, 2017 Hands-on tutorial 3 Phonon calculations with SCAN Yubo Zhang and Jianwei

More information

Molecular Dynamics. Park City June 2005 Tully

Molecular 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 information

Methods of Computer Simulation. Molecular Dynamics and Monte Carlo

Methods 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 information

Basics of Statistical Mechanics

Basics 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 information

Introduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou

Introduction 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 information

Statistical methods in atomistic computer simulations

Statistical 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 information

Thermal and Statistical Physics Department Exam Last updated November 4, L π

Thermal and Statistical Physics Department Exam Last updated November 4, L π Thermal and Statistical Physics Department Exam Last updated November 4, 013 1. a. Define the chemical potential µ. Show that two systems are in diffusive equilibrium if µ 1 =µ. You may start with F =

More information

ICONE COMPUTER SIMULATION OF DIFFUSION OF PB-BI EUTECTIC IN LIQUID SODIUM BY MOLECULAR DYNAMICS METHOD

ICONE COMPUTER SIMULATION OF DIFFUSION OF PB-BI EUTECTIC IN LIQUID SODIUM BY MOLECULAR DYNAMICS METHOD Proceedings of ICONE TH International Conference on Nuclear Engineering Arlington, VA, USA, April -, ICONE-3 COMPUTER SIMULATION OF DIFFUSION OF PB-BI EUTECTIC IN LIQUID SODIUM BY MOLECULAR DYNAMICS METHOD

More information

Molecular Simulation Background

Molecular Simulation Background Molecular Simulation Background Why Simulation? 1. Predicting properties of (new) materials 2. Understanding phenomena on a molecular scale 3. Simulating known phenomena? Example: computing the melting

More information

arxiv:cond-mat/ v1 [cond-mat.stat-mech] 8 Oct 1996

arxiv:cond-mat/ v1 [cond-mat.stat-mech] 8 Oct 1996 December 21, 2013 arxiv:cond-mat/9610066v1 [cond-mat.stat-mech] 8 Oct 1996 Some Finite Size Effects in Simulations of Glass Dynamics Jürgen Horbach, Walter Kob, Kurt Binder Institut für Physik, Johannes

More information

Javier Junquera. Statistical mechanics

Javier 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 information

Metadynamics-biased ab initio Molecular Dynamics Study of

Metadynamics-biased ab initio Molecular Dynamics Study of Metadynamics-biased ab initio Molecular Dynamics Study of Heterogeneous CO 2 Reduction via Surface Frustrated Lewis Pairs Mireille Ghoussoub 1, Shwetank Yadav 1, Kulbir Kaur Ghuman 1, Geoffrey A. Ozin

More information

Statistical Mechanics in a Nutshell

Statistical Mechanics in a Nutshell Chapter 2 Statistical Mechanics in a Nutshell Adapted from: Understanding Molecular Simulation Daan Frenkel and Berend Smit Academic Press (2001) pp. 9-22 11 2.1 Introduction In this course, we will treat

More information

Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany

Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical

More information

Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges

Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges Modifications of the Robert- Bonamy Formalism and Further Refinement Challenges Q. Ma, NASA/GISS R. H. Tipping, Univ. of Alabama C. Boulet, Univ. Paris-Sud, France Statements The presentation is not a

More information

Major Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1

Major Concepts Lecture #11 Rigoberto Hernandez. TST & Transport 1 Major Concepts Onsager s Regression Hypothesis Relaxation of a perturbation Regression of fluctuations Fluctuation-Dissipation Theorem Proof of FDT & relation to Onsager s Regression Hypothesis Response

More information

A Brief Introduction to Statistical Mechanics

A 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 information

Molecular dynamics simulation of Aquaporin-1. 4 nm

Molecular 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 information

Time reversible Born Oppenheimer molecular dynamics

Time 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 information

LECTURE 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 - 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 information

IV. Classical Molecular Dynamics

IV. 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 information

André Schleife Department of Materials Science and Engineering

André 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 information

Principles of Equilibrium Statistical Mechanics

Principles 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 information

Computational Chemistry - MD Simulations

Computational 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 information

Stress analysis of lithium-based rechargeable batteries using micro and macro scale analysis

Stress analysis of lithium-based rechargeable batteries using micro and macro scale analysis Stress analysis of lithium-based rechargeable batteries using micro and macro scale analysis Utsav Kumar Atanu K. Metya Jayant K. Singh Department of Chemical Engineering IIT Kanpur INTRODUCTION Christensen

More information

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions

Exercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,

More information

Computer simulation methods (2) Dr. Vania Calandrini

Computer simulation methods (2) Dr. Vania Calandrini Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing

More information

Rustam Z. Khaliullin University of Zürich

Rustam Z. Khaliullin University of Zürich Rustam Z. Khaliullin University of Zürich Molecular dynamics (MD) MD is a computational method for simulating time evolution of a collection of interacting atoms by numerically integrating Newton s equation

More information

Hands-on : Model Potential Molecular Dynamics

Hands-on : Model Potential Molecular Dynamics Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration

More information

Molecular and Electronic Dynamics Using the OpenAtom Software

Molecular and Electronic Dynamics Using the OpenAtom Software Molecular and Electronic Dynamics Using the OpenAtom Software Sohrab Ismail-Beigi (Yale Applied Physics) Subhasish Mandal (Yale), Minjung Kim (Yale), Raghavendra Kanakagiri (UIUC), Kavitha Chandrasekar

More information

NPTEL

NPTEL NPTEL Syllabus Nonequilibrium Statistical Mechanics - Video course COURSE OUTLINE Thermal fluctuations, Langevin dynamics, Brownian motion and diffusion, Fokker-Planck equations, linear response theory,

More information

Once Upon A Time, There Was A Certain Ludwig

Once Upon A Time, There Was A Certain Ludwig Once Upon A Time, There Was A Certain Ludwig Statistical Mechanics: Ensembles, Distributions, Entropy and Thermostatting Srinivas Mushnoori Chemical & Biochemical Engineering Rutgers, The State University

More information

Ab initio phonon calculations in mixed systems

Ab 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 information

424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393

424 Index. Eigenvalue in quantum mechanics, 174 eigenvector in quantum mechanics, 174 Einstein equation, 334, 342, 393 Index After-effect function, 368, 369 anthropic principle, 232 assumptions nature of, 242 autocorrelation function, 292 average, 18 definition of, 17 ensemble, see ensemble average ideal,23 operational,

More information

Advanced in silico drug design

Advanced in silico drug design Advanced in silico drug design RNDr. Martin Lepšík, Ph.D. Lecture: Advanced scoring Palacky University, Olomouc 2016 1 Outline 1. Scoring Definition, Types 2. Physics-based Scoring: Master Equation Terms

More information

Supporting Information for. Dynamics Study"

Supporting Information for. Dynamics Study Supporting Information for "CO 2 Adsorption and Reactivity on Rutile TiO 2 (110) in Water: An Ab Initio Molecular Dynamics Study" Konstantin Klyukin and Vitaly Alexandrov,, Department of Chemical and Biomolecular

More information

Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions

Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions Electron-phonon scattering from green s function transport combined with molecular dynamics: Applications to mobility predictions Daniele Stradi www.quantumwise.com daniele.stradi@quantumwise.com Introduction

More information

Interface Resistance and Thermal Transport in Nano-Confined Liquids

Interface Resistance and Thermal Transport in Nano-Confined Liquids 1 Interface Resistance and Thermal Transport in Nano-Confined Liquids Murat Barisik and Ali Beskok CONTENTS 1.1 Introduction...1 1.2 Onset of Continuum Behavior...2 1.3 Boundary Treatment Effects on Interface

More information

Statistical thermodynamics (mechanics)

Statistical thermodynamics (mechanics) Statistical thermodynamics mechanics) 1/15 Macroskopic quantities are a consequence of averaged behavior of many particles [tchem/simplyn.sh] 2/15 Pressure of ideal gas from kinetic theory I Molecule =

More information

Supporting Online Materials: Nature of proton. transport in a water-filled carbon nanotube and in. liquid water

Supporting Online Materials: Nature of proton. transport in a water-filled carbon nanotube and in. liquid water Supporting Online Materials: Nature of proton transport in a water-filled carbon nanotube and in liquid water Ji Chen, Xin-Zheng Li,,, Qianfan Zhang, Angelos Michaelides, and Enge Wang, ICQM and School

More information

Understanding hydrogen storage in metal organic frameworks from first principles

Understanding hydrogen storage in metal organic frameworks from first principles Understanding hydrogen storage in metal organic frameworks from first principles Understanding hydrogen storage in metal organic frameworks from first principles Sohrab Ismail-Beigi (Yale Applied Physics)

More information

Bioengineering 215. An Introduction to Molecular Dynamics for Biomolecules

Bioengineering 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 information

Molecular-dynamics study of surface segregation in liquid semiconductor alloys

Molecular-dynamics study of surface segregation in liquid semiconductor alloys PHYSICAL REVIEW B VOLUME 56, NUMBER 19 15 NOVEMBER 1997-I Molecular-dynamics study of surface segregation in liquid semiconductor alloys Wenbin Yu and D. Stroud* Department of Physics, The Ohio State University,

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statistical Mechanics Notes for Lecture 15 Consider Hamilton s equations in the form I. CLASSICAL LINEAR RESPONSE THEORY q i = H p i ṗ i = H q i We noted early in the course that an ensemble

More information

CH 240 Chemical Engineering Thermodynamics Spring 2007

CH 240 Chemical Engineering Thermodynamics Spring 2007 CH 240 Chemical Engineering Thermodynamics Spring 2007 Instructor: Nitash P. Balsara, nbalsara@berkeley.edu Graduate Assistant: Paul Albertus, albertus@berkeley.edu Course Description Covers classical

More information

Computational Physics

Computational 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 information

Chem 350: Statistical Mechanics and Chemical Kinetics. Spring Preface. Introduction 2

Chem 350: Statistical Mechanics and Chemical Kinetics. Spring Preface. Introduction 2 Preface Introduction 2 Statistical Mechanics and Chemical Kinetics: Syllabus Textbook: Thermodynamics, Statistical Thermodynamics and Kinetics 3 rd ed by Thomas Engel and Philip Reid Additional Resource

More information