Effect of Uranyl Ion Concentration on Structure and Dynamics of Aqueous Uranyl Solution: A Molecular Dynamics Simulation Study

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.

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.

Bonding trends across the series of tricarbonato-actinyl anions [(AnO2)(CO3)3]4- (An = U-Cm): the plutonium turn

Actinyl-tricarbonato anions [(AnO₂)(CO₃)₃]^4- (An = U-Cm) in various environments were investigated using theoretical approaches of quantum-mechanics, molecular-mechanics and cluster-models. Cations and solvent molecules in the 2^nd coordination sphere affect the equatorial An←O_eq bonds more than the axial An[triple bond, length as m-dash]O_ax bonds. Common actinide contraction is found for calculated and experimental axial bond lengths of ₉₂U to ₉₄Pu, though no longer for ₉₄Pu to ₉₆Cm. The tendency of U to Pu forming actinyl(vi) species dwindles away toward Cm, which already features the preferred An^III/Ln^III oxidation state of the later actinides and all lanthanides. The well known change from d-type to typical U-Pu-Cm type and then to f-type behavior is labeled as the plutonium turn, a phenomenon that is caused by f-orbital energy-decrease and f-orbital localization with increase of both nuclear charge and oxidation state, and a non-linear variation of effective f-electron population across the actinide series. Both orbital and configuration mixing and occupation of antibonding 5f type orbitals increase, weakening the AnO_ax bonds and reducing the highest possible oxidation states of the later actinides.

Determination of the Structures of Uranyl-Tri-n-butyl-Phosphate Aggregates by Coupling Experimental Results with Molecular Dynamic Simulations

The complex structure of a plutonium uranium refining by extraction (PUREX) process organic phase was characterized by combining results from experiments and molecular dynamics simulations. For the first time, the molecular interactions between tri-n-butyl phosphate (TBP) and the extracted solutes, as well as TBP aggregation after the extraction of water and/or uranyl nitrate, were described and analyzed concomitantly. Coupling molecular dynamics simulations with small- and wide-angle X-ray scattering (SWAXS) experiments can lead to simulated organic solutions that are representative of the experimental ones, even for high extractant and solute concentrations. Furthermore, this coupling is well adapted for the interpretation of SWAXS experiments without preliminary hypothesis on the size or shape of aggregates. The results link together previous literature studies obtained for each level of depiction separately (complexation or aggregation). Without uranium, or at low metal concentration, almost no aggregation was observed. At high uranium concentration, organic phases contain small [UO₂ (NO₃ )₂ (TBP)₂ ]_n polymetallic aggregates (with n=2 to 4), in which the 1:2 U/TBP stoichiometry is preserved.

Dynamics of Humic Acid and Its Interaction with Uranyl in the Presence of Hydrophobic Surface Implicated by Molecular Dynamics Simulations

This work targeted a molecular level of understanding on the dynamics of humic acid (HA) and its interaction with uranyl in the presence of hydrophobic surface mimicked by a carbon nanotube (CNT), which also represents a potential intruder in the environment accompanying with the development of nanotechnology. In aqueous phase, uranyl and HA were observed to build close contact spontaneously, driven by electrostatic interaction, leading to a more compact conformation of HA. The presence of CNT unfolds HA via π-π interactions with the aromatic rings of HA without significant perturbation on the interaction strength between HA and uranyl. These results show that the hydrophilic uranyl and the hydrophobic CNT influence the folding behavior of HA in distinct manners, which represents two fundamental mechanisms that the folding behavior of HA may be modulated in the environment, that is, uranyl enhances the folding of HA via electrostatic interactions, whereas CNT impedes its spontaneous folding via van der Waals (vdW) interactions. The work also provides molecular level of evidence on the transformation of a hydrophobic surface into a hydrophilic one via noncovalent functionalization by HA, which in turn affects the migration of HA and the cations it binds to.

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.