Theoretical studies of UO2(H2O)n(2+), NpO2(H2O)n(+), and PuO2(H2O)n(2+) complexes (n=4-6) in aqueous solution and gas...

Extensive ab initio calculations both in gas phase and solution have been carried out to study the equilibrium structure, vibrational frequencies, and bonding characteristics of various actinyl (UO2(2+), NpO2(+), and PuO2(2+)) and their hydrated forms, AnO2(H2O)n(z+) (n=4, 5, and 6). Bulk solvent effects were studied using a continuum method. The geometries were fully optimized at the coupled-cluster singles + doubles (CCSD), density-functional theory (DFT), and Møller-Plesset (MP2) level of theories. In addition vibrational frequencies have been obtained at the CCSD as well as MP2/DFT levels. The results show that both the short-range and long-range solvent effects are important. The combined discrete-continuum model, in which the ionic solute and the solvent molecules in the first and second solvation shells are treated quantum mechanically while the solvent is simulated by a continuum model, can predict accurately the bonding characteristics. Moreover, our values of solvation free energies suggest that five- and six-coordinations are equally preferred for UO2(2+), and five-coordinated species are preferred for NpO2(+) and PuO2(2+). On the basis of combined quantum-chemical and continuum treatments of the hydrated complexes, we are able to determine the optimal cavity radii for the solvation models. The coupled-cluster computations with large basis sets were employed for the vibrational spectra and equilibrium geometries both of which compare quite favorably with experiment. Our most accurate computations reveal that both five- and six-coordination complexes are important for these species.