Validated molecular-level description of NaCl dissociation, ion clustering, and solvent response in supercritical aqueous fluids across the salinity range

A June 2026 molecular-dynamics preprint (Ivanov and Alexandrovich, arXiv:2606.23349) reported the first systematic values of the degree of dissociation and the ideal dissociation constant of NaCl in supercritical H2O-NaCl fluids, from moderate salinity (~1 mol/kg) to $x_{\mathrm{NaCl}} = 0.333$ (27.8 mol/kg) at 673-1273 K and 0.1-2 GPa. It found that dissociation collapses at high salinity, so the abundance of structurally free chloride is non-monotonic in salinity and peaks near $x_{\mathrm{NaCl}} = 0.06$-$0.10$ - with direct consequences for chloride-assisted metal transport in deep crustal fluids. These values rest on a single classical force-field combination (Smith-Dang ions with SPC/fw water), and their comparison with dilute-solution conductance constants shows unresolved systematic deviations. The open problem is to establish a validated molecular-level description of NaCl speciation in supercritical aqueous fluids across the composition range up to the miscibility limit with molten NaCl: first-principles or experimental benchmarking of the concentration-, temperature-, and pressure-dependent dissociation; the size-distribution and lifetime statistics of neutral and charged multi-ion clusters, including the approach to halite crystallization; and the hydrogen-bond-network and dielectric response of water in concentrated brines that drives ion association.

Background

Water above its critical point (647.1 K, 22.1 MPa) is a supercritical fluid: a single phase whose solvent properties vary strongly with temperature and, especially, density. In the deep continental crust and upper mantle, aqueous fluids carry dissolved NaCl at concentrations from about 1 mol/kg up to tens of mol/kg, at temperatures of 673-1673 K and pressures of 0.1-2 GPa, and their thermodynamic and transport behavior controls fluid-rock interaction and the mobility of ore-forming metals. A central molecular question is how much dissolved salt is actually present as free, structurally available ions rather than as associated species. In electrolyte chemistry, a cation and an anion form an ion pair when they remain in contact (a contact ion pair) instead of being separated by solvent; the degree of dissociation, usually denoted $\alpha$, is the fraction of the electrolyte present as separated ions. Because the dielectric constant of water - its ability to screen electrostatic attraction - falls steeply as supercritical water loses density, even nominally strong electrolytes such as NaCl associate substantially there, as classical conductance experiments established for dilute solutions: Quist and Marshall (1968) measured electrical conductances of aqueous NaCl from 0 to 800 °C and to 4 kbar, and Ho, Palmer and Mesmer (1994) extended conductivity measurements to 600 °C and 300 MPa; association constants extracted from such data, however, exist only for about 0.001-0.1 mol/kg.

Molecular dynamics (MD) - numerical integration of the equations of motion of the molecular system - resolves the association state directly: the Na-Cl radial distribution function (the probability density of Na-Cl distances relative to a uniform fluid) shows a contact-pair peak near 2.75 Å separated by a minimum from the solvent-separated region, and an ion can be classified as dissociated when no ion of opposite charge lies within that first minimum (the source study uses the cutoff $q = 3.72$ Å). Classical MD of concentrated brines from ambient to supercritical conditions concluded strong association but reported no dissociation constants (Sherman and Collings 2002, 0.5-16 mol/kg). Ab initio MD (AIMD), in which interatomic forces come from density-functional theory rather than fitted classical potentials, reached similar qualitative conclusions for deep fluids at 0.5-4.5 GPa, 20-800 °C, and 1-14 mol/kg (Fowler and Sherman 2020) and for dilute salt solutions under pressure (Zhang et al. 2020), but likewise without quantitative concentration-dependent dissociation values; AIMD of NaCl dissociation in water exists only near ambient conditions (Timko, Bucher and Kuyucak 2010). Thermodynamic treatments - equilibrium models of ion association and polynuclear (multi-ion) clusters at supercritical conditions (Oelkers and Helgeson 1993) and equations of state for the H2O-NaCl system spanning the full composition range (Duan et al. 2025) - rest on activity-coefficient formalisms and PVTx data rather than molecular speciation statistics, and existing thermodynamic models of multicomponent supercritical fluids commonly treat the dissociation degree as a function of temperature and pressure but not of electrolyte concentration.

Against this backdrop, Ivanov and Alexandrovich (2026) reported the first systematic classical-MD computation of $\alpha$ and of the ideal dissociation constant $K_d = \alpha^2 x_{\mathrm{NaCl}}/(1-\alpha)$ (stoichiometric NaCl mole fraction $x_{\mathrm{NaCl}}$; all activity coefficients set to unity, so $K_d$ is a concentration-dependent ideal quantity, not a thermodynamic standard-state constant) from about 1 mol/kg ($x_{\mathrm{NaCl}} = 0.018$) to 27.8 mol/kg ($x_{\mathrm{NaCl}} = 0.333$) at 673-1273 K and 0.1-2 GPa, using the Smith-Dang ion force field (a widely used fixed-charge Na+/Cl- potential with ion-water interactions fitted to ab initio and experimental data) with the flexible simple-point-charge SPC/fw water model. They found that $\alpha$ falls from 0.3-0.7 at ~1 mol/kg to below 0.01 at the highest salinity, so that the mole fraction of structurally dissociated chloride is a non-monotonic function of salinity with a maximum near $x_{\mathrm{NaCl}} = 0.06$-0.10: in supercritical brines, unlike ordinary liquids, more salt does not mean more available chloride. Since free chloride ions form the ligand complexes that carry metals in crustal fluids, this bears directly on ore-forming transport. The authors themselves state what remains unresolved: the numerical values rest on one non-polarizable force-field combination and 'can be the subject of more precise numerical or experimental studies'; the multi-ion clusters (neutral and charged aggregates of three or more ions, whose numbers their order-of-magnitude estimates show rival the free ions in low-density fluids) were not analyzed in detail; the concentration range ends at $x_{\mathrm{NaCl}} = 0.333$ (with only two exploratory points beyond, at a single temperature and pressure), and the authors name concentrations leading to the formation of macroscopic crystals among the directions for further study - the paper noting that molten NaCl can be mixed with water in arbitrary proportions; and the water side - the hydrogen-bond network and dielectric response of concentrated brines - was not characterized.

Problem Statement

Establish a validated, quantitative molecular-level description of the state of NaCl in supercritical aqueous fluids - the H2O-NaCl system above water's critical point - as a joint function of the stoichiometric NaCl mole fraction $x_{\mathrm{NaCl}}$, temperature, and pressure, over the range 673-1273 K and 0.1-2 GPa, from moderate salinity (~1 mol/kg) up toward the limit of complete miscibility with molten NaCl ($x_{\mathrm{NaCl}} \to 1$, including the halite-saturation boundary), replacing the currently unvalidated single-force-field picture. Within this one objective, three results are jointly required:

(1) Benchmarked dissociation values: values of the degree of dissociation $\alpha(x,T,P)$ and the associated ideal dissociation constant $K_d = \alpha^2 x_{\mathrm{NaCl}}/(1-\alpha)$ at, or explicitly validated against, a first-principles (ab initio) level of molecular dynamics and/or against experiment - under a stated geometric association criterion (an ion is dissociated when no ion of opposite charge lies within the first minimum of the Na-Cl radial distribution function; the source study fixes this boundary at $q = 3.72$ Å, the Na-Cl distance of closest approach proposed by Helgeson et al. (1981) for HKF-type activity-coefficient calculations) - sufficient to confirm or correct the classical-MD values obtained with the Smith-Dang ion potential and SPC/fw water, whose systematic deviation from dilute-solution conductance constants grows with water density.

(2) Cluster statistics: the size distributions of neutral and charged multi-ion aggregates over the same composition-temperature-pressure range, their lifetimes, and their dependence on temperature, pressure, and composition, including the transition toward concentrations at which macroscopic halite crystallization becomes possible.

(3) Solvent-side mechanism: characterization of the hydrogen-bond network and dielectric response of water in moderately and highly concentrated supercritical brines, establishing how the degradation of water's screening ability controls ion association and the non-monotonic availability of structurally free chloride.

These three results jointly support the single overarching objective - a validated account of how much structurally available chloride supercritical brines actually contain as a function of salinity, temperature, and pressure. The extensions they frame separately (analogous studies of other electrolytes such as H2O-CaCl2, incorporation of the data into equations of state, and geological testing of the predicted optimal salinity $x_{\mathrm{NaCl}} \approx 0.06$-0.10) are follow-on applications outside this problem's boundary.

References

1. Ivanov, M.V., Alexandrovich, O.V., 2026. Dissociation of NaCl in supercritical aqueous fluids of moderate and high concentrations: A molecular dynamics study. arXiv:2606.23349, doi:10.48550/arXiv.2606.23349, https://arxiv.org/abs/2606.23349
2. Duan, Z., Cheng, N., Zhang, Z., Chou, I.-M., Sun, H., 2025. Molecular dynamics simulation and equation of state of the NaCl-H2O system from 573 to 1573 K, 1 to 30 kbar, and 0 to 1 m fraction of NaCl. Geochimica et Cosmochimica Acta 412, 127-141, doi:10.1016/j.gca.2025.10.022, https://doi.org/10.1016/j.gca.2025.10.022
3. Fowler, S.J., Sherman, D.M., 2020. The nature of NaCl-H2O deep fluids from ab initio molecular dynamics at 0.5-4.5 GPa, 20-800 C, and 1-14 m NaCl. Geochimica et Cosmochimica Acta 277, 243-264, doi:10.1016/j.gca.2020.03.031, https://doi.org/10.1016/j.gca.2020.03.031
4. Sherman, D.M., Collings, M.D., 2002. Ion association in concentrated NaCl brines from ambient to supercritical conditions: results from classical molecular dynamics simulations. Geochemical Transactions 3, 102, doi:10.1186/1467-4866-3-102, https://doi.org/10.1186/1467-4866-3-102
5. Zhang, C., Giberti, F., Sevgen, E., de Pablo, J.J., Gygi, F., Galli, G., 2020. Dissociation of salts in water under pressure. Nature Communications 11, 3037, doi:10.1038/s41467-020-16704-9, https://doi.org/10.1038/s41467-020-16704-9
6. Quist, A.S., Marshall, W.L., 1968. Electrical conductances of aqueous sodium chloride solutions from 0 to 800° and at pressures to 4000 bars. The Journal of Physical Chemistry 72, 684-703, doi:10.1021/j100848a050, https://doi.org/10.1021/j100848a050
7. Ho, P.C., Palmer, D.A., Mesmer, R.E., 1994. Electrical conductivity measurements of aqueous sodium chloride solutions to 600 °C and 300 MPa. Journal of Solution Chemistry 23, 997-1018, doi:10.1007/BF00974100, https://doi.org/10.1007/BF00974100
8. Elbers, M., Schmidt, C., Sternemann, C., Sahle, C.J., Jahn, S., Albers, C., Sakrowski, R., Gretarsson, H., Sundermann, M., Tolan, M., Wilke, M., 2021. Ion association in hydrothermal aqueous NaCl solutions: implications for the microscopic structure of supercritical water. Physical Chemistry Chemical Physics 23, 14845-14856, doi:10.1039/D1CP01490K, https://doi.org/10.1039/D1CP01490K
9. Oelkers, E.H., Helgeson, H.C., 1993. Calculation of dissociation constants and the relative stabilities of polynuclear clusters of 1:1 electrolytes in hydrothermal solutions at supercritical pressures and temperatures. Geochimica et Cosmochimica Acta 57, 2673-2697, doi:10.1016/0016-7037(93)90383-8, https://doi.org/10.1016/0016-7037(93)90383-8
10. Timko, J., Bucher, D., Kuyucak, S., 2010. Dissociation of NaCl in water from ab initio molecular dynamics simulations. The Journal of Chemical Physics 132, 114510, doi:10.1063/1.3360310, https://doi.org/10.1063/1.3360310