Complexation of Lanthanides and Heavy Actinides with Aqueous Sulfur-Donating Ligands

Insights into coordination and ligand trends of lanthanide complexes from the Cambridge Structural Database

Understanding lanthanide coordination chemistry can help develop new ligands for more efficient separation of lanthanides for critical materials needs. The Cambridge Structural Database (CSD) contains tens of thousands of single crystal structures of lanthanide complexes that can serve as a training ground for both fundamental chemical insights and future machine learning and generative artificial intelligence models. This work aims to understand the currently available structures of lanthanide complexes in CSD by analyzing the coordination shell, donor types, and ligand types, from the perspective of rare-earth element (REE) separations. We obtain four sets of lanthanide complexes from CSD: Subset 1, all Ln-containing complexes (49472 structures); Subset 2, mononuclear Ln complexes (27858 structures); Subset 3, mononuclear Ln complexes without cyclopentadienyl ligands (Cp) (26156 structures); Subset 4, Ln complexes with at least one 1,10-phenanthroline (phen) or its derivative as a coordinating ligand (2226 structures). The subsequent analysis of lanthanide complexes in these subsets examines the trends in coordination numbers and first shell distances as well as identifies and characterizes the ligands and donor groups. In addition, examples of Ln-complexes with commercially available complexants and phen-based ligands are interrogated in detail. This systematic investigation lays the groundwork for future data-driven ligand designs for REE separations based on the structural insights into the lanthanide coordination chemistry.

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.

Toward Metal Extraction from Regolith: Theoretical Investigation of the Solvation Structure and Dynamics of Metal Ions in Ionic Liquids

Lunar and Martian regoliths, containing feldspar, pyroxene, ilmenite, olivine, and aluminite minerals, are excellent sources of metals such as aluminum, sodium, magnesium, and iron. Ionic liquids (ILs), which are excellent solvents with extremely low vapor pressure and high electrochemical stability, can be potentially leveraged for extracting metals from regolith in an extra-terrestrial environment. A critical step in the solvation process, which determines the effectiveness of the IL solvent, is the formation of solvation shells around the metal cations. To determine the rigidity and stability of the solvation shells, which has a direct implication on the extraction of metals, we performed classical molecular dynamics simulations of dilute solutions comprising individual metal ions Na+, Mg2+, and Al3+ in two distinct ILs, [mppy][TFSI] and [mppy][HSO4]. Our results indicate that the compactness of the structure is directly related to the charge density of the metal cation and the relative size and symmetry of the IL anion. Potentials of the mean force of the metal cation with the solvating IL anion indicate the presence of energy minima with barriers that increase with the surface charge density of the cation. The increasing energy barrier leads to greater residence time of metal cations in the solvation shell, which was confirmed by evaluating corresponding autocorrelation functions. Overall, our calculations provide fundamental insights into key factors that influence the solvation of metals and can be useful in the screening of ILs for digestion of metal-containing minerals in lunar and Martian regoliths.

Defining the optimal conditions using FFNNs and NARX neural networks for modelling the extraction of Sc from aqueous solution by Cryptand-2.2.1 and Cryptand-2.1.1

The main aim of this study is to figure out how well cryptand-2.2.1 (C 2.2.1) and cryptand-2.1.1 (C 2.1.1) macrocyclic compounds (MCs) work as novel extractants for scandium (Sc) by using an artificial neural network (ANN) models in MATLAB software. Moreover, C2.2.1 and C2.1.1 have never been evaluated to recover Sc. The independent variables impacting the extraction process (concentration of MC, concentration of Sc, pH, and time), and a nonlinear autoregressive network with exogenous input (NARX) and feed-forward neural network (FFNN) models were used to estimate their optimum values. The greatest obstacle in the selective recovery process of the REEs is the similarity in their physicochemical properties, specifically their ionic radius. The recovery of Sc from the aqueous solution was experimentally evaluated, then the non-linear relationship between those parameters was predictively modeled using (NARX) and (FFNN). To confirm the extraction and stripping efficiency, an atomic absorption spectrophotometer (AAS) was employed. The results of the extraction investigations show that, for the best conditions of 0.008 mol/L MC concentration, 10 min of contact time, pH 2 of the aqueous solution, and 75 mg/L Sc initial concentration, respectively, the C 2.1.1 and C 2.2.1 extractants may reach 99 % of Sc extraction efficiency. Sc was recovered from a multi-element solution of scandium (Sc), yttrium (Y), and lanthanum (La) under these circumstances. Whereas, at a concentration of 0.3 mol/L of hydrochloric acid, the extraction of Sc was 99 %, as opposed to Y 10 % and La 7 %. The Levenberg-Marquardt training algorithm had the best training performance with an mean-squared-error, MSE, of 5.232x10-6 and 6.1387x10-5 for C 2.2.1 and C 2.1.1 respectively. The optimized FFNN architecture of 4-10-1 was constructed for modeling recovery of Sc. The extraction process was well modeled by the FFNN with an R2 of 0.999 for the two MC, indicating that the observed Sc recovery efficiency consistent with the predicted one.

Lanthanide and Actinide Ion Complexes Containing Organic Ligands Investigated by Surface-Enhanced Infrared Absorption Spectroscopy

A new technique, surface-enhanced infrared absorption (SEIRA) spectroscopy, was used for the structural investigation of lanthanide (Ln) and actinide (An) complexes containing organic ligands. We synthesized thiol derivatives of organic ligands with coordination sites similar to those of 2-[N-methyl-N-hexanethiol-amino]-2-oxoethoxy-[N',N'-diethyl]-acetamide [diglycolamide (DGA)], Cyanex-272, and N,N,N',N'-tetrakis(2-pyridinylmethyl)-1,2-ethanediamine (TPEN), which have been used for separating Ln and An through solvent extraction. These ligands were attached on a gold surface deposited on an Si prism through S-Au covalent bonds; the gold surface enhanced the IR absorption intensity of the ligands. Aqueous solutions of Ln (Eu³⁺, Gd³⁺, and Tb³⁺) and An (Am³⁺) ions were loaded onto the gold surface to form ion complexes. The IR spectra of the ion complexes were obtained using Fourier transform infrared spectroscopy in the attenuated total reflection mode. In this study, we developed a new sample preparation method for SEIRA spectroscopy that enabled us to obtain the IR spectra of the complexes with a small amount of ion solution (5 μL). This is a significant advantage for the IR measurement of radiotoxic Am³⁺ complexes. In the IR spectra of DGA, the band attributed to C═O stretching vibrations at ∼1630 cm^-1 shifted to a lower wavenumber by ∼20 cm^-1 upon complexation with Ln and An ions. Moreover, the amount of the red shift was inversely proportional to the extraction equilibrium constant reported in previous studies on solvent extraction. The coordination ability of DGA toward Ln and An ions could be assessed using the band position of the C═O band. The Cyanex-272- and TPEN-like ligands synthesized in this report also showed noticeable SEIRA signals for Ln and An complexes. This study indicates that SEIRA spectroscopy can be used for the structural investigation of ion complexes and provides a microscopic understanding of selective extraction of Ln and An.

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm2+, Eu2+, Dy2+, Yb2+, Cf2+, and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6 (DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics (MD) simulations have been performed on four divalent lanthanide halide salts (Sm2+, Eu2+, Dy2+, and Yb2+) and one divalent actinide halide salt (Cf2+) bound to three diastereoisomers of DCH18C6. Dy2+, Yb2+, Cf2+, DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm2+ and Eu2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln2+/An2+-O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln2+/An2+-N solvated in acetonitrile for systems containing iodide suggested that the N atom of acetonitrile was competitive with I- for cation interactions. Fluctuation between five and six CNs for Dy2+ and Yb2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Finally, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

Advancing the Am Extractant Design through the Interplay among Planarity, Preorganization, and Substitution Effects

Advancing the field of chemical separations is important for nearly every area of science and technology. Some of the most challenging separations are associated with the americium ion Am(III) for its extraction in the nuclear fuel cycle, ²⁴¹Am production for industrial usage, and environmental cleanup efforts. Herein, we study a series of extractants, using first-principle calculations, to identify the electronic properties that preferentially influence Am(III) binding in separations. As the most used extractant family and because it affords a high degree of functionalization, the polypyridyl family of extractants is chosen to study the effects of the planarity of the structure, preorganization of coordinating atoms, and substitution of various functional groups. The actinyl ions are used as a structurally simplified surrogate model to quickly screen the most promising candidates that can separate these metal ions. The down-selected extractants are then tested for the Am(III)/Eu(III) system. Our results show that π interactions, especially those between the central terpyridine ring and Am(III), play a crucial role in separation. Adding an electron-donating group onto the terpyridine backbone increases the binding energies to Am(III) and stabilizes Am-terpyridine coordination. Increasing the planarity of the extractant increases the binding strength as well, although this effect is found to be rather weak. Preorganizing the coordinating atoms of an extractant to their binding configuration as in the bound metal complex speeds up the binding process and significantly improves the kinetics of the separation process. This conclusion is validated by the synthesized 1,2-dihydrodipyrido[4,3-b;5,6-b]acridine (13) extractant, a preorganized derivative of the terpyridine extractant, which we experimentally showed was four times more effective than terpyridine at separating Am³⁺ from Eu³⁺ (SF_Am/Eu ∼ 23 ± 1).

In situ beam reduction of Pu(IV) and Bk(IV) as a route to trivalent transuranic coordination complexes with hydroxypyridinone chelators

The solution-state interactions of plutonium and berkelium with the octadentate chelator 3,4,3-LI(1,2-HOPO) (343-HOPO) were investigated and characterized by X-ray absorption spectroscopy, which revealed in situ reductive decomposition of the tetravalent species of both actinide metals to yield Pu(III) and Bk(III) coordination complexes. X-ray absorption near-edge structure (XANES) measurements were the first indication of in situ synchrotron redox chemistry as the Pu threshold and white-line position energies for Pu-343-HOPO were in good agreement with known diagnostic Pu(III) species, whereas Bk-343-HOPO results were found to mirror the XANES behavior of Bk(III)-DTPA. Extended X-ray absorption fine structure results revealed An-O_HOPO bond distances of 2.498 (5) and 2.415 (2) Å for Pu and Bk, respectively, which match well with bond distances obtained for trivalent actinides and 343-HOPO via density functional theory calculations. Pu(III)- and Bk(III)-343-HOPO data also provide initial insight into actinide periodicity as they can be compared with previous results with Am(III)-, Cm(III)-, Cf(III)-, and Es(III)-343-HOPO, which indicate there is likely an increase in 5f covalency and heterogeneity across the actinide series.

Uranium(IV) Thio- and Selenoether Complexes: Syntheses, Structures, and Computational Investigation of U-ER2 Interactions

Rigid thioether- and selenoether-containing pincer proligands H[AS₂ Ph 2 ] (1) and H[ASe₂ Ph 2 ] (2) were synthesized, and deprotonation provided the potassium salts [K(AS₂ Ph 2 )(dme)] (3) and [K(ASe₂ Ph 2 )(dme)₂ ] (4). Reaction of two equivalents of 3 or 4 with [UI₄ (dioxane)₂ ] afforded the uranium thioether complex [(AS₂ Ph 2 )₂ UI₂ ] (5) and the first example of a uranium-selenoether complex, [(ASe₂ Ph 2 )₂ UI₂ ] (6). X-ray structures revealed distorted square antiprismatic geometries in which the AE₂ Ph 2 ligands are κ³ -coordinated. The nature of the U-ER₂ bonding in 5 and 6, as well as methyl-free analogues of 5 and 6 and a hypothetical ether analogue, was investigated computationally (including NBO, AIM, and ELF calculations) illustrating increasing covalency from O to S to Se.