Organic and Third Phase in HNO3/TBP/n-Dodecane System: No Reverse Micelles

MDverse, shedding light on the dark matter of molecular dynamics simulations

The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation such as temperature and simulation length, and could identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore the MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and to report and standardize metadata to reuse this valuable matter.

Molecular-scale understanding of diluent effects on ligand assembly for metal ion separations

Use of metal-selective ligands in solvent extraction is instrumental in extraction of critical materials and recycling, yet, diluent effects on extraction performance are not well understood. Experimental and empirical solvent parameters have been proposed to correlate with extraction performance, but are often inadequate predictors. We follow the hypothesis that the diluents' primary influence on extraction efficiency is whether or not it hinders assembly of the bulky extracting ligands into a geometry necessary for metal complexation. This behavior is readily accessible with molecular dynamics (MD), where the atomistic description of molecules can be applied to arbitrary extractant-solvent molecules and their mixtures. Several simulated quantities are considered, from both pairwise and graph theoretical analyses, and compared to experimental distribution ratio data for americium extraction by TODGA in a series of inert, non-interacting diluents. These simple properties, especially the formation of closed triplets corresponding to the 3 : 1 ligand : metal stoichiometric solvate, suggest a potential predictive power of this approach. This methodology provides a path forward to comprehensively understand and predict diluent effects in more complex systems involving different extracting ligands and multi-component diluent mixtures.

MDverse: Shedding Light on the Dark Matter of Molecular Dynamics Simulations

The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2,000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation such as temperature and simulation length, and could identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore the MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and to report and standardize metadata to reuse this valuable matter.

Speciation and Organic Phase Structure in Nitric Acid Extraction with Trioctylamine

Understanding chemical speciation and intermolecular interactions in multicomponent liquids is essential to understanding their phase and chemical equilibria, which underpin chemical separation processes, including solvent extraction. Here we report on the extraction of nitric acid from its aqueous solutions into organic solutions of trioctylamine (TOA) in toluene, investigated with spectroscopic, X-ray scattering, and computational tools to understand molecular speciation in the organic phase and its relationship with the nanoscale structure of the organic phase. Trends in acid and water extraction clearly show two and three regimes, respectively, indicating different stoichiometric relationships, but speciation of HNO3, water, and amine in these regimes is not apparent. 1H NMR of the organic phase shows that there are at least two distinct acidic protons in the organic phase while ATR-FTIR results show that the organic phase with excess acid extraction is a mixture of trioctylammonium-nitrate ion pairs (TOAH·NO3), and undissociated HNO3 molecules. Comparison with DFT-computed IR spectra show that the chain-like configurations of TOAH·NO3·HNO3·H2O are favored over TOAH·NO3·H2O·HNO3, i.e., direct interaction between the nitrate and HNO3 molecules is more favored compared to a water-mediated interaction. SAXS of the organic phases were modeled as sums of Ornstein-Zernike (O-Z) scattering and a prepeak feature in the higher Q region that corresponds to extractant packing. The extraction of undissociated HNO3 by the ion pairs leads to an increased X-ray scattering contrast in the organic phase without any significant change in the correlation length. These results show that the organic phase nanostructure is more sensitive to the concentration of TOAH·NO3 and is relatively unaffected by excess acid extraction. These findings will enable a molecular understanding of the mechanisms behind metal extraction from acidic media with basic extractants.

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.

Molecular Forces in Liquid-Liquid Extraction

The phase transfer of ions is driven by gradients of chemical potentials rather than concentrations alone (i.e., by both the molecular forces and entropy). Extraction is a combination of high-energy interactions that correspond to short-range forces in the first solvation shell such as ion pairing or complexation forces, with supramolecular and nanoscale organization. While the latter are similar to the long-range solvent-averaged interactions in the colloidal world, in solvent extraction they are associated with lower characteristic lengths of the nanometric domain. Modeling of such complex systems is especially complicated because the two domains are coupled, whereas the resulting free energy of extraction is around k_BT to guarantee the reversibility of the practical process. Nevertheless, quantification is possible by considering a partitioning of space among the polar cores, interfacial film, and solvent. The resulting free energy of transfer can be rationalized by utilizing a combination of terms which represent strong complexation energies, counterbalanced by various entropic effects and the confinement of polar solutes in nanodomains dispersed in the diluent, together with interfacial extractant terms. We describe here this ienaics approach in the context of solvent extraction systems; it can also be applied to further complex ionic systems, such as membranes and biological interfaces.

A Telescoping View of Solute Architectures in a Complex Fluid System

Short- and long-range correlations between solutes in solvents can influence the macroscopic chemistry and physical properties of solutions in ways that are not fully understood. The class of liquids known as complex (structured) fluids-containing multiscale aggregates resulting from weak self-assembly-are especially important in energy-relevant systems employed for a variety of chemical- and biological-based purification, separation, and catalytic processes. In these, solute (mass) transfer across liquid-liquid (water, oil) phase boundaries is the core function. Oftentimes the operational success of phase transfer chemistry is dependent upon the bulk fluid structures for which a common functional motif and an archetype aggregate is the micelle. In particular, there is an emerging consensus that mass transfer and bulk organic phase behaviors-notably the critical phenomenon of phase splitting-are impacted by the effects of micellar-like aggregates in water-in-oil microemulsions. In this study, we elucidate the microscopic structures and mesoscopic architectures of metal-, water-, and acid-loaded organic phases using a combination of X-ray and neutron experimentation as well as density functional theory and molecular dynamics simulations. The key conclusion is that the transfer of metal ions between an aqueous phase and an organic one involves the formation of small mononuclear clusters typical of metal-ligand coordination chemistry, at one extreme, in the organic phase, and their aggregation to multinuclear primary clusters that self-assemble to form even larger superclusters typical of supramolecular chemistry, at the other. Our metrical results add an orthogonal perspective to the energetics-based view of phase splitting in chemical separations known as the micellar model-founded upon the interpretation of small-angle neutron scattering data-with respect to a more general phase-space (gas-liquid) model of soft matter self-assembly and particle growth. The structure hierarchy observed in the aggregation of our quinary (zirconium nitrate-nitric acid-water-tri-n-butyl phosphate-n-octane) system is relevant to understanding solution phase transitions, in general, and the function of engineered fluids with metalloamphiphiles, in particular, for mass transfer applications, such as demixing in separation and synthesis in catalysis science.

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.