Equatorial and apical solvent shells of the UO22+ ion

Enhanced uranium extraction from seawater: from the viewpoint of kinetics and thermodynamics

Uranium extraction from seawater (UES) is recognized as one of the seven pivotal chemical separations with the potential to revolutionize global paradigms. The forthcoming decade is anticipated to witness a surge in UES, driven by escalating energy demands. The oceanic reservoirs, possessing uranium quantities approximately 1000-fold higher than terrestrial mines, present a more sustainable and environmentally benign alternative. Empirical evidence from historical research indicates that adsorption emerges as the most efficacious process for uranium recovery from seawater, considering operational feasibility, cost-effectiveness, and selectivity. Over the years, scientific exploration has led to the development of a plethora of adsorbents with superior adsorption capacity. It would be efficient to design materials with a deep understanding of the adsorption from the perspective of kinetics and thermodynamics. Here, we summarize recent advancements in UES technology and the contemporary challenges encountered in this domain. Furthermore, we present our perspectives on the future trajectory of UES and finally offer our insights into this subject.

Adsorption of uranium (VI) complexes with polymer-based spherical activated carbon

Adsorption processes with carbon-based adsorbents have received substantial attention as a solution to remove uranium from drinking water. This study investigated uranium adsorption by a polymer-based spherical activated carbon (PBSAC) characterised by a uniformly smooth exterior and an extended surface of internal cavities accessible via mesopores. The static adsorption of uranium was investigated applying varying PBSAC properties and relevant solution chemistry. Spatial time-of-flight secondary ion mass spectrometry (ToF-SIMS) was employed to visualise the distribution of the different uranium species in the PBSAC. The isotherms and thermodynamics calculations revealed monolayer adsorption capacities of 28-667 mg/g and physical adsorption energies of 13-21 kJ/mol. Increasing the surface oxygen content of the PBSAC to 10 % enhanced the adsorption and reduced the equilibrium time to 2 h, while the WHO drinking water guideline of 30 µgU/L could be achieved for an initial concentration of 250 µgU/L. Uranium adsorption with PBSAC was favourable at the pH 6-8. At this pH range, uranyl carbonate complexes (UO2CO3(aq), UO2(CO3)22-, (UO2)2CO3(OH)3-) predominated in the solution, and the ToF-SIMS analysis revealed that the adsorption of these complexes occurred on the surface and inside the PBSAC due to intra-particle diffusion. For the uranyl cations (UO22+, UO2OH+) at pH 2-4, only shallow adsorption in the outermost PBSAC layers was observed. The work demonstrated the effective removal of uranium from contaminated natural water (67 µgU/L) and meeting both German (10 µgU/L) and WHO guideline concentrations. These findings also open opportunities to consider PBSAC in hybrid treatment technologies for uranium removal, for instance, from high-level radioactive waste.

Nature and coordination geometry of geologically relevant aqueous Uranium(VI) complexes up to 400 ºC: A review and new data

The structure of the uranyl aqua ion (UO₂²⁺) and a number of its inorganic complexes (specifically, UO₂Cl⁺, UO₂Cl₂⁰, UO₂SO₄⁰, [Formula: see text] , [Formula: see text] and UO₂OH₄^2-) have been characterised using X-Ray absorption spectroscopy/extended X-Ray absorption fine structure (XAS/EXAFS) at temperatures ranging from 25 to 326 ºC. Results of ab initio molecular dynamics (MD) calculations are also reported for uranyl in chloride and sulfate-bearing fluids from 25 to 400 ºC and 600 bar to 20 kilobar (kb). These results are reported alongside a comprehensive review of prior structural characterisation work with particular focus given to EXAFS works to provide a consistent and up-to-date view of the structure of these complexes under conditions relevant to U mobility in ore-forming systems and around high-grade nuclear waste repositories. Regarding reported EXAFS results, average equatorial coordination was found to decrease in uranyl and its sulfate and chloride complexes as temperature rose - the extent of this decrease differed between species and solution compositions but typically resulted in an equatorial coordination number of ∼3-4 at temperatures above 200 ºC. The [Formula: see text] complex was observed at temperatures from 25 to 247 ºC and exhibited no major structural changes over this temperature range. UO₂OH₄^2- exhibited only minor structural changes over a temperature range from 88 to 326 ºC and was suggested to manifest fivefold coordination with four hydroxyl molecules and one water molecule around its equator. Average coordination values derived from fits of the reported EXAFS data were compared to average coordination values calculated using the experimentally derived thermodynamic data for chloride complexes reported by Dargent et al. (2013) and Migdisov et al. (2018b), and for sulfate complexes reported by Alcorn et al. (2019) and Kalintsev et al. (2019). Sulfate EXAFS data were well described by available thermodynamic data, and chloride EXAFS data were described well by the thermodynamic data of Migdisov et al. (2018b), but not by the data of Dargent et al. (2013). The ab initio molecular dynamics calculations confirmed the trends in equatorial coordination observed with EXAFS and were also able to provide an insight into the effect of pressure in equatorial water coordination - for a given temperature, higher pressures appear to lead to a greater number of equatorially bound waters counteracting the temperature effect.

Plutonium complexes in water: new approach to ab initio modeling

Geometry optimization and the electronic structure calculations of Pu Z+ complexes (Z = 3–6) in water solution have been performed, within the framework of the DMol3 and Relativistic Discrete-Variational (RDV) methods. For the simulation of Pu Z+ molecular environment in aqueous solution we used 22 and 32 water molecules randomly distributed around cation. To model the effect of bulk solvent environment we used COSMO (Conductor-like Screening Model) potential for water (ε = 78.54). The obtained results showed that this approach allows the modeling of water dissociation and the formation of hydrolysis products. Our previously suggested scheme for the calculation of interaction energies between selected fragments of multi-molecular systems provides the quantitative estimation of the interaction strengths between plutonium in various oxidation states and each ligand in the first and second coordination shells in water solution.

Combining EXAFS and Computer Simulations to Refine the Structural Description of Actinyls in Water

EXAFS spectroscopy is one of the most used techniques to solve the structure of actinoid solutions. In this work a systematic analysis of the EXAFS spectra of four actinyl cations, [UO2]2+, [NpO2]2+, [NpO2]+ and [PuO2]2+ has been carried out by comparing experimental results with theoretical spectra. These were obtained by averaging individual contributions from snapshots taken from classical Molecular Dynamics simulations which employed a recently developed [AnO2]2+/+ -H2O force field based on the hydrated ion model using a quantum-mechanical (B3LYP) potential energy surface. Analysis of the complex EXAFS signal shows that both An-Oyl and An-OW single scattering paths as well as multiple scattering ones involving [AnO2]+/2+ molecular cation and first-shell water molecules are mixed up all together to produce a very complex signal. Simulated EXAFS from the B3LYP force field are in reasonable agreement for some of the cases studied, although the k= 6-8 Å-1 region is hard to be reproduced theoretically. Except uranyl, all studied actinyls are open-shell electron configurations, therefore it has been investigated how simulated EXAFS spectra are affected by minute changes of An-O bond distances produced by the inclusion of static and dynamic electron correlation in the quantum mechanical calculations. A [NpO2]+-H2O force field based on a NEVPT2 potential energy surface has been developed. The small structural changes incorporated by the electron correlation on the actinyl aqua ion geometry, typically smaller than 0.07 Å, leads to improve the simulated spectrum with respect to that obtained from the B3LYP force field. For the other open-shell actinyls, [NpO2]2+ and [PuO2]2+, a simplified strategy has been adopted to improve the simulated EXAFS spectrum. It is computed taking as reference structure the NEVPT2 optimized geometry and including the DW factors of their corresponding MD simulations employing the B3LYP force field. A better agreement between the experimental and the simulated EXAFS spectra is found, confirming the a priori guess that the inclusion of dynamic and static correlation refine the structural description of the open-shell actinyl aqua ions.

Ab initio molecular dynamics studies of hydroxide coordination of alkaline earth metals and uranyl

Ab initio molecular dynamics (AIMD) simulations of the Mg2+, Ca2+, Sr2+ and UO22+ ions in either a pure aqueous environment or an environment containing two hydroxide ions have been carried out at the density functional level of theory, employing the generalised gradient approximation via the PBE exchange-correlation functional. Calculated mean M-O bond lengths in the first solvation shell of the aquo systems compared very well to existing experimental and computational literature, with bond lengths well within values measured previously and coordination numbers in line with previously calculated values. When applied to systems containing additional hydroxide ions, the methodology revealed increased bond lengths in all systems. Proton transfer events (PTEs) were recorded and were found to be most prevalent in the strontium hydroxide systems, likely due to the low charge density of the ion and the consequent lack of hydroxide coordination. For all alkaline earths, intrashell PTEs which occurred outside of the first solvation shell were most prevalent. Only three PTEs were identified in the entire simulation data of the uranium dihydroxide system, indicating the clear impact of the increased charge density of the hexavalent uranium ion on the strength of metal-oxygen bonds in aqueous solution. Broadly, systems containing more charge dense ions were found to exhibit fewer PTEs than those containing ions of lower charge density.

Molecular Dynamics Simulation Studies on Structure, Dynamics, and Thermodynamics of Uranyl Nitrate Solution at Various Acid Concentrations

The structural and dynamical characteristics of uranyl ions in an aqueous acidic environment are of immense importance in the field of nuclear fuel reprocessing. In view of that, the structural and dynamical behavior of the uranyl ion in water has been investigated by performing molecular dynamics (MD) simulations using different force fields. All the force fields have depicted similar structural and dynamical properties except the free energy of hydration where the Guilbaud-Wipff (GW) model performs well over the others. The calculated density using MD simulations is found to be in excellent agreement with the measured experimental density, which ensures the accuracy of the adopted GW force field. The calculated surface tension and shear viscosity are seen to be increased with uranyl nitrate concentrations. At a higher concentration of about 4.0 mol/L, the supersaturation effect has been captured by an inflection in the plot of surface tension and shear viscosity against concentration because of the solution heterogeneity, which was correlated by an inflection in the scattering intensity observed by performing the dynamic light scattering experiment. The binding mode of nitrate ions with the uranyl ion is found to be concentration-dependent, and at higher concentration, it is predominantly monodentate.

UO22+ structure in solvent extraction phases resolved at molecular and supramolecular scales: a combined molecular dynamics, EXAFS and SWAXS approach

We developed a polarizable force field for unraveling the UO₂²⁺ structure in both aqueous and solvent extraction phases.

Synthesis of Microporous Covalent Phosphazene-Based Frameworks for Selective Separation of Uranium in Highly Acidic Media Based on Size-Matching Effect

On the basis of high stability of phosphorus-oxygen linkage, we constructed two microporous covalent phosphazene-based frameworks (CPFs), for the first time, by choosing hexachlorocyclotriphosphazene as a core unit and polyhydroxy aromatic compounds (hydroquinone or phloroglucinol) as monomers, named CPF-D and CPF-T, respectively. Characterization studies by using Fourier transform infrared, nuclear magnetic resonance, thermal gravimetric analysis, ⁶⁰Co γ-ray irradiation, and so forth, demonstrated that both of the CPF materials have excellent acid and radiation stability and relatively higher thermal stability. The results of batch adsorption experiments show that CPF-T is significantly more capable of sorbing uranium than CPF-D. In a pure uranium system with higher acidity (pH 1), the uranium sorption amount of CPF-T can reach up to 140 mg g^-1. Distinctively, in mixed-metal solution with 12 coexisting cations, CPF-T shows relatively stable and excellent uranium adsorption capability over a wide range of acidity (pH 4 to 3 M HNO₃), and the difference in uranium sorption amounts is less than 30% with the maximum of 0.26 mmol g^-1 at pH 4 and the minimum of 0.20 mmol g^-1 at 3 M HNO₃, which is far superior to that of the conventional solid-phase extractant (SPE) materials previously reported. The research results suggested that the sorption model based on the speculated mechanism of size-matching plus hydrogen bond network has played a dominant role in the process of uranium adsorption. The proposed strategy for the one-pot fabrication of an acid-resistant microporous framework materials by bridging the aromatic monomers via P-O bonds provides an alternative approach for the design and synthesis of new SPE materials with size-matching function desired for effective separation of uranium or other valuable metals from highly acidic environments.

Iron Vacancies Accommodate Uranyl Incorporation into Hematite

Radiotoxic uranium contamination in natural systems and nuclear waste containment can be sequestered by incorporation into naturally abundant iron (oxyhydr)oxides such as hematite (α-Fe₂O₃) during mineral growth. The stability and properties of the resulting uranium-doped material are impacted by the local coordination environment of incorporated uranium. While measurements of uranium coordination in hematite have been attempted using extended X-ray absorption fine structure (EXAFS) analysis, traditional shell-by-shell EXAFS fitting yields ambiguous results. We used hybrid functional ab initio molecular dynamics (AIMD) simulations for various defect configurations to generate synthetic EXAFS spectra which were combined with adsorbed uranyl spectra to fit experimental U L₃-edge EXAFS for U⁶⁺-doped hematite. We discovered that the hematite crystal structure accommodates a trans-dioxo uranyl-like configuration for U⁶⁺ that substitutes for structural Fe³⁺, which requires two partially protonated Fe vacancies situated at opposing corner-sharing sites. Surprisingly, the best match to experiment included significant proportions of vacancy configurations other than the minimum-energy configuration, pointing to the importance of incorporation mechanisms and kinetics in determining the state of an impurity incorporated into a host phase under low temperature hydrothermal conditions.

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and ¹⁷O NMR quadrupolar relaxation rates are investigated in the series of complexes UO₂²⁺, UO₂(CO₃)₃^4-, and (UO₂)₃(CO₃)₆^6-. Car-Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO₂²⁺ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl ¹⁷O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO₂(CO₃)₃^4- and (UO₂)₃(CO₃)₆^6- complexes show a solvent dependence of the terminal carbonate ¹⁷O and ¹³C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin-orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique ¹³C chemical shifts is correctly reproduced for (UO₂)₃(CO₃)₆^6-. NMR relaxation rate data support large ¹⁷O peak widths, but remain below those noted in the experimental literature. Comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

Trace Uranium Partitioning in a Multiphase Nano-FeOOH System

The characterization of trace elements in minerals using extended X-ray absorption fine structure (EXAFS) spectroscopy constitutes a first step toward understanding how impurities and contaminants interact with the host phase and the environment. However, limitations to EXAFS interpretation complicate the analysis of trace concentrations of impurities that are distributed across multiple phases in a heterogeneous system. Ab initio molecular dynamics (AIMD)-informed EXAFS analysis was employed to investigate the immobilization of trace uranium associated with nanophase iron (oxyhydr)oxides, a model system for the geochemical sequestration of radiotoxic actinides. The reductive transformation of ferrihydrite [Fe(OH)₃] to nanoparticulate iron oxyhydroxide minerals in the presence of uranyl (UO₂)²⁺_(aq) resulted in the preferential incorporation of U into goethite (α-FeOOH) over lepidocrocite (γ-FeOOH), even though reaction conditions favored the formation of excess lepidocrocite. This unexpected result is supported by atomically resolved transmission electron microscopy. We demonstrate how AIMD-informed EXAFS analysis lifts the strict statistical limitations and uncertainty of traditional shell-by-shell EXAFS fitting, enabling the detailed characterization of the local bonding environment, charge compensation mechanisms, and oxidation states of polyvalent impurities in complex multiphase mineral systems.

Weakly bound water structure, bond valence saturation and water dynamics at the goethite (100) surface/aqueous interface: ab initio dynamical simulations

BACKGROUND

Many important geochemical and biogeochemical reactions occur in the mineral/formation water interface of the highly abundant mineral, goethite [α-Fe(OOH)]. Ab initio molecular dynamics (AIMD) simulations of the goethite α-FeOOH (100) surface and the structure, water bond formation and dynamics of water molecules in the mineral/aqueous interface are presented. Several exchange correlation functionals were employed (PBE96, PBE96 + Grimme, and PBE0) in the simulations of a (3 × 2) goethite surface with 65 absorbed water molecules in a 3D-periodic supercell (a = 30 Å, FeOOH slab ~12 Å thick, solvation layer ~18 Å thick).

RESULTS

The lowest energy goethite (100) surface termination model was determined to have an exposed surface Fe³⁺ that was loosely capped by a water molecule and a shared hydroxide with a neighboring surface Fe³⁺. The water molecules capping surface Fe³⁺ ions were found to be loosely bound at all DFT levels with and without Grimme corrections, indicative that each surface Fe³⁺ was coordinated with only five neighbors. These long bonds were supported by bond valence theory calculations, which showed that the bond valence of the surface Fe³⁺ was saturated and surface has a neutral charge. The polarization of the water layer adjacent to the surface was found to be small and affected only the nearest water. Analysis by density difference plots and localized Boys orbitals identified three types of water molecules: those loosely bound to the surface Fe³⁺, those hydrogen bonded to the surface hydroxyl, and bulk water with tetrahedral coordination. Boys orbital analysis showed that the spin down lone pair orbital of the weakly absorbed water interact more strongly with the spin up Fe³⁺ ion. These weakly bound surface water molecules were found to rapidly exchange with the second water layer (~0.025 exchanges/ps) using a dissociative mechanism.

CONCLUSIONS

Water molecules adjacent to the surface were found to only weakly interact with the surface and as a result were readily able to exchange with the bulk water. To account for the large surface Fe-OH₂ distances in the DFT calculations it was proposed that the surface Fe³⁺ atoms, which already have their bond valence fully satisfied with only five neighbors, are under-coordinated with respect to the bulk coordination. Graphical abstract All first principle calculations, at all practically achievable levels, for the goethite 100 aqueous interface support a long bond and weak interaction between the exposed surface Fe³⁺ and water molecules capping the surface. This result is supported by bond valence theory calculations and is indicative that each surface Fe³⁺ is coordinated with only 5 neighbors.

Acidity constants and redox potentials of uranyl ions in hydrothermal solutions

We report a first principles molecular dynamics (FPMD) study of the structures, acidity constants (pK_a) and redox potentials (E⁰) of uranyl (UO₂²⁺) from ambient conditions to 573 K.

Uranyl Solvation by a Three-Dimensional Reference Interaction Site Model

We report an implementation of the three-dimensional reference interaction site model (3D RISM) that in particular addresses the treatment of the long-range Coulomb field of charged species, represented by point charges and/or a distributed charge density. A comparison of 1D and 3D results for atomic ions demonstrates a reasonable accuracy, even for a moderate size of the unit cell and a moderate grid resolution. In an application to uranyl complexes with 4-6 explicit aqua ligands and an implicit bulk solvent modeled by RISM, we show that the 3D technique is not susceptible to the deficiencies of the 1D technique exposed in our previous work [Li, Matveev, Krüger, Rösch, Comp. Theor. Chem. 2015, 1051, 151]. The 3D method eliminates the artificial superposition of explicit aqua ligands and the RISM medium and predicts essentially the same values for uranyl and uranyl-water bond lengths as a state-of-the-art polarizable continuum model. With the first solvation shell treated explicitly, the observables are nearly independent of the order of the closure relationship used when solving the set of integral equations for the various distribution functions. Furthermore, we calculated the activation barrier of water exchange with a hybrid approach that combines the 3D RISM model for the bulk aqueous solvent and a quantum mechanical description (at the level of electronic density functional theory) of uranyl interacting with explicitly represented water molecules. The calculated result agrees very well with experiment and the best theoretical estimates.

[UO2(NH3)5]Br2·NH3: synthesis, crystal structure, and speciation in liquid ammonia solution by first-principles molecular dynamics simulations

X-Ray diffraction and Car–Parrinello molecular dynamics simulations furnish insights into the speciation of uranyl(vi) in liquid ammonia, calling special attention to the effect of solvation on the U–N bond length and bond strength.

Ligand-exchange mechanism: new insight into solid-phase extraction of uranium based on a combined experimental and theoretical study

The extraction mechanism is an exchange process between the ligands on Urea-GO and the coordinated water molecules of uranyl.

Removal of fluoride and uranium by nanofiltration and reverse osmosis: a review

Inorganic contamination in drinking water, especially fluoride and uranium, has been recognized as a worldwide problem imposing a serious threat to human health. Among several treatment technologies applied for fluoride and uranium removal, nanofiltration (NF) and reverse osmosis (RO) have been studied extensively and proven to offer satisfactory results with high selectivity. In this review, a comprehensive summary and critical analysis of previous NF and RO applications on fluoride and uranium removal is presented. Fluoride retention is generally governed by size exclusion and charge interaction, while uranium retention is strongly affected by the speciation of uranium and size exclusion usually plays a predominant role for all species. Adsorption on the membrane occurs as some uranium species interact with membrane functional groups. The influence of operating conditions (pressure, crossflow velocity), water quality (concentration, solution pH), solute–solute interactions, membrane characteristics and membrane fouling on fluoride and uranium retention is critically reviewed.

Dynamics of actinyl ions in water: a molecular dynamics simulation study

The dynamics of actinyl ions (AnO₂ⁿ⁺) in aqueous solutions is important not only for the design of advanced separation processes but also for understanding the fate of actinides in the environment.

Application of Calixarenes as Macrocyclic Ligands for Uranium(VI): A Review

Calixarenes represent a well-known family of macrocyclic molecules with broad range of potential applications in chemical, analytical, and engineering materials fields. This paper covers the use of calixarenes as complexing agents for uranium(VI). The high effectiveness of calix[6]arenes in comparison to other calixarenes in uranium(VI) separation process is also presented. Processes such as liquid-liquid extraction (LLE), liquid membrane (LM) separation, and ion exchange are considered as potential fields for application of calixarenes as useful agents for bindingUO22+for effective separation from aqueous solutions containing other metal components.

Force field development for actinyl ions via quantum mechanical calculations: an approach to account for many body solvation effects

Advances in computational algorithms and methodologies make it possible to use highly accurate quantum mechanical calculations to develop force fields (pair-wise additive intermolecular potentials) for condensed phase simulations. Despite these advances, this approach faces numerous hurdles for the case of actinyl ions, AcO2(n+) (high-oxidation-state actinide dioxo cations), mainly due to the complex electronic structure resulting from an interplay of s, p, d, and f valence orbitals. Traditional methods use a pair of molecules (“dimer”) to generate a potential energy surface (PES) for force field parametrization based on the assumption that many body polarization effects are negligible. We show that this is a poor approximation for aqueous phase uranyl ions and present an alternative approach for the development of actinyl ion force fields that includes important many body solvation effects. Force fields are developed for the UO2(2+) ion with the SPC/Fw, TIP3P, TIP4P, and TIP5P water models and are validated by carrying out detailed molecular simulations on the uranyl aqua ion, one of the most characterized actinide systems. It is shown that the force fields faithfully reproduce available experimental structural data and hydration free energies. Failure to account for solvation effects when generating PES leads to overbinding between UO2(2+) and water, resulting in incorrect hydration free energies and coordination numbers. A detailed analysis of arrangement of water molecules in the first and second solvation shell of UO2(2+) is presented. The use of a simple functional form involving the sum of Lennard-Jones + Coulomb potentials makes the new force field compatible with a large number of available molecular simulation engines and common force fields.

Free energies and mechanisms of water exchange around Uranyl from first principles molecular dynamics

From density functional theory (DFT) based ab initio (Car-Parrinello) metadynamics, we compute the activation energies and mechanisms of water exchange between the first and second hydration shells of aqueous Uranyl (UO22+) using the primary hydration number of U as the reaction coordinate. The free energy and activation barrier of the water dissociation reaction [UO2(OH2)5]2+(aq) → [UO2(OH2)4]2+(aq) + H2O are 0.7 kcal and 4.7 kcal/mol respectively. The free energy is in good agreement with previous theoretical (-2.7 to +1.2 kcal/mol) and experimental (0.5 to 2.2 kcal/mol) data. The associative reaction [UO2(OH2)5]2+(aq) + H2O → [UO2(OH2)6]2+(aq) is short-lived with a free energy and activation barrier of +7.9 kcal/mol and +8.9 kca/mol respectively; it is therefore classified as associative-interchange. On the basis of the free energy differences and activation barriers, we predict that the dominant exchange mechanism between [UO2(OH2)5]2+(aq) and bulk water is dissociative.