Chemical Potentials, Activity Coefficients, and Solubility in Aqueous NaCl Solutions: Prediction by Polarizable Force Fields
Bibliographic record
Abstract
We describe a computationally efficient molecular simulation methodology for calculating the concentration dependence of the chemical potentials of both solute and solvent in aqueous electrolyte solutions, based on simulations of the salt chemical potential alone. We use our approach to study the predictions for aqueous NaCl solutions at ambient conditions of these properties by the recently developed polarizable force fields (FFs) AH/BK3 of Kiss and Baranyai (J. Chem. Phys. 2013, 138, 204507) and AH/SWM4-DP of Lamoureux and Roux (J. Phys. Chem. B 2006, 110, 3308 - 3322) and by the nonpolarizable JC FF of Joung and Cheatham tailored to SPC/E water (J. Phys. Chem. B 2008, 112, 9020 - 9041). We also consider their predictions of the concentration dependence of the electrolyte activity coefficient, the crystalline solid chemical potential, the electrolyte solubility, and the solution specific volume. We first highlight the disagreement in the literature concerning calculations of solubility by means of molecular simulation in the case of the JC FF and provide strong evidence of the correctness of our methodology based on recent independently obtained results for this important test case. We then compare the predictions of the three FFs with each other and with experiment and draw conclusions concerning their relative merits, with particular emphasis on the salt chemical potential and activity coefficient vs concentration curves and their derivatives. The latter curves have only previously been available from Kirkwood-Buff integrals, which require approximate numerical integrations over system pair correlation functions at each concentration. Unlike the case of the other FFs, the AH/BK3 curves are nearly parallel to the corresponding experimental curves at moderate and higher concentrations. This leads to an excellent prediction of the water chemical potential via the Gibbs-Duhem equation and enables the activity coefficient curve to be brought into excellent agreement with experiment by incorporating an appropriate value of the standard state chemical potential in the Henry Law convention.
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame distilled prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. Learned from the 10,348 direct Codex labels and 10,348 direct Gemma labels. Candidate is the union of thresholded teacher heads; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels or direct frontier model labels.
Codex and Gemma teacher scores by category
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.001 | 0.000 |
| Meta-epidemiology (narrow) | 0.000 | 0.000 |
| Meta-epidemiology (broad) | 0.000 | 0.000 |
| Bibliometrics | 0.000 | 0.000 |
| Science and technology studies | 0.000 | 0.000 |
| Scholarly communication | 0.000 | 0.000 |
| Open science | 0.000 | 0.000 |
| Research integrity | 0.000 | 0.000 |
| Insufficient payload (model declined to judge) | 0.000 | 0.000 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one teacher head, not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".