We report on molecular dynamics (MD) simulations for predicting the density and isothermal compressibility of an H2O-NaCl fluid as a function of temperature (673-2000 K), pressure (0.2-2.0 GPa), and salt concentration (0.0-21.9 wt%). The atomistic behavior was analyzed via the hydration number of ions and number of ion pairs. Hydration numbers of Na+ and Cl- increased with increasing pressure and decreasing temperature. Conversely, the fraction of Na-Cl ion pairs increased with decreasing pressure and increasing temperature. This hydration and association behavior is consistent with the low dielectric constant of H2O under these conditions. The presence of polynuclear clusters of Na-Cl was confirmed at high temperatures, low pressures, and high salt concentrations. We propose a purely empirical equation of state (EoS) for H2O-NaCl fluids under high temperatures and pressures that should be useful for estimating the fluid distribution in Earth's crust and upper mantle in relation to effects on earthquakes and volcanic eruptions.
ASJC Scopus subject areas