Interfacial Interaction of Titania Nanoparticles and Ligated Uranyl Species: A Relativistic DFT Investigation.

To understand interfacial behavior of actinides adsorbed onto mineral surfaces and unravel their structure-property relationship, the structures, electronic properties, and energetics of various ligated uranyl species adsorbed onto TiO2 surface nanoparticle clusters (SNCs) were examined using relativistic density functional theory. Rutile (110) and anatase (101) titania surfaces, experimentally known to be stable, were fully optimized. For the former, models studied include clean and water-free Ti27O64H20 (dry), partially hydrated (Ti27O64H20)(H2O)8 (sol) and proton-saturated [(Ti27O64H20)(H2O)8(H)2]2+ (sat), while defect-free and defected anatase SNCs involving more than 38 TiO2 units were considered. The aquouranyl sorption onto rutile SNCs is energetically preferred, with interaction energies of -8.54, -10.36, and -2.39 eV, respectively. Energy decomposition demonstrates that the sorption is dominated by orbital attractive interactions and modified by steric effects. Greater hydrogen-bonding involvement leads to increased orbital interactions (i.e., more negative energy) from dry to sol/sat complexes, while much larger steric interaction in the sat complex significantly reduces the sorption interaction (i.e., more positive energy). For dry SNC, adsorbates were varied from aquo to aquo-carbonato, to carbonato, to hydroxo uranyl species. Longer U-Osurf/U-Ti distances and more positive sorption energies were calculated upon introducing carbonato and hydroxo ligands, indicative of weaker uranyl sorption onto the substrate. This is consistent with experimental observations that the uranyl sorption rate decreases upon raising solution pH value or adding carbon dioxide. Anatase SNCs adsorbing aquouranyl are even more exothermic, because more bonds are formed than in the case of rutile. Moreover, the anatase sorption can be tuned by surface defects as well as its Ti and O stoichiometry. All the aquouranyl-SNC complexes show similar character of molecular orbitals and energetic order although differing in highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gaps and orbital energy levels, but changes can be accomplished by adding carbonato and hydroxo ligands.