Role of Solvation in the Reduction of the Uranyl Dication by Water: A Density Functional Study

Abstract

We have studied the solvation of uranyl, UO₂²⁺, and the reduced species UO(OH)²⁺ and U(OH)₂²⁺ systematically using three levels of approximation: direct application of a continuum model (M1); explicit quantum-chemical treatment of the first hydration sphere (M2); a combined quantum-chemical/continuum model approach (M3). We have optimized complexes with varying numbers of aquo ligands (n = 4−6) and compared their free energies of solvation. Models M1 and M2 have been found to recover the solvation energy only partially, underestimating it by ∼100 kcal/mol or more. With our best model M3, the calculated hydration free energy Δ_hG° of UO₂²⁺ is about −420 kcal/mol, which shifts to about −370 kcal/mol when corrected for the expected error of the model. This value agrees well with the experimentally determined interval, −437 kcal/mol < Δ_hG° < −318 kcal/mol. Complexes with 5 and 6 aquo ligands have been found to be about equally favored with models M2 and M3. The same solvation models have been applied to a two-step reduction of UO₂²⁺ by water, previously theoretically studied in the gas phase. Our results show that the solvation contribution to the reaction free energy, about 60 kcal/mol, dominates the endoergicity of the reduction.