Dissociation of NaCl in supercritical aqueous fluids of moderate and high concentrations: A molecular dynamics study

Mikhail V. Ivanov; Olga V. Alexandrovich

arxiv: 2606.23349 · v1 · pith:DIU3P3SUnew · submitted 2026-06-22 · ❄️ cond-mat.soft · cond-mat.other· physics.chem-ph· physics.geo-ph

Dissociation of NaCl in supercritical aqueous fluids of moderate and high concentrations: A molecular dynamics study

Mikhail V. Ivanov , Olga V. Alexandrovich This is my paper

Pith reviewed 2026-06-26 06:32 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.otherphysics.chem-phphysics.geo-ph

keywords molecular dynamicssupercritical fluidsNaCl dissociationion associationaqueous electrolyteshigh-temperature fluidschloride availability

0 comments

The pith

Mole fraction of dissociated Na+ and Cl- ions in supercritical NaCl solutions peaks at 6-10% concentration.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Classical molecular dynamics simulations follow NaCl association and dissociation in supercritical water from low to very high salt concentrations. The degree of dissociation drops as concentration rises because ions form pairs and larger clusters, yet the absolute number of free ions first increases and then decreases. This produces a maximum in the mole fraction of structurally dissociated ions near xNaCl = 0.06-0.10. The finding shows that raising salinity does not automatically increase the supply of free chloride ions available for reactions or transport.

Core claim

The mole fraction of structurally dissociated Na+ and Cl- ions is a non-monotonic function of the stoichiometric NaCl concentration and typically reaches a maximum at xNaCl = 0.06-0.10. This occurs because the falling degree of dissociation competes with the rising total salt amount, so that most ions remain in contact pairs and multi-ion clusters at the highest concentrations examined (up to xNaCl = 0.333).

What carries the argument

Degree of dissociation a and ideal dissociation constant Kd obtained directly from counting free ions versus contact pairs and clusters in classical molecular dynamics trajectories of H2O-NaCl mixtures.

If this is right

At moderate salinity near 1 mol/kg the degree of dissociation remains 0.3-0.7, while at xNaCl = 0.333 it falls toward zero.
Additional fixed-density runs separate the effects of temperature and density up to 1673 K.
The concentration curves supply molecular-level constraints for thermodynamic models of concentrated supercritical electrolytes.
Chloride availability in high-temperature aqueous fluids does not increase steadily with added salt.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Geochemical models of hydrothermal or volcanic fluids may need to incorporate a non-monotonic ion-availability curve rather than assuming steady increase with salinity.
Optimal salt concentrations for maximum free chloride in supercritical water could be tested in targeted conductivity or spectroscopic experiments.
The location of the maximum may shift with pressure or with the choice of force field, providing a clear target for sensitivity checks.

Load-bearing premise

The chosen classical force fields correctly reproduce the balance of ion pairing, dissociation, and clustering across the entire concentration range from roughly 1 to 28 mol/kg.

What would settle it

An experimental determination of the fraction of free ions or of electrical conductivity versus concentration at one of the studied supercritical PT points that shows either monotonic behavior or a maximum at a concentration outside 0.06-0.10.

read the original abstract

We report classical molecular dynamics simulations of NaCl association and dissociation in supercritical aqueous fluids over a wide range of salt concentrations, from moderate salinity to highly concentrated H2O-NaCl mixtures attainable at high temperatures. The degree of dissociation a and the corresponding ideal dissociation constant Kd, derived directly from a, were calculated as functions of the stoichiometric NaCl mole fraction at selected pressure-temperature (PT) conditions from 673.15 to 1273.15 K and from 0.1 to 2 GPa. At moderate salinity corresponding to a molality of approximately 1 mol/kg, NaCl remains largely dissociated a = 0.3-0.7 depending on pressure and temperature). In contrast, when the mole fraction of NaCl increases up to xNaCl = 0.333 (27.8 mol/kg), the degree of dissociation tends towards zero, and most ions form Na$^+$Cl$^-$ contact pairs and multi-ion clusters. As a result of these competing trends, the mole fraction of structurally dissociated Na$^+$ and Cl$^-$ ions is a non-monotonic function of the stoichiometric NaCl concentration and typically reaches a maximum at xNaCl = 0.06-0.10. This result shows that increasing salinity does not necessarily increase the abundance of structurally available chloride ions in supercritical aqueous fluids. Additional fixed density simulations at 1 and 7 mol/kg extend the analysis up to 1673.15 K and separate the effects of temperature and density on the associate/dissociate state of the ions. The obtained concentration dependences provide molecular-level constraints for thermodynamic descriptions of concentrated supercritical electrolytes and for evaluating chloride availability in high-temperature aqueous fluids.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

MD study finds non-monotonic free-ion mole fraction peaking near xNaCl=0.07 in supercritical NaCl solutions, but leaves force-field accuracy at high concentration untested.

read the letter

The main takeaway is that these simulations show the mole fraction of free Na+ and Cl- ions peaks at xNaCl around 0.06-0.10 and then falls as contact pairs and clusters dominate at higher salinity up to 0.333. This supplies a concrete concentration dependence that thermodynamic models of supercritical electrolytes can use as a check.

The work does a clean job of running trajectories across a range of PT conditions, counting dissociated ions directly from configurations, and adding fixed-density runs at 1 and 7 mol/kg to separate temperature and density effects. Extending the concentration range beyond most prior MD studies is the actual advance here.

The soft spot is the complete absence of force-field details, validation, system sizes, or error estimates. The non-monotonic result rests on the potentials correctly balancing ion-water and ion-ion energies across the full concentration range at supercritical conditions. Standard classical force fields are known to be sensitive to cutoffs and missing polarization in this regime, so an error in how association changes with concentration would move or erase the reported maximum. The stress-test concern holds up on the abstract alone.

This paper is for modelers in high-temperature geochemistry and chemical engineering who need molecular-level constraints on chloride speciation. A reader building or testing equations of state for concentrated brines would find the trends useful even if they rerun the systems with different potentials.

It deserves peer review. The data are new and the counting method is direct, so referees can request the missing methods information and any cross-checks against experiment or other force fields.

Referee Report

3 major / 2 minor

Summary. The manuscript reports classical molecular dynamics simulations of NaCl in supercritical water at 673–1273 K and 0.1–2 GPa, extracting the dissociation degree α directly from ion-pair configurations in the trajectories for stoichiometric mole fractions up to xNaCl = 0.333 (27.8 mol/kg). It finds that α decreases from 0.3–0.7 at ~1 mol/kg to near zero at high concentration, but the product xNaCl·α (mole fraction of structurally dissociated ions) is non-monotonic and peaks at xNaCl ≈ 0.06–0.10. Fixed-density runs at 1 and 7 mol/kg extend the temperature range to 1673 K.

Significance. If the underlying force fields are shown to be reliable, the direct-counting results supply concentration-dependent constraints on ion association that are free of fitted parameters and can be used to test thermodynamic models of supercritical electrolytes and to assess chloride availability in high-temperature fluids.

major comments (3)

[Methods] Methods section: the water model, ion parameters (e.g., Joung-Cheatham or equivalent), system sizes, cutoff schemes, and equilibration/production lengths are not specified. These choices directly control the balance between ion–water solvation and Na+–Cl– contact-pair/cluster formation energies that produce the reported non-monotonic peak; without them the central claim cannot be evaluated.
[Results] Results section (and abstract): no statistical uncertainties, block-averaging errors, or convergence checks on α(x) are reported. At the highest concentrations the cluster statistics become sparse, so the location and height of the maximum at xNaCl = 0.06–0.10 could shift within the stated range once errors are quantified.
[Results] Results section: the manuscript provides no validation of the chosen force fields against experimental dissociation constants, activity coefficients, or prior simulations at supercritical conditions and comparable salinities. Given the known sensitivity of classical pairwise potentials to polarization and long-range electrostatics at high T and high salinity, this omission leaves the non-monotonic trend as an untested prediction of the specific model.

minor comments (2)

[Abstract/Introduction] The definition of the “ideal dissociation constant Kd” derived from α should be stated explicitly (including any activity-coefficient assumptions) when first introduced.
[Figures] Figure captions should indicate the number of independent trajectories or total sampling time used to compute each α value.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive suggestions. We address each major comment below and will revise the manuscript accordingly to improve clarity, reproducibility, and context.

read point-by-point responses

Referee: [Methods] Methods section: the water model, ion parameters (e.g., Joung-Cheatham or equivalent), system sizes, cutoff schemes, and equilibration/production lengths are not specified. These choices directly control the balance between ion–water solvation and Na+–Cl– contact-pair/cluster formation energies that produce the reported non-monotonic peak; without them the central claim cannot be evaluated.

Authors: We agree that the submitted manuscript omitted a dedicated Methods section with these technical details. In the revision we will insert a complete Methods section specifying the water model, ion parameters, system sizes, cutoff and long-range electrostatic schemes, and the durations of equilibration and production runs for each state point. revision: yes
Referee: [Results] Results section (and abstract): no statistical uncertainties, block-averaging errors, or convergence checks on α(x) are reported. At the highest concentrations the cluster statistics become sparse, so the location and height of the maximum at xNaCl = 0.06–0.10 could shift within the stated range once errors are quantified.

Authors: We accept that error estimates are required. We will re-analyze the trajectories with block averaging to obtain statistical uncertainties on α(x) at each concentration and will add these uncertainties to the figures and text. This will allow readers to judge whether the reported peak location remains robust. revision: yes
Referee: [Results] Results section: the manuscript provides no validation of the chosen force fields against experimental dissociation constants, activity coefficients, or prior simulations at supercritical conditions and comparable salinities. Given the known sensitivity of classical pairwise potentials to polarization and long-range electrostatics at high T and high salinity, this omission leaves the non-monotonic trend as an untested prediction of the specific model.

Authors: We acknowledge the value of explicit validation. The revised manuscript will include a new subsection that (i) states the precise force-field combination employed, (ii) cites earlier studies that tested the same or closely related models against available supercritical dissociation data, and (iii) notes the known limitations of non-polarizable potentials at these conditions. A comprehensive new validation campaign against the full set of experimental activity coefficients and dissociation constants lies outside the scope of the present work but can be flagged as future research. revision: partial

Circularity Check

0 steps flagged

No significant circularity: dissociation fractions counted directly from MD trajectories

full rationale

The paper computes the structural dissociation fraction a by direct enumeration of ion configurations (contact pairs vs. free ions) in the MD trajectories at each concentration. The reported mole fraction of dissociated ions is obtained by the elementary product a(x) * x_NaCl. No parameters are fitted to a subset of the data and then invoked as a prediction of a related quantity; no self-citation chain supplies a uniqueness theorem or ansatz; and the non-monotonic behavior emerges from the competing concentration dependences of a(x) and x itself. The derivation chain is therefore self-contained against the simulation output and does not reduce to any of the enumerated circular patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that the chosen classical force field faithfully represents ion-water and ion-ion interactions under supercritical conditions; no free parameters are fitted to produce the reported trends, and no new entities are postulated.

axioms (1)

domain assumption Classical molecular dynamics with standard force fields can accurately capture NaCl association/dissociation and clustering in supercritical H2O-NaCl mixtures.
This assumption is invoked as the basis for interpreting all reported values of α and free-ion mole fractions from the trajectories.

pith-pipeline@v0.9.1-grok · 5856 in / 1339 out tokens · 35374 ms · 2026-06-26T06:32:46.693230+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

2 extracted references

[1]

Molecular dynamics simulatio n of aqueous NaCl solutions at high pressures and temperatures. Chem. Geology 151, 11–19. Chialvo A.A., Cummings P.T., Cochran H.D., Simonson J.M., Mesmer R.E., 1995, Na+–Cl− ion pair association in supercritical water. J. Chem. Phys. 103, 9379–9387. Chialvo A.A., Simonson J.M.,

1995
[2]

The IAPWS formulation 1995 for the thermodynamic properties of ordinary water substance for general and scientific use . J. Phys. Chem. Ref. Data 31, 387–

1995

[1] [1]

Molecular dynamics simulatio n of aqueous NaCl solutions at high pressures and temperatures. Chem. Geology 151, 11–19. Chialvo A.A., Cummings P.T., Cochran H.D., Simonson J.M., Mesmer R.E., 1995, Na+–Cl− ion pair association in supercritical water. J. Chem. Phys. 103, 9379–9387. Chialvo A.A., Simonson J.M.,

1995

[2] [2]

The IAPWS formulation 1995 for the thermodynamic properties of ordinary water substance for general and scientific use . J. Phys. Chem. Ref. Data 31, 387–

1995