Calculation of dissociation constants and the relative stabilities of polynuclear clusters of 1:1 electrolytes in hydrothermal solutions at supercritical pressures and temperatures

https://doi.org/10.1016/0016-7037(93)90383-8Get rights and content

Abstract

Monte Carlo calculations (Pitzer and Schreiber, 1987), electrostatic theory (Oelkers and Helgeson, 1990), and revised interpretation of the concentration dependence of conductance data (Oelkers and Helgeson, 1991) indicate that monatomic ions and neutral ion pairs predominate in single 1:1 electrolyte solutions only in the dilute concentration range at supercritical pressures and temperatures where the dielectric constant of H2O (ϵ) is ≲ 15. To assess geologically realistic degrees of formation of polynuclear clusters at higher concentrations, the logarithms of overall dissociation constants (log β) for triple, quadruple, quintuple, and sextuple complexes of monatomic ions in 16 single 1:1 electrolyte solutions at pressures from 500 to 4000 bars and temperatures from 400 to 800°C were generated from the results of Monte Carlo calculations reported by Gillan (1983) for the restricted primitive model using hard sphere diameters (δH) obtained by regression of the logarithms of dissociation constants (log K) for diatomic neutral ion pairs reported by Quist and Marshall (1968), Frantz and Marshall (1984), and Oelkers and Helgeson (1988). In addition, the logarithms of dissociation constants for polynuclear clusters involving more than one type of cation in the system NaCl-KCl-HCl were estimated by assuming a linear dependence of δH on cluster stoichiometry. The calculated and estimated values of log β typically decrease monotonically with increasing cluster size at a given pressure and temperature. Speciation calculations generated from the logβ and log K values using activity coefficients for both charged and neutral aqueous species (Oelkers and Helgeson, 1990, 1991) indicate that polyatomic clusters of monovalent ions are major solute species in supercritical aqueous solutions with 1:1 solute concentrations > 0.5 m at low pressures and high temperatures where the dielectric constant of H2O is low, which is in general agreement with conclusions reached by Pitzer and Schreiber (1987). It follows by analogy that polyatomic clusters involving Si, Al, Ag, Cu, Fe, S, and other elements in geologic systems may increase substantially the solubilities of rock- and ore-forming minerals in concentrated hydrothermal solutions at supercritical temperatures and pressures.

References (56)

  • H.R Corti et al.

    Cluster theory applied to aqueous (2:2) electrolytes over a wide concentration range: Osmotic coefficients, association and redissociation

    J. Chem. Soc., Faraday Trans.

    (1986)
  • H.P Eugster et al.

    Mineral solubilities and speciation in supercritical metamorphic fluids

  • J.D Frantz et al.

    Electrical conductances and ionization constants of salts, acids, and bases in supercritical aqueous fluids. I. Hydrochloric acid from 100° to 700°C and at pressures to 4000 bars

    Amer. J. Sci.

    (1984)
  • J.D Frantz et al.

    Mineralsolution equilibria: V. Solubilities of rock forming minerals in supercritical fluids

    Geochim. Cosmochim. Acta

    (1981)
  • H.L Friedman

    Lewis-Randall to McMillan-Mayer conversion for the thermodynamic excess functions of solutions. Part I. Partial free energy coefficients

    J. Soln. Chem.

    (1972)
  • H.L Friedman et al.

    Electrolyte solutions at equilibrium

  • H.L Friedman et al.

    Ion pairing and related topics

    Pure Applied Chem.

    (1979)
  • R.M Fuoss et al.

    Properties of electrolytic solutions. III. The dissociation constant

    J. Amer. Chem. Soc.

    (1933)
  • R.M Fuoss et al.

    Properties of electrolytic solutions: IV. The conductance minimum and the formation of triple ions due to the action of coulomb forces

    J. Amer. Chem. Soc.

    (1933)
  • R.M Fuoss et al.

    Properties of electrolyte solutions: XV. Thermodynamic properties of very weak electrolytes

    J. Amer. Chem. Soc.

    (1935)
  • J.A Geddes et al.

    Properties of electrolytic solutions: XVIII. Molecular polarisations and polar moments of some electrolytes in benzene solutions

    J. Faraday Soc.

    (1936)
  • M.J Gillan

    Liquid-vapor equilibrium in the restricted primitive model for ionic liquids

    Mol. Phys.

    (1983)
  • I.S Graham et al.

    A Monte Carlo study of the coexistence region of the restricted primitive model

    J. Phys. Chem.

    (1990)
  • H.C Helgeson

    Errata. Thermodynamics of minerals, reactions, and aqueous solutions at high pressures and temperatures

    Amer. J. Sci.

    (1982)
  • H.C Helgeson

    Errata II. Thermodynamics of minerals, reactions, and aqueous solutions at high pressures and temperatures

    Amer. J. Sci.

    (1985)
  • H.C Helgeson et al.

    Summary and critique of the thermodynamic properties of rock forming minerals

    Amer. J. Sci.

    (1978)
  • H.C Helgeson et al.

    Theoretical prediction of the thermodynamic behavior of aqueous electrolytes at high pressures and temperatures: IV. Calculation of activity coefficients, osmotic coefficients and apparent molal and relative partial molal properties to 600°C and 5 kb

    Amer. J. Sci.

    (1981)
  • E Hückel

    Zur Theorie konzentrierterer wässeriger Lösungen starker Elektrolyte

    Phys. Zeitschr.

    (1925)
  • Cited by (43)

    • The nature of NaCl–H<inf>2</inf>O deep fluids from ab initio molecular dynamics at 0.5–4.5 GPa, 20–800 °C, and 1–14 m NaCl

      2020, Geochimica et Cosmochimica Acta
      Citation Excerpt :

      Most Na and/or Cl species form and break down at timescales of 100–300 fs or less, but the average proportion of polyatomic entities and dissociated or undissociated water remained stable. At 14 m, all of the Na and Cl belong to polynuclear clusters and, as such, we cannot define separate NaCl, Na2Cl+, and NaCl2- species as done by Oelkers and Helgeson (1993). In summary, at high P,T, a breakdown of the water structure is associated with better agreement between ab initio-predicted densities and experimental measurements (e.g., Mantegazzi et al. (2013).

    View all citing articles on Scopus

    Present address: Laboratoire de Géochimie, Université Paul Sabatier, 38 Rue des Trente-six Ponts, 31400 Toulouse, France.

    View full text