CE 530 Molecular Simulation
|
|
- Margery Carpenter
- 6 years ago
- Views:
Transcription
1 CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo
2 2 Thermodynamic Phase Equilibria Certain thermodynamic states find that two (or more) phases may be equally stable thermodynamic phase coexistence Phases are distinguished by an order parameter; e.g. density structural order parameter magnetic or electrostatic moment Intensive-variable requirements of phase coexistence T a = T b P a = P b µ i a = µ ib, i =,,C
3 Phase Diagrams Phase behavior is described via phase diagrams two-phase regions two-phase lines 3 Field variable e.g., temperature isobar V critical point tie line L+V triple point S+V Density variable Gibbs phase rule F = 2 + C P F = number of degrees of freedom independently specified thermodynamic state variables C = number of components (chemical species) in the system P = number of phases Single component, F = 3 P L L+S S Field variable (pressure) S L critical point V Field variable (temperature)
4 Approximate Methods for Phase Equilibria by Molecular Simulation. 4 Adjust field variables until crossing of phase transition is observed e.g., lower temperature until condensation/freezing occurs Problem: metastability may obsure location of true transition isobar Temperature V Hysteresis L+V S+V L L+S S Pressure S L V Metastable region Density Temperature Try it out with the LJ NPT-MC applet
5 Approximate Methods for Phase Equilibria by Molecular Simulation 2. 5 Select density variable inside of two-phase region choose density between liquid and vapor values Problem: finite system size leads to large surface-to-volume ratios for each phase bulk behavior not modeled Metastability may still be a problem Temperature V L+V L L+S S S+V Density Try it out with the SW NVT MD applet
6 Rigorous Methods for Phase Equilibria 6 by Molecular Simulation Search for conditions that satisfy equilibrium criteria equality of temperature, pressure, chemical potentials Difficulties evaluation of chemical potential often is not easy tedious search-and-evaluation process T constant Chemical potential Simulation data Simulation data or equation of state Pressure
7 Rigorous Methods for Phase Equilibria 7 by Molecular Simulation Search for conditions that satisfy equilibrium criteria equality of temperature, pressure, chemical potentials Difficulties evaluation of chemical potential often is not easy tedious search-and-evaluation process T constant Chemical potential Simulation data or equation of state Pressure
8 Rigorous Methods for Phase Equilibria 8 by Molecular Simulation Search for conditions that satisfy equilibrium criteria equality of temperature, pressure, chemical potentials Difficulties evaluation of chemical potential often is not easy tedious search-and-evaluation process T constant Chemical potential Simulation data or equation of state Pressure
9 Rigorous Methods for Phase Equilibria 9 by Molecular Simulation Search for conditions that satisfy equilibrium criteria equality of temperature, pressure, chemical potentials Difficulties evaluation of chemical potential often is not easy tedious search-and-evaluation process T constant Chemical potential Simulation data or equation of state Pressure
10 Rigorous Methods for Phase Equilibria 0 by Molecular Simulation Search for conditions that satisfy equilibrium criteria equality of temperature, pressure, chemical potentials Difficulties evaluation of chemical potential often is not easy tedious search-and-evaluation process Coexistence T constant Chemical potential Simulation data or equation of state Pressure
11 Rigorous Methods for Phase Equilibria by Molecular Simulation Search for conditions that satisfy equilibrium criteria equality of temperature, pressure, chemical potentials Difficulties evaluation of chemical potential often is not easy tedious search-and-evaluation process Coexistence T constant Many simulations to get one point on phase diagram Chemical potential Simulation data or equation of state Pressure S L + V Pressure Temperature
12 Gibbs Ensemble. 2 Panagiotopoulos in 987 introduced a clever way to simulate coexisting phases without an interface Two simulation volumes Thermodynamic contact without physical contact How?
13 Gibbs Ensemble 2. 3 MC simulation includes moves that couple the two simulation volumes Particle exchange equilibrates chemical potential Volume exchange equilibrates pressure Incidentally, the coupled moves enforce mass and volume balance Try it out with the LJ GEMC applet
14 Gibbs Ensemble Formalism 4 Partition function GE Q ( N, V, T ) = dv Q( N, V, T ) Q( N N, V V, T ) Limiting distribution N V N = 0 0 N N N 2 V V V 2 3N N N N 2 N2 N N U( s, V) N2 N2 U( s, V2) = GE 2 N Λ π ( r, N, V ) dvds ds V ds e V ds e Q General algorithm. For each trial: with some pre-specified probability, select molecules displacement/rotation trial volume exchange trial molecule exchange trial conduct trial move and decide acceptance
15 Molecule-Exchange Trial Move Analysis of Trial Probabilities 5 Detailed specification of trial moves and probabilities Event [reverse event] Select box A [select box B] Select molecule k [select molecule k] Probability [reverse probability] /2 /2 /N A [/(N B +)] Forward-step trial probability Reverse-step trial probability ds min(, χ) N 2 A ds min(, ) 2 N + χ B Move to r new [move back to r old ] Scaled volume ds/ [ds/] Accept move [accept move] min(,χ) [min(,/χ)] χ is formulated to satisfy detailed balance
16 6 Molecule-Exchange Trial Move Analysis of Detailed Balance Forward-step trial probability ds min(, χ) N 2 A Reverse-step trial probability ds min(, ) 2 N + χ B Detailed balance π i π ij = π j π ji 3N N N Limiting N N ( A, ) ( B, ) (,, ) A NB NA NA U VA NB NB U VB A A A GE A B distribution N V dv d d Λ V d e Q V d e π s = s r s s s s
17 7 Molecule-Exchange Trial Move Analysis of Detailed Balance Forward-step trial probability ds min(, χ) N 2 A Reverse-step trial probability ds min(, ) 2 N + χ B Detailed balance π i π ij = π j π ji old new U N A N NB NB U NA NA NB+ NB+ e V min(, ) min(, ) A ds VB ds dv d A ds χ e VA ds VB ds dv s A χ 3N GE 2N = 3N GE Λ Q A Λ Q 2( NB + ) 3N N N Limiting N N ( A, ) ( B, ) (,, ) A NB NA NA U VA NB NB U VB A A A GE A B distribution N V dv d d Λ V d e Q V d e π s = s r s s s s
18 Molecule-Exchange Trial Move Analysis of Detailed Balance 8 Forward-step trial probability ds min(, χ) N 2 A Reverse-step trial probability ds min(, ) 2 N + χ B Detailed balance π i π ij = π j π ji old new U N A N NB NB U NA NA NB+ NB+ e V min(, ) min(, ) A ds VB ds dv d A ds χ e VA ds VB ds dv s A χ 3N GE 2N = 3N GE Λ Q A Λ Q 2( NB + ) e old new V χ = e N V N + U U B A A B V χ = V B A NA e N + B Δ U tot Acceptance probability
19 Volume-Exchange Trial Move Analysis of Trial Probabilities 9 Take steps in λ = ln(v A /V B ) ( ) V V V A V V dλ = + dv = dv dv A = V A B A B V V dλ A B Detailed specification of trial moves and probabilities A Step to λ new [step to λ old ] Event [reverse event] Accept move [accept move] Probability [reverse probability] dλ/δ dλ/δ min(,χ) [min(,/χ)] Forward-step trial probability Reverse-step trial probability dλ min(, χ) Δ dλ min(, ) Δ χ
20 20 Volume-Exchange Trial Move Analysis of Detailed Balance Forward-step trial probability dλ min(, χ) Δ Reverse-step trial probability dλ min(, ) Δ χ Detailed balance π i π ij = π j π ji Limiting distribution 3N N N N N ( A, ) ( B, ) (,, ) A N dλλ B NA NA U VA NB NB U V π N B A VA dλd d + s V GE A d e + s r s s = s VB ds e VQ
21 2 Volume-Exchange Trial Move Analysis of Detailed Balance Forward-step trial probability dλ min(, χ) Δ Reverse-step trial probability dλ min(, ) Δ χ old Detailed balance π i π ij = π j π ji new N A+ NA B+ NB U NA+ N NB+ B U N N e V min(, ) Ai, ds VBi, ds dλmin(, χ) e V d Aj, d VB, j d λ s s χ = 3N GE 3N GE Λ Q Δ Λ Q Δ Limiting distribution 3N N N N N ( A, ) ( B, ) (,, ) A N dλλ B NA NA U VA NB NB U V π N B A VA dλd d + s V GE A d e + s r s s = s VB ds e VQ
22 Volume-Exchange Trial Move Analysis of Detailed Balance 22 Forward-step trial probability dλ min(, χ) Δ Reverse-step trial probability dλ min(, ) Δ χ old Detailed balance π i π ij = π j π ji new N A+ NA B+ NB U NA+ N NB+ B U N N e V min(, ) Ai, ds VBi, ds dλmin(, χ) e V d Aj, d VB, j d λ s s χ = 3N GE 3N GE Λ Q Δ Λ Q Δ U old new NA+ NB+ U NA+ NB+ Ai, Bi, χ = Ai, Bi, e V V e V V N + new NB+ VB old VB new A V χ = A old VA e ΔU tot Acceptance probability
23 27 Gibbs Ensemble Limitations Limitations arise from particle-exchange requirement Molecular dynamics (polarizable models) Dense phases, or complex molecular models Solid phases N = 4n 3 (fcc)?
24 Approaching the Critical Point 28 Critical phenomena are difficult to capture by molecular simulation Density fluctuations are related to compressibility Correlations between fluctuations important to behavior in thermodynamic limit correlations have infinite range in simulation correlations limited by system size surface tension vanishes P ρ = 0 V ρ κ = = V P ρ P T Pressure V isotherm C L+V L L+S S σ 2 ρ κ S+V Density
25 29 Critical Phenomena and the Gibbs Ensemble Phase identities break down as critical point is approached Boxes swap identities (liquid vs. vapor) Both phases begin to appear in each box surface free energy goes to zero
26 30 Gibbs-Duhem Integration GD equation can be used to derive Clapeyron equation ln p σ = Δh ΔZ equation for coexistence line Treat as a nonlinear first-order ODE use simulation to evaluate right-hand side Trace an integration path the coincides with the coexistence line Field variable (pressure) S L critical point V Field variable (temperature)
27 GDI Predictor-Corrector Implementation 3 Given initial condition and slope (= Δh/ΔZ), predict new (p,t) pair. lnp Evaluate slope at new state condition lnp and use to correct estimate of new (p,t) pair lnp
28 32 GDI Simulation Algorithm Estimate pressure and temperature from predictor Begin simultaneous NpT simulation of both phases at the estimated conditions As simulation progresses, estimate new slope from ongoing simulation data and use with corrector to adjust p and T Continue simulation until the iterations and the simulation averages converge Repeat for the next state point Each simulation yields a coexistence point. Particle insertions are never attempted.
29 33 GDI Sources of Error Three primary sources of error in GDI method Stochastic error Finite step of integrator Stability lnp Incorrect initial condition? Uncertainty in slopes due to uncertainty in simulation averages 2 2 True curve 2 Quantify via propagation of error using confidence limits of simulation averages Smaller step Integrate series with every-other datum Compare predictor and corrector values Stability can be quantified Error initial condition can be corrected even after series is complete
30 GDI Applications and Extensions 34 Applied to a pressure-temperature phase equilibria in a wide variety of model systems emphasis on applications to solid-fluid coexistence Can be extended to describe coexistence in other planes variation of potential parameters inverse-power softness stiffness of polymers range of attraction in square-well and triangle-well variation of electrostatic polarizability as with free energy methods, need to identify and measure variation of free energy with potential parameter mixtures A = λ U λ
31 Semigrand Ensemble Thermodynamic formalism for mixtures d( A) = Ud PdV + µ dn+ µ 2dN2 + µ dn2 µ dn2 Rearrange d( A) = Ud PdV + µ d( N+ N2) + ( µ 2 µ ) dn2 Legendre transform = Ud PdV + µ dn + Δµ dn ( ) 2 2 d Y d( A Δµ N ) 2 2 ( ) = Ud PdV + µ dn N d Δµ Independent variables include N and Δµ 2 Dependent variables include N 2 must determine this by ensemble average ensemble includes elements differing in composition at same total N Ensemble distribution π N N N N (,, 2) 2 2 (,, N E p r N N 2 ) e e + Δµ r p = 3N Υ h N!
32 36 GDI/GE with the Semigrand Ensemble Governing differential equation (pressure-composition plane) ln p x = Δ µ Δ ΔZ 2 T, σ 2 Simple LJ Binary σ = ε = σ 2 = ε 2 = 0.75 σ 22 = ε 22 = Pressure Gibbs Ensemble Gibbs-Duhem Integration Mole Fraction of Species 2.0
33 GDI with the Semigrand Ensemble Freezing of polydisperse hard spheres appropriate model for colloids Integrate in pressure-polydispersity plane 37 Findings upper polydisersity bound to freezing re-entrant melting dp dν sat 2 Δm2 /2 ν = Δ( V / N).4.3 c = Composition (Fugacity) Density.2. solid c 2 =.0 fluid m 2 (ν) Polydispersity Diameter
34 38 Other Views of Coexistence Surface Dew and bubble lines Residue curves Azeotropes semigrand ensemble formulate appropriate differential equation integrate as in standard GDI method right-hand side involves partial molar properties Δ Δ + = Δ Δ Δ Δ Δ µ µ µ µ µ µ µ,,,,,, ) ( P P azeo sat P P azeo X Y d dp P Y P X Y X d d ) ( V L V L azeo sat v v h h d dp =
35 Tracing of Azeotropes Lennard-Jones binary, variation with intermolecular potential σ =.0; σ 2 =.05; σ 22 = variable ε =.0; ε 2 = 0.90; ε 22 =.0 T = Coexistence pressure, P* X2, Y2
36 40 Summary Thermodynamic phase behavior is interesting and useful Obvious methods for evaluating phase behavior by simulation are too approximate Rigorous thermodynamic method is tedious Gibbs ensemble revolutionized methodology two coexisting phases without an interface Gibbs-Duhem integration provides an alternative to use in situations where GE fails dense and complex fluids solids Semigrand ensemble is useful formalism for mixtures GDI can be extended in many ways
Supplemental 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 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 informationChemical Potential. Combining the First and Second Laws for a closed system, Considering (extensive properties)
Chemical Potential Combining the First and Second Laws for a closed system, Considering (extensive properties) du = TdS pdv Hence For an open system, that is, one that can gain or lose mass, U will also
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 1 David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Time/s Multi-Scale Modeling Based on SDSC Blue Horizon (SP3) 1.728 Tflops
More informationEvaluation of a Locus of Azeotropes by Molecular Simulation
Evaluation of a Locus of Azeotropes by Molecular Simulation Sandeep P. Pandit and David A. Kofke Dept. of Chemical Engineering, State University of New York at Buffalo, Buffalo, NY 14260 The technique
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 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 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 informationUB association bias algorithm applied to the simulation of hydrogen fluoride
Fluid Phase Equilibria 194 197 (2002) 249 256 UB association bias algorithm applied to the simulation of hydrogen fluoride Scott Wierzchowski, David A. Kofke Department of Chemical Engineering, University
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 informationChapter 6. Phase transitions. 6.1 Concept of phase
Chapter 6 hase transitions 6.1 Concept of phase hases are states of matter characterized by distinct macroscopic properties. ypical phases we will discuss in this chapter are liquid, solid and gas. Other
More informationThermodynamic behaviour of mixtures containing CO 2. A molecular simulation study
Thermodynamic behaviour of mixtures containing. A molecular simulation study V. Lachet, C. Nieto-Draghi, B. Creton (IFPEN) Å. Ervik, G. Skaugen, Ø. Wilhelmsen, M. Hammer (SINTEF) Introduction quality issues
More informationThe Second Law of Thermodynamics (Chapter 4)
The Second Law of Thermodynamics (Chapter 4) First Law: Energy of universe is constant: ΔE system = - ΔE surroundings Second Law: New variable, S, entropy. Changes in S, ΔS, tell us which processes made
More informationCHAPTER 4 Physical Transformations of Pure Substances.
I. Generalities. CHAPTER 4 Physical Transformations of Pure Substances. A. Definitions: 1. A phase of a substance is a form of matter that is uniform throughout in chemical composition and physical state.
More informationIntroduction. Statistical physics: microscopic foundation of thermodynamics degrees of freedom 2 3 state variables!
Introduction Thermodynamics: phenomenological description of equilibrium bulk properties of matter in terms of only a few state variables and thermodynamical laws. Statistical physics: microscopic foundation
More informationThe Gibbs Phase Rule F = 2 + C - P
The Gibbs Phase Rule The phase rule allows one to determine the number of degrees of freedom (F) or variance of a chemical system. This is useful for interpreting phase diagrams. F = 2 + C - P Where F
More informationPhase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane + Xenon and Xenon + Ethane
International Journal of ChemTech Research CODEN( USA): IJCRGG ISSN : 0974-4290 Vol.5, No.6, pp 2975-2979, Oct-Dec 2013 Phase Equilibria of binary mixtures by Molecular Simulation and PR-EOS: Methane +
More informationDIRECT EVALUATION OF VAPOUR-LIQUID EQUILIBRIA BY MOLECULAR DYNAMICS USING GIBBS-DUHEM INTEGRATION
Molecular Simulation, 1996, Vol. 17, pp. 21-39 Reprints available directly from the publisher Photocopying permitted by license only 0 1996 OPA (Overseas Publishers Association) Amsterdam B.V. Published
More informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy
More informationX B1 X B Quiz Fall points total. 1. Thermodynamics. (50 points)
3.012 Quiz 4 3.012 12.15.03 Fall 2003 100 points total 1. Thermodynamics. (50 points) a. Shown below is a portion of the nickel-titanium binary phase diagram. Using the diagram, answer the questions below.
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 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 informationPhase Diagrams. NC State University
Chemistry 433 Lecture 18 Phase Diagrams NC State University Definition of a phase diagram A phase diagram is a representation of the states of matter, solid, liquid, or gas as a function of temperature
More information1. Heterogeneous Systems and Chemical Equilibrium
1. Heterogeneous Systems and Chemical Equilibrium The preceding section involved only single phase systems. For it to be in thermodynamic equilibrium, a homogeneous system must be in thermal equilibrium
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 informationCE 530 Molecular Simulation
1 CE 530 Molecular Simulation Lecture 16 Dielectrics and Reaction Field Method David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Molecular models intramolecular
More informationCHEMICAL ENGINEERING THERMODYNAMICS. Andrew S. Rosen
CHEMICAL ENGINEERING THERMODYNAMICS Andrew S. Rosen SYMBOL DICTIONARY 1 TABLE OF CONTENTS Symbol Dictionary... 3 1. Measured Thermodynamic Properties and Other Basic Concepts... 5 1.1 Preliminary Concepts
More informationFor an incompressible β and k = 0, Equations (6.28) and (6.29) become:
Internal Energy and Entropy as Functions of T and V These are general equations relating the internal energy and entropy of homogeneous fluids of constant composition to temperature and volume. Equation
More informationI: Life and Energy. Lecture 2: Solutions and chemical potential; Osmotic pressure (B Lentz).
I: Life and Energy Lecture 1: What is life? An attempt at definition. Energy, heat, and work: Temperature and thermal equilibrium. The First Law. Thermodynamic states and state functions. Reversible and
More informationChapter 11 Solution Thermodynamics: Theory
Chapter 11 Solution Thermodynamics: Theory Chapter 6: constant composition fluids. Most commonly in chemical engineering application multi component mixtures of gases or liquids undergo (composition change).
More informationln( P vap(s) / torr) = T / K ln( P vap(l) / torr) = T / K
Chem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics Fall Semester 2017 Homework Problem Set Number 9 Solutions 1. McQuarrie and Simon, 9-4. Paraphrase: Given expressions
More informationMME 2010 METALLURGICAL THERMODYNAMICS II. Fundamentals of Thermodynamics for Systems of Constant Composition
MME 2010 METALLURGICAL THERMODYNAMICS II Fundamentals of Thermodynamics for Systems of Constant Composition Thermodynamics addresses two types of problems: 1- Computation of energy difference between two
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 informationPressure-enthalpy driven molecular dynamics for thermodynamic property calculation II: applications
Fluid Phase Equilibria 200 (2002) 93 110 Pressure-enthalpy driven molecular dynamics for thermodynamic property calculation II: applications Loukas I. Kioupis, Gaurav Arya, Edward J. Maginn Department
More informationThe Clausius-Clapeyron and the Kelvin Equations
PhD Environmental Fluid Mechanics Physics of the Atmosphere University of Trieste International Center for Theoretical Physics The Clausius-Clapeyron and the Kelvin Equations by Dario B. Giaiotti and Fulvio
More informationWEEK 6. Multiphase systems
WEEK 6 Multiphase systems Multiphase systems P 237. Processes usually deal with material being transferred from one phase (gas, liquid, or solid) to another. 6.1a Phase diagrams F = force on piston Water
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 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 informationPhase Transitions. Phys112 (S2012) 8 Phase Transitions 1
Phase Transitions cf. Kittel and Krömer chap 10 Landau Free Energy/Enthalpy Second order phase transition Ferromagnetism First order phase transition Van der Waals Clausius Clapeyron coexistence curve
More informationMS212 Thermodynamics of Materials ( 소재열역학의이해 ) Lecture Note: Chapter 7
2017 Spring Semester MS212 Thermodynamics of Materials ( 소재열역학의이해 ) Lecture Note: Chapter 7 Byungha Shin ( 신병하 ) Dept. of MSE, KAIST Largely based on lecture notes of Prof. Hyuck-Mo Lee and Prof. WooChul
More informationExam 3 Solutions. ClO g. At 200 K and a total pressure of 1.0 bar, the partial pressure ratio for the chlorine-containing compounds is p ClO2
Chemistry 360 Dr. Jean M. Standard Fall 2016 Name KEY Exam 3 Solutions 1.) (14 points) Consider the gas phase decomposition of chlorine dioxide, ClO 2, ClO 2 ( g) ClO ( g) + O ( g). At 200 K and a total
More informationAt this point, we've developed the tools and basic concepts necessary to apply
18 Lecture 0 At this point, we've developed the tools and basic concepts necessary to apply thermodynamics to a number of different systems, with the ultimate goal of describing chemically reacting systems.
More informationChapter 11. Freedom of Motion. Comparisons of the States of Matter. Liquids, Solids, and Intermolecular Forces
Liquids, Solids, and Intermolecular Forces Chapter 11 Comparisons of the States of Matter The solid and liquid states have a much higher density than the gas state The solid and liquid states have similar
More informationPhysical Chemistry Physical chemistry is the branch of chemistry that establishes and develops the principles of Chemistry in terms of the underlying concepts of Physics Physical Chemistry Main book: Atkins
More informationPhysical transformations of pure substances Boiling, freezing, and the conversion of graphite to diamond examples of phase transitions changes of
Physical transformations of pure substances Boiling, freezing, and the conversion of graphite to diamond examples of phase transitions changes of phase without change of chemical composition. In this chapter
More informationChapter 11 section 6 and Chapter 8 Sections 1-4 from Atkins
Lecture Announce: Chapter 11 section 6 and Chapter 8 Sections 1-4 from Atkins Outline: osmotic pressure electrolyte solutions phase diagrams of mixtures Gibbs phase rule liquid-vapor distillation azeotropes
More informationChapter 4 Phase Transitions. 4.1 Phenomenology Basic ideas. Partition function?!?! Thermodynamic limit Statistical Mechanics 1 Week 4
Chapter 4 Phase Transitions 4.1 Phenomenology 4.1.1 Basic ideas Partition function?!?! Thermodynamic limit 4211 Statistical Mechanics 1 Week 4 4.1.2 Phase diagrams p S S+L S+G L S+G L+G G G T p solid triple
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 informationThermodynamic condition for equilibrium between two phases a and b is G a = G b, so that during an equilibrium phase change, G ab = G a G b = 0.
CHAPTER 5 LECTURE NOTES Phases and Solutions Phase diagrams for two one component systems, CO 2 and H 2 O, are shown below. The main items to note are the following: The lines represent equilibria between
More informationNote: items marked with * you should be able to perform on a closed book exam. Chapter 10 Learning Objective Checklist
Note: items marked with * you should be able to perform on a closed book exam. Chapter 10 Learning Objective Checklist Sections 10.1-10.13 find pure component properties on a binary P-x-y or T-x-y diagram.*
More informationChapter 5. Simple Mixtures Fall Semester Physical Chemistry 1 (CHM2201)
Chapter 5. Simple Mixtures 2011 Fall Semester Physical Chemistry 1 (CHM2201) Contents The thermodynamic description of mixtures 5.1 Partial molar quantities 5.2 The thermodynamic of Mixing 5.3 The chemical
More informationChapter 4. The Physical transformations of pure substances Fall Semester Physical Chemistry 1 (CHM2201)
Chapter 4. The Physical transformations of pure substances 2011 Fall Semester Physical Chemistry 1 (CHM2201) Contents Phase Diagrams 4.1 The stabilities of phases 4.2 Phase boundaries 4.3 Three representative
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 13 I. PHASE DIAGRAMS The different phases of substances are characterized by different ranges of thermodynamic variables in which these phasesarethestablephases.
More informationMultiphase Flow and Heat Transfer
Multiphase Flow and Heat Transfer Liquid-Vapor Interface Sudheer Siddapuredddy sudheer@iitp.ac.in Department of Mechanical Engineering Indian Institution of Technology Patna Multiphase Flow and Heat Transfer
More informationPHYSICS 210A : STATISTICAL PHYSICS HW ASSIGNMENT #4 SOLUTIONS. ( p)( V) = 23kJ,
PHYSICS 210A : STATISTICAL PHYSICS HW ASSIGNMENT #4 SOLUTIONS 1) ν = 8 moles of a diatomic ideal gas are subjected to a cyclic quasistatic process, the thermodynamic path for which is an ellipse in the
More informationPhase Equilibria in a One-Component System I
5.60 spring 2005 Lecture #17 page 1 Phase Equilibria in a One-Component System I Goal: Understand the general phenomenology of phase transitions and phase coexistence conditions for a single component
More informationPressure Volume Temperature Relationship of Pure Fluids
Pressure Volume Temperature Relationship of Pure Fluids Volumetric data of substances are needed to calculate the thermodynamic properties such as internal energy and work, from which the heat requirements
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 informationEffect of adding an ideal inert gas, M
Effect of adding an ideal inert gas, M Add gas M If there is no change in volume, then the partial pressures of each of the ideal gas components remains unchanged by the addition of M. If the reaction
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 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 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 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 informationPETE 310 Lectures # 36 to 37
PETE 310 Lectures # 36 to 37 Cubic Equations of State Last Lectures Instructional Objectives Know the data needed in the EOS to evaluate fluid properties Know how to use the EOS for single and for multicomponent
More informationThermodynamics of solids 5. Unary systems. Kwangheon Park Kyung Hee University Department of Nuclear Engineering
Thermodynamics of solids 5. Unary systems Kwangheon ark Kyung Hee University Department of Nuclear Engineering 5.1. Unary heterogeneous system definition Unary system: one component system. Unary heterogeneous
More informationIntermolecular Forces and Liquids and Solids
PowerPoint Lecture Presentation by J. David Robertson University of Missouri Intermolecular Forces and Liquids and Solids Chapter 11 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction
More informationChapter 8 Phase Diagram, Relative Stability of Solid, Liquid, and Gas
Chapter 8 Phase Diagram, Relative Stability of Solid, Liquid, and Gas Three states of matter: solid, liquid, gas (plasma) At low T: Solid is most stable. At high T: liquid or gas is most stable. Ex: Most
More informationChapter 1. The Properties of Gases Fall Semester Physical Chemistry 1 (CHM2201)
Chapter 1. The Properties of Gases 2011 Fall Semester Physical Chemistry 1 (CHM2201) Contents The Perfect Gas 1.1 The states of gases 1.2 The gas laws Real Gases 1.3 Molecular interactions 1.4 The van
More informationCHEMISTRY 443, Fall, 2014 (14F) Section Number: 10 Examination 2, November 5, 2014
NAME: CHEMISTRY 443, Fall, 2014 (14F) Section Number: 10 Examination 2, November 5, 2014 Answer each question in the space provided; use back of page if extra space is needed. Answer questions so the grader
More informationPhase Separation Degree of Freedom Analysis. Binary Vapor-Liquid Systems. Azeotropic Systems. - Gibbs phase rule F C P 2 -General analysis
Lecture 5. Single Equilibrium Stages (1) Phase Separation [Ch. 4] Degree of Freedom Analysis - Gibbs phase rule F CP2 -General analysis Binary Vapor-Liquid Systems - Examples of binary system - Phase equilibrium
More informationCH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria. Dr. M. Subramanian
CH2351 Chemical Engineering Thermodynamics II Unit I, II Phase Equilibria Dr. M. Subramanian Associate Professor Department of Chemical Engineering Sri Sivasubramaniya Nadar College of Engineering Kalavakkam
More information3.012 PS 7 3.012 Issued: 11.05.04 Fall 2004 Due: 11.12.04 THERMODYNAMICS 1. single-component phase diagrams. Shown below is a hypothetical phase diagram for a single-component closed system. Answer the
More informationPhase transitions of quadrupolar fluids
Phase transitions of quadrupolar fluids Seamus F. O Shea Department of Chemistry, University of Lethbridge, Lethbridge, Alberta, Canada, T1K 3M4 Girija S. Dubey Brookhaven National Laboratory, Upton, New
More informationCH.8 Polymers: Solutions, Blends, Membranes, and Gels
CH.8 Polymers: Solutions, Blends, embranes, and Gels 8. Properties of Polymers Polymers are chain-like molecules. Linear polymer Branched polymer Cross-linked polymer Polymers show little tendency to crystallize.
More informationName: Discussion Section:
CBE 141: Chemical Engineering Thermodynamics, Spring 2017, UC Berkeley Midterm 2 FORM B March 23, 2017 Time: 80 minutes, closed-book and closed-notes, one-sided 8 ½ x 11 equation sheet allowed lease show
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 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 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 informationStatistical 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(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 informationChapter 11. Intermolecular Forces and Liquids & Solids
Chapter 11 Intermolecular Forces and Liquids & Solids The Kinetic Molecular Theory of Liquids & Solids Gases vs. Liquids & Solids difference is distance between molecules Liquids Molecules close together;
More informationChem 112 Dr. Kevin Moore
Chem 112 Dr. Kevin Moore Gas Liquid Solid Polar Covalent Bond Partial Separation of Charge Electronegativity: H 2.1 Cl 3.0 H Cl δ + δ - Dipole Moment measure of the net polarity in a molecule Q Q magnitude
More informationPhysics 408 Final Exam
Physics 408 Final Exam Name You are graded on your work (with partial credit where it is deserved) so please do not just write down answers with no explanation (or skip important steps)! Please give clear,
More informationMolecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties
12 th HLRS Results and Review Workshop Molecular modeling of hydrogen bonding fluids: vapor-liquid coexistence and interfacial properties High Performance Computing Center Stuttgart (HLRS), October 8,
More informationMolecular Simulation Study of the Effect of Pressure on the Vapor-Liquid Interface of the Square-Well Fluid
4218 Langmuir 2005, 21, 4218-4226 Molecular Simulation Study of the Effect of Pressure on the Vapor-Liquid Interface of the Square-Well Fluid Jayant K. Singh and David A. Kofke* Department of Chemical
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 informationPhase Change (State Change): A change in physical form but not the chemical identity of a substance.
CHM 123 Chapter 11 11.1-11.2 Phase change, evaporation, vapor pressure, and boiling point Phase Change (State Change): A change in physical form but not the chemical identity of a substance. Heat (Enthalpy)
More informationThermodynamics, Part 2
hermodynamics, art November 1, 014 1 hase equilibria Consider the three phase diagrams in Fig.(6.1) (a) CO, (b) water; and (c) helium. CO is a normal material with a typical and simple phase diagram. At
More informationCHEMISTRY XL-14A PHYSICAL EQUILIBRIUM. August 13, 2011 Robert Iafe
CHEMISTRY XL-14A PHYSICAL EQUILIBRIUM August 13, 2011 Robert Iafe Chapter Overview 2 Phases and Phase Transitions Solubility Colligative Properties Binary Liquid Mixtures Phases and Phase Transitions 3
More informationwater Plays dominant role in radiation All three phases emit and absorb in longwave radiation
4.,4. water Plays dominant role in radiation All three phases emit and absorb in longwave radiation Some shortwave (solar) radiation is absorbed by all phases of water Principal role in the shortwave radiation
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 informationStatistical Thermodynamics. Lecture 8: Theory of Chemical Equilibria(I)
Statistical Thermodynamics Lecture 8: Theory of Chemical Equilibria(I) Chemical Equilibria A major goal in chemistry is to predict the equilibria of chemical reactions, including the relative amounts of
More informationThermodynamic integration
Thermodynamic integration Localizing liquid-solid phase transitions Christoph Tavan Freie Universität Berlin / Technische Universität Berlin December 7, 2009 Overview Problem Theoretical basics Thermodynamic
More informationName: Discussion Section:
CBE 141: Chemical Engineering Thermodynamics, Spring 2017, UC Berkeley Midterm 2 FORM A March 23, 2017 Time: 80 minutes, closed-book and closed-notes, one-sided 8 ½ x 11 equation sheet allowed Please show
More informationLecture Phase transformations. Fys2160,
Lecture 12 01.10.2018 Phase transformations Fys2160, 2018 1 A phase transformation Discontinuous change in the properties of substance when the environent is changed infinitesimaly. Change between phases
More informationIntermolecular Forces and Liquids and Solids Chapter 11
Intermolecular Forces and Liquids and Solids Chapter 11 A phase is a homogeneous part of the system in contact with other parts of the system but separated from them by a well defined boundary. Phases
More information2D Clausius Clapeyron equation
2D Clausius Clapeyron equation 1/14 Aim: Verify the Clausius Clapeyron equation by simulations of a 2D model of matter Model: 8-4 type potential ( Lennard-Jones in 2D) u(r) = 1 r 8 1 r 4 Hard attractive
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 informationChapter 5 - Systems under pressure 62
Chapter 5 - Systems under pressure 62 CHAPTER 5 - SYSTEMS UNDER PRESSURE 5.1 Ideal gas law The quantitative study of gases goes back more than three centuries. In 1662, Robert Boyle showed that at a fixed
More informationPhysical Chemistry I CHEM 4641 Final Exam 13 questions, 30 points
Physical Chemistry I CHEM 4641 Final Exam 13 questions, 30 points Name: KEY Gas constant: R = 8.314 J mol -1 K -1 = 0.008314 kj mol -1 K -1. Boltzmann constant k = 1.381 10-23 J/K = 0.6950 cm -1 /K h =
More information