Monte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects
|
|
- Evangeline Floyd
- 6 years ago
- Views:
Transcription
1 Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects
2 Lecture 1 2
3 Lecture 1&2 3
4 Lecture 1&3 4
5 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT Isobaric-isothermal,P,T V µvt Grand-canonical µ,v,t Fluctuating (Measured) 5
6 VT Liquids Equation of State of Liquid Carbon Imposed 6
7 if force PT is difficult to calculate e.g. carbon force field Liquids Equation of State of Liquid Carbon Imposed 7
8 PT on-isotropic Systems e.g. Solids Structure and Transformation of Carbon anotube Arrays 8
9 µvt Adsorption Adsorption in Carbon anostructues Schwegler et al. 9
10 Partition function Q VT Statistical Thermodynamics 1 = Λ! dr exp β 3 Ensemble average 1 1 A = VT Q Λ! VT [ U( r )] ( ) [ ( )] r exp U r dr A β 3 Probability to find a particular configuration ( ) 1 1 r dr' ( r' r ) exp ( r' ) exp ( = δ βu βu r ) 3 Q Λ! VT Lecture 2 Free energy βf = ln ( ) Q VT 10
11 Ensemble average Lecture 2 A VT 1 = Q Λ VT = dr = dr 1 dr A( r ) exp[ βu ( r )] 3! ( ) ( dr A( r ) P( r ) r P r ) = dr P( r ) ( r ) C exp[ U ( r )] dr A = C exp[ βu ( r )] dr A A dr Generate configuration using MC: {,r,r,r,r } P MC ( ) [ ( )] MC r = C exp βu r β ( r ) exp[ βu ( r )] exp[ βu ( r )] M 1 r ! A = A( ri ) M M i= 1 Weighted Distribution with = = = dr dr dr A dr MC ( r ) P ( r ) MC P ( r ) A dr A dr MC ( r ) C exp[ βu ( r )] MC C exp[ βu ( r )] ( r ) exp[ βu ( r )] exp[ βu ( r )]
12 Monte Carlo: Detailed balance Lecture 2 o n Ko ( n) = Kn ( o) Ko ( n) = o ( ) α( o n) acc( o n) Kn ( o) = n ( ) α( n o) acc( n o) acc(o n) acc(n o) (n) α(n o) = (o) α(o n) = (n) (o)
13 Lecture 2 VT-ensemble n () exp acc( o n) ( n) = acc( n o) ( o) ( ) βu n acc( o n) acc( n o) = exp β U n ( ) U( o) 13
14 PT ensemble We control the Temperature (T) Pressure (P) umber of particles () 14
15 Scaled coordinates Intermezzo Partition function Q VT 1 = Λ! Scaled coordinates dr exp β 3 si = ri / L [ U( r )] This gives for the partition function 3 L QVT = ds exp βu 3 Λ! s ; L V = ds exp βu 3 Λ! s ; ( ) ( L) The energy depends on the real coordinates 15
16 The PT ensemble Here they are an ideal gas V 0 : total volume M : total number of particles in volume V Here they interact M- in volume V 0 -V V 0 is fixed V varies from 0 to V 0 What is the statistical thermodynamics of this ensemble?
17 The PT ensemble: Fixed V partition function V QVT = ds exp βu s ; L 3 Λ! Q MV0,V,T = Q MV0,V,T = ( V 0 V ) M 3(M ) Λ M ( )! ( V 0 V ) M 3(M ) Λ M ( )! ( )! " ( ) ds M exp βu 0 s M ; L V Λ 3! ds ds V # $ Λ 3! exp! " βu s ; L exp$ % βu s ; L ( ) ( ) & ' 17 # $
18 Q MV0,V,T = ( V 0 V ) M 3(M ) Λ M ( )! V Λ 3! ds exp$ % βu s ; L ( ) To get the Partition Function of this system, we have to integrate over all possible volumes: ( ) M Q = dv V 0 V V MV0,,T 3(M ) Λ ( M )! Λ 3! ow let us take the following limits: ds & ' exp$ % βu s ; L ( ) M M ρ = constant V0 V As the particles are an ideal gas in the big reservoir we have: & ' ρ = βp 18
19 ( ) M Q = dv V 0 V V MV0,,T 3(M ) Λ ( M )! Λ 3! We have ds exp$ % βu s ; L ( ) & ' M M V V = V 1 V V V exp M V V ( ) M ( ) M ( ) M M M ( V V) V exp[ ρv] = V exp[ βpv] This gives: To make the partition function dimensionless (see Frenkel/Smit for more info) β P QPT = dvexp PV V ds exp U s ; L 3! Λ [ ] ( β β ) 19
20 Partition function: PT Ensemble β P QPT = dvexp PV V ds exp U s ; L 3! Λ [ ] ( β β ) Probability to find a particular configuration: (, s ) exp[ β ] exp β ( s ; ) PT V V PV U L Sample a particular configuration: change of volume change of reduced coordinates Acceptance rules?? Detailed balance 20
21 Detailed balance o n Ko ( n) = Kn ( o) Ko ( n) = o ( ) α( o n) acc( o n) Kn ( o) = n ( ) α( n o) acc( n o) acc( o n) ( n) α( n o) ( n) = = acc( n o) ( o) α( o n) ( o) 21
22 acc( o n) ( n) = acc( n o) ( o) PT-ensemble (, s ) exp[ ] exp ( s ; ) PT V V βpv βu L Suppose we change the position of a randomly selected particle acc( o n) V exp PV exp U s ; L = acc( n o) [ ] ( V exp βpv exp βu s ) o ; L ( exp βu s ) n ; L = = exp β U n U o ( exp βu s ) o ; L 22 [ ] ( β β ) n { ( ) ( ) }
23 acc( o n) ( n) = acc( n o) ( o) PT-ensemble (, s ) exp[ ] exp ( s ; ) PT V V βpv βu L Suppose we change the volume of the system acc( o n) Vn exp PVn exp U s ; L n = acc( n o) [ ] ( V exp exp s ; ) o βpv o βu L o [ ] ( β β ) V n = exp βpv ( ) exp ( ) ( 0) n Vo U n U V o 23 { β }
24 Algorithm: PT Randomly change the position of a particle Randomly change the volume 24
25 25
26 26
27 27
28 PT simulations Equation of State of Lennard Jones System 28
29 Measured and Imposed Pressure Imposed pressure P Measured pressure <P> from virial P = < F dvv e βpv ds e βu (s ) # F % # & $ V % ( >,T = $ V ' dvv e βpv ds e βu (s ) exp*+ β ( F(V )+ PV ), - p(v ) = Q PT = βp Q PT dv exp*+ β ( F(V )+ PV ), - & ( ',T 29
30 # P = < F & % ( >,T = $ V ' P = P = βp Q(PT ) βp Q(PT ) # dv F & % ( $ V ' dv dvv e βpv ds e βu (s ) # F % $ V dvv e βpv ds e βu (s ),T exp[ βpv ] β exp* + β F V ( ) + PV ( ) exp*+ βf V, - V ( ), - & ( ',T 30
31 Measured and Imposed Pressure Partial integration For V=0 and V= b a fdg = [ ] b fg a exp" # β F V b a gdf ( ( ) + PV ) $ % = 0 Therefore, P = P = βp Q(PT ) βp Q(PT ) = dv exp[ βpv ] β dvp exp! " β F V ( ) exp!" βf V # $ V ( ( ) + PV ) # $ = P 31
32 Grand-canonical ensemble What are the equilibrium conditions? 32
33 Grand-canonical ensemble We impose: Temperature (T) Chemical potential (µ) Volume (V) But OT pressure 33
34 The ensemble of the total system Here they are an ideal gas Here they interact What is the statistical thermodynamics of this ensemble? 34
35 The ensemble: partition function V QVT = ds exp βu s ; L 3 Λ! Q MV0,V,T = ( V 0 V ) M 3(M ) Λ M ( )! ( ) M V0 V V MV V T M ( ) ( )! " ( ) ds M V exp βu 0 s M ; L ds V # $ Λ 3! exp! " βu s ; L ( ) Q 0,, = ds exp βu s ; L 3 3 Λ M! Λ! ( ) 35 # $
36 Q = ( V 0 V ) M V MV0 ds exp$ βu ( s ; L) &,V,T 3(M ) Λ ( M )! Λ 3! % ' To get the Partition Function of this system, we have to sum over all possible number of particles Q = ( V 0 V ) =M M V MV0 ds exp$ βu ( s ; L) &,,T 3(M ) Λ ( M )! Λ 3! % ' =0 ow let us take the following limits: M M ρ = constant V0 V As the particles are an ideal gas in the big reservoir we have: ( 3 ) µ = ktln Λ ρ B = exp( ) βµ V Qµ VT = ds exp βu s ; L 3 Λ! = 0 ( ) 36
37 Partition function: µvt Ensemble = ( ) exp βµ V Qµ VT = ds exp βu s ; L 3 = 0 Λ! Probability to find a particular configuration: ( ) ( ) ( exp βµ V, ) µ VT V s exp βu s ; L 3 Λ! Sample a particular configuration: Change of the number of particles Change of reduced coordinates Acceptance rules?? Detailed balance ( ) 37
38 Detailed balance o n Ko ( n) = Kn ( o) Ko ( n) = o ( ) α( o n) acc( o n) Kn ( o) = n ( ) α( n o) acc( n o) acc( o n) ( n) α( n o) ( n) = = acc( n o) ( o) α( o n) ( o) 38
39 acc( o n) ( n) = acc( n o) ( o) µvt-ensemble ( ) ( exp βµ V, ) µ VT V s exp βu s ; L 3 Λ! Suppose we change the position of a randomly selected particle exp ( ) acc( o n) 3 Λ! acc( n o) = exp βµ V βµ V ( ) Λ 3! ( ) n L exp βu s ; ( ) o L exp βu s ; { β U( n) U( ) } = exp 0 ( ) 39
40 acc( o n) ( n) = acc( n o) ( o) µvt-ensemble ( ) ( exp βµ V, ) µ VT V s exp βu s ; L 3 Λ! Suppose we change the number of particles of the system ( βµ ( + )) Λ ( + ) exp( βµ ) V ( ) + 1 exp 1 V exp βu s ; acc( o n) 1! = acc( n o) ( exp βu s ; L ) 3 o Λ! exp( βµ ) V = exp[ βδu ] 3 Λ + 1 ( ) ( +1 L ) n 40
41 41
42 42
43 Application: equation of state of Lennard-Jones 43
44 Application: adsorption in zeolites 44
45 Summary Ensemble Constant (Imposed) Fluctuating (Measured) Function VT,V,T P βf=-lnq(,v,t) PT,P,T V βg=-lnq(,p,t) µvt µ,v,t βω=-lnq(µ,v,t)=-βpv 45
46 Practical Points Boundaries CPU saving methods Reduced units Long ranged forces
47 Boundary effects In small systems, boundary effects are always large atoms in a simple cubic crystal 488 boundary atoms atoms in a simple cubic crystal still 6% boundary atoms.
48 Periodic boundary conditions
49 ener(x,en) energies)
50 Energy evaluation costs! The most time-consuming part of any simulation is the evaluation of all the interactions between the molecules. In general: (-1)/2 = O( 2 ) But often, intermolecular forces have a short range: Therefore, we do not have to consider interactions with far- away atoms.
51 Saving CPU Cell list Verlet List Han sur Lesse
52 Application: Lennard Jones u LJ u r *# σ & ( r) = 4ε $ % r ' ( +, " # $ ( ) = u LJ potential The Lennard-Jones potential 12 # σ $ % r ( r) r r c 0 r > r c 6 & ' ( The truncated Lennard-Jones potential The truncated and shifted Lennard-Jones potential -./ ( ) = u LJ u r # $ % ( r) u LJ ( r ) c r r c 0 r > r c
53 ener(x,en) energies)
54 Phase diagrams of Lennard Jones fluids
55 Long ranged interactions Long-ranged forces require special techniques. Coulomb interaction (1/r in 3D) Dipolar interaction (1/r 3 in 3D)...and, in a different context: Interactions through elastic stresses (1/r in 3D) Hydrodynamic interactions (1/r in 3D)
56 Reduced units Example: Particles with mass m and pair potential: Unit of length: Unit of energy: σ ε Unit of time:
57 Beyond standard MC on Boltzmann Sampling Lecture 2 Parallel tempering More to Come Tomorrow (Daan Frenkel)
58 Parallel tempering/replica Exchange Ergodicity problems can occur, especially in glassy systems: biomolecules, molecular glasses, gels, etc. The solution: go to high temperature High barriers in energy landscape: difficult to sample E low T Barriers effectively low: easy to sample E phase space high T phase space
59 Parallel tempering/replica Exchange Simulate two systems simultaneously system 1 system 2 temperature T 1 temperature T 2 e 1U 1 (r ) e 2U 2 (r ) total Boltzmann weight: e 1U 1 (r ) e 2U 2 (r )
60 Swap move Allow two systems to swap system 2 system 1 temperature T 1 temperature T 2 e 1U 2 (r ) e 2U 1 (r ) total Boltzmann weight: e 1U 2 (r ) e 2U 1 (r )
61 Acceptance rule The ratio of the new Boltzmann factor over the old one is (n) (o) = e( 2 1)[U 2 (r ) U 2 1(r )] The swap acceptance ratio is acc(1 2) = min 1,e ( 2 1)[U 2 (r ) U 1 (r )]
62 More replicas Consider M replica s in the VT ensemble at a different temperature. A swap between two systems of different temperatures (T i,t j ) is accepted if their potential energies are near. other parameters can be used: Hamiltonian exchange
63 Questions Lunch 63
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 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 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 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 informationMonte Carlo Methods in Statistical Mechanics
Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics
More informationMonte Carlo Molecular Simulation with Isobaric-Isothermal and Gibbs-NPT Ensembles. Thesis by Shouhong Du
Monte Carlo Molecular Simulation with Isobaric-Isothermal and Gibbs-NPT Ensembles Thesis by Shouhong Du In Partial Fulfillment of the Requirements For the Degree of Master of Science King Abdullah University
More informationPoint Defects in Hard-Sphere Crystals
6722 J. Phys. Chem. B 2001, 105, 6722-6727 Point Defects in Hard-Sphere Crystals Sander Pronk and Daan Frenkel* FOM Institute for Atomic and Molecular Physics, Kruislaan 407 1098 SJ Amsterdam, The Netherlands
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 informationExam TFY4230 Statistical Physics kl Wednesday 01. June 2016
TFY423 1. June 216 Side 1 av 5 Exam TFY423 Statistical Physics l 9. - 13. Wednesday 1. June 216 Problem 1. Ising ring (Points: 1+1+1 = 3) A system of Ising spins σ i = ±1 on a ring with periodic boundary
More informationThermodynamics at Small Scales. Signe Kjelstrup, Sondre Kvalvåg Schnell, Jean-Marc Simon, Thijs Vlugt, Dick Bedeaux
Thermodynamics at Small Scales Signe Kjelstrup, Sondre Kvalvåg Schnell, Jean-Marc Simon, Thijs Vlugt, Dick Bedeaux At small scales: 1. The systematic method of Hill 2. Small systems: surface energies 3.
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 informationHydrogels in charged solvents
Hydrogels in charged solvents Peter Košovan, Christian Holm, Tobias Richter 27. May 2013 Institute for Computational Physics Pfaffenwaldring 27 D-70569 Stuttgart Germany Tobias Richter (ICP, Stuttgart)
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 4 I. THE ISOTHERMAL-ISOBARIC ENSEMBLE The isothermal-isobaric ensemble is the closest mimic to the conditions under which most experiments are
More informationMolecular Modeling of Matter
Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy
More informationLecture 2+3: Simulations of Soft Matter. 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications
Lecture 2+3: Simulations of Soft Matter 1. Why Lecture 1 was irrelevant 2. Coarse graining 3. Phase equilibria 4. Applications D. Frenkel, Boulder, July 6, 2006 What distinguishes Colloids from atoms or
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 informationi=1 n i, the canonical probabilities of the micro-states [ βǫ i=1 e βǫn 1 n 1 =0 +Nk B T Nǫ 1 + e ǫ/(k BT), (IV.75) E = F + TS =
IV.G Examples The two examples of sections (IV.C and (IV.D are now reexamined in the canonical ensemble. 1. Two level systems: The impurities are described by a macro-state M (T,. Subject to the Hamiltonian
More information( ) Lectures on Hydrodynamics ( ) =Λ Λ = T x T R TR. Jean-Yve. We can always diagonalize point by point. by a Lorentz transformation
Lectures on Hydrodynamics Jean-Yve µν T µ ( x) 0 We can always diagonalize point by point by a Lorentz transformation T Λ µν x ( u ) µ and space-rotation R ( Ω) Then, in general, λ0 0 0 0 0 λx 0 0 Λ Λ
More informationPHYS 352 Homework 2 Solutions
PHYS 352 Homework 2 Solutions Aaron Mowitz (, 2, and 3) and Nachi Stern (4 and 5) Problem The purpose of doing a Legendre transform is to change a function of one or more variables into a function of variables
More information3.320 Lecture 18 (4/12/05)
3.320 Lecture 18 (4/12/05) Monte Carlo Simulation II and free energies Figure by MIT OCW. General Statistical Mechanics References D. Chandler, Introduction to Modern Statistical Mechanics D.A. McQuarrie,
More informationOn the Calculation of the Chemical Potential. Using the Particle Deletion Scheme
On the Calculation of the Chemical Potential Using the Particle Deletion Scheme Georgios C. Boulougouris,2, Ioannis G. Economou and Doros. Theodorou,3,* Molecular Modelling of Materials Laboratory, Institute
More informationMonte Carlo (MC) Simulation Methods. Elisa Fadda
Monte Carlo (MC) Simulation Methods Elisa Fadda 1011-CH328, Molecular Modelling & Drug Design 2011 Experimental Observables A system observable is a property of the system state. The system state i is
More informationFree Energy Estimation in Simulations
Physics 363/MatSE 38/ CSE 36 Free Energy Estimation in Simulations This is a brief overview of free energy estimation, much of which can be found in the article by D. Frenkel Free-Energy Computation and
More informationCopyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
More informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
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 informationComputer 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 informationPressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster
University of Rhode Island DigitalCommons@URI Chemistry Faculty Publications Chemistry 2005 Pressure Dependent Study of the Solid-Solid Phase Change in 38-Atom Lennard-Jones Cluster Dubravko Sabo University
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
More informationMonte Carlo. Lecture 15 4/9/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
Monte Carlo Lecture 15 4/9/18 1 Sampling with dynamics In Molecular Dynamics we simulate evolution of a system over time according to Newton s equations, conserving energy Averages (thermodynamic properties)
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 informationThe Partition Function Statistical Thermodynamics. NC State University
Chemistry 431 Lecture 4 The Partition Function Statistical Thermodynamics NC State University Molecular Partition Functions In general, g j is the degeneracy, ε j is the energy: = j q g e βε j We assume
More informationSupplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More informationPart II Statistical Physics
Part II Statistical Physics Theorems Based on lectures by H. S. Reall Notes taken by Dexter Chua Lent 2017 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More informationParallel Tempering Algorithm in Monte Carlo Simulation
Parallel Tempering Algorithm in Monte Carlo Simulation Tony Cheung (CUHK) Kevin Zhao (CUHK) Mentors: Ying Wai Li (ORNL) Markus Eisenbach (ORNL) Kwai Wong (UTK/ORNL) Metropolis Algorithm on Ising Model
More informationIntermolecular Forces and Monte-Carlo Integration 열역학특수연구
Intermolecular Forces and Monte-Carlo Integration 열역학특수연구 2003.3.28 Source of the lecture note. J.M.Prausnitz and others, Molecular Thermodynamics of Fluid Phase Equiliria Atkins, Physical Chemistry Lecture
More informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More information(# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble
Recall from before: Internal energy (or Entropy): &, *, - (# = %(& )(* +,(- Closed system, well-defined energy (or e.g. E± E/2): Microcanonical ensemble & = /01Ω maximized Ω: fundamental statistical quantity
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 informationMonte Carlo Simulations in Statistical Physics
Part II Monte Carlo Simulations in Statistical Physics By D.Stauffer Introduction In Statistical Physics one mostly deals with thermal motion of a system of particles at nonzero temperatures. For example,
More informationfiziks Institute for NET/JRF, GATE, IIT-JAM, JEST, TIFR and GRE in PHYSICAL SCIENCES
Content-Thermodynamics & Statistical Mechanics 1. Kinetic theory of gases..(1-13) 1.1 Basic assumption of kinetic theory 1.1.1 Pressure exerted by a gas 1.2 Gas Law for Ideal gases: 1.2.1 Boyle s Law 1.2.2
More informationChE 210B: Advanced Topics in Equilibrium Statistical Mechanics
ChE 210B: Advanced Topics in Equilibrium Statistical Mechanics Glenn Fredrickson Lecture 1 Reading: 3.1-3.5 Chandler, Chapters 1 and 2 McQuarrie This course builds on the elementary concepts of statistical
More informationStatistical 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 informationChemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model
Chemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model Mahfuzh Huda, 1 Siti Mariyah Ulfa, 1 Lukman Hakim 1 * 1 Department of Chemistry, Faculty of Mathematic
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
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 informationRecitation: 10 11/06/03
Recitation: 10 11/06/03 Ensembles and Relation to T.D. It is possible to expand both sides of the equation with F = kt lnq Q = e βe i If we expand both sides of this equation, we apparently obtain: i F
More informationMulticanonical parallel tempering
JOURNAL OF CHEMICAL PHYSICS VOLUME 116, NUMBER 13 1 APRIL 2002 Multicanonical parallel tempering Roland Faller, Qiliang Yan, and Juan J. de Pablo Department of Chemical Engineering, University of Wisconsin,
More informationThermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State
23 Bulletin of Research Center for Computing and Multimedia Studies, Hosei University, 28 (2014) Thermodynamics of Three-phase Equilibrium in Lennard Jones System with a Simplified Equation of State Yosuke
More informationChE 503 A. Z. Panagiotopoulos 1
ChE 503 A. Z. Panagiotopoulos 1 STATISTICAL MECHANICAL ENSEMLES 1 MICROSCOPIC AND MACROSCOPIC ARIALES The central question in Statistical Mechanics can be phrased as follows: If particles (atoms, molecules,
More informationThus, the volume element remains the same as required. With this transformation, the amiltonian becomes = p i m i + U(r 1 ; :::; r N ) = and the canon
G5.651: Statistical Mechanics Notes for Lecture 5 From the classical virial theorem I. TEMPERATURE AND PRESSURE ESTIMATORS hx i x j i = kt ij we arrived at the equipartition theorem: * + p i = m i NkT
More informationPhase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics
More information1 Fluctuations of the number of particles in a Bose-Einstein condensate
Exam of Quantum Fluids M1 ICFP 217-218 Alice Sinatra and Alexander Evrard The exam consists of two independant exercises. The duration is 3 hours. 1 Fluctuations of the number of particles in a Bose-Einstein
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 informationCalculation of the chemical potential in the Gibbs ensemble
MOLECULAR PHYSICS, 1989, VOL. 68, NO. 4, 951-958 Calculation of the chemical potential in the Gibbs ensemble by B. SMIT Koninklijke/Shell-Laboratorium, Amsterdam, (Shell Research B.V.) Badhuisweg 3, 1031
More informationThe Second Virial Coefficient & van der Waals Equation
V.C The Second Virial Coefficient & van der Waals Equation Let us study the second virial coefficient B, for a typical gas using eq.v.33). As discussed before, the two-body potential is characterized by
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 0 Simple Biasing Methods David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo simulation Markov chain to generate
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 informationVII.B Canonical Formulation
VII.B Canonical Formulation Using the states constructed in the previous section, we can calculate the canonical density matrix for non-interacting identical particles. In the coordinate representation
More informationAn Extended van der Waals Equation of State Based on Molecular Dynamics Simulation
J. Comput. Chem. Jpn., Vol. 8, o. 3, pp. 97 14 (9) c 9 Society of Computer Chemistry, Japan An Extended van der Waals Equation of State Based on Molecular Dynamics Simulation Yosuke KATAOKA* and Yuri YAMADA
More informationMelting line of the Lennard-Jones system, infinite size, and full potential
THE JOURNAL OF CHEMICAL PHYSICS 127, 104504 2007 Melting line of the Lennard-Jones system, infinite size, and full potential Ethan A. Mastny a and Juan J. de Pablo b Chemical and Biological Engineering
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 informationStatistical 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 informationSelf-referential Monte Carlo method for calculating the free energy of crystalline solids
Self-referential Monte Carlo method for calculating the free energy of crystalline solids M. B. Sweatman* Department of Chemical and Process Engineering, University of Strathclyde, Glasgow, G1 1XJ, United
More informationTSTC Lectures: Theoretical & Computational Chemistry
TSTC Lectures: Theoretical & Computational Chemistry Rigoberto Hernandez May 2009, Lecture #2 Statistical Mechanics: Structure S = k B log W! Boltzmann's headstone at the Vienna Zentralfriedhof.!! Photo
More informationABSTRACT. Sarvin Moghaddam, Doctor of Philosophy, 2003
ABSTRACT Title of Dissertation: FINELY-DISCRETIZED LATTICE MODELS FOR THERMODYNAMIC PROPERTIES OF FLUIDS Sarvin Moghaddam, Doctor of Philosophy, 2003 Dissertation directed by: Professor Mikhail Anisimov
More information424 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 informationalthough Boltzmann used W instead of Ω for the number of available states.
Lecture #13 1 Lecture 13 Obectives: 1. Ensembles: Be able to list the characteristics of the following: (a) icrocanonical (b) Canonical (c) Grand Canonical 2. Be able to use Lagrange s method of undetermined
More informationPhase transitions for particles in R 3
Phase transitions for particles in R 3 Sabine Jansen LMU Munich Konstanz, 29 May 208 Overview. Introduction to statistical mechanics: Partition functions and statistical ensembles 2. Phase transitions:
More informationChE 524 A. Z. Panagiotopoulos 1
ChE 524 A. Z. Panagiotopoulos 1 VIRIAL EXPANSIONS 1 As derived previously, at the limit of low densities, all classical fluids approach ideal-gas behavior: P = k B T (1) Consider the canonical partition
More informationContents. 1 Introduction and guide for this text 1. 2 Equilibrium and entropy 6. 3 Energy and how the microscopic world works 21
Preface Reference tables Table A Counting and combinatorics formulae Table B Useful integrals, expansions, and approximations Table C Extensive thermodynamic potentials Table D Intensive per-particle thermodynamic
More informationMolecular 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 informationGases and the Virial Expansion
Gases and the irial Expansion February 7, 3 First task is to examine what ensemble theory tells us about simple systems via the thermodynamic connection Calculate thermodynamic quantities: average energy,
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 informationsummary of statistical physics
summary of statistical physics Matthias Pospiech University of Hannover, Germany Contents 1 Probability moments definitions 3 2 bases of thermodynamics 4 2.1 I. law of thermodynamics..........................
More informationTHE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM
Molecular Simulation, 1987, Vol. 1, pp. 87-93 c Gordon and Breach Science Publishers S.A. THE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM M. MEZEI Department of Chemistry, Hunter
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 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 informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More information5. Systems in contact with a thermal bath
5. Systems in contact with a thermal bath So far, isolated systems (micro-canonical methods) 5.1 Constant number of particles:kittel&kroemer Chap. 3 Boltzmann factor Partition function (canonical methods)
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 informationGibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid
MOLECULAR PHYSICS, 1989, VOL. 68, No. 3, 629-635 Gibbs ensemble simulation of phase equilibrium in the hard core two-yukawa fluid model for the Lennard-Jones fluid by E. N. RUDISILL and P. T. CUMMINGS
More informationEquilibrium Molecular Thermodynamics from Kirkwood Sampling
This is an open access article published under a Creative Commons Attribution (CC-BY) License, which permits unrestricted use, distribution and reproduction in any medium, provided the author and source
More informationChemical Potential, Helmholtz Free Energy and Entropy of Argon with Kinetic Monte Carlo Simulation
Chemical Potential, Helmholtz Free Energy and Entropy of Argon with Kinetic Monte Carlo Simulation Chunyan Fan a, D. D. Do a *, D. icholson a and E. Ustinov b a School of Chemical Engineering, University
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationMetropolis Monte Carlo simulation of the Ising Model
Metropolis Monte Carlo simulation of the Ising Model Krishna Shrinivas (CH10B026) Swaroop Ramaswamy (CH10B068) May 10, 2013 Modelling and Simulation of Particulate Processes (CH5012) Introduction The Ising
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 informationElastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals
PHYSICAL REVIEW B 70, 214113 (2004) Elastic constants and the effect of strain on monovacancy concentration in fcc hard-sphere crystals Sang Kyu Kwak and David A. Kofke Department of Chemical and Biological
More informationA Monte Carlo method for chemical potential determination in single and multiple occupancy crystals
epl draft A Monte Carlo method for chemical potential determination in single and multiple occupancy crystals Nigel B. Wilding 1 and Peter Sollich 2 1 Department of Physics, University of Bath, Bath BA2
More informationMOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES
MOLECULAR DYNAMIC SIMULATION OF WATER VAPOR INTERACTION WITH VARIOUS TYPES OF PORES USING HYBRID COMPUTING STRUCTURES V.V. Korenkov 1,3, a, E.G. Nikonov 1, b, M. Popovičová 2, с 1 Joint Institute for Nuclear
More informationwhere (E) is the partition function of the uniform ensemble. Recalling that we have (E) = E (E) (E) i = ij x (E) j E = ij ln (E) E = k ij ~ S E = kt i
G25.265: Statistical Mechanics Notes for Lecture 4 I. THE CLASSICAL VIRIAL THEOREM (MICROCANONICAL DERIVATION) Consider a system with Hamiltonian H(x). Let x i and x j be specic components of the phase
More informationClassical Monte-Carlo simulations
Classical Monte-Carlo simulations Graduate Summer Institute Complex Plasmas at the Stevens Insitute of Technology Henning Baumgartner, A. Filinov, H. Kählert, P. Ludwig and M. Bonitz Christian-Albrechts-University
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 information763620SS STATISTICAL PHYSICS Solutions 5 Autumn 2012
7660SS STATISTICAL PHYSICS Solutions 5 Autumn 01 1 Classical Flow in Phase Space Determine the trajectories of classical Hamiltonian flow in -dim corresponding to a particle in a constant gravitational
More informationChapter 2 Experimental sources of intermolecular potentials
Chapter 2 Experimental sources of intermolecular potentials 2.1 Overview thermodynamical properties: heat of vaporization (Trouton s rule) crystal structures ionic crystals rare gas solids physico-chemical
More informationFor info and ordering all the 4 versions / languages of this book please visit: http://trl.lab.uic.edu/pon Contents Preface vii Chapter 1 Advances in Atomic and Molecular Nanotechnology Introduction 1
More informationBinary Hard-Sphere Mixtures Within Spherical Pores
Journal of the Korean Physical Society, Vol. 35, No. 4, October 1999, pp. 350 354 Binary Hard-Sphere Mixtures Within Spherical Pores Soon-Chul Kim Department of Physics, Andong National University, Andong
More informationPhysics 127b: Statistical Mechanics. Lecture 2: Dense Gas and the Liquid State. Mayer Cluster Expansion
Physics 27b: Statistical Mechanics Lecture 2: Dense Gas and the Liquid State Mayer Cluster Expansion This is a method to calculate the higher order terms in the virial expansion. It introduces some general
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 information