Multi-scale modelling of uranyl chloride solutions

Classical molecular dynamics simulations with explicit polarization have been successfully used to determine the structural and thermodynamic properties of binary aqueous solutions of uranyl chloride (UO2Cl2). Concentrated aqueous solutions of uranyl chloride have been studied to determine the hydration properties and the ion-ion interactions. The bond distances and the coordination number of the hydrated uranyl are in good agreement with available experimental data. Two stable positions of chloride in the second hydration shell of uranyl have been identified. The UO2(2+)-Cl(-) association constants have also been calculated using a multi-scale approach. First, the ion-ion potential averaged over the solvent configurations at infinite dilution (McMillan-Mayer potential) was calculated to establish the dissociation/association processes of UO2 (2+)-Cl(-) ion pairs in aqueous solution. Then, the association constant was calculated from this potential. The value we obtained for the association constant is in good agreement with the experimental result (KUO2Cl(+) = 1.48 l mol(-1)), but the resulting activity coefficient appears to be too low at molar concentration.