Insights into water extraction and aggregation mechanisms of malonamide-alkane mixtures.

Structure at the nanoscale in the organic phase of liquid-liquid extraction systems is often tied to separation performance. However, the weak interactions that drive extractant assembly lead to poorly defined structures that are challenging to identify. Here, we investigate the mechanism of water extraction for a malonamide extractant commonly applied to f-element separations. We measure extractant concentration fluctuations in the organic phase with small angle X-ray scattering (SAXS) before and after contact with water at fine increments of extractant concentration, finding no qualitative changes upon water uptake that might suggest significant nanoscopic reorganization of the solution. The critical composition for maximum fluctuation intensity is consistent with small water-extractant adducts. The extractant concentration dependence of water extraction is consistent with a power law close to unity in the low concentration regime, suggesting the formation of 1 : 1 water-extractant adducts as the primary extraction mechanism at low concentration. At higher extractant concentrations, the power law slope increases slightly, which we find is consistent with activity effects modeled using Flory-Huggins theory without introduction of additional extractant-water species. Molecular dynamics simulations are consistent with these findings. The decrease in interfacial tension with increasing extractant concentration shows a narrow plateau region, but it is not correlated with any change in fluctuation or water extraction trends, further suggesting no supramolecular organization such as reverse micellization. This study suggests that water extraction in this system is particularly simple: it relies on a single mechanism at all extractant concentrations, and only slightly enhances the concentration fluctuations characteristic of the dry binary extractant/diluent mixture.

Influence of Aqueous Phase Acidity on Ln(III) Coordination by N,N,N',N'-Tetraoctyldiglycolamide.

This study highlights the importance of combining distribution ratio measurements with multiple spectroscopic techniques to provide a more comprehensive understanding of organic phase Ln coordination chemistry. Solvent extraction investigations with N,N,N',N'-tetraoctyldiglycolamide (TODGA) in n-heptane reveal the sensitivity of Ln complexation to the HNO3 concentration. Distribution ratio measurements in tandem with UV-Vis demonstrated that increasing the concentration of HNO3 above 0.5 M with a constant NO3- of 1 M increases the number of coordinating TODGA molecules, from a 1:2 to a 1:3 Ln:TODGA complex. At each concentration of HNO3 considered herein (from 0.01 to 1 M), Eu lifetime analysis demonstrated no evidence of H2O coordination. Results from Fourier transform infrared investigations suggest the presence of inner-sphere NO3- under low concentrations of HNO3 when the 1:2 Ln:TODGA complex is present. Increasing the HNO3 concentration above 0.5 M increases the propensity for outer-sphere interactions by removing the coordinated NO3- and saturating the Ln coordination sphere with three TODGA molecules, resulting in the well-established cationic, trischelate homoleptic [Ln(TODGA)3]3+ complex. This work demonstrates the importance in considering the NO3- source for solvent extraction systems. In particular, for systems with an affinity for outer-sphere interactions with molar concentrations of HNO3, changing the NO3- source can change the inner-sphere coordination of the Ln complex, which, in turn, affects the separation efficacy.

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.

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.

Solvent effects on extractant conformational energetics in liquid-liquid extraction: a simulation study of molecular solvents and ionic liquids.

Extractant design in liquid-liquid extraction (LLE) is a research frontier of metal ion separations that typically focuses on the direct extractant-metal interactions. However, a more detailed understanding of energetic drivers of separations beyond primary metal coordination is often lacking, including the role of solvent in the extractant phase. In this work, we propose a new mechanism for enhancing metal-complexant energetics with nanostructured solvents. Using molecular dynamics simulations with umbrella sampling, we find that the organic solvent can reshape the energetics of the extractant's intramolecular conformational landscape. We calculate free energy profiles of different conformations of a representative bidentate extractant, n-octyl(phenyl)-N,N-diisobutyl carbamoyl methyl phosphinoxide (CMPO), in four different solvents: dodecane, tributyl phosphate (TBP), and dry and wet ionic liquid (IL) 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([EMIM][Tf2N]). By promoting reorganization of the extractant molecule into its binding conformation, our findings reveal how particular solvents can ameliorate this unfavorable step of the metal separation process. In particular, the charge alternating nanodomains formed in ILs substantially reduce the free energy penalty associated with extractant reorganization. Importantly, using alchemical free energy calculations, we find that this stabilization persists even when we explicitly include the extracted cation. These findings provide insight into the energetic drivers of metal ion separations and potentially suggest a new approach to designing effective separations using a molecular-level understanding of solvent effects.

Using Metadynamics to Reveal Extractant Conformational Free Energy Landscapes.

Understanding the impact of extractant functionalization on metal-binding energetics in liquid-liquid extraction is essential to guide the development of better separation processes. Traditionally, computational extractant design uses electronic structure calculations on metal-ligand clusters to determine the metal-binding energy of the lowest energy state. Although highly accurate, this approach does not account for all of the relevant physics encountered under experimental conditions. Such methodologies often neglect entropic contributions such as temperature effects and ligand flexibility, in addition to approximating solvent-extractant interactions with implicit solvent models. In this study, we use classical molecular dynamics simulations with an advanced sampling method, metadynamics, to map out extractant molecule conformational free energies in the condensed phase. We generate the complete conformational landscape in solution for a family of bidentate malonamide-based extractants with different functionalizations of the headgroup and the side chains. In particular, we show how such alkyl functionalization reshapes the free energy landscape, affecting the free energy penalty of organizing the extractant into the cis-like metal-binding conformation from the trans-like conformation of the free extractant in solution. Specifically, functionalizing alkyl tails to the center of the headgroup has a greater influence on increasing molecular rigidity and disfavoring the binding conformation than functionalizing side chains. These findings are consistent with trends in metal-binding energetics based on experimentally reported distribution ratios. We also consider a different bidentate extractant molecule, carbamoylmethylphosphine oxide, and show how the choice of solvent can further reshape the conformational energetic landscape. This study demonstrates the feasibility of using molecular dynamics simulations with advanced sampling techniques to investigate extractant conformational energetics in solution, which, more broadly, will enable extractant design that accounts for entropic effects and explicit solvation.

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid-liquid mixtures.

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid-liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid-liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

Critical fluctuations in liquid-liquid extraction organic phases controlled by extractant and diluent molecular structure.

Extractant aggregation in liquid-liquid extraction organic phases impacts extraction energetics and is related to the deleterious efficiency-limiting liquid-liquid phase transition known as third phase formation. Using small angle X-ray scattering, we find that structural heterogeneities across a wide range of compositions in binary mixtures of malonamide extractants and alkane diluents are well described by Ornstein-Zernike scattering. This suggests that structure in these simplified organic phases originates from the critical point associated with the liquid-liquid phase transition. To confirm this, we measure the temperature dependence of the organic phase structure, finding critical exponents consistent with the 3D Ising model. Molecular dynamics simulations were also consistent with this mechanism for extractant aggregation. Due to the absence of water or any other polar solutes required to form reverse-micellar-like nanostructures, these fluctuations are inherent to the binary extractant/diluent mixture. We also show how the molecular structure of the extractant and diluent modulate these critical concentration fluctuations by shifting the critical temperature: critical fluctuations are suppressed by increasing extractant alkyl tail lengths or decreasing diluent alkyl chain lengths. This is consistent with how extractant and diluent molecular structure are known to impact metal and acid loading capacity in many-component LLE organic phases, suggesting phase behavior of practical systems may be effectively studied in simplified organic phases. Overall, the explicit connection between molecular structure, aggregation and phase behavior demonstrated here will enable the design of more efficient separations processes.

Thiocyanate Ions Form Antiparallel Populations at the Concentrated Electrolyte/Charged Surfactant Interface.

Anions play significant roles in the separation of lanthanides and actinides. The molecular-scale details of how these anions behave at aqueous interfaces are not well understood, especially at high ionic strengths. Here, we describe the interfacial structure of thiocyanate anions at a soft charged interface up to 5 M bulk concentration with combined classical and phase-sensitive vibrational sum frequency generation (PS-VSFG) spectroscopy and molecular dynamics (MD) simulations. At low concentrations thiocyanate ions are mostly oriented with their sulfur end pointing toward the charged surfactants. The VSFG signal reaches a plateau at around 100 mM bulk concentration, followed by significant changes above 1 M. At high concentrations a new thiocyanate population emerges with their sulfur end pointing toward the bulk liquid. The -CN stretch frequency is different for up and down oriented SCN- ions, indicating different coordination environments. These results provide key molecular-level insights for the interfacial behavior of complex anions in highly concentrated solutions.

Advancing Chemical Separations: Unraveling the Structure and Dynamics of Phase Splitting in Liquid-Liquid Extraction.

Liquid-liquid extraction (LLE), the go-to process for a variety of chemical separations, is limited by spontaneous organic phase splitting upon sufficient solute loading, called third phase formation. In this study we explore the applicability of critical phenomena theory to gain insight into this deleterious phase behavior with the goal of improving separations efficiency and minimizing waste. A series of samples representative of rare earth purification were constructed to include each of one light and one heavy lanthanide (cerium and lutetium) paired with one of two common malonamide extractants (DMDOHEMA and DMDBTDMA). The resulting postextraction organic phases are chemically complex and often form rich hierarchical structures whose statics and dynamics near the critical point were probed herein with small-angle X-ray scattering and high-speed X-ray photon correlation spectroscopy. Despite their different extraction behaviors, all samples show remarkably similar critical behavior with exponents well described by classical critical point theory consistent with the 3D Ising model, where the critical behavior is characterized by fluctuations with a single diverging length scale. This unexpected result indicates a significant reduction in relevant chemical parameters at the critical point, indicating that the underlying behavior of phase transitions in LLE rely on far fewer variables than are generally assumed. The obtained scalar order parameter is attributed to the extractant fraction of the extractant/diluent mixture, revealing that other solution components and their respective concentrations simply shift the critical temperature but do not affect the nature of the critical fluctuations. These findings point to an opportunity to drastically simplify studies of liquid-liquid phase separation and phase diagram development in general while providing insights into LLE process improvement.

Mesostructuring in Liquid-Liquid Extraction Organic Phases Originating from Critical Points.

Organic phase structure plays an important role in solute extraction energetics and phase behavior of liquid-liquid extraction (LLE) systems. For a binary extractant (amphiphile)/solvent mixture of relevance to LLE, we find that the organic phase mesostructuring is consistent with extractant concentration fluctuations as the compositional isotherm traverses the Widom line above its liquid-liquid critical point. This reveals a different mechanism for the well-documented heterogeneities in LLE organic phases that are typically attributed to micellization.

Impact of Water Extraction on Malonamide Aggregation: A Molecular Dynamics and Graph Theoretic Approach.

Solution structure in liquid-liquid extraction affects the efficacy of separation; however, even for simplified organic phases, structural characterization and attribution of aggregation to intermolecular interactions are fundamental challenges. We investigate water uptake into organic phases for two malonamides commonly applied to actinide and lanthanide separations. Extracted water induces reorganization of the amphiphilic extractant molecules, although we find this rearrangement is not strongly manifested in small-angle X-ray scattering making it challenging to probe without methods such as atomistic simulation. Using a graph theoretic approach to define hydrogen bonded water/malonamide aggregates from molecular dynamics simulations, we find evidence of a characteristic aggregate size by water number that results from geometric accommodation of the surrounding malonamide molecules. This implies a degree of size selectivity inherent to these water-in-oil aggregates. Conversely, we find no evidence of a characteristic size of the aggregates with respect to their malonamide number. By defining a separate graphical representation of self-association of the amphiphilic malonamides, we quantify how water affects the local and nonlocal topology of the malonamide network, providing a basis for characterization of the structure and impact of polar solutes in increasingly complex organic phases.

Cluster Identification Using Modularity Optimization to Uncover Chemical Heterogeneity in Complex Solutions.

Structural heterogeneity is commonly manifested in solutions and liquids that feature competition of different interparticle forces. Identifying and characterizing heterogeneity across different length scales requires multimodal experimental measurement and/or the application of new techniques for the interrogation of atomistic simulation data. Within the latter, the parsing of networks of interparticle interactions (chemical networks) has been demonstrated to be a valuable tool for identifying subensembles of chemical environments. However, chemical networks can adopt a wide variety of topologies that challenge generalizable methods for identifying heterogeneous behavior, and few network analysis algorithms have been proposed for multiscale resolution. In this study, we apply a method of partitioning using the graph theoretic concept of clusters and communities. Using a modularity optimization algorithm, the cluster partition creates subgraphs based on their relative internal and external connectivities. The methodology is tested on two soft matter systems that have significantly different network topologies so as to probe its ability to identify multiple scale features and its generalizability. A binary Lennard-Jones fluid is first examined, where one component causes subgraphs that have high internal network connectivity yet are still connected to the rest of the interparticle network of interactions. The impact of connectivity and edge weighting on the cluster partition is investigated. In the second system, hierarchically organized molecular structures comprised of hydrogen bonded water molecules are identified at a liquid/liquid interface. These structures have a much more sparse network with significantly varied internal connectivity that is a challenge to differentiate from the background hydrogen bonding network of water molecules at the instantaneous interface. The organized macrostructures are effectively isolated from the background network using the cluster partition, and a time-dependent implementation allows us to reveal their reactivity. These studies indicate that cluster partitioning based upon intermolecular network connectivity patterns is broadly generalizable, depending only on user-defined intermolecular connectivity, is operable across different length scales, and is extensible to the study of dynamic phenomena.

Molecular-scale origins of solution nanostructure and excess thermodynamic properties in a water/amphiphile mixture.

The molecular and nanoscale origins of nonideality in excess thermodynamic properties are essential to understanding cosolvent mixtures, yet they remain challenging to determine. Here, we consider a binary mixture of water and an amphiphile, N,N,N',N'-tetramethylmalonamide (TMMA), which is characterized by strong hydrogen bonding between the two components and no hydrogen bonding between amphiphiles. Using molecular dynamics simulation, validated with excess volume measurements and X-ray scattering, we identify three distinct solution regimes across the composition range of the binary mixture and find that the transition between two of these regimes, marked by the water percolation threshold, is closely correlated with minima in the excess volume and excess enthalpy. Structural analysis of the simulations reveals an interplay between local interactions and solution nanostructure, determined by the relative strength of the water-water and water-amphiphile hydrogen bonding interactions. By comparison with other amphiphiles, such as linear alcohols, the relative strength of like and unlike interactions between water and amphiphile affects the relationship between thermodynamics and structural regimes. This provides insight into how molecular forces of mutual solvation interact across length scales and how they manifest in excess thermodynamic properties.

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.

Amphiphile Organization in Organic Solutions: An Alternative Explanation for Small-Angle X-ray Scattering Features in Malonamide/Alkane Mixtures.

The role of different intermolecular interactions in the aggregation of amphiphiles in an organic solvent is studied for systems of relevance to liquid-liquid extraction (LLE), a chemical process used to selectively recover metals from complex mixtures. Of specific interest is the role, or lack thereof, of hydrogen bonding, which is often assumed to be a main driver of the organic phase structural organization that has been linked to separation efficacy. Toward that end, a series of malonamide extractants in n-dodecane have been studied in the absence of any extracted aqueous solutes, including water. The series of extractants includes N,N'-dimethyl-N,N'-dibutyltetradecylmalonamide (DMDBTDMA), two of its homologs, and N,N'-dimethyl-N,N'-dioctylhexylethoxymalonamide (DMDOHEMA). This simplified model LLE system enables systematic investigation of the role of dipole-dipole and alkyl tail steric interactions in amphiphile aggregation. Small-angle X-ray scattering (SAXS) profiles computed from molecular dynamics trajectories are in good agreement with experimental SAXS data. Molecular dynamics simulations show that malonamide aggregation results from dipole-driven self-association and lacks characteristic aggregate sizes. Mid-q correlation peaks in the SAXS profiles emerge at high concentration for each malonamide. In those densely packed solutions, the correlation peaks are observed to result from alkyl tail-induced spacing between electron-rich polar head groups, with peak positions determined by the different alkyl tail lengths present in the malonamide molecule. This explanation of the SAXS correlation peaks contrasts with the prevailing literature, which attributes mesoscale features observed in small-angle scattering to the formation of microemulsions. Instead, this work finds that these features are present in the absence of water or any reverse micellar organization of the malonamides. As such, molecular-scale malonamide self-association and packing, rather than microemulsion-based colloidal-scale descriptions, is a more appropriate framework for these LLE systems.

A Geometric Measure Theory Approach to Identify Complex Structural Features on Soft Matter Surfaces.

The structural features that protrude above or below a soft matter interface are well-known to be related to interfacially mediated chemical reactivity and transport processes. It is a challenge to develop a robust algorithm for identifying these organized surface structures, as the morphology can be highly varied and they may exist on top of an interface containing significant interfacial roughness. A new algorithm that employs concepts from geometric measure theory, algebraic topology, and optimization is developed to identify candidate structures at a soft matter surface, and then, using a probabilistic approach, to rank their likelihood of being a complex structural feature. The algorithm is tested for a surfactant laden water/oil interface, where it is robust to identifying protrusions responsible for water transport against a set identified by visual inspection. To our knowledge, this is the first example of applying geometric measure theory to analyze the properties of a chemical/materials science system.

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.

Effect of Ions on the Optical Absorption Spectra of Aqueously Solvated Chromophores.

In the condensed phase, ions often create heterogeneous local environments around a solute, which may impart chemical reactivity or perturbations to physicochemical properties. Although the former has been the subject of some study, the latter-particularly as is pertains to optical absorption spectroscopy-is much less understood. In this work, the computed UV-vis absorption spectrum is examined for the aqueously solvated chromophore anion of green fluorescent protein for different local ion configurations. The strong ability of water to screen the ions from the chromophore results in little change in excitation energy compared to a purely aqueous environment. However, upon forming a contact ion pair with a sodium ion at either of the two electronegative oxygen sites of the chromophore, there is a spectral shift to either higher or lower energies. Surprisingly, our analysis suggests that the cause of the spectral shift is dominated not by the electrostatic presence of the ion but instead by ion disruption of the hydrogen bond network at the oxygen contact ion pair site.