Thermodynamic integration
|
|
- Oswin Joseph
- 5 years ago
- Views:
Transcription
1 Thermodynamic integration Localizing liquid-solid phase transitions Christoph Tavan Freie Universität Berlin / Technische Universität Berlin December 7, 2009
2 Overview Problem Theoretical basics Thermodynamic integration Localization of phase equilibria
3 60 Problem 50 P* random init. conf. fcc init. conf. bcc init. conf S. Grandner and S. H. L. Klapp: Freezing of charged colloids in slit pores. J. Chem. Phys., 129, (2008). ρ b * Simulations predict solid and fluid phase Metastable density range
4 60 Problem 50 P* ρ b * Where is the coexistence line?
5 60 Problem 50 P* ρ b * Where is the coexistence line?
6 60 Problem 50 P* ρ b * Where is the coexistence line?
7 Thermodynamics of phase equilibria Condition for equilibrium between solid (s) and liquid (l) phase: µ s (P, T) = µ l (P, T) With chemical potential µ = g = G N and Gibbs free energy G = F PV and Helmholtz free energy F.
8 Thermodynamics of phase equilibria So calculate Gibbs free energy for both phases: g = G N = F N PV N
9 Thermodynamics of phase equilibria So calculate Gibbs free energy for both phases: g = G N = F N PV N Easy: Pressure P as an ensemble average
10 Thermodynamics of phase equilibria So calculate Gibbs free energy for both phases: g = G N = F N PV N Easy: Pressure P as an ensemble average Hard: Helmholtz free energy F: F(N,V,T) = 1 ln Z(N,V,T) β
11 Thermodynamics of phase equilibria So calculate Gibbs free energy for both phases: g = G N = F N PV N Easy: Pressure P as an ensemble average Hard: Helmholtz free energy F: F(N,V,T) = 1 ln Z(N,V,T) β 1 Z(N,V,T) = Λ 3N exp [ βu(r) ] d 3N r N! with Λ = h β/2πm and β = 1/k B T.
12 Thermodynamics of phase equilibria So calculate Gibbs free energy for both phases: g = G N = F N PV N Easy: Pressure P as an ensemble average Hard: Helmholtz free energy F: F(N,V,T) = 1 ln Z(N,V,T) β 1 Z(N,V,T) = Λ 3N exp [ βu(r) ] d 3N r N! with Λ = h β/2πm and β = 1/k B T. In general impossible to evaluate directly
13 Observation: Derivatives of F Derivatives of the thermodynamic potentials are easy to clalculate: ( ) F = P(N,V,T) V NT
14 Observation: Derivatives of F Derivatives of the thermodynamic potentials are easy to clalculate: ( ) F = P(N,V,T) V NT Intensive quantities: ( ) f (ρ,t) = P(ρ,T) ρ ρ 2 T
15 Observation: Derivatives of F Derivatives of the thermodynamic potentials are easy to clalculate: ( ) F = P(N,V,T) V NT Intensive quantities: ( ) f (ρ,t) = P(ρ,T) ρ ρ 2 Calculate f = F/N by integration: T f (ρ) f (0) = ρ 0 P( ρ,t) ρ 2 d ρ
16 Observation: Derivatives of F Derivatives of the thermodynamic potentials are easy to clalculate: ( ) F = P(N,V,T) V NT Intensive quantities: ( ) f (ρ,t) = P(ρ,T) ρ ρ 2 Calculate f = F/N by integration: T f (ρ) f (0) = ρ 0 P( ρ,t) ρ 2 d ρ
17 Good for fluid phase: P* Thermodynamic integration fest (fcc) flüssig f (ρ) = ρ 0 P( ρ,t) ρ 2 d ρ ρ*
18 Good for fluid phase: Thermodynamic integration f (ρ) = ρ 0 P( ρ,t) ρ 2 d ρ For solid phase: Phase transitions break the reversibility of the integration path
19 Good for fluid phase: Thermodynamic integration f (ρ) = ρ 0 P( ρ,t) ρ 2 d ρ For solid phase: Phase transitions break the reversibility of the integration path But thermodynamic integration is still possible
20 Thermodynamic integration: solid phase Define effective interaction potential: Ũ(λ) = (1 λ)u + λu Ref λ [0,1]
21 Thermodynamic integration: solid phase Define effective interaction potential: Ũ(λ) = (1 λ)u + λu Ref λ [0,1] Calculate derivative of the free energy: ( ) F(λ) = 1 ln Z(N,V,T,λ) λ β λ NVT = 1 Z βz λ
22 Thermodynamic integration: solid phase Define effective interaction potential: Ũ(λ) = (1 λ)u + λu Ref λ [0,1] Calculate derivative of the free energy: ( ) F(λ) = 1 ln Z(N,V,T,λ) λ β λ NVT = 1 Z βz λ Ũ(λ) = λ λ
23 Thermodynamic integration: solid phase Define effective interaction potential: Ũ(λ) = (1 λ)u + λu Ref λ [0,1] Thermodynamic integration: 1 F = F Ref 0 U Ref U λ dλ Choose reference system with known free energy F Ref.
24 Einstein crystal Potential: U Ein (r) = N i=1 α ( ) ri r 2 i,0 2 Free energy: F Ein = 3 ( ) N 2π 2 β ln αβλ 2
25 Einstein crystal: choice of α Choose α such that fluctuations in the integrand become minimal: r 2 λ=1 = r 2 λ=0
26 Einstein crystal: choice of α Choose α such that fluctuations in the integrand become minimal: r 2 λ=1 = r 2 λ=0 Mean square displacement for the einstein crystal is known: r 2 λ=1 = 3 2 k BT 2 α
27 Einstein crystal: choice of α Choose α such that fluctuations in the integrand become minimal: r 2 λ=1 = r 2 λ=0 Mean square displacement for the einstein crystal is known: r 2 λ=1 = 3 2 k BT 2 α Condition: α = 3k BT r 2 λ=0
28 Constrained center of mass For λ 0 the particles are no longer bound to their lattice sites. The integrand in 1 F = F Ein becomes sharply peaked. 0 U Ein U λ dλ } {{ } F
29 Constrained center of mass For λ 0 the particles are no longer bound to their lattice sites. The integrand in 1 F = F Ein becomes sharply peaked. 0 U Ein U λ dλ } {{ } Constrain the center of mass movement. Free energy changes: F F CM = F CM Ein FCM
30 Constrained center of mass Free energy changes: F CM = F CM Ein FCM Change due to constraint for the einstein crystal: β ( ( ) F Ein F CM 3 β 2 ) Ein = 2 ln α 4π 2 m and for arbitrary interacting crystals: β ( F F CM) = ln ( ρ ) ln ( β 2πNm )
31 Constrained center of mass Free energy changes: F CM = F CM Ein FCM Change due to constraint for the einstein crystal: β ( ( ) F Ein F CM 3 β 2 ) Ein = 2 ln α 4π 2 m and for arbitrary interacting crystals: β ( F F CM) = ln ( ρ ) ln ( β 2πNm )
32 Constrained center of mass Excess free energy F ex F F id of the unconstrained crystal: βf ex N 3 2 ln 2 ln N N ( ) 2π β FCM αβ N + 3 ( ) 2π 2N ln αβ + ln ρ N ln ρ + 1 ln(2π) 2N.
33 Constrained center of mass Excess free energy F ex F F id of the unconstrained crystal: βf ex N 3 2 ln 2 ln N N ( ) 2π β FCM αβ N + 3 ( ) 2π 2N ln αβ + ln ρ N ln ρ + 1 ln(2π) 2N.
34 Constrained center of mass Excess free energy F ex F F id of the unconstrained crystal: βf ex N 3 2 ln 2 ln N N ( ) 2π β FCM αβ N + 3 ( ) 2π 2N ln αβ + ln ρ N ln(2π) ln ρ + 1 2N. System size dependence: β F CM N = const. ln N N + O ( N 1)
35 Constrained center of mass Excess free energy F ex F F id of the unconstrained crystal: βf ex N + ln N ( ) 2π N 3 2 ln const. O ( N 1) + 3 ( ) 2π αβ 2N ln αβ + ln ρ N System size dependence: β F CM N ln(2π) ln ρ + 1 2N. ln N = const. N + O ( N 1) So βf ex /N + ln N/N depends linear on 1/N.
36 System size dependence Soft sphere fcc-crystal at ρ = , using α = 132 and N = 216, 392, 810, 1728, 3600 (β F ex* )/N + ln(n)/n Diese Arbeit (3) - 5.5(2)/N Polson et al (7) (7)/N /N
37 Summary Thermodynamic integration for fluid phase: βf ex N ρ = β 0 P( ρ,t) ρ/β ρ 2 d ρ
38 Summary Thermodynamic integration for fluid phase: βf ex N ρ = β 0 P( ρ,t) ρ/β ρ 2 d ρ Thermodynamic integration for solid phase: βf ex N = 3 2 ln 2 ln N N ( ) 2π β 1 U Ein U αβ N λ dλ + 3 2N ln 0 + ln ρ N ln ρ + 1 ln(2π) 2N ( ) 2π αβ
39 Can we now find phase-equilibria? Theoretically expected form of the chemical potential vs. pressure: fest flüssig µ P f P c P s P
40 Can we now find phase-equilibria? Computer simulation results for U(r) = ε(r ) 12 soft spheres at T = 1: fest (fcc) flüssig µ* P*
41 Can we now find phase-equilibria? Difference between the previously shown curves: Polynom 2. Grades Polynom 3. Grades Polynom 4. Grades µ* s -µ* f P* Phase coexistence at P = 22.6(5)
42 Equation of state for soft spheres Coexistence at P = 22.6(5), ρ f = 1.152(7), ρ s = 1.194(7) P* P* flüssig fest ρ* ρ* fest (fcc) flüssig Phasen koexistenz
43 Summary Thermodynamic integration for liquid and solid phase. Allows calculation of Helmholtz free energy in NVT-Ensemble. Gibbs free energy follows directly. Phase coexistence can be found. Equation of state can be drawn.
44 What about the charged colloids? 60 50??? P* ρ* b Coexistence pressure still to be found... Effects due to walls?
Phase 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 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 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 informationPeter A. Monson. Department of Chemical Engineering, University of Massachusetts, Amherst, MA 01003, USA. Acknowledgements: John Edison (Utrecht)
Understanding Adsorption/Desorption Hysteresis for Fluids in Mesoporous Materials using Simple Molecular Models and Classical Density Functional Theory Peter A. Monson Department of Chemical Engineering,
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 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 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 informationGas-liquid phase separation in oppositely charged colloids: stability and interfacial tension
7 Gas-liquid phase separation in oppositely charged colloids: stability and interfacial tension We study the phase behaviour and the interfacial tension of the screened Coulomb (Yukawa) restricted primitive
More information(Crystal) Nucleation: The language
Why crystallization requires supercooling (Crystal) Nucleation: The language 2r 1. Transferring N particles from liquid to crystal yields energy. Crystal nucleus Δµ: thermodynamic driving force N is proportional
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 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 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 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 informationCritical Behavior I: Phenomenology, Universality & Scaling
Critical Behavior I: Phenomenology, Universality & Scaling H. W. Diehl Fachbereich Physik, Universität Duisburg-Essen, Campus Essen 1 Goals recall basic facts about (static equilibrium) critical behavior
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 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 informationFokker-Planck Equation with Detailed Balance
Appendix E Fokker-Planck Equation with Detailed Balance A stochastic process is simply a function of two variables, one is the time, the other is a stochastic variable X, defined by specifying: a: the
More informationFree energy calculations using molecular dynamics simulations. Anna Johansson
Free energy calculations using molecular dynamics simulations Anna Johansson 2007-03-13 Outline Introduction to concepts Why is free energy important? Calculating free energy using MD Thermodynamical Integration
More informationarxiv:cond-mat/ v1 [cond-mat.stat-mech] 9 Apr 2003
Free energy calculations of elemental sulphur crystals via molecular dynamics simulations C. Pastorino and Z. Gamba arxiv:cond-mat/0304224v1 [cond-mat.stat-mech] 9 Apr 2003 Department of Physics, Comisión
More informationComplementary approaches to high T- high p crystal structure stability and melting!
Complementary approaches to high T- high p crystal structure stability and melting! Dario ALFÈ Department of Earth Sciences & Department of Physics and Astronomy, Thomas Young Centre@UCL & London Centre
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 informationPhase Behaviour of a System of Inverse Patchy Colloids: A Simulation Study
DIPLOMARBEIT Phase Behaviour of a System of Inverse Patchy Colloids: A Simulation Study Ausgeführt am Institut für Theoretische Physik der Technischen Universität Wien unter der Anleitung von Ao. Univ.
More informationDependence of the surface tension on curvature
Dependence of the surface tension on curvature rigorously determined from the density profiles of nanodroplets Athens, st September M. T. Horsch,,, S. Eckelsbach, H. Hasse, G. Jackson, E. A. Müller, G.
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 informationDENSITY FUNCTIONAL THEORY FOR STUDIES OF MULTIPLE STATES OF INHOMOGENEOUS FLUIDS AT SOLID SURFACES AND IN PORES.
J. Smith, D. Stroud, MRS Symposium Proceedings Series, v.49, p.7-33, 998. DENSITY FUNCTIONAL THEORY FOR STUDIES OF MULTIPLE STATES OF INHOMOGENEOUS FLUIDS AT SOLID SURFACES AND IN PORES. A.. NEIMARK, and
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 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 informationChapter 6. Heat capacity, enthalpy, & entropy
Chapter 6 Heat capacity, enthalpy, & entropy 1 6.1 Introduction In this lecture, we examine the heat capacity as a function of temperature, compute the enthalpy, entropy, and Gibbs free energy, as functions
More informationGeometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M.
Geometry explains the large difference in the elastic properties of fcc and hcp crystals of hard spheres Sushko, N.; van der Schoot, P.P.A.M. Published in: Physical Review E DOI: 10.1103/PhysRevE.72.067104
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 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 informationHandout 10. Applications to Solids
ME346A Introduction to Statistical Mechanics Wei Cai Stanford University Win 2011 Handout 10. Applications to Solids February 23, 2011 Contents 1 Average kinetic and potential energy 2 2 Virial theorem
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 informationNucleation rate (m -3 s -1 ) Radius of water nano droplet (Å) 1e+00 1e-64 1e-128 1e-192 1e-256
Supplementary Figures Nucleation rate (m -3 s -1 ) 1e+00 1e-64 1e-128 1e-192 1e-256 Calculated R in bulk water Calculated R in droplet Modified CNT 20 30 40 50 60 70 Radius of water nano droplet (Å) Supplementary
More informationMOLECULAR DYNAMICS CALCULATION OF SURFACE TENSION USING THE EMBEDDED ATOM MODEL
MOLECULAR DYNAMICS CALCULATION OF SURFACE TENSION USING THE EMBEDDED ATOM MODEL David Belashchenko & N. Kravchunovskaya Moscow Steel and Alloys Institute, Russian Federation Oleg Ostrovski University of
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 informationPart1B(Advanced Physics) Statistical Physics
PartB(Advanced Physics) Statistical Physics Course Overview: 6 Lectures: uesday, hursday only 2 problem sheets, Lecture overheads + handouts. Lent erm (mainly): Brief review of Classical hermodynamics:
More informationFluid and solid phases of the Gaussian core model
J. Phys.: Condens. Matter (000) 5087 508. Printed in the UK PII: S0953-8984(00)954-6 Fluid and solid phases of the Gaussian core model A Lang, C N Likos, M Watzlawek andhlöwen Institut für Theoretische
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 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 informationII. Equilibrium Thermodynamics Lecture 7: Statistical Thermodynamics
II. Equilibrium Thermodynamics Lecture 7: Statistical Thermodynamics Notes by ChangHoon Lim (and MZB) Open circuit voltage of galvanic cell is To understand compositional effects on, we need to consider
More informationEquilibrium Phase Behavior and Maximally. Random Jammed State of Truncated Tetrahedra
Equilibrium Phase Behavior and Maximally Random Jammed State of Truncated Tetrahedra Duyu Chen,, Yang Jiao, and Salvatore Torquato,,,,, Department of Chemistry, Princeton University, Princeton, New Jersey
More information2m + U( q i), (IV.26) i=1
I.D The Ideal Gas As discussed in chapter II, micro-states of a gas of N particles correspond to points { p i, q i }, in the 6N-dimensional phase space. Ignoring the potential energy of interactions, the
More informationStudy of Phase Transition in Pure Zirconium using Monte Carlo Simulation
Study of Phase Transition in Pure Zirconium using Monte Carlo Simulation Wathid Assawasunthonnet, Abhinav Jain Department of Physics University of Illinois assawas1@illinois.edu Urbana, IL Abstract NPT
More informationFig. 3.1? Hard core potential
6 Hard Sphere Gas The interactions between the atoms or molecules of a real gas comprise a strong repulsion at short distances and a weak attraction at long distances Both of these are important in determining
More informationPre-yield non-affine fluctuations and a hidden critical point in strained crystals
Supplementary Information for: Pre-yield non-affine fluctuations and a hidden critical point in strained crystals Tamoghna Das, a,b Saswati Ganguly, b Surajit Sengupta c and Madan Rao d a Collective Interactions
More informationCH 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 information9.1 System in contact with a heat reservoir
Chapter 9 Canonical ensemble 9. System in contact with a heat reservoir We consider a small system A characterized by E, V and N in thermal interaction with a heat reservoir A 2 characterized by E 2, V
More informationChapter 2. Thermodynamics of pure metals and alloys
-1 Chapter hermodynamics of pure metals and alloys.1 Introduction Crucial to the design and processing of materials is a quantitative understanding of the phase-equilibrium information with respect to
More informationarxiv: v1 [cond-mat.soft] 3 Sep 2010
Free energies, vacancy concentrations and density distribution anisotropies in hard sphere crystals: A combined density functional and simulation study M. Oettel 1,, S. Görig 1, A. Härtel 3, H. Löwen 3,
More informationPerturbation theory calculations of model pair potential systems
Graduate Theses and Dissertations Graduate College 2016 Perturbation theory calculations of model pair potential systems Jianwu Gong Iowa State University Follow this and additional works at: http://lib.dr.iastate.edu/etd
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 informationSub -T g Relaxation in Thin Glass
Sub -T g Relaxation in Thin Glass Prabhat Gupta The Ohio State University ( Go Bucks! ) Kyoto (January 7, 2008) 2008/01/07 PK Gupta(Kyoto) 1 Outline 1. Phenomenology (Review). A. Liquid to Glass Transition
More informationStatistical 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 informationThermodynamics of phase transitions
Thermodynamics of phase transitions Katarzyna Sznajd-Weron Institute of Physics Wroc law University of Technology, Poland 7 Oct 2013, SF-MTPT Katarzyna Sznajd-Weron (WUT) Thermodynamics of phase transitions
More informationThermal 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 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 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 informationMidterm Examination April 2, 2013
CHEM-UA 652: Thermodynamics and Kinetics Professor M.E. Tuckerman Midterm Examination April 2, 203 NAME and ID NUMBER: There should be 8 pages to this exam, counting this cover sheet. Please check this
More informationSummary: Thermodynamic Potentials and Conditions of Equilibrium
Summary: Thermodynamic Potentials and Conditions of Equilibrium Isolated system: E, V, {N} controlled Entropy, S(E,V{N}) = maximum Thermal contact: T, V, {N} controlled Helmholtz free energy, F(T,V,{N})
More informationEvaluating Free Energy and Chemical Potential in Molecular Simulation
Evaluating Free Energy and Chemical Potential in Molecular Simulation Marshall T. McDonnell MSE 614 April 26, 2016 Outline Free Energy Motivation Theory Free Energy Perturbation Thermodynamic Integration
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 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 informationLiquid-vapour oscillations of water in hydrophobic nanopores
Alicante 2003 4th European Biophysics Congress Liquid-vapour oscillations of water in hydrophobic nanopores 7 th July 2003 Oliver Beckstein and Mark S. P. Sansom Department of Biochemistry, Laboratory
More informationAn explicit expression for finite-size corrections to the chemical potential
J. Phys.: Condens. Matter 1 (1989) 8659-8665. Printed in the UK An explicit expression for finite-size corrections to the chemical potential B Smitt and D Frenkelt t Koninklijke/Shell-Laboratorium, Amsterdam
More informationAnalysis 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 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 informationNENG 301 Week 8 Unary Heterogeneous Systems (DeHoff, Chap. 7, Chap )
NENG 301 Week 8 Unary Heterogeneous Systems (DeHoff, Chap. 7, Chap. 5.3-5.4) Learning objectives for Chapter 7 At the end of this chapter you will be able to: Understand the general features of a unary
More informationarxiv: v1 [cond-mat.soft] 28 Jul 2018
The pair-potential system. I. Fluid phase isotherms, isochores, and quasiuniversality Andreas Kvist Bacher, Thomas B. Schrøder, and Jeppe C. Dyre Glass and Time, IMFUFA, Department of Science and Environment,
More informationStates of matter Part 2
Physical Pharmacy Lecture 2 States of matter Part 2 Assistant Lecturer in Pharmaceutics Overview The Liquid State General properties Liquefaction of gases Vapor pressure of liquids Boiling point The Solid
More informationThermodynamics of phase transitions
Thermodynamics of phase transitions Katarzyna Sznajd-Weron Department of Theoretical of Physics Wroc law University of Science and Technology, Poland March 12, 2017 Katarzyna Sznajd-Weron (WUST) Thermodynamics
More informationFrom the ideal gas to an ideal glass: a thermodynamic route to random close packing.
Leslie V. Woodcock From the ideal gas to an ideal glass: a thermodynamic route to random close packing. Department of Chemical and Biomolecular Engineering National University of Singapore, Singapore 117576
More informationHeterogenous Nucleation in Hard Spheres Systems
University of Luxembourg, Softmatter Theory Group May, 2012 Table of contents 1 2 3 Motivation Event Driven Molecular Dynamics Time Driven MD simulation vs. Event Driven MD simulation V(r) = { if r < σ
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 informationarxiv: v1 [cond-mat.soft] 20 Jun 2008
Accurate determination of crystal structures based on averaged local bond order parameters Wolfgang Lechner and Christoph Dellago Faculty of hysics, University of Vienna, Boltzmanngasse, 19 Vienna, Austria
More informationModeling Dynamics (and Thermodynamics) of Fluids Confined in Mesoporous Materials. Escuela Giorgio Zgrablich,
Modeling Dynamics (and Thermodynamics) of Fluids Confined in Mesoporous Materials Peter A. Monson Department of Chemical Engineering, University of Massachusetts, Amherst, MA 01003, USA Escuela Giorgio
More informationThermodynamics of histories: application to systems with glassy dynamics
Thermodynamics of histories: application to systems with glassy dynamics Vivien Lecomte 1,2, Cécile Appert-Rolland 1, Estelle Pitard 3, Kristina van Duijvendijk 1,2, Frédéric van Wijland 1,2 Juan P. Garrahan
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
More informationMonte Carlo Methods. Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects
Monte Carlo Methods Ensembles (Chapter 5) Biased Sampling (Chapter 14) Practical Aspects Lecture 1 2 Lecture 1&2 3 Lecture 1&3 4 Different Ensembles Ensemble ame Constant (Imposed) VT Canonical,V,T P PT
More informationChapter 7 PHASE EQUILIBRIUM IN A ONE-COMPONENT SYSTEM
Chapter 7 PHASE EQUILIBRIUM IN A ONE-COMPONENT SYSTEM 7.1 INTRODUCTION The intensive thermodynamic properties of a system are temperature, pressure, and the chemical potentials of the various species occurring
More informationIntroduction to Phase Transitions in Statistical Physics and Field Theory
Introduction to Phase Transitions in Statistical Physics and Field Theory Motivation Basic Concepts and Facts about Phase Transitions: Phase Transitions in Fluids and Magnets Thermodynamics and Statistical
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 11 I. PHYSICAL AND CHEMICAL RELEVANCE OF FREE ENERGY In this section, we will consider some examples showing the significance of free energies.
More informationFajun Zhang, Roland Roth, Marcell Wolf, Felix Roosen-Runge, Maximilian W. A. Skoda, Robert M. J. Jacobs, Michael Stzuckie and Frank Schreiber
Soft Matter, 2012, 8, 1313 Fajun Zhang, Roland Roth, Marcell Wolf, Felix Roosen-Runge, Maximilian W. A. Skoda, Robert M. J. Jacobs, Michael Stzuckie and Frank Schreiber Universität Tübingen, Institut für
More informationLiquids and Solids. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Liquids and Solids Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. 1 Gases, Liquids and Solids Gases are compressible fluids. They have no proper volume and proper
More informationThe glass transition as a spin glass problem
The glass transition as a spin glass problem Mike Moore School of Physics and Astronomy, University of Manchester UBC 2007 Co-Authors: Joonhyun Yeo, Konkuk University Marco Tarzia, Saclay Mike Moore (Manchester)
More informationCHEM-UA 652: Thermodynamics and Kinetics
CHEM-UA 652: hermodynamics and Kinetics Notes for Lecture 6 I. SAISICAL MECHANICS OF SOLVAION: SOLVAION FREE ENERGIES We consider a solvent with coordinates r (a),...,r(a) N a in to which a solute with
More informationJanuary 6, 2015 I. INTRODUCTION FIRST ORDER PHASE TRANSITIONS. A. Basic phenomena
1 January 6, 2015 I. INTRODUCTION If I need to summarize the focus of this coming quarter, I d say it is the topic of phase transitions. Phase transitions, however, are a manifestation of interactions
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 informationthe expansion for the Helmholtz energy derived in Appendix A, part 2, the expression for the surface tension becomes: σ = ( a + ½ k(ρ) ρ 2 x ) dx
Motivation The surface tension plays a major role in interfacial phenomena. It is the fundamental quantity that determines the pressure change across a surface due to curvature. This in turn is the basis
More informationComputing phase diagrams of model liquids & self-assembly of (bio)membranes
Computing phase diagrams of model liquids & self-assembly of (bio)membranes Ulf Rørbæk Pedersen (ulf@urp.dk) Crystal structure of a simplistic model of a molecule. What is the melting temperature? A typical
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 informationHow patchy can one get and still condense? The role of dissimilar patches in the interactions of colloidal particles
Page 1 of 100 How patchy can one get and still condense? The role of dissimilar patches in the interactions of colloidal particles Paulo Teixeira Instituto Superior de Engenharia de Lisboa and Centro de
More informationPerturbation approach for equation of state for hard-sphere and Lennard Jones pure fluids
PRAMANA c Indian Academy of Sciences Vol. 76, No. 6 journal of June 2011 physics pp. 901 908 Perturbation approach for equation of state for hard-sphere and Lennard Jones pure fluids S B KHASARE and M
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 informationWetting Transitions at Fluid Interfaces and Related Topics
Wetting Transitions at Fluid Interfaces and Related Topics Kenichiro Koga Department of Chemistry, Faculty of Science, Okayama University Tsushima-Naka 3-1-1, Okayama 7-853, Japan Received April 3, 21
More information4.1 Constant (T, V, n) Experiments: The Helmholtz Free Energy
Chapter 4 Free Energies The second law allows us to determine the spontaneous direction of of a process with constant (E, V, n). Of course, there are many processes for which we cannot control (E, V, n)
More informationCRYSTAL STRUCTURE, PHASE CHANGES, AND PHASE DIAGRAMS
CRYSTAL STRUCTURE, PHASE CHANGES, AND PHASE DIAGRAMS CRYSTAL STRUCTURE CRYSTALLINE AND AMORPHOUS SOLIDS Crystalline solids have an ordered arrangement. The long range order comes about from an underlying
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 informationTHERMODYNAMICS THERMOSTATISTICS AND AN INTRODUCTION TO SECOND EDITION. University of Pennsylvania
THERMODYNAMICS AND AN INTRODUCTION TO THERMOSTATISTICS SECOND EDITION HERBERT B. University of Pennsylvania CALLEN JOHN WILEY & SONS New York Chichester Brisbane Toronto Singapore CONTENTS PART I GENERAL
More informationThe Phases of Hard Sphere Systems
The Phases of Hard Sphere Systems Tutorial for Spring 2015 ICERM workshop Crystals, Quasicrystals and Random Networks Veit Elser Cornell Department of Physics outline: order by disorder constant temperature
More information