Joy Billups Woller A THESIS. Presented to the Faculty of. The Graduate College in the University of Nebraska
|
|
- Elwin Lang
- 5 years ago
- Views:
Transcription
1 DETERMINATION OF THE THERMODYNAMIC PROPERTIES OF FLUIDS BY GIBBS ENSEMBLE MONTE CARLO COMPUTER SIMULATIONS by Joy Billups Woller A THESIS Presented to the Faculty of The Graduate College in the University of Nebraska In Partial Fulfillment of the Requirements For the Degree of Master of Science Major: Chemistry Under the Supervision of Professor Xiao C. Zeng Lincoln, Nebraska July, 1997
2 DETERMINATION OF THE THERMODYNAMIC PROPERTIES OF FLUIDS BY GIBBS ENSEMBLE MONTE CARLO COMPUTER SIMULATIONS Adviser: Xiao C. Zeng Joy Billups Woller, M.S. University of Nebraska, 1997 Gibbs ensemble Monte Carlo (GEMC) is a useful technique for studying the thermodynamic properties of fluids. The theory and methodology of GEMC are explained, and the results of its use in two different projects are discussed. The first project involves a study of the effects of size, energy, and dipole moment parameters on the vapor liquid equilibria of binary mixtures containing Lennard-Jones (LJ) and Stockmayer (SM) components. The results of this work show that the phase coexistence properties of the mixtures are sensitive to the difference between the reduced dipole moments of the two components: an increase in this difference causes the phase envelope to widen for both the LJ SM and SM SM mixtures. Additional work which compares LJ SM and SM SM results to vapor liquid equilibria data simulated using an angleaveraged approximation of long-range dipolar interactions, the Keesom potential, suggests that the Keesom potential is not suitable to approximate the dipolar interactions for mixtures when the dipolar interaction is dominant. The second project involves
3 GEMC simulations of carbon dioxide, in which each molecule is modeled using three LJ sites with a point charge centered on each atom. Parameters for the carbon dioxide interaction potentials were generated using ab initio techniques, and the vapor liquid coexistence data produced were compared to experimental and previous simulation results. The ab initio potential parameters did not produce data consistent with either experimental or previous simulation results, which used parameters calculated from fitting to the experimental data. The disparity may be caused by the fact that ab initio calculations cannot properly account for the many interactions occurring in the liquid phase, or the fact that the higher stability of some CO 2 CO 2 orientations is not considered. While further attempts to generate potential parameters using ab initio methods may produce better results, fitting of experimental data clearly provides better LJ potential parameters for use in GEMC simulations of carbon dioxide at this point.
4 To Mom and Dad with love and thanks.
5 ACKNOWLEDGMENTS I would first like to thank my adviser, Professor Xiao C. Zeng, for the opportunity to venture into the interesting world of computer simulations and for his guidance in performing research. He has been a kind and generous teacher. I must also extend sincere thanks to Dr. Guangtu Gao of Beijing, China, who spent a year in Professor Zeng s lab as a visiting scholar. Dr. Gao s patient instruction and programming skill made our collaboration a valuable learning experience for me and a productive opportunity to explore fluid equilibria. Also, thank you to Ruben Parra, who provided the ab initio parameters for our carbon dioxide work. I would also like to acknowledge the following members of our group for their support: Hong Yuan, Kwang Jin Oh, Mark Gibson, Dr. Kenichiro Koga, and Dr. David Swanson. Thank you to Professors John J. Stezowski and William H. Braunlin, who graciously agreed to be on my committee and hear my oral defense. The work contained in this thesis was supported in part by the National Science Foundation, the Office of Naval Research, and the Petroleum Research Fund administered by the American Chemical Society. I credit (or blame, depending on my mood) two people with guiding me into this fascinating, and sometimes frustrating, world of chemistry. The first is my high school teacher, Mr. John McNeilly, who did a fantastic job of laying the solid foundation upon which I ve built all I know about chemistry. The second, Dr. Richard Shoemaker, has been my mentor for the past six years, and his encouragement and friendship helped get me through my darkest hours, i.e., organic chemistry. To my parents I owe a profound debt of gratitude for the appreciation of education they instilled in me, through their words and by example. I must also thank them for
6 teaching me the value of persistence, without which I couldn t have stuck with what has been, at times, a miserable exercise in masochism. Though Mom hasn t been able to help me with homework since I started long-division, at her knee (and sometimes across it) I learned enough about character and honesty to fill a few textbooks of my own. From Dad, I have learned the value of hard work, the power of logic, and a few choice swear words to use when a cow runs the wrong way. Thank you both for being terrific parents. I d also like to acknowledge my siblings for, well, their existence. I wouldn t be who I am today without them, so it s fortunate they re around to blame. Stephanie, you re a great kid and you have a wonderful future ahead. I hope you find something you truly love to do and pursue it with reckless abandon. Wes, you are the most honorable person I know, and I m intensely proud of you. I couldn t ask for a better brother. And Brandi, though you missed out on most of the sibling squabbles that went along with childhood in our family, I hope I made up for that deficiency while you lived with us. I m proud to call you my sister. Finally, I d like to thank my husband, Kevin, whose patient understanding and kind words are my source of strength and inspiration. While we face many uncertainties about our future (will we have to live in New Jersey?), the knowledge that we face it together makes all worries seem trivial. Thank you for being a wonderful husband and true friend. I love you.
7 CONTENTS Chapter 1. An Introduction to Computer Simulations Intermolecular Forces and Potential Models Simulation Methods Gibbs Ensemble Monte Carlo...7 References...12 Chapter 2. Gibbs Ensemble Monte Carlo of Binary Mixtures Containing Dipolar Molecules Introduction Potential Models Simulation Method Results Discussion...31 References...32 Chapter 3. Gibbs Ensemble Monte Carlo of Carbon Dioxide Using Ab Initio Size and Energy Parameters Introduction Potential Model Simulation Method Results Discussion...42 References...44 Appendix A...45 Appendix B...46
8 1 CHAPTER 1 AN INTRODUCTION TO COMPUTER SIMULATIONS The use of computers to predict thermodynamic properties has a relatively short history. However, recent rapid advances in computer power and data storage capacity combined with new simulation techniques have made computer simulations increasingly practical. The information gleaned from computer simulations of molecular systems can provide insight into the interactions of molecules and the quality of the theory describing chemical systems. Once the quality of the models and techniques used is verified, computer simulations can also provide invaluable data which may be used in predicting the thermodynamic properties of fluids, data which may be used by chemical engineers and environmental scientists, among others. In the past decade, a new computer simulation method has emerged as a useful tool in computing phase coexistence behavior of fluids. That method, Gibbs ensemble Monte Carlo, is discussed in this chapter, along with the other techniques and models needed in generating phase coexistence information by computer simulations. 1.1 Intermolecular Forces and Potential Models One of the most challenging tasks in performing computer simulations of molecular systems is coming up with a mathematical description of molecules which is both realistic and computationally feasible. While complex potential models, with
9 2 numerous adjustable parameters, can provide excellent descriptions of molecular interactions, they are often too computationally demanding to perform the Monte Carlo or molecular dynamics simulations required for calculating thermodynamic properties. On the other hand, the simplest potential models, with only a few adjustable parameters, don t necessarily describe reality. A compromise option involves putting a great deal of effort into optimizing the parameters of the simple potential models until they produce a reasonable approximation of reality. The work reported in the next two chapters is concerned primarily with two potential functions: the Lennard-Jones potential, which describes non-electrostatic attractive and repulsive interactions between pairs of molecules, and the Stockmayer potential, which describes dipolar interactions between pairs of molecules. By working with these relatively simple descriptions of complex interactions, we hope to better understand how they may be applied to predictive simulations used in chemical engineering and other fields The Lennard-Jones Potential Function The Lennard-Jones (LJ) potential is one of the most widely used intermolecular potentials, valued primarily for its computational simplicity, and is written as: σ ij Ur ( ij ) = ε ij r ij ij σ r ij
10 3 where σ ij is the LJ size parameter, ε ij is the LJ energy parameter, and r ij is the separation between particles i and j. The LJ function is plotted in Figure 1.1; both the independent and dependent variables are shown as reduced variables (see Appendix A). The LJ size parameter is equal to the intermolecular separation for which the energy is zero, which is easily shown to be σ = r, where r m m is the intermolecular separation at minimum U* energy, or U ( r) = ε. The LJ size and energy parameters are usually determined by choosing some physical property and fitting values of σ ε r* that property as calculated by the potential function to experimentally determined values. Figure 1.1 The Lennard-Jones potential energy function. Asterisks indicate reduced variables (see Appendix A) The Stockmayer Potential Function For two polar molecules, the contribution of the long-range electrostatic attraction between dipoles must be added to the description of intermolecular forces. The Stockmayer potential, used in this work to describe dipolar interactions, can be written with the LJ potential function as: 12 v v v v v v v v v σ ij σ ij µ i µ j ( µ i rij)( µ j rij) Ur ( ij, µ i, µ j ) = ε ij r 6 ij r ij rij rij
11 4 where v µ i and v µ j are the vectors for the embedded point dipoles in particles i and j; the other parameters are as in the LJ function. 1 The contribution from dipolar interactions can also be represented in scalar form as: U Dipole µ iµ j = 3 ( 2cosθicosθj sinθisinθjcos φij) 1.3 r ij where the angles representing the spatial relationships between particles i and j are shown in Figure 1.2. The simplicity of working with point dipoles makes the Stockmayer potential an attractive option in modeling polar fluids. Of course, since it ignores the effects of polarization and higher-order Figure 1.2 Spatial relationships between particles i and j used to calculate Stockmayer potentials. multipoles, it oversimplifies the interactions occurring in real fluids. However, in modeling dipolar and non-polar interactions of a polar fluid, the Stockmayer potential describes the essential aspects of a polar fluid. 2 In the work described in Chapter 2, the Stockmayer potential is examined without comparing it to real fluids.
12 5 1.2 Simulation Methods Since the Lennard-Jones function falls off rapidly as the intermolecular distance increases, the energy contribution at more than a few σ has a negligible effect on individual molecules. To reduce computational demands, a cutoff radius is set, beyond which intermolecular interactions are not explicitly calculated. In other words, when r r c, where r c is the cutoff radius, the intermolecular pair potential is set to zero. At the end of a simulation, a long-range correction term may be added to the total energy to account for the accumulated effects of the truncation. To allow a simulated system, with its limited number of molecules, to model bulk fluids, periodic boundary conditions are applied. To create periodic boundary conditions, the simulation cell (that which contains the molecules being simulated) is replicated in all dimensions around the cell. In Figure 1.3, a two-dimensional representation of a ninemolecule simulation is shown. One consequence of this replication is that molecules near the sides of the simulation cell experience intermolecular interactions with the periodic images of other molecules; without periodic boundary conditions, molecules near cell boundaries would experience very different environments than those in the center of the simulation cell. Another consequence of the periodic boundary conditions is that the transfer of a molecule (using either molecular dynamics or Monte Carlo methods) across a boundary results in an image of that molecule moving back into the simulation cell on the other side.
13 6 Figure 1.3 Two-dimensional representation of periodic boundary conditions. The center cell is the simulation cell; the gray cells hold the periodic images of the simulation molecules. Typically the box length in a simulation is longer than the cutoff radius applied to the potential function for nonpolar and uncharged molecules; as a result, these molecules do not interact with their periodic images. For long-range electrostatic interactions, however, special techniques must be applied to efficiently sum the interactions between a molecule and its periodic images. The specifics of the Ewald sum, the technique used in this research, are not explained in depth here; descriptions of the method are available
14 7 from a number of sources. 3 5 Details of the Ewald sum as they apply to the specific projects in this work are discussed in the next two chapters. 1.3 Gibbs Ensemble Monte Carlo In 1987, the procedure for Gibbs ensemble Monte Carlo (GEMC) simulations was first proposed by Panagiotopoulos 6 as a means of directly calculating the phase coexistence properties of either pure fluids or mixtures. GEMC has the advantage of being able to directly calculate two-phase equilibrium properties without requiring that the free energy be calculated as in conventional methods. 7 The goal of GEMC is to create a two-phase system, without an interface, in which the temperature, pressure, and chemical potential of all components in the two phases are equal. 8 However, there is no true constant-µpt ensemble, as extensive variables would be unbounded if only intensive parameters were specified. At least one extensive variable must be fixed, and for GEMC, that variable is the number of molecules N. 5 The effect is that the simulated system has coexisting phases in which the temperature, pressure, and chemical potential (or potentials, in the case of mixtures) are fixed. Since the actual value of µ is not specified (only the difference between different phases is set so that µ = 0) GEMC does not create a true constant-µpt ensemble, but such an ensemble is effectively simulated. The algorithm for performing a GEMC simulation involves applying three types perturbations to a system of two boxes, one containing gas and one containing liquid. The
15 8 perturbations require displacing a molecule within its current simulation box, changing the volume of the two boxes, and transferring a molecule from one box to another. There are two versions of GEMC, one which couples the volume changes so that the total volume of the vapor-liquid system remains constant (NVT) and a version in which the volume changes are not coupled but are used to ensure both boxes are at the same pressure (NPT). The NPT-GEMC method, in which the coexistence pressure is specified in advance, is better for modeling mixtures and is used in the binary mixture work described in Chapter 2, while the NVT-GEMC is used to simulate pure fluids and is adopted in the carbon dioxide work described in Chapter 3. 8 The first step in performing a GEMC simulation requires setting up two boxes, one for the liquid phase and one for the gas phase. These boxes are usually set up initially in a face-centered cubic structure, and each models a region deep within the bulk fluid. Periodic boundary conditions are applied, and no interface between the two phases exists in the simulation. After the phase boxes are set up, the three different types of perturbations are applied. Appendix B contains a flow diagram showing how these perturbations are implemented in the inner loop of a GEMC computer simulation.
16 Particle Displacement In this step, a particle is moved to some random new position within its original box. As such, the box is essentially a system at constant-nvt, and particles are moved using the normal Metropolis Monte Carlo algorithm. 9 The particle displacement is accepted or rejected based on the probability given by min(1,3 move ) where 3 move is given in Equation 1.4. Figure 1.4 Particle displacement within a box. exp( βenew) 3move = = exp( β E) 1.4 exp( βe ) old and E is the change in configurational energy for the move Volume Change In NPT-GEMC, to assure the constant pressure criterion is met, the volume of each box is varied independently, and is accepted or rejected with a probability given by min(1,3 vol ), where 3 vol is given in Equation vol = I I I II I V + V exp β E + E N kbt ln I V II II II V + V I II N k ln + ( + ) BT P V V II V 1.5 where the superscript I and II designate the properties of each individual box and k B is the Boltzmann constant.
17 10 Figure 1.5 shows a single box undergoing a volume change; note that the molecules remain in the same relative positions with respect to one another and the box boundaries. Figure 1.5 Before and after a volume change in a single box. Though it is possible to vary the volumes of each box simultaneously, convergence occurs more quickly when the volume of only one box is changed at a time, so the probability equation may be simplified to I I I I V + V I 3 vol = exp β E N kbt ln + P V I 1.6 V where the superscript I designates the box undergoing a volume perturbation. Again, E is the change in configurational energy. 8 For NVT-GEMC, the volumes of both boxes are adjusted at the same time so that the total volume is conserved; an increase in volume in one box is coupled with a decrease of the same amount in the other box. The volume perturbation in this step is either accepted or rejected based on the probability given by min(1,3 vol ), where 3 vol is given in Equation The superscripts I and II refer to the different boxes. 3 vol = N ( [ exp β E I B I V k T ln + E I II + V I V N II V k T ln B II + V II V 1.7
18 Particle Transfer In this perturbation, required to achieve coexistence between the two phases, a randomly selected particle is transferred from one phase box to the other, as shown in Figure 1.6. The picture is somewhat misleading, however, as the Figure 1.6 Particle transfer between phases. transfer is actually accomplished by creating a particle in one phase (for example, box I) and simultaneously destroying a particle in the other phase (box II). During this perturbation, the coexisting phases are representative of the grand canonical ensemble (constant-µv I T and constant-µv II T), where the two boxes are at the same temperature and each component has the same chemical potential in both boxes. The probability that the transfer will be accepted is given as 3 transfer, which may be calculated using Equations 1.8, 1.9, and 1.10, in which the particle is being transferred to box I from box II. I I I I 3 = exp[ β E + ln( zv /( N + 1))] 1.8 transfer II II II II 3 = exp[ β E + ln( N / zv )] 1.9 transfer II I I II I II V ( N + 1) 3 transfer = 3transfer 3transfer = exp β E + E + kbt ln I II 1.10 V N where the activity coefficient z is exp(βµ)/λ 3, λ is the de Broglie wavelength of constituent particles, and µ is the chemical potential. 8
19 12 REFERENCES Maitland, G. C.; Rigby, M.; Smith, E. B.; Wakeham, W. A. Intermolecular Forces: Their Origin and Determination; Clarendon: Oxford, Van Leeuwen, M. E. Fluid Phase Equil. 1994, 99, 1. Ewald, P. P. Ann. Phys. 1921, 64, 253. Allen, M. P.; Tildesley, D. J. Computer Simulation of Liquids; Clarendon: Oxford, 1987; pp Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications; Academic: Boston, Panagiotopoulos, A. Z., Molec. Phys., 1987, 61, 813. Panagiotopoulos, A. Z.; Quirke, N.; Stapleton, M.; Tildesley, D. J. Mol. Phys. 1988, 63, 527. Panagiotopoulos, A. Z.; Stapleton, M. R. Fluid Phase Equil. 1989, 53, 133. Metropolis, N.; Rosenbluth, A. W.; Rosenbluth, M. N.; Teller, A. N.; Teller, E. J. Chem. Phys. 1953, 21, 1087.
Phase 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 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 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 informationGibbs Ensemble Computer Simulations of Vapor Liquid Equilibrium of Hydrazine
CHME 498 Undergraduate Research Final Report New Mexico State University Gibbs Ensemble Computer Simulations of Vapor Liquid Equilibrium of Hydrazine Summer 2017 Gustav Barraza Faculty Adviser: Dr. Martha
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 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 informationVapour-liquid equilibria for two centre Lennard-Jones diatomics and dipolar diatomics
MOLECULAR PHYSICS, 1993, VOL. 80, No. 4, 997-1007 Vapour-liquid equilibria for two centre Lennard-Jones diatomics and dipolar diatomics By GIRIJA S. DUBEYt, SEAMUS F. O'SHEA Department of Chemistry, University
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 informationMultiple time step Monte Carlo
JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Multiple time step Monte Carlo Balázs Hetényi a) Department of Chemistry, Princeton University, Princeton, NJ 08544 and Department of Chemistry
More informationChemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model
Chemical Potential of Benzene Fluid from Monte Carlo Simulation with Anisotropic United Atom Model Mahfuzh Huda, 1 Siti Mariyah Ulfa, 1 Lukman Hakim 1 * 1 Department of Chemistry, Faculty of Mathematic
More informationMolecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid
Molecular dynamics simulation of the liquid vapor interface: The Lennard-Jones fluid Matthias Mecke and Jochen Winkelmann a) Institut für Physikalische Chemie, Universität Halle-Wittenberg, Geusaer Str.,
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 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 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 informationTHE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM
Molecular Simulation, 1987, Vol. 1, pp. 87-93 c Gordon and Breach Science Publishers S.A. THE DETAILED BALANCE ENERGY-SCALED DISPLACEMENT MONTE CARLO ALGORITHM M. MEZEI Department of Chemistry, Hunter
More informationIntroduction Statistical Thermodynamics. Monday, January 6, 14
Introduction Statistical Thermodynamics 1 Molecular Simulations Molecular dynamics: solve equations of motion Monte Carlo: importance sampling r 1 r 2 r n MD MC r 1 r 2 2 r n 2 3 3 4 4 Questions How can
More informationCE 530 Molecular Simulation
CE 530 Molecular Simulation Lecture 20 Phase Equilibria David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Thermodynamic Phase Equilibria Certain thermodynamic states
More 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 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 informationThermodynamics and structural properties of the dipolar Yukawa fluid. Citation Journal Of Chemical Physics, 1999, v. 111 n. 1, p.
Title Thermodynamics and structural properties of the dipolar Yukawa fluid Author(s) Szalai, I; Henderson, D; Boda, D; Chan, KY Citation Journal Of Chemical Physics, 1999, v. 111 n. 1, p. 337-344 Issued
More informationIntermolecular Potentials and The Second Virial Coefficient
Intermolecular Potentials and The Second Virial Coefficient Patrick L. Holt Department of Chemistry & Physics Bellarmine University Louisville, Kentucky 40205 pholt@bellarmine.edu Copyright 2004 by the
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 informationCopyright 2001 University of Cambridge. Not to be quoted or copied without permission.
Course MP3 Lecture 4 13/11/2006 Monte Carlo method I An introduction to the use of the Monte Carlo method in materials modelling Dr James Elliott 4.1 Why Monte Carlo? The name derives from the association
More informationEquation of State for Free Energy of Homogeneous Nucleation in Supersaturated Lennard-Jones Vapor Phase Derived by Monte Carlo Simulations
Equation of State for Free Energy of Homogeneous Nucleation in Supersaturated Lennard-Jones Vapor Phase Derived by Monte Carlo Simulations YAMADA Yuri March, 2003 Hosei University ABSTRACT Monte Carlo
More informationSIMCON - Computer Simulation of Condensed Matter
Coordinating unit: 230 - ETSETB - Barcelona School of Telecommunications Engineering Teaching unit: 748 - FIS - Department of Physics Academic year: Degree: 2017 BACHELOR'S DEGREE IN ENGINEERING PHYSICS
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 informationComputer simulation of vapor-liquid equilibria of linear dipolar fluids: Departures from the principle of corresponding states
Computer simulation of vapor-liquid equilibria of linear dipolar fluids: Departures from the principle of corresponding states Benito Garzón, Santiago Lago, a) and Carlos Vega Departamento de Química Física,
More informationEffect of a density-dependent potential on the phase behaviour of fluids
MOLECULAR PHYSICS, 1992, VOL. 77, No. 6, 1021-1031 Effect of a density-dependent potential on the phase behaviour of fluids By B. SMIT Koninklijke/Shell-Laboratorium, Amsterdam (Shell Research B.V.), PO
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 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 informationMONTE CARLO METHOD. Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002.
MONTE CARLO METHOD Reference1: Smit Frenkel, Understanding molecular simulation, second edition, Academic press, 2002. Reference 2: David P. Landau., Kurt Binder., A Guide to Monte Carlo Simulations in
More informationMIHALY MEZEI. Department of Physiology and Biophysics, Mount Sinai School of Medicine, CUNY. New York, NY 10029
Molecular Simulation, 1992, Vol. 9, pp, 257-261 c Gordon and Breach THEORETICAL CALCULATION OF THE LIQUID-VAPOR COEXISTENCE CURVE OF WATER, CHLOROFORM AND METHANOL WITH THECAVITY-BIASED MONTE CARLO METHOD
More informationABSTRACT. Sarvin Moghaddam, Doctor of Philosophy, 2003
ABSTRACT Title of Dissertation: FINELY-DISCRETIZED LATTICE MODELS FOR THERMODYNAMIC PROPERTIES OF FLUIDS Sarvin Moghaddam, Doctor of Philosophy, 2003 Dissertation directed by: Professor Mikhail Anisimov
More informationCalculation of the chemical potential in the Gibbs ensemble
MOLECULAR PHYSICS, 1989, VOL. 68, NO. 4, 951-958 Calculation of the chemical potential in the Gibbs ensemble by B. SMIT Koninklijke/Shell-Laboratorium, Amsterdam, (Shell Research B.V.) Badhuisweg 3, 1031
More informationMonte Carlo simulations of vapor-liquid-liquid equilibrium of ternary
1 Monte Carlo simulations of vapor-- equilibrium of ternary petrochemical mixtures Suren Moodley a, Kim Bolton b, and Deresh Ramjugernath a a Thermodynamics Research Unit, School of Chemical Engineering,
More informationIntroduction to Simulation - Lectures 17, 18. Molecular Dynamics. Nicolas Hadjiconstantinou
Introduction to Simulation - Lectures 17, 18 Molecular Dynamics Nicolas Hadjiconstantinou Molecular Dynamics Molecular dynamics is a technique for computing the equilibrium and non-equilibrium properties
More informationREVIEW: Derivation of the Mean Field Result
Lecture 18: Mean Field and Metropolis Ising Model Solutions 1 REVIEW: Derivation of the Mean Field Result The Critical Temperature Dependence The Mean Field Approximation assumes that there will be an
More informationTheoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes
J. At. Mol. Sci. doi: 10.4208/jams.121011.011412a Vol. 3, No. 4, pp. 367-374 November 2012 Theoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes Xiu-Ying Liu a,, Li-Ying Zhang
More informationWang-Landau Monte Carlo simulation. Aleš Vítek IT4I, VP3
Wang-Landau Monte Carlo simulation Aleš Vítek IT4I, VP3 PART 1 Classical Monte Carlo Methods Outline 1. Statistical thermodynamics, ensembles 2. Numerical evaluation of integrals, crude Monte Carlo (MC)
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 informationMonte Carlo Methods in Statistical Mechanics
Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics
More 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 informationMonte Carlo Molecular Simulation with Isobaric-Isothermal and Gibbs-NPT Ensembles. Thesis by Shouhong Du
Monte Carlo Molecular Simulation with Isobaric-Isothermal and Gibbs-NPT Ensembles Thesis by Shouhong Du In Partial Fulfillment of the Requirements For the Degree of Master of Science King Abdullah University
More 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 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 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 informationSchool of Chemical Engineering, Universidad del Valle, Ciudad Universitaria Melendez, Building 336, Apartado 25360, Cali, Colombia
Vapor-phase chemical equilibrium and combined chemical and vapor-liquid equilibrium for the ternary system ethylene + water + ethanol from reaction-ensemble and reactive Gibbs-ensemble molecular simulations
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 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 informationEpisode 1: Phis wants to be a physicist
Illustration: Xia Hong Script: Xia Hong 12/2012 Phis is 14 years old. Like most kids of her age, she s a bit anxious about finding the person inside her She tried singing. Hmm You are very talented, Phis!
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 informationPORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES
PORE SIZE DISTRIBUTION OF CARBON WITH DIFFERENT PROBE MOLECULES Atichat Wongkoblap*, Worapot Intomya, Warangkhana Somrup, Sorod Charoensuk, Supunnee Junpirom and Chaiyot Tangsathitkulchai School of Chemical
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 informationKINETIC MONTE CARLO SIMULATION OF MIXTURES: PHASE EQUILIBRIA
KINETIC MONTE CARLO SIMULATION OF MIXTURES: PHASE EQUILIBRIA SHILIANG JOHNATHAN TAN BEng (Chem) A thesis submitted for the degree of Master of Philosophy at The University of Queensland in 2015 School
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 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 informationMolecular Simulation towards Efficient and Representative Subsurface Reservoirs Modeling. Dissertation by Ahmad Kadoura
Molecular Simulation towards Efficient and Representative Subsurface Reservoirs Modeling Dissertation by Ahmad Kadoura In Partial Fulfillment of the Requirements For the Degree of Doctor of Philosophy
More informationRemoving the mystery of entropy and thermodynamics. Part 3
Removing the mystery of entropy and thermodynamics. Part 3 arvey S. Leff a,b Physics Department Reed College, Portland, Oregon USA August 3, 20 Introduction In Part 3 of this five-part article, [, 2] simple
More informationModeling & Simulation of Glass Structure
Modeling & Simulation of Glass Structure VCG Lecture 19 John Kieffer Department of Materials Science and Engineering University of Michigan 1 Overview Historical perspective Simulation methodologies Theoretical
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 informationSTATISTICAL AND THERMAL PHYSICS
Phys 362 Spring 2015 STATISTICAL AND THERMAL PHYSICS Phys 362 Spring 2015 Instructor: Office: Professor David Collins WS 228B Phone: 248-1787 email: Office Hours: dacollin@coloradomesa.edu M 9 10 am, 11
More informationThe Intersection of Chemistry and Biology: An Interview with Professor W. E. Moerner
The Intersection of Chemistry and Biology: An Interview with Professor W. E. Moerner Joseph Nicolls Stanford University Professor W.E Moerner earned two B.S. degrees, in Physics and Electrical Engineering,
More informationMultiple time step Monte Carlo simulations: Application to charged systems with Ewald summation
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 1 1 JULY 004 Multiple time step Monte Carlo simulations: Application to charged systems with Ewald summation Katarzyna Bernacki a) Department of Chemistry
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 informationA New Intermolecular Potential Model for the n- Alkane Homologous Series
A New Intermolecular Potential Model for the n- Alkane Homologous Series Jeffrey R. Errington and Athanassios Z. Panagiotopoulos School of Chemical Engineering Cornell University Ithaca, NY 14853 and Institute
More informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
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 informationMultiscale Materials Modeling
Multiscale Materials Modeling Lecture 02 Capabilities of Classical Molecular Simulation These notes created by David Keffer, University of Tennessee, Knoxville, 2009. Outline Capabilities of Classical
More informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 4 I. THE ISOTHERMAL-ISOBARIC ENSEMBLE The isothermal-isobaric ensemble is the closest mimic to the conditions under which most experiments are
More informationIntermolecular Model Potentials and Virial Coefficients from Acoustic Data
JASEM ISSN 1119-8362 All rights reserved Full-text Available Online at https://www.ajol.info/index.php/jasem http://www.bioline.org.br/ja J. Appl. Sci. Environ. Manage. Vol.22 (2) 246-251. February 2018
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 informationICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below
ICCP Project 2 - Advanced Monte Carlo Methods Choose one of the three options below Introduction In statistical physics Monte Carlo methods are considered to have started in the Manhattan project (1940
More informationHow Different is the Dynamics of a Lennard-Jones Binary Fluid from One-Component Lennard-Jones Fluid? 1
How Different is the Dynamics of a Lennard-Jones Binary Fluid from One-Component Lennard-Jones Fluid? Takayuki NARUMI and Michio TOKUYAMA Summary We investigate the dynamics of liquids and supercooled
More informationChemistry Physical Chemistry I Fall 2018
Chemistry 309 - Physical Chemistry I Fall 2018 Instructor: Office Hours: Dr. Samuel A. Abrash C-208 Gottwald Science Center Work: 289-8248 Home: 323-7363 Cell: 363-2597 sabrash@richmond.edu www.richmond.edu/~sabrash
More informationTitle Theory of solutions in the energy r of the molecular flexibility Author(s) Matubayasi, N; Nakahara, M Citation JOURNAL OF CHEMICAL PHYSICS (2003), 9702 Issue Date 2003-11-08 URL http://hdl.handle.net/2433/50354
More informationarxiv: v1 [physics.comp-ph] 22 Apr 2009
Comprehensive study of the vapour-liquid equilibria of the pure two-centre Lennard-Jones plus pointdipole fluid Jürgen Stoll, Jadran Vrabec, Hans Hasse arxiv:0904.3413v1 [physics.comp-ph] 22 Apr 2009 Institut
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 informationMonte Carlo Study of the Crystalline and Amorphous NaK Alloy
Available online at www.sciencedirect.com ScienceDirect Procedia Computer Science 08C (207) 25 22 International Conference on Computational Science, ICCS 207, 2-4 June 207, Zurich, Switzerland Monte Carlo
More informationFaculty: Andrew Carr, Ryan Felix, Stephanie Gould, James Hebda, Karla McCain, John Richardson, Lindsay Zack
CHEMISTRY Chair: Bradley Smucker (Fall 2017) and Andrew Carr (Spring 2018) Faculty: Andrew Carr, Ryan Felix, Stephanie Gould, James Hebda, Karla McCain, John Richardson, Lindsay Zack Adjunct Instructor:
More informationCHEM-UP! D A Y 3 S P R I N G The Academic Support Daytona State College (Chem-Up 3, Page 1 of 36)
CHEM-UP! D A Y 3 S P R I N G 2 0 1 5 The Academic Support Center @ Daytona State College (Chem-Up 3, Page 1 of 36) Welcome to Chem Up! An Introduction to Basic Chemistry Concepts Please sign in. The Academic
More informationComparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures
Fluid Phase Equilibria 178 (2001) 87 95 Comparison of different mixing rules for prediction of density and residual internal energy of binary and ternary Lennard Jones mixtures Jian Chen a,, Jian-Guo Mi
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 informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More informationMD simulation of methane in nanochannels
MD simulation of methane in nanochannels COCIM, Arica, Chile M. Horsch, M. Heitzig, and J. Vrabec University of Stuttgart November 6, 2008 Scope and structure Molecular model for graphite and the fluid-wall
More informationAn Introduction to Two Phase Molecular Dynamics Simulation
An Introduction to Two Phase Molecular Dynamics Simulation David Keffer Department of Materials Science & Engineering University of Tennessee, Knoxville date begun: April 19, 2016 date last updated: April
More informationMonte Carlo methods for phase equilibria of uids 2 1. Introduction The phase behavior of uids and their mixtures is of central importance to many tech
TOPICAL REVIEW Monte Carlo methods for phase equilibria of uids Athanassios Z. Panagiotopoulos z Institute for Physical Science and Technology and Department of Chemical Engineering, University of Maryland,
More informationChemistry Physical Chemistry I Fall 2017
Chemistry 309 - Physical Chemistry I Fall 2017 Instructor: Office Hours: Dr. Samuel A. Abrash C208 Gottwald Science Center Work: 289-8248 Home: 323-7363 Cell: 363-2597 sabrash@richmond.edu www.richmond.edu/~sabrash
More informationProcedure for Setting Goals for an Introductory Physics Class
Procedure for Setting Goals for an Introductory Physics Class Pat Heller, Ken Heller, Vince Kuo University of Minnesota Important Contributions from Tom Foster, Francis Lawrenz Details at http://groups.physics.umn.edu/physed
More informationCHEMISTRY, B.S. (SCIENCE)
Chemistry, B.S. (Science) 1 CHEMISTRY, B.S. (SCIENCE) Begin Campus: Any Penn State Campus End Campus: University Park Program Description This major provides a strong foundation in the theory and practice
More informationChem 4501 Introduction to Thermodynamics, 3 Credits Kinetics, and Statistical Mechanics
Chem 4501 Introduction to hermodynamics, 3 Credits Kinetics, and Statistical Mechanics Module Number 2 Active Learning Answers and Optional Problems/Solutions 1. McQuarrie and Simon, 2-6. Paraphrase: How
More informationComputer simulations as concrete models for student reasoning
Computer simulations as concrete models for student reasoning Jan Tobochnik Department of Physics Kalamazoo College Kalamazoo MI 49006 In many thermal physics courses, students become preoccupied with
More informationCHE 371: Kinetics and Thermodynamics Fall 2008
CHE 371: Kinetics and Thermodynamics Fall 2008 Class Meetings: Lecture: M, T, W, F 9:00 AM, Olin 103 Laboratory: T, W 1:30-5:20 PM, R 2:30-6:20 PM Instructor: Prof. Amanda Nienow, Nobel 106C, 933-7327,
More informationINTERMOLECULAR FORCES
INTERMOLECULAR FORCES Their Origin and Determination By GEOFFREY C. MAITLAND Senior Research Scientist Schlumberger Cambridge Research, Cambridge MAURICE RIGBY Lecturer in the Department of Chemistry King's
More informationIntroduction to Thermodynamic States Gases
Chapter 1 Introduction to Thermodynamic States Gases We begin our study in thermodynamics with a survey of the properties of gases. Gases are one of the first things students study in general chemistry.
More informationWelcome to Chem Up! Agenda for Day 3. How Does Learning Work? I went to class, I ve read my textbook, now what? 8/14/2014
Welcome to Chem Up! An Introduction to Basic Chemistry Concepts Please sign in. Agenda for Day 3 How Does Learning Work? Memorization Techniques Compounds and Mixtures Bonding Chemical Reactions Hydrogen
More informationPHY 6500 Thermal and Statistical Physics - Fall 2017
PHY 6500 Thermal and Statistical Physics - Fall 2017 Time: M, F 12:30 PM 2:10 PM. From 08/30/17 to 12/19/17 Place: Room 185 Physics Research Building Lecturer: Boris Nadgorny E-mail: nadgorny@physics.wayne.edu
More informationGreen Chemistry Commitment
Info Kit Green Chemistry Commitment What is the Green Chemistry Commitment? The Green Chemistry Commitment (GCC) is a consortium program that unites the green chemistry community around shared goals and
More informationThermodynamic and transport properties of carbon dioxide from molecular simulation
Thermodynamic and transport properties of carbon dioxide from molecular simulation Carlos Nieto-Draghi, Theodorus de Bruin, Javier Pérez-Pellitero, Josep Bonet Avalos, and Allan D. Mackie Citation: The
More informationAdsorption of Lennard-Jones Fluids in Carbon Slit Pores of a Finite Length. AComputer Simulation Study
1 Invited Contribution Adsorption of Lennard-Jones Fluids in Carbon Slit Pores of a Finite Length. AComputer Simulation Study A. Wongkoblap 1, S. Junpirom 2 and D.D. Do 1 * (1) Department of Chemical Engineering,
More informationThe Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity
Draft: September 29, 1999 The Effect of Model Internal Flexibility Upon NEMD Simulations of Viscosity N. G. Fuller 1 and R. L. Rowley 1,2 Abstract The influence of model flexibility upon simulated viscosity
More information