Extensive atomistic molecular dynamics simulation results suggest that probability of counter ion occupying first solvation shell of uranyl ion depends both on its nature and concentration. In general, uranyl-counter ion complexes with pentagonal bi-pyramidal structure with five ligands in the equatorial plane of the linear UO22+ ion are observed. In case of nitrate ion, pure aqua complexes at lower concentrations and mixed mono-nitro aqua complexes are observed at higher concentrations, whereas in case of sulphate and carbonate ions, no pure aqua complexes are observed. The NO3− and SO42− ions act as unidentate, but CO32− acts both as uni- and bi-dentate ligand. In addition, polynuclear uranium complexes with bridging SO42− and CO32− ligands are observed. Relative strength of binding of counter ions with uranyl ion from PMF calculations follows the order CO32− > SO42− > NO3− with the contact pair free energies of about −29.0, −14.0 and −1.1 kcal/mol respectively.
Read full abstract