TBP at the Water−Oil Interface: The Effect of TBP Concentration and Water Acidity Investigated by Molecular Dynamics Simulations

Mineral Acid Co-Extraction in Reactive Extraction of Lactic Acid Using a Thymol-Menthol Deep Eutectic Solvent as a Green Modifier

Carboxylic acids can be isolated from fermentation broths using reactive liquid-liquid extraction, offering an alternative to the environmentally harmful state-of-the-art process of precipitating calcium lactate. To enhance the sustainability of liquid-liquid extraction processes, greener solvents, such as natural deep eutectic solvents, are investigated. However, fermentation broths often exhibit pH values unsuitable for carboxylic acid extraction, which can be adjusted using mineral acids, though mineral acids may be co-extracted. In this study, we systematically examine the co-extraction of hydrochloric, nitric, sulfuric, and phosphoric acid during extraction and back-extraction of lactic acid. The solvent phase consisted of tri-n-octylamine, trioctylphosphine oxide, or tributyl phosphate diluted in a thymol-menthol deep eutectic solvent. The back-extraction was conducted using a diluent swing with p-cymene as the antisolvent and water as the receiving phase. Tri-n-octylamine showed the highest efficiency for lactic acid (up to 29.8%) but also the highest co-extraction of mineral acids (up to 50.9%). In contrast, trioctylphosphine oxide exhibited a lower but more selective lactic acid extraction (5.94%) with low mineral acids co-extraction (0.135%). Overall, the highest co-extraction was observed for phosphoric acid and the lowest for nitric acid. In conclusion, the selected solvent phase composition and mineral acid influence the co-extraction and, thus, final product purity. The successful application of the natural deep eutectic solvent as the modifier enhances the sustainability of liquid-liquid extraction processes.

Structure and dynamics of dissociated and undissociated forms of nitric acid and their implications in interfacial mass transfer: insights from molecular dynamics simulations

Nitric acid (HNO3) is widely used in various chemical and nuclear industries. Therefore, it is important to develop an understanding of the different forms of nitric acid for its practical applications. Molecular dynamics (MD) simulation is one of the best tools to investigate the behavior of concentrated nitric acid in aqueous solution with various forms together with pure nitric acid to identify a suitable model of nitric acid for use in simulations of biphasic systems for interfacial mass transfer. The Mulliken partial charge embedded OPLS-AA force field was used to model the neutral nitric acid, hydronium ion and nitrate ion, and it was found that the Mulliken partial charge embedded force field works quite well. The computed density of the dissociated and mixed-form acid was in good agreement with the experimental values. In water, the HNO3 molecule was seen to be coordinated with three water molecules in the first sphere of coordination. The distribution of water surrounding the HNO3 molecule and nitrate ion was corroborated by the DFT-optimized hydrated cluster. The calculated diffusivity values of the neutral acid and ions were significantly higher in the mixed form of nitric acid, which is an important dynamic quantity controlling the kinetics of the liquid-liquid interfacial extraction. The structural analysis revealed that the local aggregation is minimized when both forms of acid are present together in the solution. The water-ion and water-neutral acid interactions were predicted to be enhanced, as confirmed by H-bond studies. The shear viscosity of the mixed acid exhibited excellent agreement with the experimental values, which again confirms the consideration of the mixed form of nitric acid. The simulated value of surface tension for the mixed form of acid also appeared to be quite accurate based on the surface tension of water. The mixed form of nitric acid comprising both forms of acid is the best representation for nitric acid to be considered for MD simulations of biphasic systems. The mixed form of nitric acid established that the concentrated nitric acid may not be present either in the fully dissociated form or fully undissociated form in the solution.

Aqueous Interfaces in Chemical Separations

Chemical separations play a vital role in refinery and reprocessing of critical materials, such as platinum group metals, rare earths, and actinides. The choice of separation system─whether it is liquid-liquid extraction (LLE), sorbents, or membranes─depends on specific needs and applications. In almost all separation processes, the desired metal ions adsorb or transfer across an aqueous interface, such as the solid/liquid interface in sorbents or oil/water interfaces in LLE. Despite these separation technologies being extensively used for decades, our understanding of the molecular-scale mechanisms governing ion adsorption and transport at interfaces remains limited. This knowledge gap presents a significant challenge in meeting the increasing demands for these critical materials due to their growing use in advanced technologies. Fortunately, recent advancements in surface-specific experimental and computational techniques offer promising avenues to bridge this gap and facilitate the development of next-generation separation systems. Interestingly, unanswered questions regarding interfacial phenomena in chemical separations hold great relevance to various fields, including energy storage, geochemistry, and atmospheric chemistry. Therefore, the model interfacial systems developed for studying chemical separations, such as amphiphilic molecules assembled at a solid/water, air/water, or oil/water interface, may have far-reaching implications, extending beyond separations and opening doors to addressing a wide range of scientific inquiries. This perspective discusses recent interfacial studies elucidating amphiphile-ion interactions in chemical separations of metal ions. These studies provide direct, molecular-scale information about solute and solvent behavior at aqueous interfaces, including multivalent and complex ions in highly concentrated solutions, which play key roles in LLE of critical materials.

Uracil/H+ Symport by FurE Refines Aspects of the Rocking-bundle Mechanism of APC-type Transporters

Transporters mediate the uptake of solutes, metabolites and drugs across the cell membrane. The eukaryotic FurE nucleobase/H+ symporter of Aspergillus nidulans has been used as a model protein to address structure-function relationships in the APC transporter superfamily, members of which are characterized by the LeuT-fold and seem to operate by the so-called 'rocking-bundle' mechanism. In this study, we reveal the binding mode, translocation and release pathway of uracil/H+ by FurE using path collective variable, funnel metadynamics and rational mutational analysis. Our study reveals a stepwise, induced-fit, mechanism of ordered sequential transport of proton and uracil, which in turn suggests that FurE, functions as a multi-step gated pore, rather than employing 'rocking' of compact domains, as often proposed for APC transporters. Finally, our work supports that specific residues of the cytoplasmic N-tail are involved in substrate translocation, in line with their essentiality for FurE function.

Recent Advances in the Study of Trivalent Lanthanides and Actinides by Phosphinic and Thiophosphinic Ligands in Condensed Phases

The separation of trivalent actinides and lanthanides is a key step in the sustainable development of nuclear energy, and it is currently mainly realized via liquid-liquid extraction techniques. The underlying mechanism is complicated and remains ambiguous, which hinders the further development of extraction. Herein, to better understand the mechanism of the extraction, the contributing factors for the extraction are discussed (specifically, the sulfur-donating ligand, Cyanex301) by combing molecular dynamics simulations and experiments. This work is expected to contribute to improve our systematic understanding on a molecular scale of the extraction of lanthanides and actinides, and to assist in the extensive studies on the design and optimization of novel ligands with improved performance.

Molecular dynamics simulation study of adsorption of anionic–nonionic surfactants at oil/water interfaces

Four anionic–nonionic surfactants with the same headgroups and different units of oxygen ethyl (EO) and oxygen propyl (PO) were adopted to investigate the influence on oil/water interfacial tensions in this article. Molecular dynamics (MD) simulations were conducted to study the interfacial property of the four surfactants. Four parameters were proposed to reveal the effecting mechanism of molecular structure on interfacial tension, which included the interfacial thickness, order parameter of the hydrophobic chain, radial distribution function, and the solvent accessible surface area. In addition, the electrostatic potential of the four surfactants was calculated. The research results indicated that the interface facial mask formed by the surfactants, which contained three EO or three PO units was more stable, and it was easier for the surfactants of six EO or six PO units to form a microemulsion at higher concentrations. The adsorption mechanism of the anionic–nonionic surfactant systems at the oil/water interfaces was supplemented at a molecular level, which provided fundamental guidance for an in-depth understanding of the optimal selection of the surfactants in enhancing oil recovery.

Four anionic–nonionic surfactants with the same headgroups and different units of oxygen ethyl (EO) and oxygen propyl (PO) were adopted to investigate the influence on oil/water interfacial tensions in this article.

Unexpected inverse correlations and cooperativity in ion-pair phase transfer

Liquid/liquid extraction is one of the most widely used separation and purification methods, where a forefront of research is the study of transport mechanisms for solute partitioning and the relationships that these have to solution structure at the phase boundary. To date, organized surface features that include protrusions, water-fingers, and molecular hinges have been reported. Many of these equilibrium studies have focused upon small-molecule transport – yet the extent to which the complexity of the solute, and the competition between different solutes, influence transport mechanisms have not been explored. Here we report molecular dynamics simulations that demonstrate that a metal salt (LiNO₃) can be transported via a protrusion mechanism that is remarkably similar to that reported for H₂O by tri-butyl phosphate (TBP), a process that involves dimeric assemblies. Yet the LiNO₃ out-competes H₂O for a bridging position between the extracting TBP dimer, which in-turn changes the preferred transport pathway of H₂O. Examining the electrolyte concentration dependence on ion-pair transport unexpectedly reveals an inverse correlation with the extracting surfactant concentration. As [LiNO₃] increases, surface adsorbed TBP becomes a limiting reactant in correlation with an increased negative surface charge induced by excess interfacial NO₃⁻, however the rate of transport is enhanced. Within the highly dynamic interfacial environment, we hypothesize that this unique cooperative effect may be due to perturbed surface organization that either decreases the energy of formation of transporting protrusion motifs or makes it easier for these self-assembled species to disengage from the surface.

A forefront of research in separations science (specifically liquid–liquid extraction) is the study of transport mechanisms for solute partitioning, and the relationships that these have to solution structure at the phase boundary.

Using Microemulsions: Formulation Based on Knowledge of Their Mesostructure

Microemulsions, as thermodynamically stable mixtures of oil, water, and surfactant, are known and have been studied for more than 70 years. However, even today there are still quite a number of unclear aspects, and more recent research work has modified and extended our picture. This review gives a short overview of how the understanding of microemulsions has developed, the current view on their properties and structural features, and in particular, how they are related to applications. We also discuss more recent developments regarding nonclassical microemulsions such as surfactant-free (ultraflexible) microemulsions or ones containing uncommon solvents or amphiphiles (like antagonistic salts). These new findings challenge to some extent our previous understanding of microemulsions, which therefore has to be extended to look at the different types of microemulsions in a unified way. In particular, the flexibility of the amphiphilic film is the key property to classify different microemulsion types and their properties in this review. Such a classification of microemulsions requires a thorough determination of their structural properties, and therefore, the experimental methods to determine microemulsion structure and dynamics are reviewed briefly, with a particular emphasis on recent developments in the field of direct imaging by means of electron microscopy. Based on this classification of microemulsions, we then discuss their applications, where the application demands have to be met by the properties of the microemulsion, which in turn are controlled by the flexibility of their amphiphilic interface. Another frequently important aspect for applications is the control of the rheological properties. Normally, microemulsions are low viscous and therefore enhancing viscosity has to be achieved by either having high concentrations (often not wished for) or additives, which do not significantly interfere with the microemulsion. Accordingly, this review gives a comprehensive account of the properties of microemulsions, including most recent developments and bringing them together from a united viewpoint, with an emphasis on how this affects the way of formulating microemulsions for a given application with desired properties.

Interaction of epoxy-based hydrogels and water: A molecular dynamics simulation study

Biomaterials play a crucial role in tissue engineering as a functional replacement, regenerative medicines, supportive scaffold for guided tissue growth, and drug delivery devices. The term biomaterial refers to metals, ceramics, and polymers account for the vast majority. In the case of polymers, hydrogels have emerged as active materials for an immense variety of applications. Epoxy-based hydrogels possess a unique network structure that enables very high levels of hydrophilicity and biocompatibility. Hydrogel such as Medipacs Epoxy Polymers (MEPs) models were constructed to understand water's behavior at the water/hydrogel interface and hydrogel network. We computed the Gibbs dividing surface (GDS) to define the MEP/water interface, and all the physicochemical properties were computed based on GDS. We calculated the radial distribution function (RDF), the 2D surface roughness of the immersed MEPs. RDF analysis confirmed that the first hydration shell is at a distance of 1.86 Å, and most of the water molecules are near the hydroxyl group of the MEPs network. Hydrogen bonds (H-bonds) analysis was performed, and the observation suggested that the disruption of the H-bonds between MEP chains leads to an increase in the polymer matrix's void spaces. These void spaces are filled with diffused water molecules, leading to swelling of the MEP hydrogel. The swelling parameter was estimated from the fitted curve of the yz-lattice of the simulation cell. The MEP/water interface simulation results provide insightful information regarding the design strategy of epoxy-based hydrogel and other hydrogels vital for biomedical applications.

An octanol hinge opens the door to water transport

Despite their prevalent use as a surrogate for partitioning of pharmacologically active solutes across lipid membranes, the mechanism of transport across water/octanol phase boundaries has remained unexplored. Using molecular dynamics, graph theoretical, cluster analysis, and Langevin dynamics, we reveal an elegant mechanism for the simplest solute, water. Self-assembled octanol at the interface reversibly binds water and swings like the hinge of a door to bring water into a semi-organized second interfacial layer (a “bilayer island”). This mechanism is distinct from well-known lipid flipping and water transport processes in protein-free membranes, highlighting important limitations in the water/octanol proxy. Interestingly, the collective and reversible behavior is well-described by a double well potential energy function, with the two stable states being the water bound to the hinge on either side of the interface. The function of the hinge for transport, coupled with the underlying double well energy landscape, is akin to a molecular switch or shuttle that functions under equilibrium and is driven by the differential free energies of solvation of H₂O across the interface. This example successfully operates within the dynamic motion of instantaneous surface fluctuations, a feature that expands upon traditional approaches toward controlled solute transport that act to avoid or circumvent the dynamic nature of the interface.

Despite their pharmacological relevance, the mechanism of transport across water/octanol phase boundaries has remained unexplored. Octanol molecular assemblies are demonstrated to reversibly bind water and swing like the hinge of a door.

Understanding the Complex Surface Interplay for Extraction: A Molecular Dynamics Study

By means of classical molecular dynamics simulation the interfacial properties of methanol and n‐dodecane, which are two potential candidate solvents for use in non‐aqueous liquid–liquid extraction, were assessed. The question of how the interface changes depending on the concentration of extractant (tri‐n‐butyl phosphate) and salt (LiCl) is addressed. Two different models to represent systems were used to evaluate how LiCl and tri‐n‐butyl phosphate affect mutual miscibility, and how the last‐named behaves depending on the chemical environment. Tri‐n‐butyl phosphate increases the mutual solubility of the solvents, whereas LiCl counteracts it. The extractant was found to be mostly adsorbed on the interface between the solvents, and therefore the structural features of the adsorption were investigated. Adsorption of tri‐n‐butyl phosphate changes depending on its concentration and the presence of LiCl. It exhibits a preferential orientation in which the butyl chains point at the n‐dodecane phase and the phosphate group points at the methanol phase. For high concentrations of tri‐n‐butyl phosphate, its molecular orientation is preserved by diffusion of the excess molecules into both the methanol and n‐dodecane phases. However, LiCl hinders the diffusion into the methanol phase, and thus increases the concentration of tri‐n‐butyl phosphate at the interface and forces a rearrangement with subsequent loss of orientation.

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.

Influence of Light and Temperature on the Extractability of Cerium(IV) as a Surrogate of Plutonium(IV) and its Effect on the Simulation of an Accidental Fire in the PUREX Process

Modeling of plutonium(IV) behavior during an accidental fire in a reprocessing plant was considered using various non-radioactive metallic surrogates. Among those elements, cerium(IV) was supposed to be a suitable candidate due to possible formation of a complex with TPB, but its extractability and stability have not been studied previously under representative plutonium uranium reduction extraction (PUREX) conditions. In this work, we investigated the chemical analogy between cerium(IV) and plutonium(IV) in this extractive process and combustion thereof. Distribution ratios are reported for acidities of 1-4 mol L^-1 in equal volumes of nitric acid and a 30:70 mixture of tributylphosphate and hydrogenated tetrapropylene. The influences of light, temperature, and extraction time were studied by UV-vis spectroscopy. The results showed that cerium(IV) is extracted quantitatively but is reduced over time to cerium(III) in the organic mixture. Spectrophotometric investigations of this reaction kinetics revealed an apparent rate constant k of 0.021 ± 0.002 mol^0.5 L^0.5 min^-1 at 298 K and an apparent fractional reaction order of 0.5. The activation energy of this reduction was found to be around 82 ± 2 kJ mol^-1 by the Arrhenius plot method. The combustion of mono- and biphasic solutions prepared with a cerium(IV) concentration of 10 g L^-1 revealed that the extracted complexes, Ce₂O·6NO₃·3TBP_(org) or Ce₄O₄·8NO₃·6TBP_(org), are reduced during the combustion. Compositions of the resulting ashes and soot were analyzed and highlighted the presence of pyrophosphates and polycyclic aromatic hydrocarbons, with some traces of cerium. Ce(IV) is not suitable to represent Pu(IV) from a chemical point of view in HNO₃/TBP-HTP solutions.

Mechanisms of phase separation in temperature-responsive acidic aqueous biphasic systems

The thermal and acid responsive behaviour of bulky phosphonium-based ILs is elucidated using a mixed experimental and computational approach.

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.

A High-Efficiency Approach for the Synthesis of N235-Impregnated Resins and the Application in Enhanced Adsorption and Separation of Vanadium(V)

Trialkylamine (N235)-tributyl phosphate (TBP) impregnated resins (N-TIRs) were prepared, so as to evaluate the effects of the addition of TBP on the preparation and adsorption performance of N235-impregnated resins (NIRs). The results show that TBP can obviously increase the impregnation ratio and shorten the impregnation equilibrium time of the N-TIRs when compared to that of the NIRs (57.73% versus 36.95% and 5 min versus 240 min). It is confirmed that TBP can interact with N235 during the impregnation process, which shorten the adsorption equilibrium time and increases the adsorption capacity of the N-TIRs for V(V) when compared to that of the NIRs (6 h versus 10 h and 50.95 mg·g−1 versus 46.73 mg·g−1). The kinetics fitting results demonstrate that the adsorption of V(V) onto N-TIRs and NIRs all conform to pseudo-second order kinetic model and chemical reaction is the rate-limiting step of the whole adsorption process. In the meanwhile, the reaction constant (Ks) implies that the chemical reaction rate of V(V) with the impregnated extractants in N-TIRs is faster than that in NIRs. The N-TIRs present higher stability and selectivity than NIRs. This study manifests that the addition of a secondary reagent may be a potential and novel technique on the preparation of SIRs and the enhancement of adsorption and separation for ions.

Comprehensive Peptide Ion Structure Studies Using Ion Mobility Techniques: Part 3. Relating Solution-Phase to Gas-Phase Structures

Molecular dynamics (MD) simulations have been utilized to study peptide ion conformer establishment during the electrospray process. An explicit water model is used for nanodroplets containing a model peptide and hydronium ions. Simulations are conducted at 300 K for two different peptide ion charge configurations and for droplets containing varying numbers of hydronium ions. For all conditions, modeling has been performed until production of the gas-phase ions and the resultant conformers have been compared to proposed gas-phase structures. The latter species were obtained from previous studies in which in silico candidate structures were filtered according to ion mobility and hydrogen-deuterium exchange (HDX) reactivity matches. Results from the present study present three key findings namely (1) the evidence from ion production modeling supports previous structure refinement studies based on mobility and HDX reactivity matching, (2) the modeling of the electrospray process is significantly improved by utilizing initial droplets existing below but close to the calculated Rayleigh limit, and (3) peptide ions in the nanodroplets sample significantly different conformers than those in the bulk solution due to altered physicochemical properties of the solvent. Graphical Abstract ᅟ.

Self-assembled highly ordered acid layers in precisely sulfonated polyethylene produce efficient proton transport

Recent advances in polymer synthesis have allowed remarkable control over chain microstructure and conformation. Capitalizing on such developments, here we create well-controlled chain folding in sulfonated polyethylene, leading to highly uniform hydrated acid layers of subnanometre thickness with high proton conductivity. The linear polyethylene contains sulfonic acid groups pendant to precisely every twenty-first carbon atom that induce tight chain folds to form the hydrated layers, while the methylene segments crystallize. The proton conductivity is on par with Nafion 117, the benchmark for fuel cell membranes. We demonstrate that well-controlled hairpin chain folding can be utilized for proton conductivity within a crystalline polymer structure, and we project that this structure could be adapted for ion transport. This layered polyethylene-based structure is an innovative and versatile design paradigm for functional polymer membranes, opening doors to efficient and selective transport of other ions and small molecules on appropriate selection of functional groups.

Influences of Cation Ratio, Anion Type, and Water Content on Polytypism of Layered Double Hydroxides

Layered double hydroxides (LDHs) are a significant sink of anions (CO₃^2-, SO₄^2-, NO₃^-, Cl^-, etc.) and divalent transition-metal cations in soil. The anion exchange capacity gives rise to functional materials. The stability of LDHs is determined by the interaction between cation-bearing layers and intercalated water and anions, which is correlated with polytypism and coordination structure. A systematic investigation is performed to show the influence of cation ratio, anion type, and water content on polytypism, swelling behavior, and interlayer structure of Mg-Al-LDHs using molecular dynamics simulations. LDHs intercalated with NO₃^- ions exhibit a polytype transition from 3 R₁ (three-layer rhombohedral polytype) to 1 T (one-layer trigonal polytype) with increasing water content. NO₃^- ions exhibit a D_{3 h} point group symmetry at low water contents. The polytype transition coincides with the complete transformation into tilted NO₃^- ion with a C_{2 v} point group symmetry. The transition appears at a lower water content when the Mg/Al ratio is lower. LDHs with SO₄^2- ions exhibit a three-stage polytypism. The first and last stages are 3 R₁. The intermediate stage could be 1 T or a mixture of different O(octahedra)-type interlayers, which depends on the cation ratio. The relative popularity of SO₄^2- ions with a C _s point group symmetry is characteristic for the intermediate stage, while mostly SO₄^2- ions exhibit a C_{3 v} symmetry. There is no clear relevance between cation ratio and water content at which a polytype transition happens. The configurational adjustments of NO₃^- and SO₄^2- ions facilitate the swelling behavior of LDHs. LDHs with CO₃^2- or Cl^- ions always maintain a 3 R₁ polytype irrespective of water content and hardly swell. The configurations of anions and water reflect local coordination structure due to hydrogen bonds. The layer-stacking way influences long-ranged Coulombic interactions. Hydrogen-bonding structure and long-ranged Coulombic interactions collectively determine polytypism and stability of LDHs.

Molecular dynamics simulation for the test of calibrated OPLS-AA force field for binary liquid mixture of tri-iso-amyl phosphate and n-dodecane

Tri-isoamyl phosphate (TiAP) has been proposed to be an alternative for tri-butyl phosphate (TBP) in the Plutonium Uranium Extraction (PUREX) process. Recently, we have successfully calibrated and tested all-atom optimized potentials for liquid simulations using Mulliken partial charges for pure TiAP, TBP, and dodecane by performing molecular dynamics (MD) simulation. It is of immense importance to extend this potential for the various molecular properties of TiAP and TiAP/n-dodecane binary mixtures using MD simulation. Earlier, efforts were devoted to find out a suitable force field which can explain both structural and dynamical properties by empirical parameterization. Therefore, the present MD study reports the structural, dynamical, and thermodynamical properties with different mole fractions of TiAP-dodecane mixtures at the entire range of mole fraction of 0-1 employing our calibrated Mulliken embedded optimized potentials for liquid simulation (OPLS) force field. The calculated electric dipole moment of TiAP was seen to be almost unaffected by the TiAP concentration in the dodecane diluent. The calculated liquid densities of the TiAP-dodecane mixture are in good agreement with the experimental data. The mixture densities at different temperatures are also studied which was found to be reduced with temperature as expected. The plot of diffusivities for TiAP and dodecane against mole fraction in the binary mixture intersects at a composition in the range of 25%-30% of TiAP in dodecane, which is very much closer to the TBP/n-dodecane composition used in the PUREX process. The excess volume of mixing was found to be positive for the entire range of mole fraction and the excess enthalpy of mixing was shown to be endothermic for the TBP/n-dodecane mixture as well as TiAP/n-dodecane mixture as reported experimentally. The spatial pair correlation functions are evaluated between TiAP-TiAP and TiAP-dodecane molecules. Further, shear viscosity has been computed by performing the non-equilibrium molecular dynamics employing the periodic perturbation method. The calculated shear viscosity of the binary mixture is found to be in excellent agreement with the experimental values. The use of the newly calibrated OPLS force field embedding Mulliken charges is shown to be equally reliable in predicting the structural and dynamical properties for the mixture without incorporating any arbitrary scaling in the force field or Lennard-Jones parameters. Further, the present MD simulation results demonstrate that the Stokes-Einstein relation breaks down at the molecular level. The present methodology might be adopted to evaluate the liquid state properties of an aqueous-organic biphasic system, which is of great significance in the interfacial science and technology.

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.

Ion Transport Mechanisms in Liquid-Liquid Interface

Interfacial liquid-liquid ion transport is of crucial importance to biotechnology and industrial separation processes including nuclear elements and rare earths. A water-in-oil microemulsion is formulated here with density and dimensions amenable to atomistic molecular dynamics simulation, facilitating convergent theoretical and experimental approaches to elucidate interfacial ion transport mechanisms. Lutetium(III) cations are transported from the 5 nm diameter water pools into the surrounding oil using an extractant (a lipophilic ligand). Changes in ion coordination sphere and interactions between the interfacial components are studied using a combination of synchrotron X-ray scattering, spectroscopy, and atomistic molecular dynamics simulations. Contrary to existing hypotheses, our model system shows no evidence of interfacial extractant monolayers, but rather ions are exchanged through water channels that penetrate the surfactant monolayer and connect to the oil-based extractant. Our results highlight the dynamic nature of the oil-water interface and show that lipophilic ion shuttles need not form flat monolayer structures to facilitate ion transport across the liquid-liquid interface.

TBP Assisted Uranyl Extraction in Water-Dodecane Biphasic System: Insights from Molecular Dynamics Simulation

In the PUREX (Plutonium Uranium Recovery by Extraction Process) process, the extraction of uranyl ion from dissolver solution to the organic phase is influenced by co extraction of the other species, such as water and nitric acid and it is assumed that the presence of water or acid droplets in the organic phase intensifies the coordination mechanism of TBP. The present study illustrates the uranyl extraction from the aqueous phase to the organic phase using molecular dynamics (MD) simulation. Here, we consider the biphasic systems to gain insights into the characteristics of the interface and humidity of the organic phase under different acidic and neutral conditions. MD being a force field method, can’t satisfactorily model the bond making and breaking process therefore a priori choice has been made concerning the different status of proton for the acidic phase. Further, the importance of charge species transferability during uranyl-TBP complexation have been investigated considering two different models of uranyl nitrate; united UO2(NO3)2 complex and separate UO2 2+ and NO3 – ions. From the results, it is recommended to use the ionic uranyl model with separate UO2 2+ and NO3 – to study the structural and dynamical properties of extracted uranyl ions in the organic phase. Also, it was noticed that extracted uranyl ions in the organic phase are not completely dehydrated but are surrounded by water molecules. In other words the results show co extraction of other species such as water and acid molecules to the organic phase. Most remarkably, the present study evident that the neutral HNO3 effectively represents the acidity effect for the receiving phase in terms of acid/water extraction and their aggregation to form water droplet, especially when ionic model of uranyl nitrate is considered.

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.