Experimental Coordination Environment of Uranyl(VI) in Aqueous Solution

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.

A Coupled EXAFS-Molecular Dynamics Study on PuO2+ and NpO2+ Hydration: The Importance of Electron Correlation in Force-Field Building

The physicochemical properties of the monovalent actinyl cations, PuO₂⁺ and NpO₂⁺, in water have been studied by means of classical molecular dynamic simulations. A specific set of cation-water intermolecular potentials based on ab initio potential energy surfaces has been built on the basis of the hydrated ion concept. The TIP4P water model was adopted. Given the paramagnetic character of these actinyls, the cation-water interaction energies were computed from highly correlated wave functions using the NEVPT2 method. It is shown that the multideterminantal character of the wave function has a relevant effect on the main distances of the hydrated molecular cations. Several structural, dynamical, and energetic properties of the aqueous solutions have been obtained and analyzed. Structural RDF analysis gives An-O_yl distances of 1.82 and 1.84 Å and An-O(water) distances of 2.51 and 2.53 Å for PuO₂⁺ and NpO₂⁺ in water, respectively. Experimental EXAFS spectra from dilute aqueous solutions of PuO₂⁺ and NpO₂⁺ are revisited and analyzed, assuming tetra- and pentahydration of the actinyl cations. Simulated EXAFS spectra have been computed from the snapshots of the MD simulations. Good agreement with the experimental information available is found. The global analysis leads us to conclude that both PuO₂⁺ and NpO₂⁺ cations in water are stable pentahydrated aqua ions.

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.

Aquatic interaction of uranium with two naturally ubiquitous pyrazine compounds: Speciation studies by experiment and theory

The present studies interpret the speciation of uranyl (UO₂²⁺) with the most ubiquitous class of natural species named pyrazines in terms of stability, speciation and its identification, thermodynamics, spectral properties determined by a range of experimental techniques and further evidenced by theoretical insights. UO₂²⁺ forms ML and ML₂ kind of species with a qualitative detection of ML₃ species, while the ESI-MS identified the formation of all the complexes including ML₃. Both the ligands act as bidentate chelators with a difference in ring size and coordinating atoms in the complex formed. The ML₃ complexes involve the third ligand participation as monodentate via carboxylate only due to the restricted coordination number and space around the UO₂²⁺ ion to accommodate three ligand molecules in its primary coordination sphere. All the complexes are found to be endothermic and purely entropy driven formations. The complex formations showed redshift in the absorption spectra and the shift was further enhanced from ML to ML₂ formation. The UO₂²⁺ ion redox properties are used to explore the redox potential and heterogeneous electron-transfer kinetic parameters as a function of pH and concentration of UO₂²⁺ in presence of pyrazine carboxylates. Interestingly, the cyclic voltammograms identified the ligands also as redox sensitive. The theoretical calculation gave inputs to understand the complex formation at the molecular level with major emphasis on geometry optimization, energetics, bonding parameters, molecular orbital diagrams and bond critical point analyses. The experimental observations in combination with theoretical addendum provided detailed knowledge on the interaction of UO₂²⁺ with pyrazine-2-carboxylate and pyrazine-2,3-dicarboxylates.

Theoretical Insight into Coordination Chemistry of Am(VI) and Am(V) with Phenanthroline Ligand: Implications for High Oxidation State Based Minor Actinide Separation

Despite continuing and burgeoning interest in americium (Am) coordination chemistry in recent years, investigations of the electronic structures and bonding chemistry of high oxidation state americium complexes and their implications for minor actinide separation remain relatively less explored to date. Here, we used density functional theory (DFT) to create high oxidation states of americium but experimentally feasible models of Am(V) and Am(VI) complexes of phenanthroline ligand (DAPhen) as [AmO₂(L)]^1+/2+ and [AmO₃(L)]¹⁺ (L = 2,9-bis[(N,N-dimethyl)-carbonyl]-1,10-phenanthroline (oxo-DAPhen, L_O) and 2,9-bis[(N,N-dimethyl)-thio-carbonyl]-1,10-phenanthroline (thio-DAPhen, L_S)), meanwhile comparing these with [UO₂(L)]²⁺. On the basis of the calculations, the Am(V) and Am(VI) oxidation state are thermodynamically feasible and can be stabilized by DAPhen ligands. From a comparative study, the strength of thio-DAPhen in the separation of high oxidation state Am emerges better than does oxo-DAPhen, which relates to the nature, energy level, and spatial arrangement of their frontier orbitals. This study provides fundamental knowledge toward understanding the transuranic separations processes, which has implications in designing new, more selective extraction processes for the separation of Am from curium (Cm) as well as lanthanide.

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.

Competition/Cooperation between Humic Acid and Graphene Oxide in Uranyl Adsorption Implicated by Molecular Dynamics Simulations

Molecular dynamics (MD) simulations were performed to investigate the influence of curvature and backbone rigidity of an oxygenated surface, here graphene oxide (GO), on its adsorption of uranyl in collaboration with humic acid (HA). The planar curvature of GO was found to be beneficial in impeding the folding of HA. This, together with its rigidity that helps stabilize the extended conformation of HA, offered rich binding sites to interact with uranyl with only marginal loss of binding strength. According to our simulations, the interaction between uranyl and GO was mainly driven by electrostatic interactions. The presence of HA not only provided multiple sites to compete/cooperate with GO for adsorption of free uranyl but also interacted with GO acting as a "bridge" to connect uranyl and GO. The potential of mean force (PMF) profiles implied that HA significantly enhanced the interaction strength between uranyl and GO and stabilized the uranyl-GO complex. Meanwhile, GO could reduce the diffusion coefficients of uranyl and HA and retard their migrations in aqueous solution. This work provides theoretical hints on the GO-based remediation strategies for the sites contaminated by uranium or other heavy metal ions and oxygenated organic pollutants.

Density functional theory (DFT) calculations of VI/V reduction potentials of uranyl coordination complexes in non-aqueous solutions

Of particular interest within the +6 uranium complexes is the linear uranyl(vi) cation and it forms numerous coordination complexes in solution and exhibits incongruent redox behavior depending on coordinating ligands. In this study, to determine the reduction potentials of uranyl complexes in non-aqueous solutions, a hybrid density functional theory (DFT) approach was used in which two different DFT functionals, B3LYP and M06, were applied. Bulk solvent effects were invoked through the conductor-like polarizable continuum model. The solute cavities were described with the united-atom Kohn-Sham (UAKS) cavity definition. Inside the cavity the dielectric constant matches the value of a vacuum and outside the cavity the dielectric constant value is the same as that of the solvent of interest, for example, dimethyl sulfoxide (DMSO), dimethylformamide (DMF), dichloromethane (DCM), acetonitrile and pyridine. With the help of the Nernst equation the calculated reduction potentials with respect to the ferrocene (Fc) reference electrode are converted into reduction free energies (RFEs). Uranyl complexes of organic ligands which range from mono- to hexa-dentate coordination modes were investigated in non-aqueous solutions of DMSO, DMF, DCM, acetonitrile and pyridine solutions. The effect of the spin-orbit correction and the reference electrode correction on the RFEs and various methods such as the direct method and the isodesmic reaction model were explored. Overall, our computational determination of RFEs of uranyl complexes in various non-aqueous solutions demonstrates that the RFEs can be obtained within ∼0.2 eV of experimental values.

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.

Tris-amidoximate uranyl complexes via η2 binding mode coordinated in aqueous solution shown by X-ray absorption spectroscopy and density functional theory methods

The present study sheds some light on the long-standing debate concerning the coordination properties between uranyl ions and the amidoxime ligand, which is a key ingredient for achieving efficient extraction of uranium. Using X-ray absorption fine structure combined with theoretical simulation methods, the binding mode and bonding nature of a uranyl-amidoxime complex in aqueous solution were determined for the first time. The results show that in a highly concentrated amidoxime solution the preferred binding mode between UO₂²⁺ and the amidoxime ligand is η² coordination with tris-amidoximate species. In such a uranyl-amidoximate complex with η² binding motif, strong covalent interaction and orbital hybridization between U 5f/6d and (N, O) 2p should be responsible for the excellent binding ability of the amidoximate ligand to uranyl. The study was performed directly in aqueous solution to avoid the possible binding mode differences caused by crystallization of a single-crystal sample. This work also is an example of the simultaneous study of local structure and electronic structure in solution systems using combined diagnostic tools.

Structural and magnetic susceptibility characterization of Pu(v) aqua ion using sonochemistry as a facile synthesis method

The facile sonochemical preparation of pure, stable and concentrated Pu(v) aqueous solutions allowed to investigate its solvation environment and magnetic properties.

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.

Crown ether complexes of actinyls: a computational assessment of AnO2(15-crown-5)2+ (An = U, Np, Pu, Am, Cm)

Computational characterization of AnO₂²⁺–(15-crown-5) complexes (An = U, Np, Pu, Am, and Cm) reveals actinyl insertion coordination to crown ether.

Theoretical study on complexation of U(vi) with ODA, IDA and TDA based on density functional theory

Precise structures of U(vi) complexes with ODA, IDA and TDA in different binding modes.

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.

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.