Actinide-lanthanide separation (ALSEP) has been a topic of interest in recent years as it has been shown to selectively extract problematic metals from spent nuclear fuel. However, the process suffers from slow kinetics, prohibiting it from being applied to nuclear facilities. In an effort to improve the process, many fundamental studies have been performed, but the majority have only focused on the thermodynamics of separation. Therefore, to understand the mechanism behind the ALSEP process, molecular dynamics (MD) simulations were utilized to obtain the dynamics and solvation characteristics for an organic extractant, 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEHEHP). Simulations were conducted with both pure and biphasic solvent systems to evaluate the complex solvent interactions within the ALSEP extraction method. The MD simulations revealed solvation and dynamical behaviors that are consistent with the experimentally observed chemical properties of HEHEHP for the pure solvent systems (e.g., hydrophobic/hydrophilic behaviors of the polar head group and alkyl chains and dimer formation between the ligands within an organic solvent). When present in a biphasic solvent system, interfacial behaviors of the ligand revealed that, at low concentrations, the alkyl side chains of HEHEHP were parallel to the interfacial plane. Upon increasing the concentration to 0.75 M, tendency for the parallel orientation decreased and a more perpendicular-like orientation was observed. Analysis of ligand solvation energies in different solvents through the thermodynamic integration method demonstrated favorability toward n-dodecane and biphasic solvents, which is in agreement with the previous experimental findings.

The Actinide-Lanthanide Separation Process (ALSEP) is a solvent extraction approach for separating relevant trivalent minor actinides (e.g., americium and curium) from used nuclear fuel. However, relatively slow kinetics in the stripping step of the process restricts process throughput when scaled for industrial implementation. To assist in identifying specific kinetic barriers associated with the separation, the solvation and dynamic behaviors of the two organic extractants in the current ALSEP implementation, N,N,N′,N′-tetra(2ethylhexyl)diglycolamide (T2EHDGA) and 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEHEHP), were probed through molecular dynamics (MD) simulations. The simulations examined the effects of extractant and nitric acid concentration on the interfacial behavior of the extractants in three solvent systems (n-dodecane, water, and n-dodecane + water). Solvation analyses of T2EHDGA revealed expected amphiphilic behavior in pure solvent systems. In a nitric-acid-free biphasic solvent, it was found that T2EHDGA expressed similar interfacial conformations as HEHEHP, suggesting that a parallel-like configuration, relative to the interface, is adopted at low concentrations. When HNO 3 was introduced to biphasic systems containing a single molecule of extractant, HEHEHP was observed to retain a relatively parallel orientation while the T2EHDGA orientation was no longer affected by the presence of the interface. At bulk extractant concentrations, representative of the ALSEP process, the presence of nitric acid had minimal impact on the ligand orientation. Calculated diffusion constants showed that only some systems involving T2EHDGA were affected by the presence of acid.

Table S1. Calculated relativistic charges and coordinates for T2EHDGA. atom Coordinates charges (a.u.) X Y Z C25 3.54