Liquid-Liquid Extraction of Uranyl by TBP: The TBP and Ions Models and Related Interfacial Features Revisited by MD and PMF Simulations

We report a molecular dynamics (MD) study of biphasic systems involved in the liquid-liquid extraction of uranyl nitrate by tri-n-butylphosphate (TBP) to hexane, from "pH neutral" or acidic (3 M nitric acid) aqueous solutions, to assess the model dependence of the surface activity and partitioning of TBP alone, of its UO2(NO3)(2)(TBP)(2) complex, and of UO2(NO3)(2) or UO22+ uncomplexed. For this purpose, we first compare several electrostatic representations of TBP with regards to its polarity and conformational properties, its interactions with H2O, HNO3, and UO2(NO3)(2) species, its relative free energies of solvation in water or oil environments, the properties of the pure TBP liquid and of the pure-TBP/water interface. The free energies of transfer of TBP, UO2(NO3)(2), UO22+, and the UO2(NO3)(2)(TBP)(2) complex across the water/oil interface are then investigated by potential of mean force (PMF) calculations, comparing different TBP models and two charge models of uranyl nitrate. Describing uranyl and nitrate ions with integer charges (+2 and -1, respectively) is shown to exaggerate the hydrophilicity and surface activity of the UO2(NO3)(2)(TBP)(2) complex. With more appropriate ESP charges, mimicking charge transfer and polarization effects in the UO2(NO3)(2) moiety or in the whole complex, the latter is no more surface active. This feature is confirmed by MD, PMF, and mixing-demixing simulations with or without polarization. Furthermore, with ESP charges, pulling the UO2(NO3)(2) species to the TBP phase affords the formation of UO2(NO3)(2)(TBP)(2) at the interface, followed by its energetically favorable extraction. The neutral complexes should therefore not accumulate at the interface during the extraction process, but diffuse. to the oil phase. A similar feature is found for an UO2(NO3)(2)(Amide)(2) neutral complex with fatty amide extracting ligands, calling for further simulations and experimental studies (e.g., time evolution of the nonlinear spectroscopic signature and of surface tension) on the interfacial landscape upon ion extraction.