- Open Access
Ion association in concentrated NaCl brines from ambient to supercritical conditions: results from classical molecular dynamics simulations
© The Royal Society of Chemistry and the Division of Geochemistry of the American Chemical Society 2002
- Received: 05 September 2002
- Accepted: 18 November 2002
- Published: 28 November 2002
Highly concentrated NaCl brines are important geothermal fluids; chloride complexation of metals in such brines increases the solubility of minerals and plays a fundamental role in the genesis of hydrothermal ore deposits. There is experimental evidence that the molecular nature of the NaCl–water system changes over the pressure–temperature range of the Earth's crust. A transition of concentrated NaCl–H2O brines to a "hydrous molten salt" at high P and T has been argued to stabilize an aqueous fluid phase in the deep crust.
In this work, we have done molecular dynamic simulations using classical potentials to determine the nature of concentrated (0.5–16 m) NaCl–water mixtures under ambient (25°C, 1 bar), hydrothermal (325°C, 1 kbar) and deep crustal (625°C, 15 kbar) conditions. We used the well-established SPCE model for water together with the Smith and Dang Lennard-Jones potentials for the ions (J. Chem. Phys., 1994, 100, 3757). With increasing temperature at 1 kbar, the dielectric constant of water decreases to give extensive ion-association and the formation of polyatomic (Na n Cl m ) n-m clusters in addition to simple NaCl ion pairs. Large polyatomic (Na n Cl m ) n-m clusters resemble what would be expected in a hydrous NaCl melt in which water and NaCl were completely miscible. Although ion association decreases with pressure, temperatures of 625°C are not enough to overcome pressures of 15 kbar; consequently, there is still enhanced Na–Cl association in brines under deep crustal conditions.
- Pair Distribution Function
- Concentrate Brine
- Polynuclear Cluster
- Restricted Primitive Model
- Polyatomic Cluster
The compositions of fluid inclusions indicate that extremely concentrated NaCl brines can be involved in the transport of metals and the formation of hydrothermal ore deposits. There is a variety of evidence that Na and Cl ions become highly associated into NaCl ion pairs in high-temperature aqueous solutions at pressures below 2 kbar. Under such conditions, water and NaCl appear to behave like ideal mixtures[2, 3] so that the activity of water, equals its mole fraction ( ). Above 4 kbar, however, there is evidence that NaCl ion pairs become dissociated in NaCl–water mixtures. For example, Quist and Marshall find an increase in the electrical conductance of NaCl–H2Q mixtures between 2 and 4 kbar at 800°C. More recently, Aranovich and Newton measured the activity of water in concentrated NaCl solutions at 600–900°C and 2–15 kbar using the brucite-periclase equilibrium. At high pressure, the activity of water is well approximated by . This is the athermal activity expected if the NaCl–water system behaved as a hydrous salt melt consisting of water and completely dissociated NaCl. The much smaller activity of water under these conditions means that a free solution phase could coexist with an H2O-saturated granitic melt in the deep crust. Such a chloride-rich solution phase must play a fundamental role in scavenging metals prior to the formation of hydrothermal ore deposits.
Atomistic simulations of the NaCl–water system could confirm the interpretation of the experimental data and explore PT regimes that are not easily accessible in the laboratory. The simplest atomistic approach is the "restricted primitive model" where the solvent is replaced by a dielectric continuum and the ions are approximated as hard spheres with equal radii. Restricted primitive model simulations done by Pitzer and Schrieber and Gillan were used by Oelkers and Helgeson[8, 9] to derive association constants of polynuclear clusters in 1 : 1 electrolytes in supercritical aqueous solutions up to 4 kbar. These calculations predicted considerable ion pairing and clustering of Na–Cl in supercritical aqueous fluids. However, the use of a simple constant dielectric to model the solvent means we cannot understand any of the short-range physical interactions. It would be better to use an explicit molecular representation of the solvent. Koneshan and Rasaiah have shown that the degree of ion association in 1 m NaCl (at 683 K with a solvent density of 0.175 g cm-3) is signficantly different for continuum and molecular models simulations of the solvent.
A molecular treatment of the solvent and ion interactions is practical for large systems if we can express the interatomic interactions in terms of two-body potential functions. We can use molecular dynamics (MD) or Monte Carlo methods to sample the statistical mechanics of the system and determine not only thermodynamic quantities but also structural information and ionic speciation.
The SPC/E (Extended Simple Point Charge) water model of Berendsen et al. gives a very simple parameterization for the effective charges and short-range potential for water. In the SPC/E model we represent water molecules by three point charges in a rigid geometry with an HOH bond angle slightly different (109.5°) from that in the gas-phase H2O molecule (104.5°). Although the model does not allow any polarization or dissociation of the water molecules, it accurately reproduces the thermodynamic and dielectric properties of water, at least along the liquid–vapor coexistence curve. In particular, the SPC/E model predicts the critical point of water to be at 362–374°C with a density of 0.29–0.326 gcm-3. For real water, the critical point is at 374°C with a density of 0.322 g cm-3. The SPC/E model also gives a good prediction of the dielectric constant of water and its temperature dependence: the predicted values of ε(SPCE) = 81.0 at 25°C and ε(SPCE) = 6. at T(c) vs. 78.0 and 5.3, respectively, in real water. It is important that the dielectric properties are reproduced as accurately as possible if we wish to predict ion clustering NaCl electrolytes.
Lennard-Jones parameters from Smith and Dang, 199414
Formal ionic charges q i of +1 and -1 are used for Na+ and Cl-. The parameters were used in molecular dynamical simulations of a NaCl ion pair in 216 water molecules. This was a fairly small simulation but it gave very accurate predictions of experimental hydration energies and complex geometries. Smith and Dang predict that at 25°C, Na+ cations will be surrounded by 5.8 water molecules with a Na-O bond length of 2.33 Å. The Cl- anions are solvated by 6.9 waters with a Cl-H distance of 2.22 Å. These are in good agreement with those observed experimentally.
We are assuming that the interatomic potentials we are using can describe the NaCl–water system at the densities of interest at high pressure and temperature. Physically this means that we are assuming that the bonding and charge distributions in the water molecules at 15 kbar, 625°C(p = 1–1.3 g cm-3) are reasonably similar to those at ambient conditions. More sophisticated calculations based on quantum mechanical (ah initio) molecular dynamics are needed to explore this problem. At present, however, such calculations are limited to very small systems for very short run times and cannot address the problem of NaCl ion association as done here. Classical calculations for the equation of state of water at high P and T give good predictions of the water equations of state. This suggests that the potentials used here should be reasonably valid at density conditions of deep crustal fluids.
The formation of ion-pairs and clusters in NaCl solutions has been studied by atomistic MD simulations of dilute to 1 molar NaCl solutions[10, 19–22] Few simulations, however, have addressed the nature of concentrated NaCl solutions at elevated pressure and temperature. Brodholt obtained pair distribution functions and compared the densities of the simulated solutions to those obtained from the equations of state of Archer and of Anderko and Pitzer and found very good agreement. This was done for solutions up to 5.22 m in concentration and between temperatures of 300 to 2573 K, and under pressures of 1 bar up to 5 kbar. Although Brodholt noted the existence of polynuclear clusters, their existence was not directly tested in his study. The objective of the work presented here is to understand the nature of extremely concentrated (up to 16 m) NaCl solutions at conditions from ambient to those of the lower continental crust (625°C, 15 kbar).
Molecular dynamics simulations were done on systems consisting of 975–1100 H2O molecules and 10–160 NaCl pairs. We believe our systems are large enough to describe the system insofar as simulations on much smaller systems (200–300 H2O, 5–20 NaCl pairs) gave essentially the same speciation and density. The simulations were conducted using the MOLDY molecular dynamics program. MOLDY implements periodic boundary conditions and the link-cell method, and uses quaternions to describe molecular orientations. It also uses the Ewald summation in the evaluation of coulombic interactions, and the link-cell algorithm for evaluation of the interactions. The real (rC) and k-space (kC) interaction cut-offs and Ewald parameter (α) were determined using eqn. (2)-(4), such that an accuracy of E = exp(-p) = 10-5 in the coulombic potential energy was achieved:
where N is the number of molecules, V is the volume of the cell and tR/tF has been empirically determined to be 5.5 for MOLDY. For a typical box length of 33 Å (e.g., in a simulation of 1100 H2O molecules and 80 NaCl pairs at 1 kbar), this gave rC = 12.5 Å. We used a time step of 0.25 fs as larger steps caused the system to be unstable. Each run was equilibrated for 25 ps (for temperature scaling, discussed below) and statistics were accumulated for 100–150 ps.
Statistics were collected on a constant NPT ensemble. The pressure was kept constant using the method of Parrinello and Rahman, and fixing all the off-diagonal components of the stress matrix to zero. During the equilibration phase, the temperature of the simulation was controlled by re-scaling the velocities every 20 time steps by the factor
where g is the number of degrees of freedom, kB is the Boltzmann constant, T the desired temperature, and <k> is the rolling average of the system kinetic energy over the previous 20 time steps. The scaling was carried out separately for translational and rotational components. The re-scaling was turned off once the temperature of the system has reached equilibrium about the desired value, in order to generate a realisation of the constant NPT ensemble. During the unsealed runs, temperatures were constant to within ~ 3%.
3.2 Structure of water in NaCl brines
3.2 Na and Cl solvation and ion-association
The coordination numbers and bond lengths for dilute NaCl solutions at 25°C, 1 bar are in good agreement with those from previous calculations and experiment. For dilute NaCl, we predict a coordination number of 5.3 with only waters in the solvation sphere for Na+ at 25°C, I bar. With increasing NaCl molality, we find increased Cl complexation of Na. At 25°C and 1 bar the Na+ ions in an 8 m solution of NaCl have 4.3 waters and 1.0 Cl ligands in their first coordination shell. There is some evidence of solvent-separated Na–Cl ion pairs from the peak in the Na–Cl pair distribution function at 5 Å. However, this disappears with increasing temperature.
3.3 Dynamics of Na–Cl association
Since conductance measurements have been used to measure ion association in NaCl brines, it is of interest to determine the mobility of ions and clusters as a function of pressure, temperature and concentration. In the dilute limit, the limiting conductance (Λ0) of a monovalent electrolyte is
where e is the electronic charge, F is the Faraday constant, kB is Boltzmann's constant and D+,- are the diffusion constants for the individual ions. The self-diffusion constant (D) of a species can be determined from the Einstein relation
The simulations show that with increasing temperature there is increased association of Na+ and Cl- at the expense of solvation by water molecules. This is expected in terms of the decreased dielectric constant of water. The effect of pressure up to 15 kbar is not enough to overcome the effect of temperature to 625°C. Strong association between Na+ and Cl- ions persist to extreme (625°C, 15 kbar) conditions. There is little fundamental difference between the structures of the NaCl–water mixtures at 1 kbar, 325°C and 15 kbar, 625°C. The molecular picture of fully dissociated ions in high pressure NaCl–water mixtures that has been put forward by previous authors, therefore, is not correct. However, the large clusters have short lifetimes and the residence time of a Na or Cl ion in a polyatomic cluster is comparable to the time between atomic collisions. Because of this dynamic nature of the ion association, the transition from a conventional electrolyte to a hydrous molten salt is continuous with pressure and temperature.
MDC's studentship was supported by NERC and a CASE award from Rio Tinto Technology. Computer time on the T3E at CSAR Manchester was provided by NERC grant GR9/03550 (administered by John Brodholt and David Price, UCL). Computer time was also provided by the Laboratory for Advanced Computation in the Mathematical Sciences, University of Bristol.
- Roedder E: Fluid inclusions as samples of ore fluids. Geochemistry of hydrothermal ore deposits. Edited by: HL Barnes. 1967, Holt, Rinehart, and Winston, Inc., New York, 515-574.Google Scholar
- Quist AS, Marshall WL: Electrical conductances of aqueous sodium chloride solutions from 0 to 800°C at at pressures up to 4000 bars. J Phys Chem. 1968, 72: 684-703. 10.1021/j100848a050.View ArticleGoogle Scholar
- Franz G: The brucite-periclase equilibrium at reduced H2O activities: some information about the system H2O–NaCl. Am J Sci. 1982, 282: 1325-1339.View ArticleGoogle Scholar
- Aranovich LYa, Newton RC: H2O activity in concentrated NaCl solutions at high pressures and temperatures measured by the brucite-periclase equilibrium. Contrib Mineral Petrol. 1996, 125: 200-212. 10.1007/s004100050216.View ArticleGoogle Scholar
- Pitzer KS, Simonson JM: Thermodynamics of multi-component, miscible ionic systems: theory and equations. J Phys Chem. 1986, 90: 3005-3009. 10.1021/j100404a042.View ArticleGoogle Scholar
- Pitzer KS, Schrieber DR: The restricted primitive model of ionic fluids. Properties of the vapor and the critical region. Mol Phys. 1987, 60: 1067-1078. 10.1080/00268978700100711.View ArticleGoogle Scholar
- Gillan MJ: Liquid-vapour equilibrium in the restricted primitive model for ionic liquids. Mol Phys. 1983, 49: 421-442. 10.1080/00268978300101271.View ArticleGoogle Scholar
- Oelkers EH, Helgeson HC: Calculation of dissociation constants and the relative stabilities of polynuclear clusters of 1 : 1 electrolytes in hydrothermal solutions at supercritical pressures and temperatures. Geochim Cosmochim Acta. 1993, 57: 2673-2697. 10.1016/0016-7037(93)90383-8.View ArticleGoogle Scholar
- Oelkers EH, Helgeson HC: Multiple ion association in supercritical aqueous solutions of single electrolytes. Science. 1993, 261: 888-891.View ArticleGoogle Scholar
- Koneshan S, Rasaiah JC: Computer simulation studies of aqueous sodium chloride solutions at 298 K and 683 K. J Chem Phys. 2000, 113: 8125-8137. 10.1063/1.1314341.View ArticleGoogle Scholar
- Berendsen HJC, Grigera JR, Straatsma TP: The missing term in effective pair potentials. J Phys Chem. 1987, 91: 6269-6271. 10.1021/j100308a038.View ArticleGoogle Scholar
- Guissani Y, Guillot B: A computer simulation study of the liquid-vapour coexistence curve of water. J Chem Phys. 1993, 98: 8221-8235. 10.1063/1.464527.View ArticleGoogle Scholar
- Haar L, Gallagher JS, Kell GS: NBS/NRC Steam Tables. 1984, Hemisphere, Bristol, PAGoogle Scholar
- Smith DE, Dang LX: Computer simulations of NaCl association in polarizable water. J Chem Phys. 1994, 100: 3757-3766. 10.1063/1.466363.View ArticleGoogle Scholar
- Powell DH, Neilson GW, Enderby GW: The structure of Cl- in aqueous solution-an experimental determination of G(C1H)(R) and G(C1O)(R). J Phys: Condens Matter. 1993, 5: 5723-5730. 10.1088/0953-8984/5/32/003.Google Scholar
- Liu W, Sakane S, Wood RH, Doren DJ: The Hydration Free Energy of Aqueous Na+ and Cl- at High Temperatures Predicted by Ab Initiol/Classical Free Energy Perturbation: 973 K with 0.535 g cm-3 and 573 K with 0.725 g cm-3. J Phys Chem A. 2002, 106: 1409-1418. 10.1021/jp013568i.View ArticleGoogle Scholar
- Brodholt JP, Wood BJ: Simulations of the structure and thermodynamic properties of water at high pressures and temperatures. J Geophys Res. 1993, 98: 519-536.View ArticleGoogle Scholar
- Brodholt JP, Wood BJ: Measurements of the PVT properties of water to 25 kbars and 1600°C from synthetic fluid inclusions in corundum. Geochim Cosmochim Acta. 1994, 58: 2143-2148. 10.1016/0016-7037(94)90292-5.View ArticleGoogle Scholar
- Chialvo AA, Cummings PT, Simonson JM, Mesmer RE: Temperature and density effects on the high temperature ionic speciation in dilute Na+/Cl- aqueous solutions. J Chem Phys. 1996, 105: 9248-9257. 10.1063/1.472756.View ArticleGoogle Scholar
- Chialvo AA, Cummings PT, Simonson JM, Mesmer RE: Molecular simulation study of speciation in supercritical aqueous NaCl solutions. J Mol Liq. 1997, 73–4: 361-372. 10.1016/S0167-7322(97)00079-2.View ArticleGoogle Scholar
- Driesner T, Seward TM, Tironi LG: Molecular dynamics simulation study of ionic hydration and ion association in dilute and 1 molal aqueous sodium chloride solutions from ambient to supercritical conditions. Geochim Cosmochim Acta. 1998, 62: 3095-3107. 10.1016/S0016-7037(98)00207-5.View ArticleGoogle Scholar
- Degreve L, da Silva FLB: Detailed microscopic study of 1 M aqueous NaCl solution by computer simulations. J Mol Liq. 2000, 87: 217-232. 10.1016/S0167-7322(00)00122-7.View ArticleGoogle Scholar
- Brodholt JP: Molecular dynamics simulations of aqueous NaCl solutions at high pressures and temperatures. Chem Geol. 1998, 151: 11-19. 10.1016/S0009-2541(98)00066-7.View ArticleGoogle Scholar
- Archer DG: Thermodynamic properties of the NaCl + H2O system: 2. Thermodynamic properties of NaCl(aq), NaCl-2H2O(cr), and phase equilibria. J Phys Chem Ref Data. 1992, 21: 793-829.View ArticleGoogle Scholar
- Anderko A, Pitzer KS: Equation-of-state representation of phase equilibria and volumetric properties of the system NaCl–H2O above 573 K. Geochim Cosmochim Acta. 1993, 57: 1657-1680. 10.1016/0016-7037(93)90105-6.View ArticleGoogle Scholar
- Allen MP, Tildesley DJ: Computer Simulation of Liquids. 1987, Clarendon Press, OxfordGoogle Scholar
- Refson K: Moldy : a portable molecular dynamics simulation program for serial and parallel computers. Comput Phys Commun. 2000, 126: 310-329. 10.1016/S0010-4655(99)00496-8.View ArticleGoogle Scholar
- Parrinello M, Rahman A: Polymorphic transitions in single crystals: a new molecular dynamics method. J Appl Phys. 1981, 52: 7182-7190. 10.1063/1.328693.View ArticleGoogle Scholar
- Soper AK, Bruni F, Ricci MA: Site-site pair correlation functions of water from 25 to 400°C: Analysis of new and old diffraction data. J Chem Phys. 1997, 106: 247-254. 10.1063/1.473030.View ArticleGoogle Scholar
- Zhu S, Robinson GW: Molecular-dynamics computer simulation of an aqueous NaCl solution: Structure. J Chem Phys. 1992, 97: 4336-4348. 10.1063/1.463903.View ArticleGoogle Scholar
- Oelkers EH, Helgeson HC: Triple-ion anions and polynuclear complexing in supercritical electrolyte solutions. Geochim Cosmochim Acta. 1990, 54: 727-738. 10.1016/0016-7037(90)90368-U.View ArticleGoogle Scholar
- Oelkers EH, Helgeson EHHC: Calculation of activity coefficients and degrees of formation of neutral ion pairs in supercritical electrolyte solutions. Geochim Cosmochim Acta. 1991, 55: 1235-1251. 10.1016/0016-7037(91)90303-M.View ArticleGoogle Scholar
- Lyubartsev AP, Laaksonen A: Concentration effects in aqueous NaCl solutions. A molecular dynamics simulation. J Phys Chem. 1996, 100: 16410-16418. 10.1021/jp961317h.View ArticleGoogle Scholar
- Sharygin AV, Wood RH, Zimmerman GH, Balashov VN: Multiple ion association versus redissociation in aqueous NaCl and KC1 at high temperatures. J Phys Chem B. 2002, 106: 7121-7134. 10.1021/jp013647u.View ArticleGoogle Scholar
- Atkins PW: Physical Chemistry. Edited by: WH Freeman. 1978, San Francisco, CA, 1032-Google Scholar