Simulations of dissociation constants in low pressure supercritical water

This article reports molecular dynamics simulations of the dissociation of hydrochloric acid and sodium hydroxide in water from ambient to supercritical temperatures at a fixed pressure of 250 atm. Corrosion of reaction vessels is known to be a serious problem of supercritical water, and acid/base dissociation can be a significant contributing factor to this. The SPC/e model was used in conjunction with solute models determined from density functional calculations and OPLSAA Lennard–Jones parameters. Radial distribution functions were calculated, and these show a significant increase in solute–solvent ordering upon forming the product ions at all temperatures. For both dissociations, rapidly decreasing entropy of reaction was found to be the controlling thermodynamic factor, and this is thought to arise due to the ions produced from dissociation maintaining a relatively high density and ordered solvation shell compared to the reactants. The change in entropy of reaction reaches a minimum at the critical temperature. The values of pKa and pKb were calculated and both increased with temperature, in qualitative agreement with other work, until a maximum value at 748 K, after which there was a slight decrease