hydrate systems Gránásy Research Institute for Solid State Physics & Optics H-1525 Budapest, POB 49, Hungary László
|
|
- Penelope Hood
- 5 years ago
- Views:
Transcription
1 Phase field theory of crystal nucleation: Application to the hard-sphere and CO 2 Bjørn Kvamme Department of Physics, University of Bergen Allégaten 55, N-5007 N Bergen, Norway László Gránásy Research Institute for Solid State Physics & Optics H-1525 Budapest, POB 49, Hungary hydrate systems
2
3 Crystalline solidification: nucleation & growth Phase field theory
4 Nucleation Heterophase Heterophase fluctuations: fluctuations: Wolde & Frenkel,, 1996 Gasser et al., 2001 Auer & Frenkel,, 2001 Yonezawa,, 1991 Classical Classical (sharp (sharp intf.) intf.) picture: picture: cluster cluster free free energy energy (W)( (W) ) = = volumetric volumetric + + interfacial interfacial contributions contributions Crystal-liquid Crystal-liquid liquid interface interface Crit. Crit. fluct. fluct.. typically typically ~ ~ molecules molecules Diffuse interface model is is needed!!!!!! (e.g. (e.g. phase phase field field theory) theory) Davidchack & Laird, 1998
5 II. Phase-field theory of nucleation in hard sphere system Free energy functional: (standard phase field theory, see e.g. Warren & Boettinger 1995, adopted to HS) F bt 2 = d r ( m ) + f ( m, φ ) 2 where f(m,f) ) = w T g(m) ) + [1 - p(m)] f S (f) + p(m) f L (f) g(m) ) = ¼ m 2 (1 - m) 2, p(m) ) = m 3 (10-15m + 6m6 2 ), f L,S =polynomials reproducing the MD results, Model parameters: w & b intf. intf.. free free energy energy & intf. intf.. thickness they they can can be be fixed fixed in in equilibrium & predictions be be made made for for nucleation without adjustable parameters
6 Critical fluctuations: Extremum of grand potential functional (unstable equilibrium): W = F - m dr r Euler-Lagrange eqs. for phase field m (non-conserved) conserved): δω ω ω ω 0 = = = bt δm m m m 2 m and for volume fraction f (conserved): 0 = δω δφ = ω φ ω φ = ω φ ω φ
7
8 Critical fluctuations: spherical symmetry Phase-field: m + 2/r m = Dm(m)/( /(bt) Volume fraction: f = f(m) implicit equation Boundary cond.: r fi : m fi1, 1, f fi f ; r = 0: = 0: m, f = 0 W*
9 Hard sphere system: - Thermodynamics known accurately (Hall, 1970) - Interface profiles for density available 10% - 90% interface thickness, d (Davidchack & Laird, 1998) - Interface free energy from MD (Davidchack & Laird, 2000) orientation g s 2 /kt (111) (100) (110) g & d b & w
10 Free energy surface: For σ = 890 nm & (111) interface
11 (X-X min )/(X max -X min ) Equilibrium profiles: z/σ (111) Normalized distributions: (X-X( min )/(X max red: MD peak height blue: MD filtered density black solid: PFT [1-m(z)] black dashed: PFT volume fraction max -X min ) good agreement for (111) & (100) less satisfactory for (110) (X-X min )/(X max -X min ) (X-X min )/(X max -X min ) Compare solid line with solid, dashed with dashed z/σ (100) z/σ (110)
12 Nucleation : : Radial order parameter profiles: solid line: 1-m dashed line: (φ-φ )/(φ S -φ ) 0.6 X r/σ calculated with average interface properties
13 Nucleation barrier: No adjustable parameters Classical droplet model overestimates J by 3 to 5 orders of magnitude PFT upper curve: (110) central curve: average lower curve: (111) W*/kT CNT Φ calculated assuming spherical shape triangles: MC (Auer( & Frenkel,, 2001)
14 Tolman length: Curvature correction by Tolman (1949): g(r)» g / (1+ 2d2 T /R) δ T = Tolman length Route (i): δ T = R e -R p (exact) Route(ii): δ T = (R( p /2)[γ /γ 1] Route(ii): δ T = (R e /2)[γ /γ 1] R p & R e : radius of surf. of tension & of equimolar surface. solid line: route (i) dashed line: route (ii) dotted line: route (iii) δ T /σ Φ calculated assuming spherical shape
15 III. Phase-field theory of nucleation of hydrate nucleation Free energy functional: (Warren & Boettinger 1995) F bt 2 = d r ( φ ) + f ( φ, c) 2 where f(f,c) = w T g(f) ) + [1 - p(f)] f + p(f) f S L g(f) ) = ¼ f 2 (1 - f) 2, p(f) ) = f 3 (10-15f + 6f6 2 ), f L,S = polynomials for aqueous solution of CO 2 and CO 2 -hydrate Model parameters: w & b intf. intf.. free free energy energy & intf. intf.. thickness calculation of of nucl. nucl.. barrier barrier can can be be made made without adjustable parameters
16
17 Thermodynamics A thermodynamic description of the hydrate phase according to the van der Waal Platteuw approach requires the chemical potential of CO2 as well as the impact of CO2 on the inclusion in the cavity.
18 Note that: - the use of the DFT (Density Functional Theory) or PFT (Phase Field Theory) requires continuous description of the thermodynamics in the transitions from one phase to the other (no sharp interface). - consistency is required for all components in all phases - Stochiometry of the hydrate is an implicit function that is proportional to the exponent of the chemical potential of the guest molecule through the cavity partition function I.e.: the most proper way to describe the guest molecule chemical potential is so that it is consistent with the same model that describes the interactions with the cavity water molecules
19 a g = e x l n p[ ( β h ( µ ω ) h ( ω g ) d ω u )] β = 1 / k B T a= exp ([ βµ gu )] g= ln( hhd βω )() ω β= / 1kTB Extended adsorbtion Theory 0 µ = µ ν RT ln( 1 + aij ) w w i i j Single cavity integration (small guest molecules) Harmonic oscillator approach (large guest molecules) π a b V cage [ ( 1 β µ β ln b) ] = exp + = ( m / β 2πh ) [ βw v ] exp ( ) 2 3/ 2 dv ( ) [ β µ ] a = exp g u g = ln( βhω ) h( ω ) dω β = 1 / k T B f µ + x µ S Hydrate Hydrate = x CO 2 CO 2 Hydrate w Hydrate w
20 θ ij a 1 + ij a ij where θ ij is the filling fraction of type i guest in cavity type j Assuming that CO 2 only fills large cavity: x Hydrate CO 2 = 3 23 θ ij I.e.: The stochiometry of the hydrate Is an implicit function of the progress of The nucleation and has to be treated as such
21 Lquid solution description, CO2 µ aqueous CO 2 r aqueous r ( T, P, x) = µ CO 2 ( T, P, x) + RT ln( x CO 2γ CO 2 ) where superscript aqueous denotes liquid water phase and infinity symbol denote infinite dilution as reference state (activity coefficient of CO 2 goes to unity as the molefraction of CO 2 approaches zero). Infinite dilution chemical potential for model system From Molecular Dynamics simulations using Thermodynamic integration. Activity coefficients fitted to equilibrium solubility data.
22 Lquid solution description, H2O aqueous r aqueous r µ H O ( T, P, x) = µ H 2O pure ( T, P, x) + RT ln( xh 2Oγ H 2 2 O ) Activity coefficient for water calculated using Gibbs-Duhem f L = x CO aqueous r aqueous r 2µ CO2 ( T, P, x) + xh2 Oµ H2O ( T, P, x)
23 Free energy surface: b & w γ ice = 29.1 mj/m 2 & d = 1 nm 10-90
24 Radial profiles: Nucleation : : W* = J W* CNT = J φ, c CO red: 1-φ blue: c CO2 dashed: CNT Initial CO 2 content: 3.3% (saturated) r (A) b & w γ ice = 29.1 mj/m 2 & d = 1 nm 10-90
25 Plans: (1) Interface profiles & interface free energy from Molecular Dynamics (2) Full dynamics phase field theory for nucleation and growth to model hydrate formation under realistic conditions
26 Work in progress Interface free energy using a thermodynamic cycle and thermodynamic integration (Mezei s approach). Cycle sismilar to the approach of Broughton and Gilmer (1986), as also applied by several other groups. Direct calculation of absolute chemical potential of CO2 using thermodynamic integrations (Overheads describing the approaches is available)
27 Summary: (1) The phase field theory: - predicts the height of the nucleation barrier reasonably for the hard sphere liquid (without adjustable parameters); - the Tolman length decreases with increasing size of the fluctuations. (2) We formulated a phase field theory for CO 2 hydrate nucleation in aqueous CO 2 solution Financed by ESA Prodex No /00/NL/SFe Forms part of ESA MAP Project AO
28
29 Hydrate melting Surface tension? γ=70 ±9 mn/m Interface thickness?
30 INTERFACIAL TENSION Interface Run, ns T, K P t, P n, atm h z, Å γ, mn/m γ, mn/m atm Sim liquidliquid Sim liquidliquid Sim liquidliquid Exp a liquidliquid 275 < Sim liquid-gas NVT sat NVT sat 39 6 Exp a liquid-gas 275 < Exp b liquid-gas Exp b liquid-gas Sim liquid-gas NVT half NVT half Exp b liquid-gas Sim this work; NVT sat stands for NVT simulation of liquid-gas interface corresponding to saturated CO 2 density; NVT half, for simulation at half saturated density; a experimental data [13]; b experimental data [12].
31 CO2 hydrate Water System: Characteristic Density Profiles at 22 and 48 ps
32 Trajectories of CH 4 guests Liquid-phase CO 2 molecule trapped near the interface
33 TI Cleaving aproach to interfacial free energy estimation: Reversible path Broughton & Gilmer: Molecular Dynamics of Interface. VI J. Chem. Phys., Vol. 84, No. 10, 15 May 1986
34 POLYNOMIAL PATH FOR THERMODYNAMIC INTEGRATION Mezei MJ. Comp Chem, 13, 651 (1992) E(λ, X N ) = f 1 (λ)*e 1 (X N ) + f 0 (λ)*e 0 (X N ) f 0 (0) =0 and f 0 (1) = 1 f 1 (0) =1 and f 1 (1) = 0 A = A = 1 < E ( λ λ > d λ λ 1 A ) / 0 0
35 E = ( ulj + ues ) = ij i E i Polynomial path, Mezei (1992), introduces different exponents for potential terms differing in inverse distance powers: A = 1 [ k 1 < ] i k i i > k i 0 λ E d λ λ
36 Chemical potential of CO 2 in structure I hydrate 1. Conversion of CO 2 to CH 4 via TI a) scale down charges b) convert LJ parameters 2. Single-guest scaling in CH 4 hydrate 3. Steps 1+2 CO 2 chemical potential
Nucleation 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 informationSurface analysis algorithms in the mardyn program and the ls1 project
Surface analysis algorithms in the mardyn program and the ls1 project Stuttgart, 15 th December 1 M. T. Horsch Surface tension The virial route Bakker-Buff equation: γ R 2 out in dz z Normal pressure decays
More informationof Nebraska - Lincoln
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry 10-1-2006 Homogeneous nucleation at high supersaturation
More informationCrystal nucleation for a model of globular proteins
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 17 1 MAY 2004 Crystal nucleation for a model of globular proteins Andrey Shiryayev and James D. Gunton Department of Physics, Lehigh University, Bethlehem,
More informationInterfaces and interfacial energy
Interfaces and interfacial energy 1/14 kinds: l/g }{{ l/l } mobile s/g s/l s/s Example. Estimate the percetage of water molecules on the surface of a fog droplet of diameter (i) 0.1 mm (naked eye visibility
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 informationSurface property corrected modification of the classical nucleation theory
CCP5 Annual Meeting Surface property corrected modification of the classical nucleation theory Sheffield Hallam University, September 15, 2010 Martin Horsch, Hans Hasse, and Jadran Vrabec The critical
More informationMOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE
MOLECULAR DYNAMICS SIMULATION OF HETEROGENEOUS NUCLEATION OF LIQUID DROPLET ON SOLID SURFACE Tatsuto Kimura* and Shigeo Maruyama** *Department of Mechanical Engineering, The University of Tokyo 7-- Hongo,
More informationVIII. Phase Transformations. Lecture 38: Nucleation and Spinodal Decomposition
VIII. Phase Transformations Lecture 38: Nucleation and Spinodal Decomposition MIT Student In this lecture we will study the onset of phase transformation for phases that differ only in their equilibrium
More informationAtomic Transport & Phase Transformations Lecture III-2
Atomic Transport & Phase Transformations Lecture III-2 PD Dr. Nikolay Zotov zotov@imw.uni-stuttgart.de Atomic Transport & Phase Transformations Part III Lectures Solid State Reactions Short Description
More informationwww.mrs.org/publications/bulletin Nucleation and the Solid Liquid Interfacial Free Energy David T.Wu, László Gránásy, and Frans Spaepen Abstract This article reviews the current understanding of the fundamentals
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 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 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 informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More 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 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 informationA MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE
A MOLECULAR DYNAMICS SIMULATION OF A BUBBLE NUCLEATION ON SOLID SURFACE Shigeo Maruyama and Tatsuto Kimura Department of Mechanical Engineering The University of Tokyo 7-- Hongo, Bunkyo-ku, Tokyo -866,
More informationSurface Tension of the Vapor Liquid Interface with Finite Curvature
Colloid Journal, Vol. 65, No. 4, 00, pp. 44045. Translated from Kolloidnyi Zhurnal, Vol. 65, No. 4, 00, pp. 480494. Original Russian Text Copyright 00 by Zhukhovitskii. Surface Tension of the VaporLiquid
More informationA Thermodynamic Study of Methane Hydrates Formation In Glass Beads
AJChE 2016, Vol. 16, No. 1, 15 22 A Thermodynamic Study of Methane Hydrates Formation In Glass Beads Tintin Mutiara *,1,2 Budhijanto 1 I Made Bendiyasa 1 Imam Prasetyo 1 1 Department of Chemical Engineering,
More informationDirect determination of the Tolman length from the bulk pressures of liquid drops via molecular dynamics simulations
Direct determination of the Tolman length from the bulk pressures of liquid drops via molecular dynamics simulations Alan E. van Giessen, and Edgar M. Blokhuis Citation: The Journal of Chemical Physics
More informationHill s nano-thermodynamics is equivalent with Gibbs thermodynamics for curved surfaces
Hill s nano-thermodynamics is equivalent with Gibbs thermodynamics for curved surfaces arxiv:1806.03444v1 [physics.chem-ph] 9 Jun 2018 Dick Bedeaux and Signe Kjelstrup Porelab, Department of Chemistry,
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 informationStep 1. Step 2. g l = g v. dg = 0 We have shown that over a plane surface of water. g v g l = ρ v R v T ln e/e sat. this can be rewritten
The basic question is what makes the existence of a droplet thermodynamically preferable to the existence only of water vapor. We have already derived an expression for the saturation vapor pressure over
More informationSpeeding up path integral simulations
Speeding up path integral simulations Thomas Markland and David Manolopoulos Department of Chemistry University of Oxford Funded by the US Office of Naval Research and the UK EPSRC Outline 1. Ring polymer
More informationHeterogeneous nucleation on mesoscopic wettable particles: A hybrid thermodynamic/densityfunctional
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published esearch - Department of Chemistry 7--00 Heterogeneous nucleation on mesoscopic wettable
More informationNumerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry
Numerical Solution of Integral Equations in Solidification and Melting with Spherical Symmetry V. S. Ajaev and J. Tausch 2 Southern Methodist University ajaev@smu.edu 2 Southern Methodist University tausch@smu.edu
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 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 informationPhases of matter and phase diagrams
Phases of matter and phase diagrams Transition to Supercritical CO2 Water Ice Vapor Pressure and Boiling Point Liquids boil when the external pressure equals the vapor pressure. Temperature of boiling
More informationCalculation of the crystal-melt interfacial free energy of succinonitrile from molecular simulation
THE JOURNAL OF CHEMICAL PHYSICS 124, 044707 2006 Calculation of the crystal-melt interfacial free energy of succinonitrile from molecular simulation Xiaobing Feng and Brian B. Laird a Department of Chemistry,
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 informationA patching model for surface tension of spherical droplet and Tolman length. II
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Xiao Cheng Zeng Publications Published Research - Department of Chemistry -5-999 A patching model for surface tension of
More informationNUCLETION OF GAS HYDRATES
NUCLETION OF GAS HYDRATES Dimo Kashchiev Institute of Physical Chemistry, Bulgarian Academy of Sciences, Sofia Detailed presentation of the results in: D.Kashchiev, A.Firoozabadi, J. Crystal Growth 241
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 informationVapor-to-Droplet Transition in a Lennard-Jones Fluid: Simulation Study of Nucleation Barriers Using the Ghost Field Method
5962 J. Phys. Chem. B 2005, 109, 5962-5976 Vapor-to-Droplet Transition in a Lennard-Jones Fluid: Simulation Study of Nucleation Barriers Using the Ghost Field Method Alexander V. Neimark* and Aleksey Vishnyakov
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More 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 informationAtomic Transport & Phase Transformations Lecture III-1
Atomic Transport & Phase Transformations Lecture III-1 PD Dr. Nikolay Zotov zotov@imw.uni-stuttgart.de Atomic Transport & Phase Transformations Part III Lectures Solid State Reactions Short Description
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 informationFundamentals of Hydrates, Climate Perspectives, and Energy Potentials
CCUS Student Week 2018 Fundamentals of Hydrates, Climate Perspectives, and Energy Potentials Luis Zerpa Center for Hydrate Research Colorado School of Mines October 18, 2018 Golden, CO What are Gas Hydrates?
More informationSurface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water
Surface thermodynamics of planar, cylindrical, and spherical vapour-liquid interfaces of water Gabriel V. Lau, Ian J. Ford, Patricia A. Hunt, Erich A. Müller, and George Jackson Citation: The Journal of
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 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 informationSupplemental Material for Curvature-dependence of the liquid-vapor surface tension beyond the Tolman approximation
Supplemental Material for Curvature-dependence of the liquid-vapor surface tension beyond the Tolman approximation Nicolas Bruot and Frédéric Caupin Institut Lumière Matière, UMR536 Université Claude Bernard
More informationLecture 3. Phenomena at Liquid-gas and Liquid-Liquid interfaces. I
Lecture 3 Phenomena at Liquid-gas and Liquid-Liquid interfaces. I Adsorption at Gas-Liquid interface Measurements of equilibrium adsorption surface tension measurements (Wilhelmy plate) surface analysis
More informationChapter 17. Fundamentals of Atmospheric Modeling
Overhead Slides for Chapter 17 of Fundamentals of Atmospheric Modeling by Mark Z. Jacobson Department of Civil & Environmental Engineering Stanford University Stanford, CA 94305-4020 August 21, 1998 Mass
More informationEquilibria in Materials
009 fall Advanced Physical Metallurgy Phase Equilibria in Materials 10. 13. 009 Eun Soo Park Office: 33-316 Telephone: 880-71 Email: espark@snu.ac.kr Office hours: by an appointment 1 Contents for previous
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 informationTHE MODIFIED YOUNG S EQUATION FOR THE CONTACT ANGLE OF A SMALL SESSILE DROP FROM AN INTERFACE DISPLACEMENT MODEL
International Journal of Modern Physics B, Vol. 13, No. 7 (1999) 355 359 c World Scientific Publishing Company THE MODIFIED YOUNG S EQUATION FOR THE CONTACT ANGLE OF A SMALL SESSILE DROP FROM AN INTERFACE
More informationWhere do ions solvate?
PRAMANA c Indian Academy of Sciences Vol. 64, No. 6 journal of June 25 physics pp. 957 961 YAN LEVIN Instituto de Física, Universidade Federal do Rio Grande do Sul, Caixa Postal 1551, CEP 9151-97, Porto
More informationHyeyoung Shin a, Tod A. Pascal ab, William A. Goddard III abc*, and Hyungjun Kim a* Korea
The Scaled Effective Solvent Method for Predicting the Equilibrium Ensemble of Structures with Analysis of Thermodynamic Properties of Amorphous Polyethylene Glycol-Water Mixtures Hyeyoung Shin a, Tod
More information6 Hydrophobic interactions
The Physics and Chemistry of Water 6 Hydrophobic interactions A non-polar molecule in water disrupts the H- bond structure by forcing some water molecules to give up their hydrogen bonds. As a result,
More informationThe transport of mass at the nano-scale during. evaporation of droplets: the Hertz-Knudsen equation. at the nano-scale
Page of 0 0 0 0 0 The transport of mass at the nano-scale during evaporation of droplets: the Hertz-Knudsen equation at the nano-scale Marcin Zientara,, Daniel Jakubczyk, Marek Litniewski, and Robert Hołyst
More informationChapter 9 Generation of (Nano)Particles by Growth
Chapter 9 Generation of (Nano)Particles by Growth 9.1 Nucleation (1) Supersaturation Thermodynamics assumes a phase change takes place when there reaches Saturation of vapor in a gas, Saturation of solute
More informationToward a Quantitative Theory of Ultrasmall Liquid Droplets and VaporsLiquid Nucleation
4988 Ind. Eng. Chem. Res. 2008, 47, 4988 4995 Toward a Quantitative Theory of Ultrasmall Liquid Droplets and VaporsLiquid Nucleation Zhidong Li and Jianzhong Wu* Department of Chemical and EnVironmental
More informationINTRODUCTION TO MODERN THERMODYNAMICS
INTRODUCTION TO MODERN THERMODYNAMICS Dilip Kondepudi Thurman D Kitchin Professor of Chemistry Wake Forest University John Wiley & Sons, Ltd CONTENTS Preface xiii PART I THE FORMALIS1VI OF MODERN THER1VIODYNAMICS
More informationarxiv: v1 [cond-mat.stat-mech] 27 Nov 2015
The solid-liquid interfacial free energy out of equilibrium arxiv:1511.8668v1 [cond-mat.stat-mech] 27 Nov 215 Bingqing Cheng, 1 Gareth A. Tribello, 2 and Michele Ceriotti 1, 1 Laboratory of Computational
More informationOrder-parameter-based Monte Carlo simulation of crystallization
THE JOURNAL OF CHEMICAL PHYSICS 124, 134102 2006 Order-parameter-based Monte Carlo simulation of crystallization Manan Chopra Department of Chemical Engineering, University of Wisconsin, Madison, Wisconsin
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 informationPHASE-FIELD ANALYSIS OF NUCLEATION PHENOMENA Ph.D. thesis. Gyula TÓTH
PHASE-FIELD ANALYSIS OF NUCLEATION PHENOMENA Ph.D. thesis Gyula TÓTH Supervisors: Prof. László Gránásy Research Institute for Solid State Physics and Optics, Budapest, Hungary Brunel Centre for Advanced
More informationA coherent picture for water at extreme negative pressure
A coherent picture for water at extreme negative pressure Mouna El Mekki Azouzi 1, Claire Ramboz 2, Jean-François Lenain 3 and Frédéric Caupin 1 1 Laboratoire de Physique de la Matière Condensée et Nanostructures,
More informationThe role of specific cations and water entropy on the stability of branched DNA motif structures
The role of specific cations and water entropy on the stability of branched DNA motif structures Tod A Pascal, William A Goddard III, Prabal K Maiti and Nagarajan Vaidehi 1. Figures Supporting Information
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
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 informationDendritic solidification of binary alloys
University of Ljubljana Faculty of Mathematics and Physics Department of Physics Seminar I b - 3rd year, 2nd cycle Dendritic solidification of binary alloys Author: Tadej Dobravec Mentor: prof. dr. Božidar
More informationMetropolis, 2D Ising model
Metropolis, 2D Ising model You can get visual understanding from the java applets available, like: http://physics.ucsc.edu/~peter/ising/ising.html Average value of spin is magnetization. Abs of this as
More informationThermodynamic expansion of nucleation free-energy barrier and size of critical nucleus near the vapor-liquid coexistence
JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 FEBRUARY 1999 hermodynamic expansion of nucleation free-energy barrier and size of critical nucleus near the vapor-liquid coexistence Kenichiro Koga
More informationNucleation: Theory and Applications to Protein Solutions and Colloidal Suspensions
Nucleation: Theory and Applications to Protein Solutions and Colloidal Suspensions Richard P. Sear Department of Physics, University of Surrey Guildford, Surrey GU2 7XH, United Kingdom r.sear@surrey.ac.uk
More informationGFD 2006 Lecture 2: Diffusion-controlled solidification
GFD 2006 Lecture 2: Diffusion-controlled solidification Grae Worster; notes by Victor Tsai and Dan Goldberg March 15, 2007 1 Finishing off Lecture 1 As shown in Lecture 1, an approximation for the diffusion
More informationTransition Theory Abbreviated Derivation [ A - B - C] # E o. Reaction Coordinate. [ ] # æ Æ
Transition Theory Abbreviated Derivation A + BC æ Æ AB + C [ A - B - C] # E A BC D E o AB, C Reaction Coordinate A + BC æ æ Æ æ A - B - C [ ] # æ Æ æ A - B + C The rate of reaction is the frequency of
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 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 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 informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: hermodynamics and Kinetics Notes for Lecture 14 I. HE CLAEYRON EQUAION he Clapeyron attempts to answer the question of what the shape of a two-phase coexistence line is. In the - plane,
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 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 informationMonte Carlo Calculations of Effective Surface Tension for Small Clusters
Monte Carlo Calculations of Effective Surface Tension for Small Clusters Barbara N. Hale Physics Department and Cloud and Aerosol Science Laboratory, University of Missouri- Rolla, Rolla, MO 65401, USA
More informationTemperature and Pressure Controls
Ensembles Temperature and Pressure Controls 1. (E, V, N) microcanonical (constant energy) 2. (T, V, N) canonical, constant volume 3. (T, P N) constant pressure 4. (T, V, µ) grand canonical #2, 3 or 4 are
More 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 informationSupporting Information: On Localized Vapor Pressure Gradients Governing Condensation and Frost Phenomena
Supporting Information: On Localized Vapor Pressure Gradients Governing Condensation and Frost Phenomena Saurabh Nath and Jonathan B. Boreyko Department of Biomedical Engineering and Mechanics, Virginia
More informationMolecular Dynamics Simulation Study of Transport Properties of Diatomic Gases
MD Simulation of Diatomic Gases Bull. Korean Chem. Soc. 14, Vol. 35, No. 1 357 http://dx.doi.org/1.51/bkcs.14.35.1.357 Molecular Dynamics Simulation Study of Transport Properties of Diatomic Gases Song
More informationPhase field approach to heterogeneous crystal nucleation in alloys
Phase field approach to heterogeneous crystal nucleation in alloys James A. Warren, 1 Tamás Pusztai, 2 László Környei, 2 and László Gránásy 3 1 Metallurgy Division, National Institute of Standards and
More informationA New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform Phase Fluid
Commun. Theor. Phys. (Beijing, China) 39 (2003) pp. 231 237 c International Academic Publishers Vol. 39, No. 2, February 15, 2003 A New Uniform Phase Bridge Functional: Test and Its Application to Non-uniform
More informationAtomistic molecular simulations for engineering applications: methods, tools and results. Jadran Vrabec
Atomistic molecular simulations for engineering applications: methods, tools and results Jadran Vrabec Motivation Simulation methods vary in their level of detail The more detail, the more predictive power
More informationL = I ω = const. I = 2 3 MR2. When the balloon shrinks (because the air inside it cools down), the moment of inertia decreases, R = 1. L = I ω = I ω.
PHY 30 K. Solutions for mid-term test #3. Problem 1: Out in space, there are no forces acting on the balloon except gravity and hence no torques (with respect to an axis through the center of mass). Consequently,
More informationUnusual Entropy of Adsorbed Methane on Zeolite Templated Carbon. Supporting Information. Part 2: Statistical Mechanical Model
Unusual Entropy of Adsorbed Methane on Zeolite Templated Carbon Supporting Information Part 2: Statistical Mechanical Model Nicholas P. Stadie*, Maxwell Murialdo, Channing C. Ahn, and Brent Fultz W. M.
More informationCLUSTERING AND WATER ACTIVITY IN CONCENTRATED SUCROSE SOLUTIONS
EURO FOOD S WATER CLUSTERING AND WATER ACTIVITY IN CONCENTRATED SUCROSE SOLUTIONS Jean GENOTELLE 1, Mohamed MATHLOUTHI 2,BarbaraROGE 2, Maciej STARZAK 3 1 Former Director of GTS Sèvres, France 1 Université
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 informationPhase-field modeling of binary alloy solidification with coupled heat and solute diffusion
PHYSICAL REVIEW E 69, 051607 (004) Phase-field modeling of binary alloy solidification with coupled heat and solute diffusion J. C. Ramirez, 1 C. Beckermann, 1, * A. Karma, and H.-J. Diepers,3 1 Department
More informationTwo-Dimensional Spin-Polarized Hydrogen at Zero
J Low Temp Phys (2013) 171:685 692 DOI 10.1007/s10909-012-0756-7 Two-Dimensional Spin-Polarized Hydrogen at Zero Temperature L. Vranješ Markić J. Boronat Received: 13 July 2012 / Accepted: 16 September
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 informationSpecific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate
Specific Ion Solvtion in Ethylene Carbonate and Propylene Carbonate A. Arslanargin, A. Powers, S. Rick, T. Pollard, T. Beck Univ Cincinnati Chemistry Support: NSF, OSC TSRC 2016 November 2, 2016 A. Arslanargin,
More informationFree energy concept Free energy approach LBM implementation Parameters
BINARY LIQUID MODEL A. Kuzmin J. Derksen Department of Chemical and Materials Engineering University of Alberta Canada August 22,2011 / LBM Workshop OUTLINE 1 FREE ENERGY CONCEPT 2 FREE ENERGY APPROACH
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 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 informationNUCLEATION IN FINITE SYSTEMS: THEORY AND COMPUTER SIMULATION*t. M. RAO and B. J. BERNE. Chemistry Dept, Columbia University, New York, U.S.A.
NUCLEATION IN FINITE SYSTEMS: THEORY AND COMPUTER SIMULATION*t M. RAO and B. J. BERNE Chemistry Dept, Columbia University, New York, U.S.A. Abstract. The free energy of formation of a droplet in a finite
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 informationFernanda C Bononi 1, Ted Hullar 2, Cort Anastasio 2, Davide Donadio 1
Fernanda C Bononi 1, Ted Hullar 2, Cort Anastasio 2, Davide Donadio 1 1 Department of Chemistry 2 Department of Air, Land and Water Resources UC Davis Ice surface in climate modeling: sea water ice polar
More informationPrecipitation. Size! Shape! Size distribution! Agglomeration!
Precipitation Size! Shape! Size distribution! Agglomeration! Precipitation Four major questions: 1. Why do molecules/ions precipitate? 2. What determines the size? 3. What determines the size distribution?
More information