Structure, Kinetics, and Thermodynamics of the Aqueous Uranyl(VI) Cation

Analyzing (3-Aminopropyl)triethoxysilane-Functionalized Porous Silica for Aqueous Uranium Removal: A Study on the Adsorption Behavior

This study synthesized (3-aminopropyl)triethoxysilane-functionalized porous silica (AP@MPS) to adsorb aqueous uranium (U(VI)). To comprehensively analyze the surface properties of the AP@MPS materials, a combination of SEM, BET, XPS, NMR, and zeta potential tests were conducted. The adsorption experiments for U(VI) revealed the rapid and efficient adsorption capacity of AP@MPS, with the solution condition of a constant solution pH = 6.5, an initial U(VI) concentration of 600 mg × L-1, a maximum U(VI) capacity of AP@MPS reaching 381.44 mg-U per gram of adsorbent, and a removal rate = 63.6%. Among the four types of AP@MPS with different average pore sizes tested, the one with an average pore size of 2.7 nm exhibited the highest U(VI) capacity, particularly at a pH of 6.5. The adsorption data exhibited a strong fit with the Langmuir model, and the calculated adsorption energy aligned closely with the findings from the Potential of Mean Force (PMF) analysis. The outcomes obtained using the Surface Complex Formation Model (SCFM) highlight the dominance of the coulombic force ΔG0coul as the principal component of the adsorption energy (ΔG0ads). This work garnered insights into the adsorption mechanism by meticulously examining the ΔG0ads across a pH ranging from 4 to 8. In essence, this study's findings furnish crucial insights for the future design of analogous adsorbents, thereby advancing the realm of uranium(VI) removal methodologies.

Solvent effects on extractant conformational energetics in liquid-liquid extraction: a simulation study of molecular solvents and ionic liquids

Extractant design in liquid-liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant-metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. We calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf2N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

An electrochemical technique for sensing uranium adsorption and desorption

Uranium is a toxic, heavy metal that can pose elevated health risks if leached into the environment. Uranium mobility is dependent on many factors, including speciation, solution pH, ligands in solution, and presence of surfaces. Surface adsorption is one phenomenon that inhibits uranium mobility in the environment and is studied as a naturally occurring phenomenon as well an intentional tool for environmental remediation. This work presents and validates a potentiometric, electrochemical technique for sensing uranium adsorption on and desorption from an electrochemically active surface solely through changes of the electrode potential. This novel electrochemical technique presents a new lens to study adsorption that complements external techniques (e.g., spectroscopy). Indication of adsorption and desorption via the electrochemical technique are gravimetrically validated using an electrochemical quartz crystal microbalance. This work contributes a unique, complementary technique that corroborates the adsorption of uranium on an electrode surface.

In Silico Simulations Reveal Molecular Mechanism of Uranyl Ion Toxicity towards DNA-Binding Domain of PARP-1 Protein

The molecular toxicity of the uranyl ion (UO22+) in living cells is primarily determined by its high affinity to both native and potential metal-binding sites that commonly occur in the structure of biomolecules. Recent advances in computational and experimental research have shed light on the structural properties and functional impacts of uranyl binding to proteins, organic ligands, nucleic acids, and their complexes. In the present work, we report the results of the computational investigation of the uranyl-mediated loss of DNA-binding activity of PARP-1, a eukaryotic enzyme that participates in DNA repair, cell differentiation, and the induction of inflammation. The latest experimental studies have shown that the uranyl ion directly interacts with its DNA-binding subdomains, zinc fingers Zn1 and Zn2, and alters their tertiary structure. Here, we propose an atomistic mechanism underlying this process and compute the free energy change along the suggested pathway. Our Quantum Mechanics/Molecular Mechanics (QM/MM) simulations of the Zn2-UO22+ complex indicate that the uranyl ion replaces zinc in its native binding site. However, the resulting state is destroyed due to the spontaneous internal hydrolysis of the U-Cys162 coordination bond. Despite the enthalpy of hydrolysis being +2.8 kcal/mol, the overall reaction free energy change is -0.6 kcal/mol, which is attributed to the loss of domain's native tertiary structure originally maintained by a zinc ion. The subsequent reorganization of the binding site includes the association of the uranyl ion with the Glu190/Asp191 acidic cluster and significant perturbations in the domain's tertiary structure driven by a further decrease in the free energy by 6.8 kcal/mol. The disruption of the DNA-binding interface revealed in our study is consistent with previous experimental findings and explains the loss of PARP-like zinc fingers' affinity for nucleic acids.

Force Field Parameterization of Actinyl Molecular Cations Using the 12-6-4 Model

In this work, a set of 12-6-4 force fields (FFs) parameters were developed for the actinyl molecular cations, AnO2n+ (n = 1, 2), from uranium to plutonium for classical molecular dynamics (MD) for four water models: TIP3P, SPC/E, OPC3, and TIP4Pew. Such a non-bonded potential model taking into account the induced dipole between the metallic center and the surrounding molecules has shown better performances for various cations than the classic 12-6 non-bonded potentials. The parametrization method proposed elsewhere for metallic cations has been extended to these molecular cations. In contrast to the actinyl 12-6 FFs from the literature, the new models reproduce correctly both solvation and thermodynamic properties, thanks to the inclusion of the induced dipole term (C4). The transferability of such force fields was assessed by performing MD simulations of carbonato actinyl species, which are highly implicated in actinide migration or actinide extraction from seawater. A highly satisfying agreement was found when comparing the EXAFS signals computed from our MD simulation to the experimental ones. The set of FFs developed here opens new possibilities for the study of actinide chemistry.

Stochastic Model of Solvent Exchange in the First Coordination Shell of Aqua Ions

Ion microsolvation is a basic, yet fundamental, process of ionic solutions underlying many relevant phenomena in either biological or nanotechnological applications, such as solvent reorganization energy, ion transport, catalytic activity, and so on. As a consequence, it is a topic of extensive investigations by various experimental techniques, ranging from X-ray diffraction to NMR relaxation and from calorimetry to vibrational spectroscopy, and theoretical approaches, especially those based on molecular dynamics (MD) simulations. The conventional microscopic view of ion solvation is usually provided by a "static" cluster model representing the first ion-solvent coordination shell. Despite the merits of such a simple model, however, ion coordination in solution should be better regarded as a complex population of dynamically interchanging molecular configurations. Such a more comprehensive view is more subtle to characterize and often elusive to standard approaches. In this work, we report on an effective computational strategy aiming at providing a detailed picture of solvent coordination and exchange around aqua ions, thus including the main structural, thermodynamic, and dynamic properties of ion microsolvation, such as the most probable first-shell complex structures, the corresponding free energies, the interchanging energy barriers, and the solvent-exchange rates. Assuming the solvent coordination number as an effective reaction coordinate and combining MD simulations with enhanced sampling and master-equation approaches, we propose a stochastic model suitable for properly describing, at the same time, the thermodynamics and kinetics of ion-water coordination. The model is successfully tested toward various divalent ions (Ca²⁺, Zn²⁺, Hg²⁺, and Cd²⁺) in aqueous solution, considering also the case of a high ionic concentration. Results show a very good agreement with those issuing from brute-force MD simulations, when available, and support the reliable prediction of rare ion-water complexes and slow water exchange rates not easily accessible to usual computational methods.

A computational investigation of orbital overlap versus energy degeneracy covalency in [UE2]2+ (E = O, S, Se, Te) complexes

Covalency in analogues of uranyl with heavy chalcogens is explored using DFT, and traced to increased energy-degeneracy as the group is descended.

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.

Uranyl Arsenate Complexes in Aqueous Solution: Insights from First-Principles Molecular Dynamics Simulations

In this study, the structures and acidity constants (p K_a's) of uranyl arsenate complexes in solutions have been revealed by using the first principle molecular dynamics technique. The results show that uranyl and arsenate form stable complexes with the U/As ratios of 1:1 and 1:2, and the bidentate complexation between U and As is highly favored. Speciation-pH distributions are derived based on free energy and p K_a calculations, which indicate that for the 1:1 species, UO₂(H₂AsO₄)(H₂O)₃⁺ is the major species at pH < 7, while UO₂(HAsO₄)(H₂O)₃⁰ and UO₂(AsO₄)(H₂O)₃^- dominate in acid-to-alkaline and extreme alkaline pH ranges. For the 1:2 species, UO₂(H₂AsO₄)₂(H₂O)⁰ is dominant under acid-to-neutral pH conditions, while UO₂(HAsO₄)(H₂AsO₄)(H₂O)^-, UO₂(HAsO₄)(HAsO₄)(H₂O)^2-, and UO₂(AsO₄)(HAsO₄)(H₂O)^3- become the major forms in the pH range of 7.2-10.7, 10.7-12.1, and >12.1, respectively.

A molecular dynamics investigation of actinyl–ligand speciation in aqueous solution

Actinyl ions (AnO₂ⁿ⁺), the form in which actinides are commonly found in aqueous solution, are important species in the nuclear fuel cycle.

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.

First-principles molecular dynamics simulation of the Ca2UO2(CO3)3 complex in water

First principles molecular dynamics simulation reveals the structure and solvation of the Ca₂UO₂(CO₃)₃ complex in water and the hydrogen bonding network that differentiates the two Ca ions.

Theoretical insight into the binding affinity enhancement of serine with the uranyl ion through phosphorylation

Various type complexes of uranyl ion with serine and phosphoserine were investigated which showed phosphorylation induced binding affinity enhancement of serine with the uranyl ion.

Mechanism of adsorption affinity and capacity of Mg(OH)2 to uranyl revealed by molecular dynamics simulation

The redistribution of surface OH groups results in a negative charge site which facilitates the uranyl adsorption.

Passage of TBP–uranyl complexes from aqueous–organic interface to the organic phase: insights from molecular dynamics simulation

Water/organic interface representing TBP orientation for neutral versus acidic interface and occurrence of UO₂²⁺–TBP–NO₃⁻ species in various stoichiometry.

Molecular dynamics simulation study of distribution and dynamics of aqueous solutions of uranyl ions: the effect of varying temperature and concentration

Investigating the characteristics of actinyl ions has been of great interest due to their direct relevance in the nuclear fuel cycle. All-atom molecular dynamics simulations have been employed to study the orientational structure and dynamics of aqueous solutions of uranyl ions of various concentrations. The orientational structure of water around a uranyl ion has been thoroughly investigated by calculating different orientational probability distributions corresponding to different molecular axes of water. The orientational distribution of water molecules in the first coordination shell of a uranyl ion is found to be markedly different from that in bulk water. Analysis of counterion distribution around the uranyl ion reveals the presence of nitrate ions along with water molecules in the first solvation shell. From the comparison of the number of coordinated water and nitrate ions at various uranyl nitrate concentrations, it is evident that these two species compete for occupying the first solvation shell of the uranyl ion. Orientational dynamics of water molecules about different molecular axes of water in the vicinity of uranyl ions have also been investigated and decreasing orientational mobility of water with increasing uranyl concentration has been found. However, it is observed that the orientational dynamics remains more or less the same whether we consider all the water molecules in the aqueous solution or only the solvation shell water molecules. The effect of temperature on the translational and orientational characteristics of the aqueous uranyl solutions has also been studied in detail.

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.

Solvation structure and transport properties of alkali cations in dimethyl sulfoxide under exogenous static electric fields

A combination of molecular dynamics simulations and pulsed field gradient nuclear magnetic resonance spectroscopy is used to investigate the role of exogenous electric fields on the solvation structure and dynamics of alkali ions in dimethyl sulfoxide (DMSO) and as a function of temperature. Good agreement was obtained, for select alkali ions in the absence of an electric field, between calculated and experimentally determined diffusion coefficients normalized to that of pure DMSO. Our results indicate that temperatures of up to 400 K and external electric fields of up to 1 V nm(-1) have minimal effects on the solvation structure of the smaller alkali cations (Li(+) and Na(+)) due to their relatively strong ion-solvent interactions, whereas the solvation structures of the larger alkali cations (K(+), Rb(+), and Cs(+)) are significantly affected. In addition, although the DMSO exchange dynamics in the first solvation shell differ markedly for the two groups, the drift velocities and mobilities are not significantly affected by the nature of the alkali ion. Overall, although exogenous electric fields induce a drift displacement, their presence does not significantly affect the random diffusive displacement of the alkali ions in DMSO. System temperature is found to have generally a stronger influence on dynamical properties, such as the DMSO exchange dynamics and the ion mobilities, than the presence of electric fields.

Investigations on preferential Pu(IV) extraction over U(VI) by N,N-dihexyloctanamide versus tri-n-butyl phosphate: evidence through small angle neutron scattering and DFT studies

Straight chain amide N,N-dihexyloctanamide (DHOA) has been found to be a promising alternative extractant to tri-n-butyl phosphate (TBP) for the reprocessing of irradiated uranium- and thorium-based fuels. Unlike TBP, DHOA displays preferential extraction of Pu(IV) over U(VI) at higher acidities (≥3 M HNO3) and poor extraction at lower acidities. Density functional theory (DFT) based calculations have been carried out on the structures and relative binding energies of U(VI) and Pu(IV) with the extractant molecules. These calculations suggest that the differential hardness of the two extractants is responsible for the preferential binding/complexation of TBP to uranyl, whereas the softer DHOA and the bulky nature of the extractant lead to stronger binding/complexation of DHOA to Pu(IV). In conjunction with quantum chemical calculations, small angle neutron scattering (SANS) measurements have also been performed for understanding the stoichiometry of the complex formed that leads to relatively lower extraction of Th(IV) (a model for Pu(IV)) as compared to U(VI) using DHOA and TBP as the extractants. The combined experimental and theoretical studies helped us to understand the superior complexation/extraction behavior of Pu(IV) over U(VI) with DHOA.

Molecular dynamics simulations of uranyl and uranyl carbonate adsorption at aluminosilicate surfaces

Adsorption at mineral surfaces is a critical factor controlling the mobility of uranium(VI) in aqueous environments. Therefore, molecular dynamics (MD) simulations were performed to investigate uranyl(VI) adsorption onto two neutral aluminosilicate surfaces, namely, the orthoclase (001) surface and the octahedral aluminum sheet of the kaolinite (001) surface. Although uranyl preferentially adsorbs as a bidentate inner-sphere complex on both surfaces, the free energy of adsorption on the orthoclase surface (-15 kcal mol(-1)) is significantly more favorable than that at the kaolinite surface (-3 kcal mol(-1)), which is attributed to differences in surface functional groups and the ability of the orthoclase surface to release a surface potassium ion upon uranyl adsorption. The structures of the adsorbed complexes compare favorably with X-ray absorption spectroscopy results. Simulations of the adsorption of uranyl complexes with up to three carbonate ligands revealed that uranyl complexes coordinated to up to two carbonate ions are stable on the orthoclase surface whereas uranyl carbonate surface complexes are unfavored at the kaolinite surface. Combining the MD-derived equilibrium adsorption constants for orthoclase with aqueous equilibrium constants for uranyl carbonate species indicates the presence of adsorbed uranium complexes with one or two carbonates under alkaline conditions, in support of current uranium(VI) surface complexation models.

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.