Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
|
|
- Kristin Bailey
- 5 years ago
- Views:
Transcription
1 Phase Equilibria and Molecular Solutions Jan G. Korvink and Evgenii Rudnyi IMTEK Albert Ludwig University Freiburg, Germany
2 Preliminaries Learning Goals Phase Equilibria Phase diagrams and classical thermodynamics Simulations of phase equilibria Molecular solutions Quantum chemistry: molecule in the solution Simulating molecule in the solution References Leach, A.R., Molecular modelling: principles and applications, 2001
3 On-line resources Reviews of A.Z. Panagiotopoulos Monte Carlo Methods for Phase Equilibria of Fluids Gibbs Ensemble Techniques Molecular simulation of phase equilibria Free Energy Simulations Introduction to Continuum Electrostatics, with Molecular Applications (Mike Gilson)
4 Phase diagrams Introduction to phase diagrams Equilibria criteria Phase molar Gibbs energy as a thermodynamic model Thermodynamics simulation CALPHAD Problem to compute a partition function
5 Introduction to phase diagrams One-component phase diagram: S - solid, L - liquid, V - vapor Temperature-density and temperature-pressure T Τ ρ p
6 Equilibria criteria U, V const -> max S T, V const -> min F T, p const -> min G Phase 1: Phase 2: T 1 S 1 = du 1 + p 1 dv 1 µ 1 dn 1 T 2 S 2 = du 2 + p 2 dv 2 µ 2 dn 2 U - internal energy, V - volume, n - number of moles, T - temperature, p - pressure, µ - chemical potential Panagiotopolus, Maximize S = S 1 + S 2 at U 1 + U 2 = const and V 1 + V 2 = const Result: thermal T 1 = T 2, mechanical p 1 = p 2, and phase µ 1 = µ 2 equilibria criteria
7 Phase molar Gibbs energy If the molar Gibbs energy function G m ( T, p, x) is known for all the phases in question, it is possible to compute all the equilibrium properties {mole fraction x i = n i ( n i )} G G Gibbs energy G = ng m, entropy S =, volume V =, T p enthalpy H = G + TS, internal energy U = H - pv, Helmholtz energy F = G - pv Chemical potential G µ i = n i Phase equilibria by employing the equilibrium criterion
8 Thermodynamics simulation Two component systems Two phases: solid solution + liquid solution Regular solution model mix G = RT( x A ln x A + x B ln x B ) + Ωx A x B Each phase characterized by an interaction parameter - Ω. If it is zero we have an ideal solution Continuous change in Ω leads to different topology reference?
9 CALPHAD Community to produce the databases of molar Gibbs energies to compute the chemical and phase equilibria in the multicomponent systems Phenomenological approach: to choose an analytical form for the molar Gibbs energy to determine unknowns parameters from all available experimental results
10 Simulating partition function 1 Partition function: Z = , h dn exp{ βe} dp N dr N β = 1 k B T N! Helmholtz energy: F = U o k B TlnZ Impossible to compute, Metropolis algorithm can estimate only D = exp{ βe} D( p N, r N ) dp N dr N exp{ βe} dp N dr N Trick 1 = exp{ βe} exp{ βe}, then we can write 1 exp{ βe} exp{ βe} dp N dr N -- = where D = exp{ βe} Z exp{ βe} dp N dr N Does not help: it is necessary to sample high energy regions
11 Simulations of phase equilibria Computing difference between two free energies Particle insertion Direct simulation Gibbs ensemble method More tricks Example: water
12 Free Energy Difference, I Two systems: reference U 1 and in question U 2 Coupling parameter between two systems U( λ) = ( 1 λ)u 1 + λu 2 Partition function depend on λ Z = a exp{ βu( λ) } d r N F λ NVT = exp{ βe}( U( λ) λ) dr N = exp{ βe} dr N U λ, <>-ensemble average Then, thermodynamic integration: F 2 F 1 = U dλ λ Run a series simulation at different λ, then estimate integral numerically
13 Free Energy Difference, II Two systems: reference U 1 and in question U 2 F 2 F 1 = k B Tln( Z 2 Z 1 ) Write expression for Z and use 1 = exp{ β 1 U 1 } exp{ βu 1 } Thermodynamic perturbation Monte-Carlo over the system 1: Monte-Carlo over the system 2: F 2 F 1 = k B T F 2 F 1 = k B T U 2 U ln exp k B T U 1 U ln exp k B T 1 2 In principle, this is equivalent, but in practice, this is not.
14 Particle insertion: Widom approach First state - N particles, second state - N+1 particles µ F = = F n 2 F 1 = k B Tln( Z N + 1 Z N ) TV Because of change in the number of particles µ ex = µ µ id = k B T U test ln exp k B T 1 It is possible to derive similar for the particle deletion, but the performance is worse Problem with crystals: no place to insert a particle
15 Direct simulation The main problem is introducing interface Table 1: Percentage of particles in the interface of a cubic domain (from Frencel&Smit) N P int 78% 49% 14% 6% The equilibration time is too long The metastable states (overcooled liquid or oversaturated vapor) Computationally expensive Example: two-phase three component system (Panagiotopoulus, J. Phys.: Condens. Matter 12 (2000) R25)
16 Gibbs ensemble method Simulating two-phase equilibria without interface Two different boxes at given temperature Panagiotopoulus, Chemical potential in the implicit form
17 More tricks Non-Boltzman distributions Sampling high energy states: Umbrella sampling Configuration bias Monte-Carlo Cluster moves Histogram-reweighting method: Grand Canonical Monte-Carlo
18 Example: water It is necessary to derive specialized force fields Panagiotopoulos, 2000
19 QC: Molecule in solution Molecules in the electrical field Extreme case Macroscopic continuum models Electrostatic energy Discrete modeling of the solvation shell
20 Molecules in the electrical field Electrical field leads to the molecule polarization (non-zero net dipole moment) The field of these dipoles runs against the inducing field Macroscopic result - dielectric q constant U 1 q πεε o r 12 hexane - 1.9, benzene - 2.3, water - 78 The electronic distribution in the molecule in the solution is different =
21 Extreme case Solution chemistry Na 2 SO 4 = Na + + SO However, SO 4 2- is not stable in the gas phase Second electron affinity is negative SO 2-4 = SO e Quantum chemistry (or force field) of individual SO 2-4 is of no use
22 Macroscopic continuum models E = E o + E s E s = E es + E disp + E cav, total = isolated molecule + solvation energy, solvation = electrostatic + dispersion + cavity If we can not use E o, then let us solve Schrödinger equation with modified Hamiltonian H = H o + U Problems: Whether we can use macroscopic ε near cavity We take solvent at some temperature (ε depend on T), but we describe ground molecule energy E s - solvation energy or Gibbs solvation energy
23 Electrostatic interaction Kirkwood, molecule - N point charges { Q i, r i }, spherical cavity 1 ( E es --Σ, Σ n + 1) ( 1 ε) n n+ ( n+ 1)ε Q Q ( r r ) n = P i j ij i j a 2n + 1 n ( cosθ ij ) n - multipole expansion, a - cavity radius, P - Legendre polynomial n = 0, monopole, Born solvation energy for the ions, Q 2 E es = a ε n = 1, dipole, Onsager s reaction field 1 E es = -- 2( ε 1) µ2 2ε + 1 a 3 For non-spheric cavities, solve Poisson equation
24 Discrete solvation shell To model solvation shell, add point charges around the molecule in question Advantage More real interactions No cavity geometry No mascroscopic parameters Problems How to construct the solvation shell Solve Schrödinger equation for a cluster of the solute molecule and several solvent molecules chemphys/faeder/gallery/ Optimize geometry
25 Simulating molecule in solution Continuum models: Poisson-Boltzman equation Example: trypsin Stochastic mechanics
26 Continuum models Take protein and make a cavity based on van der Waals radii of atoms Write Poisson equation ε φ = 4πρ Use ε for solvent out of the cavity, and ε=2-4 within the cavity to take into account the polarization of the atoms Model a concentration of the mobile ions by Boltzman equation n = n bulk exp( Qφ ( k B T )) Combine and, to make life simple, linearize k is related to Debye-Hückel radius ε φ kφ = 4πρ Solve by the boundary element method or the finite difference method
27 Results From Leach s book 3D Electrostatic isopotential contours around trypsin [Marquart et al 1983]. Contours are drawn at -1kT (red) and +1kT (blue). The trypsin inhibitor is also shown with its electrostatic potential mapped onto the molecular surface.
28 Stochastic mechanics One solute molecule + a number of solvent molecules Non-periodic, stochastic boundary conditions Langevin dynamics Mẋ + Cẋ + V( x) = DW () t Kinetic energy: mass matrix and acceleration Dissipation energy to surrounding: dumping matrix and velocity Potential function RHS: normalized white noise, W(t) - Wiener process Fluctuation-dissipation theorem DD T = 2k B TC Typically dumping matrix is modeled as C = γm
29 Summary Phase Equilibria Phase diagrams and classical thermodynamics Simulations of phase equilibria Molecular solutions Quantum chemistry: molecule in the solution Simulating molecule in the solution
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Statistical Thermodynamics and Monte-Carlo Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals From Micro to Macro Statistical Mechanics (Statistical
More 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 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 informationThe change in free energy on transferring an ion from a medium of low dielectric constantε1 to one of high dielectric constant ε2:
The Born Energy of an Ion The free energy density of an electric field E arising from a charge is ½(ε 0 ε E 2 ) per unit volume Integrating the energy density of an ion over all of space = Born energy:
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 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 information510 Subject Index. Hamiltonian 33, 86, 88, 89 Hamilton operator 34, 164, 166
Subject Index Ab-initio calculation 24, 122, 161. 165 Acentric factor 279, 338 Activity absolute 258, 295 coefficient 7 definition 7 Atom 23 Atomic units 93 Avogadro number 5, 92 Axilrod-Teller-forces
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 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 information2. Derive ideal mixing and the Flory-Huggins models from the van der Waals mixture partition function.
Lecture #5 1 Lecture 5 Objectives: 1. Identify athermal and residual terms from the van der Waals mixture partition function.. Derive ideal mixing and the Flory-Huggins models from the van der Waals mixture
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 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 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 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 informationCOSMO-RS Theory. The Basics
Theory The Basics From µ to properties Property µ 1 µ 2 activity coefficient vapor pressure Infinite dilution Gas phase Pure compound Pure bulk compound Partition coefficient Phase 1 Phase 2 Liquid-liquid
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 informationPart II: Statistical Physics
Chapter 6: Boltzmann Statistics SDSMT, Physics Fall Semester: Oct. - Dec., 2014 1 Introduction: Very brief 2 Boltzmann Factor Isolated System and System of Interest Boltzmann Factor The Partition Function
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 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 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 information1. Electricity and Magnetism (Fall 1995, Part 1) A metal sphere has a radius R and a charge Q.
1. Electricity and Magnetism (Fall 1995, Part 1) A metal sphere has a radius R and a charge Q. (a) Compute the electric part of the Maxwell stress tensor T ij (r) = 1 {E i E j 12 } 4π E2 δ ij both inside
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 informationPart II: Statistical Physics
Chapter 6: Boltzmann Statistics SDSMT, Physics Fall Semester: Oct. - Dec., 2013 1 Introduction: Very brief 2 Boltzmann Factor Isolated System and System of Interest Boltzmann Factor The Partition Function
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 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 informationLecture 6. NONELECTROLYTE SOLUTONS
Lecture 6. NONELECTROLYTE SOLUTONS NONELECTROLYTE SOLUTIONS SOLUTIONS single phase homogeneous mixture of two or more components NONELECTROLYTES do not contain ionic species. CONCENTRATION UNITS percent
More informationModeling in Solution. Important Observations. Energy Concepts
Modeling in Solution 1 Important Observations Various electrostatic effects for a solvated molecule are often less important than for an isolated gaseous molecule when the molecule is dissolved in a solvent
More information3.320: Lecture 19 (4/14/05) Free Energies and physical Coarse-graining. ,T) + < σ > dµ
3.320: Lecture 19 (4/14/05) F(µ,T) = F(µ ref,t) + < σ > dµ µ µ ref Free Energies and physical Coarse-graining T S(T) = S(T ref ) + T T ref C V T dt Non-Boltzmann sampling and Umbrella sampling Simple
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
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 informationStatistical thermodynamics for MD and MC simulations
Statistical thermodynamics for MD and MC simulations knowing 2 atoms and wishing to know 10 23 of them Marcus Elstner and Tomáš Kubař 22 June 2016 Introduction Thermodynamic properties of molecular systems
More informationMOLECULAR INTERACTIONS NOTES
- 1 - MOLECULAR INTERACTIONS NOTES Summary of Fundamental Molecular Interactions Q1Q Ion-Ion U ( r) = 4 r πεε o µ Q cosθ U ( r) = 4πεε o r µ 1µ U ( r) 3 r µ 1 µ U ( r) 6 r µ 1 α U ( r) 6 r Ion-Dipole Dipole-Dipole
More informationMolecular Driving Forces
Molecular Driving Forces Statistical Thermodynamics in Chemistry and Biology SUBGfittingen 7 At 216 513 073 / / Ken A. Dill Sarina Bromberg With the assistance of Dirk Stigter on the Electrostatics chapters
More informationSome properties of the Helmholtz free energy
Some properties of the Helmholtz free energy Energy slope is T U(S, ) From the properties of U vs S, it is clear that the Helmholtz free energy is always algebraically less than the internal energy U.
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 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 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 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 informationto satisfy the large number approximations, W W sys can be small.
Chapter 12. The canonical ensemble To discuss systems at constant T, we need to embed them with a diathermal wall in a heat bath. Note that only the system and bath need to be large for W tot and W bath
More informationPhysics is time symmetric Nature is not
Fundamental theories of physics don t depend on the direction of time Newtonian Physics Electromagnetism Relativity Quantum Mechanics Physics is time symmetric Nature is not II law of thermodynamics -
More informationChapter 2. Dielectric Theories
Chapter Dielectric Theories . Dielectric Theories 1.1. Introduction Measurements of dielectric properties of materials is very important because it provide vital information regarding the material characteristics,
More informationClassical Thermodynamics. Dr. Massimo Mella School of Chemistry Cardiff University
Classical Thermodynamics Dr. Massimo Mella School of Chemistry Cardiff University E-mail:MellaM@cardiff.ac.uk The background The field of Thermodynamics emerged as a consequence of the necessity to understand
More informationSuggestions for Further Reading
Contents Preface viii 1 From Microscopic to Macroscopic Behavior 1 1.1 Introduction........................................ 1 1.2 Some Qualitative Observations............................. 2 1.3 Doing
More informationCH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions
CH.7 Fugacities in Liquid Mixtures: Models and Theories of Solutions The aim of solution theory is to express the properties of liquid mixture in terms of intermolecular forces and liquid structure. The
More informationNon equilibrium thermodynamics: foundations, scope, and extension to the meso scale. Miguel Rubi
Non equilibrium thermodynamics: foundations, scope, and extension to the meso scale Miguel Rubi References S.R. de Groot and P. Mazur, Non equilibrium Thermodynamics, Dover, New York, 1984 J.M. Vilar and
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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 21 May 2001 with E. A. Spiegel
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 informationMolecular Modeling -- Lecture 15 Surfaces and electrostatics
Molecular Modeling -- Lecture 15 Surfaces and electrostatics Molecular surfaces The Hydrophobic Effect Electrostatics Poisson-Boltzmann Equation Electrostatic maps Electrostatic surfaces in MOE 15.1 The
More informationToday s Lecture: Atmosphere finish primitive equations, mostly thermodynamics
Today s Lecture: Atmosphere finish primitive equations, mostly thermodynamics Reference Peixoto and Oort, Sec. 3.1, 3.2, 3.4, 3.5 (but skip the discussion of oceans until next week); Ch. 10 Thermodynamic
More informationLecture 6 Free Energy
Lecture 6 Free Energy James Chou BCMP21 Spring 28 A quick review of the last lecture I. Principle of Maximum Entropy Equilibrium = A system reaching a state of maximum entropy. Equilibrium = All microstates
More informationFinal Exam for Physics 176. Professor Greenside Wednesday, April 29, 2009
Print your name clearly: Signature: I agree to neither give nor receive aid during this exam Final Exam for Physics 76 Professor Greenside Wednesday, April 29, 2009 This exam is closed book and will last
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 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 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 informationMolecular Dynamics for Microsystems. Introduction Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany
Molecular Dynamics for Microsystems. Introduction Evgenii B. Rudnyi and Jan G. Korvink IMTEK Albert Ludwig University Freiburg, Germany Preliminaries Learning Goals Why Molecular Simulation? Molecular
More informationINTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition
INTRODUCTION TO о JLXJLA Из А lv-/xvj_y JrJrl Y üv_>l3 Second Edition Kerson Huang CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is an imprint of the Taylor & Francis Group an Informa
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 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 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 informationTitle Super- and subcritical hydration of Thermodynamics of hydration Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2000), 8109 Issue Date 2000-05-08 URL http://hdl.handle.net/2433/50350
More informationElectrolytes. Chapter Basics = = 131 2[ ]. (c) From both of the above = = 120 8[
Chapter 1 Electrolytes 1.1 Basics Here we consider species that dissociate into positively and negatively charged species in solution. 1. Consider: 1 H (g) + 1 Cl (g) + ()+ () = { } = (+ )+ ( ) = 167[
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 informationUNIVERSITY OF OSLO FACULTY OF MATHEMATICS AND NATURAL SCIENCES
UNIVERSITY OF OSLO FCULTY OF MTHEMTICS ND NTURL SCIENCES Exam in: FYS430, Statistical Mechanics Day of exam: Jun.6. 203 Problem :. The relative fluctuations in an extensive quantity, like the energy, depends
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 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 informationIntro. Each particle has energy that we assume to be an integer E i. Any single-particle energy is equally probable for exchange, except zero, assume
Intro Take N particles 5 5 5 5 5 5 Each particle has energy that we assume to be an integer E i (above, all have 5) Particle pairs can exchange energy E i! E i +1andE j! E j 1 5 4 5 6 5 5 Total number
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 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 informationIMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER
IMPROVED METHOD FOR CALCULATING SURFACE TENSION AND APPLICATION TO WATER ABSTRACT Hong Peng 1, Anh V Nguyen 2 and Greg R Birkett* School of Chemical Engineering, The University of Queensland Brisbane,
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 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 informationANSWERS TO PROBLEM SET 1
CHM1485: Molecular Dynamics and Chemical Dynamics in Liquids ANSWERS TO PROBLEM SET 1 1 Properties of the Liouville Operator a. Give the equations of motion in terms of the Hamiltonian. Show that they
More informationElementary Lectures in Statistical Mechanics
George DJ. Phillies Elementary Lectures in Statistical Mechanics With 51 Illustrations Springer Contents Preface References v vii I Fundamentals: Separable Classical Systems 1 Lecture 1. Introduction 3
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 informationThermal & Statistical Physics Study Questions for the Spring 2018 Department Exam December 6, 2017
Thermal & Statistical Physics Study Questions for the Spring 018 Department Exam December 6, 017 1. a. Define the chemical potential. Show that two systems are in diffusive equilibrium if 1. You may start
More informationCan a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water?
Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? Ruhong Zhou 1 and Bruce J. Berne 2 1 IBM Thomas J. Watson Research Center; and 2 Department of Chemistry,
More informationMolecular Aggregation
Molecular Aggregation Structure Analysis and Molecular Simulation of Crystals and Liquids ANGELO GAVEZZOTTI University of Milano OXFORD UNIVERSITY PRESS Contents PART I FUNDAMENTALS 1 The molecule: structure,
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 informationConcentrating on the system
Concentrating on the system Entropy is the basic concept for discussing the direction of natural change, but to use it we have to analyze changes in both the system and its surroundings. We have seen that
More information(i) T, p, N Gibbs free energy G (ii) T, p, µ no thermodynamic potential, since T, p, µ are not independent of each other (iii) S, p, N Enthalpy H
Solutions exam 2 roblem 1 a Which of those quantities defines a thermodynamic potential Why? 2 points i T, p, N Gibbs free energy G ii T, p, µ no thermodynamic potential, since T, p, µ are not independent
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 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 informationFluid Equations for Rarefied Gases
1 Fluid Equations for Rarefied Gases Jean-Luc Thiffeault Department of Applied Physics and Applied Mathematics Columbia University http://plasma.ap.columbia.edu/~jeanluc 23 March 2001 with E. A. Spiegel
More information1. Thermodynamics 1.1. A macroscopic view of matter
1. Thermodynamics 1.1. A macroscopic view of matter Intensive: independent of the amount of substance, e.g. temperature,pressure. Extensive: depends on the amount of substance, e.g. internal energy, enthalpy.
More informationPhysics 4230 Final Examination 10 May 2007
Physics 43 Final Examination May 7 In each problem, be sure to give the reasoning for your answer and define any variables you create. If you use a general formula, state that formula clearly before manipulating
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 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 informationChapter 6 Thermodynamic Properties of Fluids
Chapter 6 Thermodynamic Properties of Fluids Initial purpose in this chapter is to develop from the first and second laws the fundamental property relations which underlie the mathematical structure of
More informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
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 informationFree energy, electrostatics, and the hydrophobic effect
Protein Physics 2016 Lecture 3, January 26 Free energy, electrostatics, and the hydrophobic effect Magnus Andersson magnus.andersson@scilifelab.se Theoretical & Computational Biophysics Recap Protein structure
More informationAtoms, electrons and Solids
Atoms, electrons and Solids Shell model of an atom negative electron orbiting a positive nucleus QM tells that to minimize total energy the electrons fill up shells. Each orbit in a shell has a specific
More informationINTERMOLECULAR AND SURFACE FORCES
INTERMOLECULAR AND SURFACE FORCES SECOND EDITION JACOB N. ISRAELACHVILI Department of Chemical & Nuclear Engineering and Materials Department University of California, Santa Barbara California, USA ACADEMIC
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 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 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 informationList of Comprehensive Exams Topics
List of Comprehensive Exams Topics Mechanics 1. Basic Mechanics Newton s laws and conservation laws, the virial theorem 2. The Lagrangian and Hamiltonian Formalism The Lagrange formalism and the principle
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 information