Using Reaktoro for mineral and gas solubility calculations with the Extended UNIQUAC model

The extended UNIQUAC model is a thermodynamic model able to estimate thermodynamic properties of aqueous electrolyte solutions under a wide range of temperature, pressure, and composition conditions. Thermodynamic properties include species activity coefficients, excess molar Gibbs energy, excess molar enthalpy, excess molar heat capacity. These properties are important for aqueous speciation calculations, mineral and gas solubility computations, chemical kinetic modeling of mineral dissolution and precipitation, and in reactive transport simulations considering chemically complex aqueous electrolyte solutions. In this paper we present a brief literature review on the extended UNIQUAC model, we report on its im-plementation in C++ in the Reaktoro framework for modeling chemically reactive systems, and we show its use from Python for computing mineral and gas solubilities in aqueous solutions at a wide range of temperature, pressure, and salinity conditions. We validated the calculations against experimental data and against results obtained through the software ScaleCERE implementing the extended UNIQUAC model. Our conclusion is that the extended UNIQUAC model has been successfully implemented into the Reaktoro framework, thereby providing a suitable activity model for geochemical and reactive transport modeling.