Reverse Aggregates as Adaptive Self‐Assembled Systems for Selective Liquid‐Liquid Cation Extraction

Hierarchical Aggregation in a Complex Fluid─The Role of Isomeric Interconversion

There is an ever-increasing body of evidence that metallic complexes involving amphiliphic ligands do not form normal solutions in organic solvents. Instead, they form complex fluids with intricate structures. For example, the metallic complexes may aggregate into clusters, and these clusters themselves may aggregate into superclusters. To gain a deeper insight into the mechanisms at play, we have used an improved force field to conduct extensive molecular dynamics simulations of a system composed of zirconium nitrate, water, nitric acid, tri-n-butyl phosphate, and n-octane. The important new finding is that a dynamic equilibrium between the cis and trans isomers of the metal complex is likely to play a key role in the aggregation behavior. The isolated cis and trans isomers have similar energies, but simulation indicates that the clusters consist predominantly of cis isomers. With increasing metal concentration, we hypothesize that more clustering occurs and the chemical equilibrium shifts toward the cis isomer. It is possible that such isomeric effects play a role in the liquid-liquid extraction of other species and the inclusion of such effects in flow sheet modeling may lead to a better description of the process.

Molecular Dynamics Study of the Aggregation Behavior of N,N,N',N'-Tetraoctyl Diglycolamide

Liquid–liquid extraction is a commonly used technique to separate metals and is a process that has particular relevance to the nuclear industry. There has been a drive to use environmentally friendly ligands composed only of carbon, hydrogen, nitrogen, and oxygen. One example is the i-SANEX process that has been developed to separate minor actinides from spent nuclear fuel. The underlying science of such processes, is, however, both complex and intriguing. Recent research indicates that the liquid phases involved are frequently structured fluids with a hierarchical organization of aggregates. Effective flow-sheet modeling of such processes is likely to benefit from the knowledge of the fundamental properties of these phases. As a stepping stone toward this, we have performed molecular dynamics simulations on a metal free i-SANEX system composed of the ligand N,N,N′,N′-tetraoctyl diglycolamide (TODGA), diluent hydrogenated tetrapropylene (TPH), and polar species water and nitric acid. We have also studied the effects of adding n-octanol and swapping TPH for n-dodecane. It would seem sensible to understand this simpler system before introducing metal complexes. Such an understanding would ideally arise from studying the system’s properties over a wide range of compositions. The large number of components, however, precludes a comprehensive scan of compositions, so we have chosen to study a fixed concentration of TODGA while varying the concentrations of water and nitric acid over a substantial range. Reverse aggregates are observed, with polar species in the interior in contact with the polar portions of the TODGA molecules and the organic diluent on the exterior in contact with the TODGA alkyl chains. These aggregates are irregular in shape and grow in size as the amount of water and nitric acid increases. At a sufficiently high polar content, a single extended cluster forms corresponding to the third phase formation. No well-defined bonding motifs were observed between the polar species and TODGA. The cluster size distribution fits an isodesmic model, where the Gibbs energy change of adding a TODGA molecule to a cluster ranges between 4.5 and 7.0 kJ mol^–1, depending on the system composition. The addition of n-octanol was found to reduce the degree of aggregation, with n-octanol acting as a co-surfactant. Exchanging the diluent TPH for n-dodecane also decreased the aggregation. We present evidence that this is due to the greater penetration of n-dodecane into the reverse aggregates. It is known, however, that the propensity for the third phase formation is greater with n-dodecane as the diluent than is the case with TPH, but we argue that these two results are not contradictory. This research casts light on the driving forces for aggregation, informs process engineers as to what species are present, and indicates that flow-sheet liquid–liquid extraction modeling might benefit by incorporating an isodesmic aggregation approach.

Liquid/liquid interface in periodic boundary condition

We study how surface phenomena can change the interface geometry in liquid-liquid two-phase systems with periodic boundary conditions. Without any curvature effect on surface tension, planar (slab), cylindrical, and spherical structures are successively obtained as a function of the total composition and elongation of the box, in accordance with molecular dynamics simulations for a water/heptane system. The curvature effects described by Tolman relationship desymmetrize the phase diagram by stabilizing a concavity but it leads to inconsistencies with high curvature. Helfrich model partially resolves this and predicts the possible presence of shells reflecting a frustrated system.

A microfluidic study of synergic liquid–liquid extraction of rare earth elements

A membrane based liquid–liquid extraction microfluidic device coupled with X-ray fluorescence enables the first quantification of free energies of transfer dependence for a complex mixture of rare earth elements and iron using synergic extractants.

Colloidal Model for the Prediction of the Extraction of Rare Earths Assisted by the Acidic Extractant

We propose the statistical thermodynamic model for the prediction of the liquid-liquid extraction efficiency in the case of rare-earth metal cations using the common bis(2-ethyl-hexyl)phosphoric acid (HDEHP) extractant. In this soft matter-based approach, the solutes are modeled as colloids. The leading terms in free-energy representation account for: the complexation, the formation of a highly curved extractant film, lateral interactions between the different extractant head groups in the film, configurational entropy of ions and water molecules, the dimerization, and the acidity of the HDEHP extractant. We provided a full framework for the multicomponent study of extraction systems. By taking into account these different contributions, we are able to establish the relation between the extraction and general complexation at any pH in the system. This further allowed us to rationalize the well-defined optimum in the extraction engineering design. Calculations show that there are multiple extraction regimes even in the case of lanthanide/acid system only. Each of these regimes is controlled by the formation of different species in the solvent phase, ranging from multiple metal cation-filled aggregates (at the low acid concentrations in the aqueous phase), to the pure acid-filled aggregates (at the high acid concentrations in the aqueous phase). These results are contrary to a long-standing opinion that liquid-liquid extraction can be modeled with only a few species. Therefore, a traditional multiple equilibria approach is abandoned in favor of polydisperse spherical aggregate formations, which are in dynamic equilibrium.

Multicomponent Model for the Prediction of Nuclear Waste/Rare-Earth Extraction Processes

We develop a minimal model for the prediction of solvent extraction. We consider a rare earth extraction system for which the solvent phase is similar to water-poor microemulsions. All physical molecular quantities used in the calculation can be measured separately. The model takes into account competition complexation, mixing entropy of complexed species, differences of salt concentrations between the two phases, and the surfactant nature of extractant molecules. We consider the practical case where rare earths are extracted from iron nitrates in the presence of acids with a common neutral complexing extractant. The solvent wetting of the reverse aggregates is taken into account via the spontaneous packing. All the water-in-oil reverse aggregates are supposed to be spherical on average. The minimal model captures several features observed in practice: reverse aggregates with different water and extractant content coexist dynamically with monomeric extractant molecules at and above a critical aggregate concentration (CAC). The CAC decreases upon the addition of electrolytes in the aqueous phase. The free energy of transfer of an ion to the organic phase is lower than the driving complexation. The commonly observed log-log relation used to determine the apparent stoichiometry of complexation is valid as a guideline but should be used with care. The results point to the fact that stoichiometry, as well as the probabilities of a particular aggregate, is dependent on the composition of the entire system, namely the extractant and the target solutes' concentrations. Moreover, the experimentally observed dependence of the extraction efficiency on branching of the extractant chains in a given solvent can be quantified. The evolution of the distribution coefficient of particular rare earth, acid, or other different metallic cations can be studied as a function of initial extractant concentration through the whole region that is typically used by chemical engineers. For every chemical species involved in the calculation, the model is able to predict the exact equilibrium concentration in both the aqueous and the solvent phases at a given thermodynamic temperature.

Thermodynamic Description of Synergy in Solvent Extraction: II Thermodynamic Balance of Driving Forces Implied in Synergistic Extraction

In the second part of this study, we analyze the free energy of transfer in the case of synergistic solvent extraction. This free energy of the transfer of an ion in dynamic equilibrium between two coexisting phases is decomposed into four driving forces combining long-range interactions with the classical complexation free energy associated with the nearest neighbors. We demonstrate how the organometallic complexation is counterbalanced by the cost in free energy related to structural change on the colloidal scale in the solvent phase. These molecular forces of synergistic extraction are driven not only by the entropic term associated with the tight packing of electrolytes in the solvent and by the free energy cost of coextracting water toward the hydrophilic core of the reverse aggregates present but also by the entropic costs in the formation of the reverse aggregate and by the interfacial bending energy of the extractant molecules packed around the extracted species. Considering the sum of the terms, we can rationalize the synergy observed, which cannot be explained by classical extraction modeling. We show an industrial synergistic mixture combining an amide and a phosphate complexing site, where the most efficient/selective mixture is observed for a minimal bending energy and maximal complexation energy.

Synergy in Extraction System Chemistry: Combining Configurational Entropy, Film Bending, and Perturbation of Complexation

Iron-uranium selectivity in liquid-liquid extraction depends not only on the mole fraction of extractants, but also on the nature of the diluent used, even if the diluent has no complexation interaction with the extracted ions. Modeling strong nonlinearity is difficult to parametrize without a large number of parameters, interpreted as "apparent constants". We determine in this paper the synergy curve versus mole fraction of HDEHP-TOPO (di(2-ethylexyl) phosphoric acid/tri-n-octyl phosphine oxide) and compare the free energy of aggregation to the free energy of extraction in various diluents. There is always a concomitant maximum of the two quantities, but with a gradual influence on intensity. The diluent is wetting the chains of the reverse aggregates responsible of the extraction. We show here that the intensity of the unexplained synergy peak is strongly dependent on the "penetrating" or "nonpenetrating" nature of the diluent. This experimental determination allows us to attribute the synergy to a combination of entropic effects favoring extraction, opposed to perturbation of the first coordination sphere by penetration as well as surfactant film bending energy.

Complexation-induced supramolecular assembly drives metal-ion extraction

Combining experiment with theory reveals the role of self-assembly and complexation in metal-ion transfer through the water-oil interface. The coordinating metal salt Eu(NO3)3 was extracted from water into oil by a lipophilic neutral amphiphile. Molecular dynamics simulations were coupled to experimental spectroscopic and X-ray scattering techniques to investigate how local coordination interactions between the metal ion and ligands in the organic phase combine with long-range interactions to produce spontaneous changes in the solvent microstructure. Extraction of the Eu(3+)-3(NO3(-)) ion pairs involves incorporation of the "hard" metal complex into the core of "soft" aggregates. This seeds the formation of reverse micelles that draw the water and "free" amphiphile into nanoscale hydrophilic domains. The reverse micelles interact through attractive van der Waals interactions and coalesce into rod-shaped polynuclear Eu(III) -containing aggregates with metal centers bridged by nitrate. These preorganized hydrophilic domains, containing high densities of O-donor ligands and anions, provide improved Eu(III) solvation environments that help drive interfacial transfer, as is reflected by the increasing Eu(III) partitioning ratios (oil/aqueous) despite the organic phase approaching saturation. For the first time, this multiscale approach links metal-ion coordination with nanoscale structure to reveal the free-energy balance that drives the phase transfer of neutral metal salts.