We have studied the solvation of uranyl, UO22+, and the reduced species UO(OH)(2+) and U(OH)(2)(2+) 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 similar to100 kcal/mol or more. With our best model M3, the calculated hydration free energy Delta(h)Gdegrees of UO22+ 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 < Delta(h)Gdegrees < -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 UO22+ 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.