Journal
JOURNAL OF PHYSICAL CHEMISTRY B
Volume 116, Issue 18, Pages 5468-5478Publisher
AMER CHEMICAL SOC
DOI: 10.1021/jp301447z
Keywords
-
Categories
Funding
- Natural Sciences and Engineering Research Council of Canada [OGP1041]
- Ontario Research Fund
- SHARCNET grid computing facility
- Grant Agency of the Czech Republic [203/08/0094]
- Ministry of Education, Youth and Sports of the Czech Republic [KONTAKT II LH1020]
- European Community [COST TD0802]
Ask authors/readers for more resources
We extend the osmotic ensemble Monte Carlo (OEMC) molecular simulation method (Moucka et al. J. Phys Chem. B 2011, 115, 7849-7861) for directly calculating the aqueous solubility of electrolytes and for calculating their chemical potentials as functions of concentration to cases involving electrolyte hydrates and mixed electrolytes, including invariant points involving simultaneous precipitation of several solutes. The method utilizes a particular semigrand canonical ensemble, which performs simulations of the solution at a fixed number of solvent molecules, pressure, temperature, and specified overall electrolyte chemical potential. It avoids calculations for the solid phase, incorporating available solid chemical potential data from thermochemical tables, which are based on well-defined reference states, or from other sources. We apply the method to a range of alkali halides in water and to selected examples involving LiCl monohydrate, mixed electrolyte solutions involving water and hydrochloric acid, and invariant points in these solvents. The method uses several existing force-field models from the literature, and the results are compared with experiment. The calculated results agree qualitatively well with the experimental trends and are of reasonable accuracy. The accuracy of the calculated solubility is highly dependent on the solid chemical potential value and also on the force-field model used. Our results indicate that pairwise additive effective force-field models developed for the solution phase are unlikely to also be good models for the corresponding crystalline solid. We find that, in our OEMC simulations, each ionic force-field model is characterized by a limiting value of the total solution chemical potential and a corresponding aqueous concentration. For higher values of the imposed chemical potential, the solid phase in the simulation grows in size without limit.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available