Standard thermodynamic properties of solutes in supercritical solvents: simulation and theory
|
|
- Briana Morrison
- 6 years ago
- Views:
Transcription
1 Chemical Physics Letters 381 (2003) Standard thermodynamic properties of solutes in supercritical solvents: simulation and theory J.L. Alvarez a,c, R. Fernandez-Prini a,b, *, E. Marceca b a Unidad Actividad Quımica, Comision Nacional Energıa Atomica, Av. Libertador 8250, uenos Aires, Argentina b INQUIMAE, Facultad Ciencias Exactas y Naturales, Universidad de uenos Aires, Ciudad Universitaria, Pabellon II, uenos Aires, Argentina c Dto. Ingenierıa Quımica, FRA, Universidad Tecnologica Nacional, Medrano 951, uenos Aires, Argentina Received 11 September 2003; in final form 15 October 2003 Published online: 4 November 2003 Abstract We report a modified simulation procedure to calculate partial molar quantities of dilute solutions in supercritical fluids through the Krichevskii function. This procedure circumvents the difficulties posed by the very large solvent compressibility and expansivity in the near-critical region of the solutions where other simulation procedures cannot be used. Ó 2003 Elsevier.V. All rights reserved. 1. Introduction In 1987 Ciccotti, Frenkel and coworkers [1] proposed a procedure to calculate partial molar quantities of binary fluid mixtures by simulation using a modified Widom insertion method [2]. eing partial molar quantities derivatives of extensive properties, the procedure used up to then to calculate their values required two simulation experiments carried out in the NpT Gibbs ensemble for two different compositions of the binary mixture. That procedure duplicated the error and the number of runs, because of the need to use two * Corresponding author. Fax: address: rfprini@cnea.gov.ar (R. Fernandez-Prini). simulation experiments. The procedure proposed by Ciccotti, Frenkel and coworkers involves the (virtual) exchange of a particle of component A by one molecule of compound in the Gibbs ensemble yielding the difference between the partial molar quantities of the two components by means of a single computational experiment. The equations for the partial molar volumes and enthalpies were derived [3] and these expressions were used to calculate the partial molar properties of Lennard Jones Ar Kr mixtures with the exchange procedure (EXP). The calculations were successful, but the equations derived for the EXP have an intrinsic limitation, the authors considered that the volume and the enthalpy of the mixtures were independent of pressure and of temperature, respectively. Since /$ - see front matter Ó 2003 Elsevier.V. All rights reserved. doi: /j.cplett
2 772 J.L. Alvarez et al. / Chemical Physics Letters 381 (2003) there is a special interest in the determination of standard (infinite dilution) partial molar quantities which are the corner-stone properties to describe the behaviour of solutes, we shall illustrate the limitations of the EXP equations with the simpler expressions for the partial molar volume of a solute at infinite dilution. The difference between the partial molar volumes V A and V of the two components, is given by [1,3] ex V expð bdua! V A V ¼ Þ N A ;N expð bdua! ex Þ hvi NA ;N ; NA ;N ð1þ where DUA! ex is the change in the energy of the system of N A molecules A and N molecules when one A is (virtually) exchanged for one molecule; h i denotes an average taken over the ensemble of ðn A þ N Þ particles. If the change of volume with pressure is taken into account, ðv A V Þ becomes V A V Dh i E V þ p ov expð bdu ex op A! Þ T N ¼ A ;N expð bdua! ex Þ N A ;N V þ p ov : ð2þ op T N A ;N When the dilute solutions are close to the solventõs critical point, ðov =opþ T is very large becoming infinity at the critical point for solutions at infinite dilution. Taking into account the difficulties of the experimental determination of the properties of solutes dissolved in supercritical fluids, which have become of great importance in chemical processes and in studying the physical chemistry of solvation phenomena, it would be very valuable to use simulation methods to determine, or at least estimate, the values of those quantities. The method devised by Ciccotti, Frenkel and coworkers will intrinsically fail for supercritical solutions which are characterized by large compressibilities. For attractive solute A (that interacts more strongly with the solvent than two solvent molecules do) VA 1 is large and negative for wide ranges of temperature and density. A similar situation occurs for partial molar enthalpies at infinite dilution since the heat capacities also become very large in the near-critical region, so the dependence of the enthalpy on temperature cannot be neglected. On the other hand, it has been established that partial molar quantities are not convenient to base the description of the near-critical behaviour because these quantities diverge at the solventõs critical point [4]. Using the volume as the archetypical property, the soluteõs partial molar VA 1 is given by " # 1 V 1 A ¼ V 1 op þ j T ¼ V ox ð1 þ j T JÞ; ð3þ A V ;T where x A is the mole fraction of A, J is called the Krichevskii function, a well behaved quantity at the solventõs critical point and j T is the isothermal compressibility of solvent, a strongly diverging quantity at the critical point of. This is the reason why the favoured description of the thermodynamic properties of near-critical and supercritical solutes is based upon J [5,6]. From the definition of J it is natural to use a modified EXP for the calculation of J by simulation experiments in the NVT ensemble. The main goal of this preliminary communication is to test the feasibility of this approach and evaluate its performance for simple attractive solutes and different reduced temperature T red ¼ ðt =T C Þ, where T C is the solventõs critical temperature. Consequently we have restricted this study to binary Lennard Jones systems, because it is possible to compare the results of the simulation procedure hereby proposed for J with an equation of state that fits many simulation data for Lennard Jones fluids [7] and also with the results of an inhomogeneous integral equation [8], the latter equation has recently been proved successful to describe the behaviour of near-critical systems [8,9]. 2. asis of the calculation procedure The application of the EXP to a function f ðr i Þ in the NVT ensemble leads to the following general expression:
3 J.L. Alvarez et al. / Chemical Physics Letters 381 (2003) hf ðr i Þi NA þ1;n 1 ¼ f ðr iþ expð bdua! ex Table 1 N A ;N A lres ÞŠ : Lennard Jones parameters used r (nm) (e=k) (K) Ref. ð4þ Ne [14] The average value of the function for the system Kr [15] consisting of N A þ 1; N 1 particles f ðr i Þ in the NVT ensemble, may be calculated by averaging the function f ðr i Þ expð bdua! ex Þ for the ensemble of N A ; N particles and dividing that quantity by Xe [15] A lres Š, where lres I is the residual chemical potential of component I. Eq. (4) was derived by a straightforward application of the published equations [1,3]. In our case f ðr i Þ is the pressure of the system when one A (solute) is added replacing a (solvent) molecule. For our purpose we need to calculate the difference in pressure p of the system having ðn A þ 1; N 1Þ particles p NA þ1;n 1 N A þ1;n and p 1 NA ;N N A ;N. According to Eq. (4) this difference is obtained averaging the two quantities in the N A ; N system. Using Eq. (4) we get J ¼ p N A þ1;n 1 expð bdua! ex Þ, A lres ÞŠ p NA ;N x A : N A ;N ð5þ The pressure was calculated in each Monte Carlo simulation step with the virial equation [10]. In every step a solvent molecule was chosen at random and exchanged by a solute particle and DU ex A! and the pressure were calculated for the virtual N Aþ1 ; N 1 system.the quantity inside the angular brackets of Eq. (5) was calculated in every step of the simulation and then averaged. With the same program the value of the difference of residual chemical potentials was calculated with the expression [1,3] A lres ÞŠ ¼ expð bdu ex A! Þ N A ;N : ð6þ The systems studied in this work consisted of 125 Lennard Jones solvent molecules having the intermolecular parameters of neon. One of the solvent molecules was exchanged in every Monte Carlo step by a Lennard Jones Kr or Xe atom. Runs were made for fluid densities within the range of interest for applications of supercritical solvents and at two values of T red > 1 for each ÔsoluteÕ atom. The size of the simulation box was adjusted in every run to accommodate the 125 particles in the box at the chosen density, and the interactions were truncated at a distance equal to half of the box length. The number of steps of each run was between 0.2 and depending on T red and strength of solute solvent interaction. The contribution to J of particles beyond the cutoff distance was calculated using for the tail contribution to pressure the equation given by Frenkel and Smit [11]. Applying this expression for the exchange of one molecule by that of the A solute molecule, we obtain J tail Dptail A! x A ¼ 32pq 2 e r 3 " e A e r A r 3 T A T #; ð7þ where T IJ gives the contribution of the integral of the distance dependent function between the cutoff reduced distance z cut IJ ¼ r=r IJ and infinity T IJ ¼ 1 2 : ð8þ 3z 3 IJ 9z 9 IJ The parameters of the Lennard Jones potential for the three atoms are given in Table 1. The parameters for the interaction of solvent molecules with the solute where obtained from those of the pure components using the simple Lorentz erthelot combining rule. 3. Equation of state and integral equation used for comparison There are two other tools that we have used to check the results obtained with the EXP simulation experiments. The first one is the JZG equation
4 774 J.L. Alvarez et al. / Chemical Physics Letters 381 (2003) Table 2 Values of J (MPa) calculated with the Monte Carlo exchange procedure Solute T red qr Xe 1.30 )45.1, )51.5 )84.1 )128.2, )126.2 Xe 1.10 )44.6 )46.5 )103.3 Kr 1.10 )26.7, )32.7 )71.6 Kr 1.03 )30.5 )72.2, )80.7 qr Xe 1.30 )172.6 Xe )187.6 )186.6 Kr 1.10 )128.3 )138.4 Kr 1.03 )86.4, )92.6 )109.4 )146.7 of state developed by Johnson et al. [7] from results of many MD and MC simulations for Lennard Jones pure fluids. The equation of state fitted to the simulation results was based upon a modified enedict Webb Rubin formulation, thus it is able to describe also the critical region of the fluid. Johnson et al. [7] remark that their equation of state may be also applied to conformal binary mixtures, hence it should be applicable to our Lennard Jones atomic binary mixtures. For this purpose we used the JZG equation employing the van der Waals one-fluid theory, as recommended by the authors, the pressure was then differentiated with respect to composition (cf. Eq. (3)) and finally the limit to infinite dilution was taken. The second tool used in the present work to compare with simulation results was the integral equation for inhomogeneous fluids, known as hydrostatic hypernetted chain (HHNC) equation [12]. We have shown recently [13] that HHNC can be successfully employed to account for many observations of equilibrium and structural properties of near-critical solutions. In this case J was calculated using the solvent solvent and the solvent solute direct correlation functions according to the equation bj ¼ q 2 ½^c ð0þ ^c A ð0þš; ð9þ where q is the density of the solvent and ^c IJ ð0þis the Fourier transform of the direct correlation function for the wave vector equal to zero. The fact noted previously [13] that to get a more reliable result for J with HHNC it was convenient to truncate the integral ^c A ð0þ after three molecular diameters from the solute, introduced some uncertainty about the correctness of the values Fig. 1. Values of J for Xe dissolved in Ne. (}) MC results; solid curve, JZG equation of state; dashed curve, HHNC equation. (a) T red ¼ 1.30, (b) T red ¼ 1.10.
5 J.L. Alvarez et al. / Chemical Physics Letters 381 (2003) obtained for J. The results of the present work show conclusively that the truncation procedure does not introduce significant changes. 4. Results and discussion The values of J obtained for the two binary mixtures are reported on Table 2, the thermodynamic state points studied correspond to densities typical of those used in studies involving supercritical solvents. The reproducibility of the calculated values is 10% in the worst case and more frequenty close to 5%. It should be stressed that the reproducibility decreases substantially when T red approaches unity and also when the solute solvent interaction is larger [8,16], this is why we did not calculate J for Xe in Ne at T red < 1:10 with the EXP. Hence, it would be misleading to give an overall value for the error, we preferred to make duplicate runs for five state points, shown in Table 2, which are representative of data reproducibility. For Kr in Ne even at T red ¼ 1:03 and reduced density 0.25, which is very close to the critical density of the solvent, the value of J is very reasonable. Figs. 1 and 2 are plots of the calculated J with EXP and with the JZG equation [7] and with HHNC; the results calculated with the Monte Carlo EXP are in very statisfactory agreement with the others. The two other procedures yield values of J that are very close indeed with the exception of the high density region (cf. Figs. 1b and 2a,b), since that density region is outside of our main interest to cover conditions typical of supercritical solutions, we shall not discuss the possible causes for this discrepancy. 5. Conclusion Our work shows that the EXP is a valuable tool to calculate J and that this route also enables the calculation of VA 1 using the value of the Krichevskii function, the compressibility of the pure solvent and Eq. (3). The precision of the results of EXP is limited by the closeness to the solventõs critical point and by the strength of the solute solvent interaction. However the use of a larger simulation box, that will imply the use of more particles, will allow approaching closer to the solventõs critical point, and/or enable to calculate properties for systems with larger solute solvent interactions [13]. This will be our next goal which will be used to calculate values of the Kritchevskii function at the solventõs critical point for dilute aqueous solutions. Acknowledgements Fig. 2. Values of J for Kr dissolved in Ne. (}) Monte Carlo results; solid curve JZG equation of state; dashed curve HHNC equation. (a) T red ¼ 1.10, (b) T red ¼ The authors are grateful to ANPCyT for partial economic support and to Dr. D.H. Larıa for making available to us the basic Monte Carlo program.
6 776 J.L. Alvarez et al. / Chemical Physics Letters 381 (2003) References [1] P. Sindzingre, G. Ciccotti, C. Massobrio, D. Frenkel, Chem. Phys. Lett. 136 (1987) 35. [2]. Widom, J. Chem. Phys. 39 (1963) 2008;. Widom, J. Phys. Chem. 86 (1982) [3] P. Sindzingre, C. Massobrio, G. Ciccotti, D. Frenkel, Chem. Phys. 129 (1989) 213. [4] R. Fernandez-Prini, M.L. Japas, Rev. Chem. Soc. 23 (1994) 155. [5] J. Alvarez, R. Fernandez-Prini, M.L. Japas, Ind. Eng. Chem. Res. 39 (2000) [6] J.P. OÕConnell, A.V. Sharygin, R.H. Wood, Ind. Eng. Chem. Res. 35 (1996) [7] J.K. Johnson, J.A. Zollweg, K.E. Gubbins, Mol. Phys. 78 (1993) 591. [8] R. Fernandez-Pini, J. Phys. Chem. 106 (2002) [9] L. ronstein, D.P. Fernandez, R. Fernandez-Prini, J. Chem. Phys. 117 (2002) 220. [10] J-P. Hansen, I.R. McDonald, Theory of Simple Liquids, Academic Press, New York, [11] D. Frenkel,. Smit, Understanding Molecular Simulation, Academic Press, New York, [12] Y. Zhou, G. Stell, J. Chem. Phys. 92 (1990) [13] G. Sciaini, E. Marceca, R. Fernandez-Prini, Phys. Chem. Chem. Phys. 4 (2002) [14] S.A. Egorov, A. Yethiraj, J.L. Skinner, Chem. Phys. Lett. 317 (2000) 558. [15] L.S. Tee, S. Gotoh, W.E. Stewart, Ind. Eng. Chem. Fundam. 5 (1966) 356. [16] S.A. Egorov, J. Chem. Phys. 112 (2000) 7138.
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces
Supplemental Material for Temperature-sensitive colloidal phase behavior induced by critical Casimir forces Minh Triet Dang, 1 Ana Vila Verde, 2 Van Duc Nguyen, 1 Peter G. Bolhuis, 3 and Peter Schall 1
More 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 informationHenry s Law Constants of Methane and Acid Gases at Pressures above the Saturation Line of Water
Henry s Law Constants of Methane and Acid Gases at Pressures above the Saturation Line of Water Josef Sedlbauer and Vladimir Majer 2* Department of Chemistry, Technical University of Liberec, 46 7 Liberec,
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 informationAn explicit expression for finite-size corrections to the chemical potential
J. Phys.: Condens. Matter 1 (1989) 8659-8665. Printed in the UK An explicit expression for finite-size corrections to the chemical potential B Smitt and D Frenkelt t Koninklijke/Shell-Laboratorium, Amsterdam
More informationTheory of infinite dilution chemical potential
Fluid Phase Equilibria Journal Volume 85, 141-151, 1993 141 Esam Z. Hamada and G.Ali Mansoorib, * a Department of Chemical Engineering, King Fahd University of Petroleum and Minerals, Dhahran 31261 (Saudi
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 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 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 informationAn accurate expression for radial distribution function of the Lennard-Jones fluid
CHEMICAL PHYSICS Volume 3, Pages -5, 5 DOI:.6/j.chemphys.4.9.7 An accurate expression for radial distribution function of the Lennard-Jones fluid Ali Morsali a, *, Elaheh K. Goharshadi a, G.Ali Mansoori
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 informationHydrophobicity in Lennard-Jones solutions
PAPER www.rsc.org/pccp Physical Chemistry Chemical Physics Hydrophobicity in Lennard-Jones solutions Mario Ishizai, Hidei Tanaa and Kenichiro Koga* Received 9th September 2010, Accepted 12th October 2010
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 informationSome comments on the double retrograde vaporization
J. Chem. Thermodynamics 35 (2003) 583 589 www.elsevier.com/locate/jct Some comments on the double retrograde vaporization Ulrich K. Deiters * Institute of Physical Chemistry, University of Cologne, Luxemburger
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 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 informationEQUATION OF STATE DEVELOPMENT
EQUATION OF STATE DEVELOPMENT I. Nieuwoudt* & M du Rand Institute for Thermal Separation Technology, Department of Chemical Engineering, University of Stellenbosch, Private bag X1, Matieland, 760, South
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 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 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 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 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 informationCritical temperature of in nitely long chains from Wertheim s perturbation theory
MOLECULAR PHYSICS, 2000, VOL. 98, NO. 17, 1295 ± 1308 Critical temperature of in nitely long chains from Wertheim s perturbation theory C. VEGA* and L. G. MACDOWELL Departamento de QuõÂ mica Fõ Â sica,
More informationSITARAM K. CHAVAN * and MADHURI N. HEMADE ABSTRACT INTRODUCTION
Int. J. Chem. Sci.: 11(1), 013, 619-67 ISSN 097-768X www.sadgurupublications.com DENSITIES, VISCOSITIES AND EXCESS THERMODYNAMIC PROPERTIES OF MONOMETHYL AMMONIUM CHLORIDE IN TETRAHYDROFURAN AND WATER
More informationCalculating thermodynamic properties from perturbation theory I. An analytic representation of square-well potential hard-sphere perturbation theory
Ž. Fluid Phase Equilibria 154 1999 1 1 Calculating thermodynamic properties from perturbation theory I. An analytic representation of square-well potential hard-sphere perturbation theory Bing-Jian Zhang
More informationChemistry. Lecture 10 Maxwell Relations. NC State University
Chemistry Lecture 10 Maxwell Relations NC State University Thermodynamic state functions expressed in differential form We have seen that the internal energy is conserved and depends on mechanical (dw)
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 informationCritical Properties of Isobaric Processes of Lennard-Jones Gases
Critical Properties of Isobaric Processes of Lennard-Jones Gases Akira Matsumoto Department of Material Sciences, College of Integrated Arts Sciences, Osaka Prefecture University, Sakai, Osaka, 599-8531,
More informationNumerical Aspects of the SAFT Equation of State
University of Rhode Island DigitalCommons@URI Senior Honors Projects Honors Program at the University of Rhode Island 006 Numerical Aspects of the SAFT Equation of State Leah M. Octavio University of Rhode
More informationKinetics of thermal decoloration of a photomerocyanine in mixtures of protic and nonpolar solvents
Journal of Photochemistry and Photobiology A: Chemistry 153 (2002) 25 31 Kinetics of thermal decoloration of a photomerocyanine in mixtures of protic and nonpolar solvents Germán Sciaini a, Diana E. Wetzler
More informationV.E Mean Field Theory of Condensation
V.E Mean Field heory of Condensation In principle, all properties of the interacting system, including phase separation, are contained within the thermodynamic potentials that can be obtained by evaluating
More informationEquation of state of additive hard-disk fluid mixtures: A critical analysis of two recent proposals
PHYSICAL REVIEW E 66, 0310 00 Equation of state of additive hard-disk fluid mixtures: A critical analysis of two recent proposals M. López de Haro* Centro de Investigación en Energía, UNAM, Temixco, Morelos
More informationVariational Approach to the Equilibrium Thermodynamic Properties of Simple Liquids. I*
THE JOURNAL OF CHEMICAL PHYSICS VOLUME 51, NUMBER 11 1 DECEMBER 1969 Variational Approach to the Equilibrium Thermodynamic Properties of Simple Liquids. I* G.Ali MANSOORi (1) AND Frank B. CANFIELD (2)
More informationTheory of Interfacial Tension of Partially Miscible Liquids
Theory of Interfacial Tension of Partially Miscible Liquids M.-E. BOUDH-HIR and G.A. MANSOORI * University of Illinois at Chicago (M/C 063) Chicago, Illinois USA 60607-7052 Abstract The aim of this work
More informationPartial molar volumes at infinite dilution in aqueous solutions of NaCl, LiCl, NaBr, and CsBr at temperatures from 550 K to 725 K
J. Chem. Thermodynamics 1998, 3, 312 Partial molar volumes at infinite dilution in aqueous solutions of NaCl, LiCl, NaBr, and CsBr at temperatures from 55 K to 725 K Josef Sedlbauer, Department of Chemistry,
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 informationVolumetric Study of the Binary Mixtures Containing a Branched Hexane and an Isomeric Chlorobutane
78 Journal of Applied Solution Chemistry and Modeling, 205, 4, 78-84 Volumetric Study of the Binary Mixtures Containing a Branched Hexane and an Isomeric Chlorobutane Hernando Guerrero, Félix M. Royo and
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 informationChapter 3 PROPERTIES OF PURE SUBSTANCES
Thermodynamics: An Engineering Approach Seventh Edition Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2011 Chapter 3 PROPERTIES OF PURE SUBSTANCES Copyright The McGraw-Hill Companies, Inc. Permission
More informationInfluence of thermodynamic variables and molecular parameters in the solubility of gases in liquids
Pure & Appl. Chem., Vol. 6, No. 11, pp. 083-090,1990. Printed in Great Britain. @ 1990 IUPAC Influence of thermodynamic variables and molecular parameters in the solubility of gases in liquids Roberto
More informationOn the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models
JOURNAL OF CHEMICAL PHYSICS VOLUME 109, NUMBER 12 22 SEPTEMBER 1998 On the local and nonlocal components of solvation thermodynamics and their relation to solvation shell models Nobuyuki Matubayasi Institute
More informationDETERMINATION OF THE POTENTIAL ENERGY SURFACES OF REFRIGERANT MIXTURES AND THEIR GAS TRANSPORT COEFFICIENTS
THERMAL SCIENCE: Year 07, Vo., No. 6B, pp. 85-858 85 DETERMINATION OF THE POTENTIAL ENERGY SURFACES OF REFRIGERANT MIXTURES AND THEIR GAS TRANSPORT COEFFICIENTS Introduction by Bo SONG, Xiaopo WANG *,
More informationExam 1 Solutions 100 points
Chemistry 360 Fall 018 Dr. Jean M. Standard September 19, 018 Name KEY Exam 1 Solutions 100 points 1.) (14 points) A chunk of gold metal weighing 100.0 g at 800 K is dropped into 100.0 g of liquid water
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 informationTHERMODYNAMIC CONSISTENCY TESTS FOR PHASE EQUILIBRIUM IN LIQUID SOLUTE+SUPERCRITICAL SOLVENT MIXTURES
THERMODYNAMIC CONSISTENCY TESTS FOR PHASE EQUILIBRIUM IN LIQUID SOLUTE+SUPERCRITICAL SOLVENT MIXTURES José O. Valderrama 1, and Víctor H. Alvarez 1 Fac. de Ingeniería, Depto. Ing. Mecánica, Univ. de la
More informationThe International Association for the Properties of Water and Steam
IAPWS G7-04 The International Association for the Properties of Water and Steam Kyoto, Japan September 2004 Guideline on the Henry s Constant and Vapor-Liquid Distribution Constant for Gases in H 2 O and
More informationMolecular simulation of adsorption from dilute solutions
Vol. 52 No. 3/2005 685 689 on-line at: www.actabp.pl Molecular simulation of adsorption from dilute solutions Werner Billes Rupert Tscheliessnig and Johann Fischer Institut für Verfahrens- und Energietechnik
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 informationRESEARCH NOTE. DOUGLAS HENDERSON? Department of Chemistry and Biochemistry, Brigham Young University, Provo, Utah USA
MOLECULAR PHYSICS, 1999, VOL. 96, No. 7, 1145-1149 RESEARCH NOTE A simple theory for the partial molar volumes of a binary mixture DOUGLAS HENDERSON? Department of Chemistry and Biochemistry, Brigham Young
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 informationThermodynamic Functions at Isobaric Process of van der Waals Gases
Thermodynamic Functions at Isobaric Process of van der Waals Gases Akira Matsumoto Department of Material Sciences, College of Integrated Arts Sciences, Osaka Prefecture University, Sakai, Osaka, 599-853,
More informationEvaluation of the pressure tensor and surface tension for molecular fluids with discontinuous potentials using the volume perturbation method
Evaluation of the pressure tensor and surface tension for molecular fluids with discontinuous potentials using the volume perturbation method Guadalupe Jiménez-Serratos, Carlos Vega, and Alejandro Gil-Villegas
More informationPhase equilibria properties of binary and ternary systems containing isopropyl ether + isobutanol + benzene at K.
Phase equilibria properties of binary and ternary systems containing isopropyl ether + isobutanol + benzene at 313.15 K. R.M. Villamañán 1, M.C. Martín 2, C.R. Chamorro 2, M.A. Villamañán 2, J.J. Segovia
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 informationChE 524 A. Z. Panagiotopoulos 1
ChE 524 A. Z. Panagiotopoulos 1 VIRIAL EXPANSIONS 1 As derived previously, at the limit of low densities, all classical fluids approach ideal-gas behavior: P = k B T (1) Consider the canonical partition
More informationIdeal Gas Behavior. NC State University
Chemistry 331 Lecture 6 Ideal Gas Behavior NC State University Macroscopic variables P, T Pressure is a force per unit area (P= F/A) The force arises from the change in momentum as particles hit an object
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 informationA method to obtain thermodynamic fundamental equations. André Serrenho, Tânia Sousa, Tiago Domingos
A method to obtain thermodynamic fundamental equations. André Serrenho, Tânia Sousa, Tiago Domingos Environmental and Energy Section, DEM, Instituto Superior Técnico Av. Rovisco Pais, 1, 1049-001 Lisboa,
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 informationEquations of State. Equations of State (EoS)
Equations of State (EoS) Equations of State From molecular considerations, identify which intermolecular interactions are significant (including estimating relative strengths of dipole moments, polarizability,
More informationPermeation of Hexane Isomers across ZSM-5 Zeolite Membranes
2618 Ind. Eng. Chem. Res. 2000, 39, 2618-2622 Permeation of Hexane Isomers across ZSM-5 Zeolite Membranes Rajamani Krishna* and Dietmar Paschek Department of Chemical Engineering, University of Amsterdam,
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 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 informationCHEM-UA 652: Thermodynamics and Kinetics
1 CHEM-UA 652: Thermodynamics and Kinetics Notes for Lecture 4 I. THE ISOTHERMAL-ISOBARIC ENSEMBLE The isothermal-isobaric ensemble is the closest mimic to the conditions under which most experiments are
More informationChapter 3 PROPERTIES OF PURE SUBSTANCES
Thermodynamics: An Engineering Approach Seventh Edition in SI Units Yunus A. Cengel, Michael A. Boles McGraw-Hill, 2011 Chapter 3 PROPERTIES OF PURE SUBSTANCES Copyright The McGraw-Hill Companies, Inc.
More informationA new algorithm for Reverse Monte Carlo simulations
A new algorithm for Reverse Monte Carlo simulations Fernando Lus B. da Silva, Bo Svensson, Torbjörn Åkesson, and Bo Jönsson Citation: The Journal of Chemical Physics 109, 2624 (1998); doi: 10.1063/1.476861
More informationSolutions and Ions. Pure Substances
Class #4 Solutions and Ions CHEM 107 L.S. Brown Texas A&M University Pure Substances Pure substance: described completely by a single chemical formula Fixed composition 1 Mixtures Combination of 2 or more
More informationError analysis in Barker s method: extension to ternary systems
Fluid hase quilibria 154 1999 05 11 rror analysis in Barer s method: extension to ternary systems I.M.A. Fonseca, L.Q. Lobo ) Departamento de ngenharia Quımica, UniÕersidade de Coimbra, 3000 Coimbra, ortugal
More informationThermodynamics I. Properties of Pure Substances
Thermodynamics I Properties of Pure Substances Dr.-Eng. Zayed Al-Hamamre 1 Content Pure substance Phases of a pure substance Phase-change processes of pure substances o Compressed liquid, Saturated liquid,
More informationChemical Potential, Helmholtz Free Energy and Entropy of Argon with Kinetic Monte Carlo Simulation
Chemical Potential, Helmholtz Free Energy and Entropy of Argon with Kinetic Monte Carlo Simulation Chunyan Fan a, D. D. Do a *, D. icholson a and E. Ustinov b a School of Chemical Engineering, University
More informationThermophysical Properties of Ethane from Cubic Equations of State
Thermophysical Properties of Ethane from Cubic Equations of State MIHAELA NOUR, DANIELA DUNA, MIRELA IONITA, VIOREL FEROIU *, DAN GEANA Politehnica University Bucharest, Department of Inorganic Chemistry,
More informationMeasurement and Calculation of Physico-Chemical Properties of Binary Mixtures Containing Xylene and 1- Alkanol
Chemical Methodologies 2(2018) 308-314 Chemical Methodologies Journal homepage: http://chemmethod.com Original Research article Measurement and Calculation of Physico-Chemical Properties of Binary Mixtures
More informationKirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies
Supporting Information for Kirkwood-Buff Integrals for Aqueous Urea Solutions Based upon the Quantum Chemical Electrostatic Potential and Interaction Energies Shuntaro Chiba, 1* Tadaomi Furuta, 2 and Seishi
More informationCHEM 254 EXPERIMENT 2 Critical point determination for SF 6
CHEM 254 EXPERIMENT 2 Critical point determination for SF 6 The equation of state of a gas defines the relationship between the pressure, temperature and volume of the gas. For ideal gases the equation
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 informationChapter 1. The Properties of Gases Fall Semester Physical Chemistry 1 (CHM2201)
Chapter 1. The Properties of Gases 2011 Fall Semester Physical Chemistry 1 (CHM2201) Contents The Perfect Gas 1.1 The states of gases 1.2 The gas laws Real Gases 1.3 Molecular interactions 1.4 The van
More 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 informationChapter 4 Phase Transitions. 4.1 Phenomenology Basic ideas. Partition function?!?! Thermodynamic limit Statistical Mechanics 1 Week 4
Chapter 4 Phase Transitions 4.1 Phenomenology 4.1.1 Basic ideas Partition function?!?! Thermodynamic limit 4211 Statistical Mechanics 1 Week 4 4.1.2 Phase diagrams p S S+L S+G L S+G L+G G G T p solid triple
More informationPreliminary Evaluation of the SPUNG Equation of State for Modelling the Thermodynamic Properties of CO 2 Water Mixtures
Available online at www.sciencedirect.com Energy Procedia 26 (2012 ) 90 97 2 nd Trondheim Gas Technology Conference Preliminary Evaluation of the SPUNG Equation of State for Modelling the Thermodynamic
More informationEffect of entropy on the dynamics of supercooled liquids: new results from high pressure data
Downloaded from http://polymerphysics.net Philosophical Magazine, Vol. 87, Nos. 3 5, 21 January 11 February 27, 459 467 Downloaded By: [Naval Research Laboratory] At: 21:43 26 March 27 Effect of entropy
More informationVol. 114 (2008) ACTA PHYSICA POLONICA A No. 6 A
Vol. 114 (2008) ACTA PHYSICA POLONICA A No. 6 A Optical and Acoustical Methods in Science and Technology Effects of Solvation of 2-Methylpyridine and 2,6-Dimethylpyridine in Dilute Solutions in Water and
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 informationReal Gases. Sections (Atkins 6th Ed.), (Atkins 7-9th Eds.)
Real Gases Sections 1.4-1.6 (Atkins 6th Ed.), 1.3-1.5 (Atkins 7-9th Eds.) Molecular Interactions Compression factor Virial coefficients Condensation Critical Constants Van der Waals Equation Corresponding
More informationProperties of real fluids in critical region: third virial coefficient
Indian J hys (February 2014) 88(2):185 191 DOI 10.1007/s12648-013-0402-5 ORIGINAL AER roperties of real fluids in critical region: third virial coefficient R Khordad*, B Mirhosseini and M M Mirhosseini
More informationReview of differential and integral calculus and introduction to multivariate differential calculus.
Chemistry 2301 Introduction: Review of terminology used in thermodynamics Review of differential and integral calculus and introduction to multivariate differential calculus. The properties of real gases:
More informationChapter 3 PROPERTIES OF PURE SUBSTANCES SUMMARY
Chapter 3 PROPERTIES OF PURE SUBSTANCES SUMMARY PURE SUBSTANCE Pure substance: A substance that has a fixed chemical composition throughout. Compressed liquid (sub-cooled liquid): A substance that it is
More informationPressure-enthalpy driven molecular dynamics for thermodynamic property calculation II: applications
Fluid Phase Equilibria 200 (2002) 93 110 Pressure-enthalpy driven molecular dynamics for thermodynamic property calculation II: applications Loukas I. Kioupis, Gaurav Arya, Edward J. Maginn Department
More informationChapter 11 Solution Thermodynamics: Theory
Chapter 11 Solution Thermodynamics: Theory Chapter 6: constant composition fluids. Most commonly in chemical engineering application multi component mixtures of gases or liquids undergo (composition change).
More informationChain length effects on aqueous alkane solubility near the solvent s critical point
Fluid Phase Equilibria 183 184 (2001) 289 294 Chain length effects on aqueous alkane solubility near the solvent s critical point Eric M. Yezdimer a,b,, Ariel A. Chialvo c,d, Peter T. Cummings b,e a Department
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 informationModelling the Solubility of Solid Aromatic Compounds in Supercritical Fluids
Modelling the Solubility of Solid Aromatic Compounds in Supercritical Fluids VIOREL FEROIU*, OCTAVIAN PARTENIE, DAN GEANA Politehnica University of Bucharest, Department of Applied Physical Chemistry and
More informationImperfect Gases. NC State University
Chemistry 431 Lecture 3 Imperfect Gases NC State University The Compression Factor One way to represent the relationship between ideal and real gases is to plot the deviation from ideality as the gas is
More informationPhase equilibria for the oxygen water system up to elevated temperatures and pressures
Fluid Phase Equilibria 222 223 (2004) 39 47 Phase equilibria for the oxygen water system up to elevated temperatures and pressures Xiaoyan Ji a,b, Xiaohua Lu b, Jinyue Yan a,c, a Department of Chemical
More informationr sat,l T sr sat,l T rf rh Ž 4.
Fluid Phase Equilibria 150 151 1998 215 223 Extended corresponding states for pure polar and non-polar fluids: an improved method for component shape factor prediction Isabel M. Marrucho a, James F. Ely
More informationMolecular Modeling of Matter
Molecular Modeling of Matter Keith E. Gubbins Lecture 1: Introduction to Statistical Mechanics and Molecular Simulation Common Assumptions Can treat kinetic energy of molecular motion and potential energy
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 informationMean spherical model-structure of liquid argon
Prami0a, Vol. 6, No 5, 1976, pp. 284-290. Printed in ndia. Mean spherical model-structure of liquid argon R V GOPALA RAO and T NAMMALVAR Department of Physical Chemistry, Jadavpur University, Calcutta
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 informationA Corresponding State Theory for the Viscosity of Liquids Bull. Korean Chem. Soc. 2008, Vol. 29, No Articles
A Corresponding State Theory for the Viscosity of Liquids Bull. Korean Chem. Soc. 2008, Vol. 29, No. 1 33 Articles A Corresponding State Theory for the Viscosity of Liquids Wonsoo Kim * and Sukbae Lee
More informationKEMS448 Physical Chemistry Advanced Laboratory Work. Freezing Point Depression
KEMS448 Physical Chemistry Advanced Laboratory Work Freezing Point Depression 1 Introduction Colligative properties are properties of liquids that depend only on the amount of dissolved matter (concentration),
More information