A Molecular Dynamics Study of Tributyl Phosphate and Diamyl Amyl Phosphonate Self-Aggregation in Dodecane and Octane

The pervasive impact of critical fluctuations in liquid-liquid extraction organic phases

Liquid-liquid extraction is an essential chemical separation technique where polar solutes are extracted from an aqueous phase into a nonpolar organic solvent by amphiphilic extractant molecules. A fundamental limitation to the efficiency of this important technology is third phase formation, wherein the organic phase splits upon sufficient loading of polar solutes. The nanoscale drivers of phase splitting are challenging to understand in the complex hierarchically structured organic phases. In this study, we demonstrate that the organic phase structure and phase behavior are fundamentally connected in a way than can be understood with critical phenomena theory. For a series of binary mixtures of trialkyl phosphate extractants with linear alkane diluents, we combine small angle x-ray scattering and molecular dynamics simulations to demonstrate how the organic phase mesostructure over a wide range of compositions is dominated by critical concentration fluctuations associated with the critical point of the third phase formation phase transition. These findings reconcile many longstanding inconsistencies in the literature where small angle scattering features, also consistent with such critical fluctuations, were interpreted as reverse micellar-like particles. Overall, this study shows how the organic phase mesostructure and phase behavior are intrinsically linked, deepening our understanding of both and providing a new framework for using molecular structure and thermodynamic variables to control mesostructure and phase behavior in liquid-liquid extraction.

Hydrogen-Bond-Driven Chemical Separations: Elucidating the Interfacial Steps of Self-Assembly in Solvent Extraction

Chemical separations, particularly liquid extractions, are pervasive in academic and industrial laboratories, yet a mechanistic understanding of the events governing their function are obscured by interfacial phenomena that are notoriously difficult to measure. In this work, we investigate the fundamental steps of ligand self-assembly as driven by changes in the interfacial H-bonding network using vibrational sum frequency generation. Our results show how the bulk pH modulates the interfacial structure of extractants at the buried oil/aqueous interface via the formation of unique H-bonding networks that order and bridge ligands to produce self-assembled aggregates. These extended H-bonded structures are key to the subsequent extraction of Co²⁺ from the aqueous phase in promoting micelle formation and subsequent ejection of the said micelle into the oil phase. The combination of static and time-resolved measurements reveals the events underlying complexities of liquid extractions at high [Co²⁺]:[ligand] ratios by showing an evolution of interfacially assembled structures that are readily tuned on a chemical basis by altering the compositions of the aqueous phase. The results of this work point to new principles to design-applied separations through the manipulation of surface charge, electrostatic screening, and the associated H-bonding networks that arise at the interface to facilitate organization and subsequent extraction.

Hierarchical phenomena in multicomponent liquids: simulation methods, analysis, chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

Molecular Insights into Water Clusters Formed in Tributylphosphate-Di-(2-ethylhexyl)phosphoric Acid Extractant Systems from Experiments and Molecular Dynamics Simulations

Di-(2-ethylhexyl)phosphoric acid (D2EHPA) and tributylphosphate (TBP) are two of the most studied and researched organophosphorous extractants. D2EHPA is an acidic extractant, offering both hydrogen bond donor and acceptor sites while TBP, a neutral extractant, only offers a single acceptor site per molecule. In spite of this, it is observed that 1 M D2EHPA in dodecane is a poorer extractant for water than 1 M TBP in dodecane. The objective of present work is to look into the molecular interactions that cause such behavior. Experiments were carried out with varying molar ratios of TBP and D2EHPA in the organic dodecane phase. Total extractant concentration was kept constant at 1 M with dodecane as diluent. Water extraction was quantified by measuring the moisture content of the organic phase after equilibration. ¹H and ³¹P NMR spectra of the organic phase samples were recorded to study the change in the chemical environment upon extraction. Small angle X-ray scattering data of water saturated extractant phases were analyzed for the possibility of a reverse micellar aggregate formation. Molecular dynamics simulations could calculate free energies in quantitative agreement with experiments. Experimental and simulation studies showed that aggregation in the organic phase was promoted by the presence of water. This combined approach, of experiments and simulation, has shown that water is indispensable for the formation of ordered aggregates of extractants in nonpolar organic solvents. It is seen that, in the organic phase, around 80% of water's hydrogen bonds are with extractant molecules rather than with itself. The analysis clearly indicates that, rather than forming an aqueous core surrounded by extractant, water acts as a bridge between extractant molecules.

Surfactant-enhanced heterogeneity of the aqueous interface drives water extraction into organic solvents

Liquid/liquid extraction (LLE) is one of the most industrially relevant separations methods. Adsorbed surfactant is demonstrated to enhance interfacial heterogeneity and lead to water protrusions that form the basis for transport into the organic phase.

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.

Network analysis and percolation transition in hydrogen bonded clusters: nitric acid and water extracted by tributyl phosphate

Network analysis of hydrogen bonded clusters formed in simulation by extraction of nitric acid and water by TBP interprets cluster topologies and identifies the mechanism for third phase formation.

Quantifying Dimer and Trimer Formation by Tri-n-butyl Phosphates in n-Dodecane: Molecular Dynamics Simulations

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.