UNIVERSITY of CALIFORNIA RIVERSIDE Density Functional

Total Page:16

File Type:pdf, Size:1020Kb

UNIVERSITY of CALIFORNIA RIVERSIDE Density Functional UNIVERSITY OF CALIFORNIA RIVERSIDE Density Functional Theory for High-Throughput Screening of Chemicals and Materials A Dissertation submitted in partial satisfaction of the requirements for the degree of Doctor of Philosophy in Chemical and Environmental Engineering by Jia Fu June 2015 Dissertation Committee: Dr. Jianzhong Wu, Chairperson Dr. David Kisailus Dr. Chia-En Chang Copyright by Jia Fu 2015 The Dissertation of Jia Fu is approved: Committee Chairperson University of California, Riverside Acknowledgements I feel grateful to many people for their great support over the years. First and foremost, I would like to thank my supervisor Prof. Jianzhong Wu for his superb and insightful guidance. Learned from him, I start to know how to think as a theoretician, to solve problems in an elegant way and to enjoy the beauty of statistical mechanics. His strong enthusiasm motivate me and shows an excellent example on how to be a successful chemical engineer and professor. I sincerely appreciate my PhD committee, Prof. David Kisailus, Prof. Chia-en A. Chang for their great advice and help in my qualifying exam and final defense. I also wish to express my special thanks to Prof. De-en Jiang, Prof. Shuangliang Zhao, Dr. Zhehui Jin, Dr. Ke Wang, Dr. Yu Liu and Dr. Jehoon Kim for their valuable support during my PhD studying. Also all the good friends that I have made in Riverside for their encouragement throughout this special period. What’s more, I greatly appreciate the endless support and patience of my parents during this long journey, who support and encourage me more than anyone else to pursue an academic career. iv In this dissertation, some parts are the reprint of the materials appeared in Journal of Physical Chemistry Letter (Volume 4, Issue 21, Page 3687-3691, October 2013), Journal Computer-Aided Molecular Design (Volume 28, Issue 3, Page 299-304, March 2014), Chemical Engineering Science (Volume 126, Page 370-382, April 2015), Langmuir (Volume 29, Issue 42, Page 12997-13002, 2013), Journal of Physical Chemistry C (Volume 119, Issue 10, Page 5374-5385, 2015), Journal of Colloid and Interface Science (Volume 438, Page 191-195, January 2015). This work is financially sponsored by the US Department of Energy (DOE) (DE-FG02-06ER46296 and DE-AC02-05CH11231), the National Science Foundation (NSF-CBET-0852353, NSF-CBET-1404046). This work utilizes supercomputer from the National Energy Research Scientific Computing Center (NERSC). v ABSTRACT OF THE DISSERTATION Density Functional Theory for High-Throughput Screening of Chemicals and Materials by Jia Fu Doctor of Philosophy Graduate Program in Chemical and Environmental Engineering University of California, Riverside, June 2015 Professor Jianzhong Wu, Chairperson Multiscale modelling is an essential part in chemical engineering field, recently with the developments in physical and computer sciences, computer modelling plays a more and more important role in quantitative predictions of chemical reactions and thermodynamic data, covering broad applications to problems of practical concern. Whereas classical simulations and quantum density functional theory methods have become routine tools for engineering applications, the power of classical density functional theory has not been widely recognized, which has already shown great success in many model systems. Considering its efficiency and accuracy, it’s necessary to establish a systematic way to utilize this method applied in realistic applications. During this PhD research, I mainly focus on developing classical density functional theory for large-scale screening of chemicals and materials, especially for organic compounds solvation free energy and gas storage and transport properties in nanoporous materials. vi For the hydration free energy calculation part, we try to build a reliable multiscale procedure to predict the water solvation free energies of organic compounds at ambient conditions, by combination of quantum mechanics, statistical mechanics and molecular mechanics force fields. Using the experimental data of hydration free energies as the benchmark, we show that the new method, taking the explicit solvent molecules into consideration, shows a comparable accuracy as full atomistic molecular simulations, but over one magnitude faster calculation speed. We also find that the theoretical results are sensitive to the input generation methods, i.e. selection of quantum-mechanical methods for determining atomic charges and solute configurations, the assignment of the van der Waals parameters and the solvent models. In generally speaking, HF/6- 31G/Vac/OPLS/ChelpG gives the best prediction, in particular those less hydrophilic ones, the overall average unsiged error is 1.35 kcal/mol compared to 700 common organic molecular hydration free energy database, close to the experimental measurement chemical accuracy limit 1 kcal/mol. On the other hand, we have investigated the four representative versions of non- local density functionals as efficient methods to predict gas adsorption in both simple slit pore model system and a large library of real metal-organic frameworks (MOFs) under a broad range of temperatures and pressures, compared with grand canonical Monte Carlo simulation data. Overall all the four methods are reasonably accurate in comparison with the simulation results, but our results show for each specific gas and interested condition, there is a best candidate: for DOE’s hydrogen storage target condition, it’s FMSA, while for ARPA-E’s methane storage target condition, it’s MFA. Besides, combing with entropy vii scaling method, we can also predict the gas self-diffusion coefficient in the same one time calculation with an average calculating time 30 seconds, far beyond the classical molecular simulation ways. From a computational perspective, the classical density functional theory is shown as a versatile and promising tool for large-scale screening of thousands of chemicals and nanostructured materials with a modern desktop computer. viii Table of Contents Title Page Copyright Page Approval Page Acknowledgements ............................................................................................................ iv Abstract of the dissertation ............................................................................................... vi Chapter 1. Introduction .......................................................................................................1 1.1 Background ..............................................................................................................1 1.1.1 Quantum Mechanics .........................................................................................1 I. Wavefunction Methods..................................................................................1 II. Quantum Density Functional Theory ...........................................................2 1.1.2 Classical Simulation Methods ..........................................................................6 I. Molecular Mechanic Force Fields .................................................................6 II. Molecular Dynamics ..................................................................................10 III. Monte Carlo Method .................................................................................11 1.1.4 Hydration Free Energy ..................................................................................13 1.1.3 Gas Adsorption in Nanoporous Materials .....................................................16 1. 2 Dissertation organization .......................................................................................18 Chapter 2. Basic Concepts of Classical Density Functional Theory .................................29 2. 1 Density Profile ........................................................................................................29 2. 2 Density Functional Theory .....................................................................................31 ix 2. 3 Intrinsic Helmholtz Free Energy ............................................................................32 2. 4 Euler-Lagrange Equation .......................................................................................33 2. 5 Excess Helmholtz Free Energy ..............................................................................33 Chapter 3. Computer-Aided Materials Design .................................................................40 3.1 Density functional methods for fast screening of metal-organic frameworks for hydrogen storage ...........................................................................................................40 3.2 Extension of excess-entropy scaling for the self-diffusivity of confined gases ......69 3.3 Classical Density Functional Theory for Systematic Predictions of Thermodynamic and Kinetic Data for CH4 Adsorption in Metal-Organic Frameworks ..........................90 3.4 Seeking Metal-Organic Frameworks (MOFs) for Methane Storage in Natural Gas Vehicles .......................................................................................................................129 3.5 Nitrogen-Doped Porous Aromatic Frameworks for Enhanced CO2 Adsorption ..155 Chapter 4. Fast Screening Method for Predicting Hydration Free Energy .....................194 4.1 High-throughput prediction of the hydration free energies of small molecules from a classical density functional theory ..............................................................................194 4.2
Recommended publications
  • Chapter 6 - Engineering Equations of State II
    Introductory Chemical Engineering Thermodynamics By J.R. Elliott and C.T. Lira Chapter 6 - Engineering Equations of State II. Generalized Fluid Properties The principle of two-parameter corresponding states 50 50 40 40 30 30 T=705K 20 20 ρ T=286K L 10 10 T=470 T=191K T=329K 0 0 0 0.1 0.2 0.3 0.4 0 0.1 0.2 0.3 0.4 0.5 -10 T=133K -10 ρV Density (g/cc) Density (g/cc) VdW Pressure (bars) in Methane VdW Pressure (bars) in Pentane Critical Definitions: Tc - critical temperature - the temperature above which no liquid can exist. Pc - critical pressure - the pressure above which no vapor can exist. ω - acentric factor - a third parameter which helps to specify the vapor pressure curve which, in turn, affects the rest of the thermodynamic variables. ∂ ∂ 2 P = P = 0 and 0 at Tc and Pc Note: at the critical point, ∂ρ ∂ρ 2 T T Chapter 6 - Engineering Equations of State Slide 1 II. Generalized Fluid Properties The van der Waals (1873) Equation Of State (vdW-EOS) Based on some semi-empirical reasoning about the ways that temperature and density affect the pressure, van der Waals (1873) developed the equation below, which he considered to be fairly crude. We will discuss the reasoning at the end of the chapter, but it is useful to see what the equations are and how we use them before deriving the details. The vdW-EOS is: bρ aρ 1 aρ Z =+1 −= − ()1− bρ RT()1− bρ RT van der Waals’ trick for characterizing the difference between subcritical and supercritical fluids was to recognize that, at the critical point, ∂P ∂ 2 P = 0 and = 0 at T , P ∂ρ ∂ρ 2 c c T T Since there are only two “undetermined parameters” in his EOS (a and b), he has reduced the problem to one of two equations and two unknowns.
    [Show full text]
  • On the Van Der Waals Gas, Contact Geometry and the Toda Chain
    entropy Article On the van der Waals Gas, Contact Geometry and the Toda Chain Diego Alarcón, P. Fernández de Córdoba ID , J. M. Isidro * ID and Carlos Orea Instituto Universitario de Matemática Pura y Aplicada, Universidad Politécnica de Valencia, 46022 Valencia, Spain; [email protected] (D.A.); [email protected] (P.F.d.C.); [email protected] (C.O.) * Correspondence: [email protected] Received: 7 June 2018; Accepted: 16 July 2018; Published: 26 July 2018 Abstract: A Toda–chain symmetry is shown to underlie the van der Waals gas and its close cousin, the ideal gas. Links to contact geometry are explored. Keywords: van der Waals gas; contact geometry; Toda chain 1. Introduction The contact geometry of the classical van der Waals gas [1] is described geometrically using a five-dimensional contact manifold M [2] that can be endowed with the local coordinates U (internal energy), S (entropy), V (volume), T (temperature) and p (pressure). This description corresponds to a choice of the fundamental equation, in the energy representation, in which U depends on the two extensive variables S and V. One defines the corresponding momenta T = ¶U/¶S and −p = ¶U/¶V. Then, the standard contact form on M reads [3,4] a = dU + TdS − pdV. (1) One can introduce Poisson brackets on the four-dimensional Poisson manifold P (a submanifold of M) spanned by the coordinates S, V and their conjugate variables T, −p, the nonvanishing brackets being fS, Tg = 1, fV, −pg = 1. (2) Given now an equation of state f (p, T,...) = 0, (3) one can make the replacements T = ¶U/¶S, −p = ¶U/¶V in order to obtain ¶U ¶U f − , ,..
    [Show full text]
  • Comments on Failures of Van Der Waals' Equation at the Gas–Liquid
    Comments on Failures of van der Waals’ Equation at the Gas–Liquid Critical Point, L. V. Woodcock, International Journal of Thermophysics (2018) 39:120 I.H. Umirzakov Institute of Thermophysics, Novosibirsk, Russia [email protected] Abstract These comments are a response to the discussion presented in the above paper concerning the “New comment on Gibbs Density Surface of Fluid Argon: Revised Critical Parameters” by Umirzakov. Here we show that: Woodcock’s results obtained for the dependencies for the isochoric heat capacity, excess Gibbs energy and coexisting difference functional of argon, and coexisting densities of liquid and vapor of the van der Waals fluid and presented in all Figures are incorrect; his Table includes incorrect values of coexisting difference functional; his paper includes many incorrect equations, mathematical and logical errors and physically incorrect assertions concerning the temperature dependences of the isochoric heat capacity and entropy of real fluids; most of the his conclusions are based on the above errors, incorrect data, incorrect comparisons and incorrect dependencies; and most of his conclusions are invalid. We also show that the van der Waals equation of state quantitatively describes the dependencies of saturation pressure on vapor density and temperature near critical point, and the equation of state can describe qualitatively the reduced excess Gibbs energy, rigidity and densities of coexisting liquid and vapor of argon, including the region near critical point. Keywords Coexistence · Critical point · First-order phase transition · Liquid · Phase equilibrium · Vapor 1. Introduction Our comments are a response to a discussion of the article “New comment on Gibbs Density Surface of Fluid Argon: Revised Critical Parameters” by Umirzakov [1] held in the paper [2].
    [Show full text]
  • Convolutional Neural Networks for Atomistic Systems
    Convolutional neural networks for atomistic systems Kevin Ryczko Department of Physics, University of Ottawa Kyle Mills Department of Physics, University of Ontario Institute of Technology Iryna Luchak Department of Electrical & Computer Engineering, University of British Columbia Christa Homenick National Research Council of Canada Isaac Tamblyn Department of Physics, University of Ontario Institute of Technology Department of Physics, University of Ottawa National Research Council of Canada Abstract We introduce a new method, called CNNAS (convolutional neural networks for atom- istic systems), for calculating the total energy of atomic systems which rivals the com- putational cost of empirical potentials while maintaining the accuracy of ab initio calculations. This method uses deep convolutional neural networks (CNNs), where the input to these networks are simple representations of the atomic structure. We use this approach to predict energies obtained using density functional theory (DFT) for 2D hexagonal lattices of various types. Using a dataset consisting of graphene, hexagonal boron nitride (hBN), and graphene-hBN heterostructures, with and with- out defects, we trained a deep CNN that is capable of predicting DFT energies to an extremely high accuracy, with a mean absolute error (MAE) of 0.198 meV / atom (maximum absolute error of 16.1 meV / atom). To explore our new methodology, arXiv:1706.09496v2 [cond-mat.mtrl-sci] 20 Mar 2018 Email addresses: [email protected] (Kevin Ryczko), [email protected] (Isaac Tamblyn) Preprint submitted to Computational Materials Science March 21, 2018 we investigate the ability of a deep neural network (DNN) in predicting a Lennard- Jones energy and separation distance for a dataset of dimer molecules in both two and three dimensions.
    [Show full text]
  • Chemical Engineering Thermodynamics
    CHEMICAL ENGINEERING THERMODYNAMICS Andrew S. Rosen SYMBOL DICTIONARY | 1 TABLE OF CONTENTS Symbol Dictionary ........................................................................................................................ 3 1. Measured Thermodynamic Properties and Other Basic Concepts .................................. 5 1.1 Preliminary Concepts – The Language of Thermodynamics ........................................................ 5 1.2 Measured Thermodynamic Properties .......................................................................................... 5 1.2.1 Volume .................................................................................................................................................... 5 1.2.2 Temperature ............................................................................................................................................. 5 1.2.3 Pressure .................................................................................................................................................... 6 1.3 Equilibrium ................................................................................................................................... 7 1.3.1 Fundamental Definitions .......................................................................................................................... 7 1.3.2 Independent and Dependent Thermodynamic Properties ........................................................................ 7 1.3.3 Phases .....................................................................................................................................................
    [Show full text]
  • The Van Der Waals Equation of State
    PhD – Environmental Fluid Mechanics – Physics of the Atmosphere University of Trieste – International Center for Theoretical Physics The van der Waals equation of state by Dario B. Giaiotti and Fulvio Stel (1) Regional Meteorological Observatory, via Oberdan, 18/A I-33040 Visco (UD) - ITALY Abstract This lecture deals with a more general form of the equation of state, called van der Waals equation, which gives a better description of reality both under the conceptual and numerical point of view. The distinction from gas and vapour naturally springs out from this equation as well as the existence of an unstable region of the P-V diagram, nor foreseen neither described by the classical ideal gas law. Even if this unstable region is not correctly described by the van der Waals equation under the analytical point of view, with added assumptions (conceptual patch) it is possible to give at least a qualitative description of reality. Last but not least the van der Waals equation foreseen the existence of negative pressures. These negative pressures exist in nature (even if they are unstable) and are fundamental for some biotic mechanisms How to obtain the van der Waals (1873) An equation of state is a relationship between the thermodynamic variables that fully define the state of a system. If the system is air (that is essentially a mixture of nitrogen and oxygen) or water vapour, four thermodynamic variables are enough to determine the state of the system and they can be a combination of pressure, volume, mass (or number of moles) and temperature. In the atmospheric physics the equation of state usually used is the well known pV =n RT where p is pressure n is the number of moles, T is the absolute temperature and R is the ideal gas constant (R = 8.3143 J mol-1 K-1).
    [Show full text]
  • Equilibrium Critical Phenomena in Fluids and Mixtures
    : wil Phenomena I Fluids and Mixtures: 'w'^m^ Bibliography \ I i "Word Descriptors National of ac \oo cop 1^ UNITED STATES DEPARTMENT OF COMMERCE • Maurice H. Stans, Secretary NATIONAL BUREAU OF STANDARDS • Lewis M. Branscomb, Director Equilibrium Critical Phenomena In Fluids and Mixtures: A Comprehensive Bibliography With Key-Word Descriptors Stella Michaels, Melville S. Green, and Sigurd Y. Larsen Institute for Basic Standards National Bureau of Standards Washington, D. C. 20234 4. S . National Bureau of Standards, Special Publication 327 Nat. Bur. Stand. (U.S.), Spec. Publ. 327, 235 pages (June 1970) CODEN: XNBSA Issued June 1970 For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, D.C. 20402 (Order by SD Catalog No. C 13.10:327), Price $4.00. NATtONAL BUREAU OF STAHOAROS AUG 3 1970 1^8106 Contents 1. Introduction i±i^^ ^ 2. Bibliography 1 3. Bibliographic References 182 4. Abbreviations 183 5. Author Index 191 6. Subject Index 207 Library of Congress Catalog Card Number 7O-6O632O ii Equilibrium Critical Phenomena in Fluids and Mixtures: A Comprehensive Bibliography with Key-Word Descriptors Stella Michaels*, Melville S. Green*, and Sigurd Y. Larsen* This bibliography of 1088 citations comprehensively covers relevant research conducted throughout the world between January 1, 1950 through December 31, 1967. Each entry is charac- terized by specific key word descriptors, of which there are approximately 1500, and is indexed both by subject and by author. In the case of foreign language publications, effort was made to find translations which are also cited. Key words: Binary liquid mixtures; critical opalescence; critical phenomena; critical point; critical region; equilibrium critical phenomena; gases; liquid-vapor systems; liquids; phase transitions; ternairy liquid mixtures; thermodynamics 1.
    [Show full text]
  • Understanding and Minimizing Density Functional Failures Using Dispersion Corrections
    Understanding and Minimizing Density Functional Failures Using Dispersion Corrections THÈSE NO 5535 (2012) PRÉSENTÉE LE 23 octobre 2012 À LA FACULTÉ DES SCIENCES DE BASE LABORATOIRE DE DESIGN MOLÉCULAIRE COMPUTATIONNEL PROGRAMME DOCTORAL EN CHIMIE ET GÉNIE CHIMIQUE ÉCOLE POLYTECHNIQUE FÉDÉRALE DE LAUSANNE POUR L'OBTENTION DU GRADE DE DOCTEUR ÈS SCIENCES PAR Stephan Niklaus STEINMANN acceptée sur proposition du jury: Prof. J. Zhu, président du jury Prof. A.-C. Corminboeuf, directrice de thèse Dr I. Tavernelli, rapporteur Dr A. Tkatchenko, rapporteur Dr T. A. Wesolowski, rapporteur Suisse 2012 To Schrödinger’s cat Il n’est pas de destin qui ne se surmonte par le mépris. — Albert Camus, Le mythe de Sisyphe But it ain’t about how hard ya hit. It’s about how hard you can get it and keep moving forward. — Rocky Balboa Acknowledgements First and foremost, my advisor Prof. Clémence Corminboeuf has fostered my scientific curios- ity and provided the framework of the four years I spent happily pursuing the doctoral studies here in Lausanne. I am immensely grateful to her for having accepted me in her group and for believing more in me than I ever did myself. In particular, I have profited tremendously from her constant availability to share her opinion and to discuss various chemical questions. If this thesis is written well, it is because of all the great help and advices Clémence has given me. The members of the “Laboratory for Computational Molecular Design” have never failed to support me. Since my early days, I admire Matthew Wodrich’s love for alkanes and Fabrice Avaltroni’s passion for Bash-scripting and Daniel Jana’s remarkable aptitude to make beautiful pictures.
    [Show full text]
  • Thermodynamic 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 and Sciences, Osaka Prefecture University, Sakai, Osaka, 599-8531, Japan Reprint requests to Prof. A. Matsumoto; Fax. 81-0722-54-9927; E-mail: [email protected] Z. Naturforsch. 55 a, 851-855 (2000); received June 29, 2000 The thermodynamic functions for the van der Waals equation are investigated at isobaric process. The Gibbs free energy is expressed as the sum of the Helmholtz free energy and PV, and the volume in this case is described as the implicit function of the cubic equation for V in the van der Waals equation. Furthermore, the Gibbs free energy is given as a function of the reduced temperature, pressure and volume, introducing a reduced equation of state. Volume, enthalpy, entropy, heat capacity, thermal expansivity, and isothermal compressibility are given as functions of the reduced temperature, pressure and volume, respectively. Some thermodynamic quantities are calculated numerically and drawn graphically. The heat capacity, thermal expansivity, and isothermal compressibility diverge to infinity at the critical point. This suggests that a second-order phase transition may occur at the critical point. Key words: Van der Waals Gases; Isobaric Process; Second-order Phase Transition. 1. Introduction the properties of the determinant. In the case of these functions, a singularity will occur thermodynamically There is a continuing problem concerning the ther- at a critical point, while the volume V(T, P) becomes modynamic behaviours in a supercritical state. The a complicated expression as a root of the cubic equa- van der Waals equation as one of many empirical tion.
    [Show full text]
  • University of Illinois at Chicago Department of Physics
    University of Illinois at Chicago Department of Physics Thermodynamics & Statistical Mechanics Qualifying Examination January 8, 2010 9.00 am - 12.00 pm Full credit can be achieved from completely correct answers to 4 questions. If the student attempts all 5 questions, all of the an- swers will be graded, and the top 4 scores will be counted towards the exam's total score. Various equations, constants, etc. are provided on the last page of the exam. 1. Van der Waals Fluid A mole of fluid obeys the van der Waals equation of state: RT a p = V −b − V 2 , where R, a and b are positive constants, and V > b. (a) Show that the specific heat at constant volume (CV ) is a function of T only. (b) Is the results still true if the fluid obeys the Dieterici equation of state RT 1 p = V −b exp −V RT ? Here, R, a and b are positive constants, and V > b. (c) Find the entropy S(T;V ) of the van der Waals fluid in the case when CV is independent of T . 1 2. Simple paramagnet 1 A simple paramagnet consists of N spin- 2 particles immersed in an external magnetic field 1 Hext pointing in the z-direction. Each spin- 2 particle can only be in either an \up" or \down" state, and the energy of a particle with spin pointing down is +µB, while the energy of particles with spin up is −µB. (a) Show that the thermodynamic potential is Hext G(T;Hext) = −aT log 2 cosh b T , where T is the temperature.
    [Show full text]
  • Ideal Gas Equation of State.Pdf
    EQUATION OF STATE The Ideal-Gas Equation of State Any equation that relates the pressure, temperature, and specific volume of a substance is called an equation of state. The simplest and best known equation of state for substances in the gas phase is the ideal-gas equation of state. Gas and vapor are often used as synonymous words. The vapor phase of a substance is called a gas when it is above the critical temperature. Vapor usually implies a gas that is not far from a state of condensation. It is experimentally observed that at a low pressure the volume of a gas is proportional to its temperature: 1 or Pv = RT where R is the gas constant. The above equation is called the ideal-gas equation of state (ideal gas relation). Since R is a constant for a gas, one can write: P v P v R = 1 1 = 2 2 T T 1 2 where subscripts 1 and 2 denote two states of an ideal gas. Note that the constant R is different for each gas; see Tables A1 and A2 in Cengel book. The ideal gas equation of state can be written in several forms: where / is the molar specific volume. That is the volume per unit mole. Ru = 8.314 kJ / (kmol. K) is the universal gas constant, R = Ru /M. The molar mass M: is defined as the mass of one mole of a substance (in gmole or kgmol). The mass of a system is equal to the product of its molar mass M and the mole number N: m = MN (kg) See Table A-1 for R and M for several substances.
    [Show full text]
  • Van Der Waals Equation of State and PVT – Properties of Real Fluid
    Van der Waals equation of state and PVT – properties of real fluid Umirzakov I. H. Institute of Thermophysics, Pr. Lavrenteva St., 1, Novosibirsk, Russia, 630090 e-mail: [email protected] Abstract It is shown that: in the case when two parameters of the Van der Waals equation of state are defined from the critical temperature and pressure the exact parametrical solution of the equations of the liquid-vapor phase equilibrium of the Van der Waals fluid quantitatively describes the experimental dependencies of the saturated pressure of argon on the temperature and reduced vapor density; it can describe qualitatively the temperature dependencies of the reduced vapor and liquid densities of argon; and it gives the quantitative description of the temperature dependencies of the reduced densities near critical point. When the parameters are defined from the critical pressure and density the parametric solution describes quantitatively the experimental dependencies of the saturated pressure of argon on the density and reduced temperature, it can describe qualitatively the dependencies of the vapor and liquid densities on the reduced temperature, and it gives the quantitative description of the dependencies of the densities on the reduced temperature near critical point. If the parameters are defined from the critical temperature and density then the exact solution describes quantitatively the experimental dependencies of the reduced saturated pressure of argon on the density and temperature, it describes qualitatively the temperature dependencies
    [Show full text]