Aqueous electrolyte solution modelling: Some limitations of the Pitzer equations

Abstract Despite intense efforts, general thermodynamic modelling of aqueous electrolyte solutions still presents a difficult challenge, with no obvious method of choice. Even though the Pitzer equations seemingly provide a well-established theoretical framework applicable to many chemical systems over a wide range of temperatures and pressures, they are not as widely adopted as their early promise might have suggested. This is strikingly illustrated by the simultaneous appearance in the literature of numerous, different (and potentially incompatible) Pitzer models alongside a proliferation of alternative theoretical approaches with inferior capabilities. To better understand this problem, the ability of the Pitzer equations to represent the physicochemical properties of aqueous solutions has been systematically investigated for exemplar electrolyte systems. Pitzer ion-interaction parameters have been calculated for selected systems by least-squares regression analysis of published solution data for activity coefficients, osmotic coefficients, relative enthalpies, heat capacities, volumes and densities to high temperatures and pressures. Although satisfactory fits can be achieved when the ranges of conditions are carefully chosen and when sufficient data are available to constrain the regression, the fits obtained tend otherwise to be unsatisfactory. The Pitzer equations do not cope well with gaps and other deficiencies in the regressed data. Profound difficulties, poorly recognized hitherto, can also arise because of variation in the sensitivity of the Pitzer functions to values for different physicochemical properties when these are combined. Given the dimensionality of numerous related thermodynamic properties, all changing as functions of composition, temperature and pressure, these problems are difficult to detect, let alone address, especially in multicomponent systems. The growing practice of improving fits simply by adding basis functions (thereby increasing the number of adjustable parameters) should be deprecated because it increases the likelihood of error propagation, introduces subjectivity, makes independent verification difficult and has deleterious implications for both automated data processing and for consistency between thermodynamic models.

[1]  W. J. Frederick,et al.  Evaluation of Pitzer ion interaction parameters of aqueous electrolytes at 25.degree.C. 1. Single salt parameters , 1988 .

[2]  D. G. Archer,et al.  Thermodynamic Properties of the NaCl+H2O System. II. Thermodynamic Properties of NaCl(aq), NaCl⋅2H2(cr), and Phase Equilibria , 1992 .

[3]  J. R. Grigera,et al.  Simple electrolytes in the mean spherical approximation , 1976 .

[4]  R. W. Carter,et al.  Thermodynamic Properties of the NaCl + H2O System. 4. Heat Capacities of H2O and NaCl(aq) in Cold-Stable and Supercooled States , 2000 .

[5]  Kenneth S. Pitzer,et al.  Thermodynamics of multicomponent, miscible, ionic systems: theory and equations , 1986 .

[6]  I. Grenthe,et al.  TDB-4 TEMPERATURE CORRECTIONS TO THERMODYNAMIC DATA AND ENTHALPY CALCULATIONS , 1999 .

[7]  E. Königsberger,et al.  FactSage and ChemApp: Two tools for the prediction of multiphase chemical equilibria in solutions , 2008 .

[8]  D. Rowland,et al.  JESS, a joint expert speciation system--V: approaching thermodynamic property prediction for multicomponent concentrated aqueous electrolyte solutions. , 2010, Talanta.

[9]  J. A. Rard,et al.  Thermodynamics of Electrolytes. 13. Ionic Strength Dependence of Higher-Order Terms; Equations for CaCl2 and MgCl2 , 1999 .

[10]  J. A. Rard,et al.  Conversion of parameters between different variants of Pitzer's ion-interaction model, both with and without ionic strength dependent higher-order terms , 2003 .

[11]  N. Møller,et al.  The prediction of mineral solubilities in natural waters: The Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O system to high ionic strengths at 25°C , 1984 .

[12]  K. Pitzer,et al.  Volumetric Properties of Single Aqueous Electrolytes from Zero to Saturation Concentration at 298.15 °K Represented by Pitzer’s Ion‐Interaction Equations , 1996 .

[13]  Kenneth S. Pitzer,et al.  Thermodynamics of concentrated electrolyte mixtures and the prediction of mineral solubilities to high temperatures for mixtures in the system Na-K-Mg-Cl-SO4-OH-H2O , 1987 .

[14]  D. Rowland,et al.  An investigation of Zdanovskii's rule for predicting the water activity of multicomponent aqueous strong electrolyte solutions , 2012 .

[15]  Kevin Murray,et al.  JESS - a Joint Expert Speciation System , 1989 .

[16]  D. Grant-Taylor.  Partial molar volumes of sodium chloride solutions at 200 bar, and temperatures from 175 to 350°C , 1981 .

[17]  P. Saluja,et al.  Apparent molar heat capacities and volumes of aqueous solutions of magnesium chloride, calcium chloride, and strontium chloride at elevated temperatures , 1987 .

[18]  E. Königsberger,et al.  JESS, a Joint Expert Speciation System--IV: a large database of aqueous solution physicochemical properties with an automatic means of achieving thermodynamic consistency. , 2010, Talanta.

[19]  Hongbin Cao,et al.  Determination and Modeling of Solubility for CaSO4·2H2O–NH4+–Cl––SO42––NO3––H2O System , 2012 .

[20]  J. Simonson,et al.  Vapor pressures and isopiestic molalities of concentrated CaCl2(aq), CaBr2(aq), and NaCl(aq) to T = 523 K , 2005 .

[21]  J. A. Rard,et al.  Conversion and optimization of the parameters from an extended form of the ion-interaction model for Ca(NO3)2(aq) and NaNO3(aq) to those of the standard Pitzer model, and an assessment of the accuracy of the parameter temperature representations , 2005 .

[22]  G. Hefter.  Some highs and lows (and in-betweens) of solubility measurements of solid electrolytes , 2013 .

[23]  E. Königsberger,et al.  Densities of Concentrated Alkaline Aluminate Solutions at Temperatures from (323 to 573) K and 10 MPa Pressure , 2010 .

[24]  G. Iglesias-Silva,et al.  Osmotic and Activity Coefficients Using a Modified Pitzer Equation for Strong Electrolytes 1:1 and 1:2 at 298.15 K , 2002 .

[25]  D. E. White,et al.  Heat capacity of aqueous magnesium chloride from 349 to 598 K , 1988 .

[26]  S. Lamperski,et al.  Individual activity coefficients of a solvent primitive model electrolyte calculated from the inverse grand-canonical Monte Carlo simulation and MSA theory , 2011 .

[27]  J. A. Rard,et al.  Isopiestic determination of the osmotic and activity coefficients of aqueous magnesium chloride solutions at 25 .degree.C , 1981 .

[28]  James S. Johnson,et al.  Osmotic and activity coefficients for binary mixtures of sodium chloride, sodium sulfate, magnesium sulfate, and magnesium chloride in water at 25.deg.. III. Treatment with the ions as components , 1970 .

[29]  Ji-guang Li,et al.  Extension of the Simple Equations for Prediction of the Properties of Mixed Electrolyte Solutions to the Mixed Ionic Liquid Solutions , 2011 .

[30]  G. Hefter,et al.  Volumes of MgCl2(aq) at temperatures from 298 K to 623 K and pressures up to 30 MPa , 1997 .

[31]  R. Neueder,et al.  An Attempt of a General Overview , 2009 .

[32]  T. Deng,et al.  Solid–liquid metastable equilibria in the quaternary system Li2SO4 + MgSO4 + Na2SO4 + H2O at T = 263.15 K , 2010 .

[33]  J. N. Brönsted.  CALCULATION OF THE OSMOTIC AND ACTIVITY FUNCTIONS IN SOLUTIONS OF UNI-UNIVALENT SALTS , 1922 .

[34]  N. Møller,et al.  The prediction of mineral solubilities in natural waters: A chemical equilibrium model for the Na-Ca-Cl-SO4-H2O system, to high temperature and concentration , 1988 .

[35]  K. Pitzer,et al.  Thermodynamic Properties of Aqueous Potassium Fluoride under Superambient Conditions , 1999 .

[36]  Kenneth S. Pitzer,et al.  Volumetric Properties of Aqueous Sodium Chloride Solutions , 1982 .

[37]  J. I. Partanen.  Re-evaluation of the Thermodynamic Activity Quantities in Pure Aqueous Solutions of Chlorate, Perchlorate, and Bromate Salts of Lithium, Sodium or Potassium Ions at 298.15 K , 2012, Journal of Solution Chemistry.

[38]  Robert N. Goldberg,et al.  Evaluated activity and osmotic coefficients for aqueous solutions: The alkaline earth metal halides , 1978 .

[39]  Thermodynamics of impurities in hydrometallurgical processes , 2011 .

[40]  J. A. Rard,et al.  Isopiestic determination of the osmotic and activity coefficients of aqueous mixtures of sodium chloride and magnesium chloride at 25.degree.C , 1982 .

[41]  Eric W. Lemmon,et al.  A Formulation for the Static Permittivity of Water and Steam at Temperatures from 238 K to 873 K at Pressures up to 1200 MPa, Including Derivatives and Debye–Hückel Coefficients , 1997 .

[42]  K. Pitzer,et al.  Thermodynamic Properties of Aqueous Potassium Sulfate under Superambient Conditions , 1999 .

[43]  Jun Ji,et al.  Solubility in the ternary system LiCl + MgCl2 + H2O at 60 and 75°C , 2010 .

[44]  John H. Weare,et al.  The prediction of mineral solubilities in natural waters: the NaKMgCaClSO4H2O system from zero to high concentration at 25° C , 1980 .

[45]  R. F. Platford.  Isopiestic measurements on the system water-sodium chloride-magnesium chloride at 25.deg. , 1968 .

[46]  P. Spitzer,et al.  Measurement of pH. Definition, standards, and procedures (IUPAC Recommendations 2002) , 2002 .

[47]  A. Felmy,et al.  Thermodynamic properties of aqueous calcium nitrate {Ca ( NO3)2}to the temperature 373 K including new enthalpy of dilution data , 2000 .

[48]  D. Fraenkel.  Monoprotic mineral acids analyzed by the smaller-ion shell model of strong electrolyte solutions. , 2011, The journal of physical chemistry. B.

[49]  G. Marion,et al.  Mineral solubilities in the Na-K-Mg-Ca-Cl-SO4-H2O system: a re-evaluation of the sulfate chemistry in the Spencer-Møller-Weare model , 1999 .

[50]  E. Königsberger,et al.  Heat capacities of aqueous solutions of sodium hydroxide and water ionization up to 300 °C at 10 MPa , 2008 .

[51]  Meng-Ting Hsieh,et al.  A predictive model for the excess gibbs free energy of fully dissociated electrolyte solutions , 2011 .

[52]  Erich Königsberger,et al.  A new optimization routine for chemsage , 1995 .

[53]  V. Afanas’ev.  Solvation of electrolytes and nonelectrolytes in aqueous solutions. , 2011, The journal of physical chemistry. B.

[54]  J. Wood.  Thermodynamics of brine-salt equilibria — II. The system NaCl-KCl-H2O from 0 to 200°C , 1976 .

[55]  K. Pitzer.  Theoretical considerations of solubility with emphasis on mixed aqueous electrolytes , 1986 .

[56]  K. Pitzer,et al.  Thermodynamic properties of Aqueous Magnesium Chloride Solutions from 250 to 600 K and to 100 MPa , 1998 .

[57]  K. Pitzer,et al.  Thermodynamics of electrolytes. II. Activity and osmotic coefficients for strong electrolytes with one or both ions univalent , 1973 .

[58]  M. Donohue,et al.  Recent Advances in Modeling Thermodynamic Properties of Aqueous Strong Electrolyte Systems , 1997 .

[59]  K. Pitzer,et al.  Thermodynamics of aqueous KOH over the full range to saturation and to 573 K , 1996 .

[60]  E. Königsberger,et al.  Densities of NaOH(aq) at Temperatures from (323 to 573) K and 10 MPa Pressure , 2007 .

[61]  E. Königsberger,et al.  A generic and updatable Pitzer characterization of aqueous binary electrolyte solutions at 1 bar and 25 °C , 2011 .

[62]  S. Altmann,et al.  Modelling coupled chemico-osmotic and advective–diffusive transport of nitrate salts in the Callovo-Oxfordian Clay , 2012 .

[63]  Rui Sun,et al.  Prediction of vapor-liquid equilibrium and PVTx properties of geological fluid system with SAFT-LJ EOS including multi-polar contribution. Part I: Application to H2O-CO2 system , 2010 .

[64]  R. Mesmer,et al.  Thermodynamic properties of aqueous solutions of the alkali metal chlorides to 250.degree.C , 1983 .

[65]  Christopher M.A. Brett,et al.  Measurement of pH. Definition, standards, and procedures , 2002 .

[66]  J. N. Brönsted.  STUDIES ON SOLUBILITY. IV. THE PRINCIPLE OF THE SPECIFIC INTERACTION OF IONS , 1922 .

[67]  J. Simonson,et al.  CaCl2(aq) at elevated temperatures. Enthalpies of dilution, isopiestic molalities, and thermodynamic properties , 1994 .

[68]  J. A. Rard,et al.  Isopiestic Investigation of the Osmotic and Activity Coefficients of Aqueous NaBr and the Solubility of NaBr·2H2O(cr) at 298.15 K: Thermodynamic Properties of the NaBr + H2O System over Wide Ranges of Temperature and Pressure , 1995 .

[69]  Gunnar Eriksson,et al.  FactSage thermochemical software and databases - recent developments , 2009 .

[70]  R. D. Lanier.  Activity coefficients of NaCl in aqueous three-component solutions by cation-sensitive glass electrodes. [25/sup 0/C] , 1965 .

[71]  Gordon Atkinson,et al.  Thermodynamics of concentrated electrolyte mixtures. 5. A review of the thermodynamic properties of aqueous calcium chloride in the temperature range 273.15-373.15 K , 1985 .

[72]  R. Mesmer,et al.  The enthalpy of dilution of aqueous sodium chloride to 673 K using a new heat-flow and liquid-flow microcalorimeter. Excess thermodynamic properties and their pressure coefficients , 1984 .

[73]  G. Perron,et al.  Heat capacities and volumes of NaCl, MgCl2, CaCl2, and NiCl2 up to 6 molal in water , 1981 .

[74]  C. Christov.  Chemical Equilibrium Model of Solution Behavior and Bishofite (MgCl2•6H2O(cr)) and Hydrogen―Carnallite (HCl•MgCl2•7H2O(cr)) Solubility in the MgCl2 + H2O and HCl•MgCl2 + H2O Systems to High Acid Concentration at (0 to 100) °C , 2009 .

[75]  K. Thomsen,et al.  Comparison of activity coefficient models for electrolyte systems , 2009 .

[76]  J. Gates,et al.  Volumetric properties of aqueous NaCl solutions from 0.0025 to 5.0 mol · kg−1, 323 to 600 K, and 0.1 to 40 MPa , 1988 .

[77]  W. Kunz,et al.  MSA-NRTL model for the description of the thermodynamic properties of electrolyte solutions , 2002 .

[78]  O. Redlich,et al.  The Molal Volumes of Electrolytes , 1964 .

[79]  Kenneth S. Pitzer,et al.  Activity Coefficients in Electrolyte Solutions , 2017 .

[80]  J. A. Rard,et al.  Conversion of Parameters Among Variants of Scatchard’s Neutral-Electrolyte Model for Electrolyte Mixtures that Have Different Numbers of Mixing Terms , 2008 .

[81]  D. G. Archer,et al.  Modeling of the thermodynamics of electrolyte solutions to high temperatures including ion association: application to hydrochloric acid , 1990 .

[82]  Chau‐Chyun Chen,et al.  Symmetric Electrolyte Nonrandom Two-Liquid Activity Coefficient Model , 2009 .

[83]  D. Stoica,et al.  pH buffer assessment and Pitzer’s equations , 2011 .

[84]  Kenneth S. Pitzer,et al.  Thermodynamic Properties of Aqueous Sodium Chloride Solutions , 1984 .

[85]  E. Königsberger,et al.  Heat capacities of aqueous sodium hydroxide/aluminate mixtures and prediction of the solubility constant of boehmite up to 300 °C , 2010 .

[86]  C. Perrin.  Numerical Recipes in Fortran 90: The Art of Scientific Computing, Second Edition, Volume 2 (3 CD-ROMs and Manual) By William H. Press, Saul A. Teukolsky, William T. Vetterling, and Brian P. Flannery. Cambridge University Press: New York, 1996. , 1997 .

[87]  K. Pitzer,et al.  Thermodynamics of NaOH(aq) in hydrothermal solutions , 1987 .

[88]  J. Simonson,et al.  The enthalpy of dilution and apparent molar heat capacity of NaOH(aq) to 523 K and 40 MPa , 1989 .

[89]  Kenneth S. Pitzer,et al.  Thermodynamics of aqueous magnesium chloride, calcium chloride, and strontium chloride at elevated temperatures , 1987 .

[90]  Kenneth S. Pitzer,et al.  Thermodynamics of electrolytes. V. effects of higher-order electrostatic terms , 1975 .

[91]  Kenneth S. Pitzer,et al.  Thermodynamics of electrolytes. I. Theoretical basis and general equations , 1973 .

[92]  G. Scatchard,et al.  Osmotic and activity coefficients for binary mixtures of sodium chloride, sodium sulfate, magnesium sulfate, and magnesium chloride in water at 25.deg.. I. Isopiestic measurements on the four systems with common ions , 1968 .

[93]  R. Mesmer,et al.  Thermodynamics of aqueous solutions of the alkali metal sulfates , 1986 .

[94]  Kenneth S. Pitzer,et al.  Densities and apparent molar volumes of aqueous magnesium sulfate and sodium sulfate to 473 K and 100 bar , 1986 .

[95]  E. Königsberger,et al.  Comprehensive Model of Synthetic Bayer Liquors. Part 1. Overview , 2005 .

[96]  P. M. May,et al.  Database of chemical reactions designed to achieve thermodynamic consistency automatically , 2001 .

[97]  J. A. Rard,et al.  Isopiestic Determination of the Osmotic and Activity Coefficients of K2SO4(aq) at the Temperatures 298.15 and 323.15 K, and Revision of the Thermodynamic Properties of the K2SO4 + H2O System , 2002 .

[98]  J. S. Rowlinson,et al.  Molecular Thermodynamics of Fluid-Phase Equilibria , 1969 .

[99]  D. Rowland,et al.  A Pitzer-based characterization of aqueous magnesium chloride, calcium chloride and potassium iodide solution densities to high temperature and pressure , 2013 .

[100]  W. Voigt.  Chemistry of salts in aqueous solutions: Applications, experiments, and theory , 2011 .

[101]  R. Mesmer,et al.  Aqueous solutions of the alkaline-earth metal chlorides at elevated temperatures. Isopiestic molalities and thermodynamic properties , 1996 .

[102]  T. Deng,et al.  Solid-liquid metastable equilibria in the quaternary system (NaCl-KCl-CaCl2-H2O) at 288.15K , 2008 .

[103]  Klaus Hack,et al.  The thermochemistry library ChemApp and its applications , 2007 .

[104]  G. Eriksson,et al.  ChemSage—A computer program for the calculation of complex chemical equilibria , 1990 .

[105]  P. Saluja,et al.  Apparent Molar Heat Capacities and Volumes of Mixed Electrolytes: [NaCl(aq) + CaCl2(aq)], [NaCl(aq) + MgCl2(aq)], and [CaCl2(aq) + MgCl2(aq)] , 1995 .

[106]  J. A. Rard,et al.  Isopiestic determination of the osmotic coefficients of aqueous CaCl/sub 2/ solutions at 25/sup 0/C , 1977 .

[107]  Xueyu Song,et al.  A molecular Debye-Hückel theory and its applications to electrolyte solutions. , 2011, The Journal of chemical physics.

[108]  D. L. Parkhurst,et al.  User's guide to PHREEQC (Version 2)-a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations , 1999 .

[109]  Freeman Dyson,et al.  A meeting with Enrico Fermi , 2004, Nature.

[110]  R. Young,et al.  Modeling phase equilibria and speciation in mixed-solvent electrolyte systems , 2004 .

[111]  K. Pitzer,et al.  Ion interaction approach to calculations of volumetric properties of aqueous multiple-solute electrolyte solutions , 1995 .

[112]  J. A. Rard,et al.  Isopiestic determination of the osmotic and activity coefficients of aqueous cesium chloride, strontium chloride, and mixtures of sodium chloride and cesium chloride at 25.degree.C , 1982 .

[113]  D. Chakravarti,et al.  Measurement of pH , 2008 .

[114]  D. G. Archer.  Thermodynamic Properties of the NaBr+H2O System , 1991 .

[115]  K. Pitzer,et al.  Thermodynamics of aqueous calcium chloride , 1983 .

[116]  Kenneth S. Pitzer,et al.  Volumetric Ion Interaction Parameters for Single-Solute Aqueous Electrolyte Solutions at Various Temperatures , 2000 .

[117]  P. V. Gaans.  Extracting robust and physically meaningful pitzer parameters for thermodynamic properties of electrolytes from experimental data. Ridge regression as an aid to the problem of intercorrelation , 1991 .

[118]  P. M. May,et al.  JESS, A joint expert speciation system-I. Raison d'être. , 1991, Talanta.

[119]  A. Haghtalab,et al.  Nonelectrolyte NRTL-NRF model to study thermodynamics of strong and weak electrolyte solutions , 2011 .

[120]  William H. Press,et al.  Numerical Recipes: FORTRAN , 1988 .

[121]  Yahong Li,et al.  Solubility prediction for the HCl–MgCl2–H2O system at 40 ∘C and solubility equilibrium constant calculation for HCl ⋅ MgCl2 ⋅ 7H2O at 40 ∘C , 2006 .

[122]  S. Marshall,et al.  Least-Squares Analysis of Osmotic Coefficient Data at 25 .degree.C According to Pitzer's Equation. 1. 1:1 Electrolytes , 1995 .

[123]  G. Hefter,et al.  The Experimental Determination of Solubilities: Hefter/Determination of Solubilities , 2003 .

[124]  J. Gmehling,et al.  Thermodynamics of Phase Equilibria in Aqueous Strong Electrolyte Systems , 2011 .

[125]  E. A. Guggenheim.  L. The specific thermodynamic properties of aqueous solutions of strong electrolytes , 1935 .

[126]  Kenneth S. Pitzer,et al.  Apparent molar heat capacity and other thermodynamic properties of aqueous potassium chloride solutions to high temperatures and pressures , 1988 .

[127]  H. Helgeson,et al.  Theoretical prediction of the thermodynamic behavior of aqueous electrolytes at high pressures and temperatures , 1974 .

[128]  H. Helgeson,et al.  Theoretical prediction of the thermodynamic behavior of aqueous electrolytes by high pressures and temperatures; IV, Calculation of activity coefficients, osmotic coefficients, and apparent molal and standard and relative partial molal properties to 600 degrees C and 5kb , 1981 .

[129]  J. A. Rard,et al.  Enthalpies of Dilution of NdCl3(aq) at Temperatures from 297.89 K to 372.08 K and an Extended Pitzer Ion-Interaction Model for the NdCl3 + H2O System , 2004 .

[130]  J. Ananthaswamy,et al.  Thermodynamics of concentrated electrolyte mixtures. I. Activity coefficients in aqueous NaCl−CaCl2 at 25°C , 1982 .

[131]  D. G. Archer.  Thermodynamic Properties of the KCl+H2O System , 1999 .

[132]  F. Millero,et al.  A Chemical Equilibrium Model for Natural Waters , 1998 .

[133]  J. F. Waters,et al.  The free proton concentration scale for seawater pH , 2013 .

[134]  J. A. Rard,et al.  Isopiestic determination of the osmotic and activity coefficients of aqueous mixtures of NaCl and MgCl/sub 2/ at 25/sup 0/C , 1987 .

[135]  D. Rowland,et al.  Thermodynamics of strong aqueous electrolyte solutions at t = 25 °c described by the Hückel equations , 2014 .

[136]  F. Millero,et al.  Apparent Molal Volumes and Heat Capacities of Alkaline Earth Chlorides in Water at 25 "C , 1974 .

[137]  Andrzej Anderko,et al.  A speciation-based model for mixed-solvent electrolyte systems , 2002 .

[138]  Kaj Thomsen,et al.  Modeling of heavy metal salt solubility using the extended UNIQUAC model , 2002 .

[139]  D. G. Archer.  Thermodynamic Properties of the NaNO3+H2O System , 2000 .

[140]  K. Ballerat-Busserolles,et al.  Apparent molar volumes and heat capacities of aqueous magnesium chloride and cadmium chloride at temperatures from 278.15 K to 393.15 K at the pressure 0.35 MPa: a comparison of ion–ion interactions , 2000 .