Tri-n-butyl phosphate (TBP) is an important extractant used in the solvent extraction process for recovering uranium and plutonium from used nuclear fuel. An atomistic molecular dynamics study was used to understand the fundamental molecular-level behavior of extracting agents in solution. Atomistic parametrization was carried out using the AMBER force field to model the TBP molecule and n-dodecane molecule, a commonly used organic solvent. Validation of the optimized force field was accomplished through various thermophysical properties of pure TBP and pure n-dodecane in the bulk liquid phase. The mass density, dipole moment, self-diffusion coefficient, and heat of vaporization were calculated from our simulations and compared favorably with experimental values. The molecular structure of TBPs in n-dodecane at a dilute TBP concentration was examined based on radial distribution functions. 1D and 2D potential mean force studies were carried out to establish the criteria for identifying TBP aggregates. The dimerization constant of TBP in the TBP/n-dodecane mixture was also obtained and matched the experimental value.

Tri-n-butyl phosphate (TBP), a representative of neutral organophosphorous ligands, is an important extractant used in the solvent extraction process for the recovery of uranium and plutonium from spent nuclear fuel. Microscopic pictures of TBP isomerism and its behavior in n-dodecane diluent were investigated utilizing MD simulations with previously optimized force field parameters for TBP and n-dodecane. Potential mean force (PMF) calculations on a single TBP molecule show seven probable TBP isomers. Radial distribution functions (RDFs) of TBP suggest the existence of TBP trimers at high TBP concentrations in addition to dimers. 2D PMF calculations were performed to determine the angle and distance criteria for TBP trimers. The dimerization and trimerization constants of TBP in n-dodecane were obtained and match our own experimental values using the FTIR technique. The new insights into the conformational behaviors of the TBP molecule as a monomer and as part of an aggregate could greatly aid in the understanding of the complexation between TBP and metal ions in a solvent extraction system.

Tri-n-butyl phosphate (TBP), a representative of neutral organophosphorous metal-ion-extracting reagents, is an important ligand used in solvent extraction processes for the recovery of uranium and plutonium from spent nuclear fuel, as well as other non-nuclear applications. Ligand-ligand and organic solvent-ligand interactions play an important role in these processes. The self-association behavior of TBP in various alkane diluents of different chain lengths (8, 12, and 16 carbons) and a branched alkane (iso-octane) was investigated by Fourier transform infrared spectroscopic measurements. By careful deconvolution of the spectra into multiple peaks, our results indicate that TBP self-associates to form not only dimers, as previous studies showed, but also trimers in the practical concentration range. Using a mathematical fitting procedure, the dimerization and trimerization constants were determined. As expected, these equilibrium constants are dependent on the solvent used. As the alkane chain for linear hydrocarbon solvents becomes longer, dimerization decreases whereas trimerization increases. For the more branched hydrocarbon, we observe a significantly higher dimerization constant. These effects are most likely due to the intermolecular van der Waals interactions between the butyl tails of each TBP molecule and the diluent hydrocarbon chain as all solvents in this study are relatively nonpolar.

We report a classical rate theory approach to predict the exchange mechanism that occurs between water and aqueous uranyl ion. Using our water and ion-water polarizable force field and molecular dynamics techniques, we computed the potentials of mean force for the uranyl ionwater pair as a function of different pressures at ambient temperature. These potentials of mean force were used to calculate rate constants using transition rate theory; the transmission coefficients also were examined using the reactive flux method and Grote-Hynes approach. The computed activation volumes are positive; thus, the mechanism of this particular water-exchange is a dissociative process.

In solvent extraction processes for recovering metal ions from used nuclear fuel, as well as other industrial applications, a better understanding of the metal complex phase transfer phenomenon would greatly aid ligand design and process optimization. We have approached this challenge by utilizing the classical molecular dynamics simulations technique to gain visual appreciation of the vapor/liquid and liquid/liquid interface between tri-n-butyl phosphate (TBP) and n-dodecane with air and water. In this study, we successfully reparameterized polarizable force fields for TBP and n-dodecane that accurately reproduced several of their thermophysical properties such as density, heat of vaporization, and dipole moment. Our models were able to predict the surface and interfacial tension of different systems when compared to experimental results that were also performed by us. Through this study, we gained atomistic understanding of the behaviors of TBP and n-dodecane at the interface against air and water, useful in further computational studies of such systems. Finally, our studies indicate that the initial configuration of a simulation may have a large effect on the final result.