Insights into the Coordination Behavior of Methyl-Substituted Phosphinic Acids with Actinides

The electronic structure and complexation behavior of methyl-substituted phosphinic acids with U(VI) and Pu(IV) were explored by applying quantum chemical methods. In contrast to Ingold's classification, our results indicate that the methyl group is electron-withdrawing, reducing the phosphoryl group electron density in substituted phosphinic acids. The magnitude of the computed complexation energy values increases along with the series, PA → MPA → DMPA, and MP → MMP → MDMP, implying an increasing complexation tendency upon methyl group substitution for both U(VI) and Pu(IV) complexes. One of the nitrate groups in UO₂(NO₃)₂•2L complexes (L = PA, MPA, and DMPA) is in monodentate coordination mode due to the additional stability gained from O₂N-O···H hydrogen bonding interactions with the acidic H atoms of respective ligands. The calculation indicates marginally stronger metal-ligand interactions in Pu(IV) complexes compared to that in U(VI), which is supported by the computed complexation energies, M-O_P bond lengths, ν(P═O), the extent of metal-ligand charge transfer, and properties of M-O_P bond critical points. The energy landscape of substituted phosphinic acid ligands is further analyzed within the framework of the activation strain model to explain the energetic preference of certain conformers.

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.

Influence of uranyl complexation on the reaction kinetics of the dodecane radical cation with used nuclear fuel extraction ligands (TBP, DEHBA, and DEHiBA)

Specialized extractant ligands - such as tri-butyl phosphate (TBP), N,N-di-(2-ethylhexyl)butyramide (DEHBA), and N,N-di-2-ethylhexylisobutryamide (DEHiBA) - have been developed for the recovery of uranium from used nuclear fuel by reprocessing solvent extraction technologies. These ligands must function in the presence of an intense multi-component radiation field, and thus it is critical that their radiolytic behaviour be thoroughly evaluated. This is especially true for their metal complexes, where there is negligible information on the influence of complexation on radiolytic reactivity, despite the prevalence of metal complexes in used nuclear fuel reprocessing solvent systems. Here we present a kinetic investigation into the effect of uranyl (UO₂²⁺) complexation on the reaction kinetics of the dodecane radical cation (RH˙⁺) with TBP, DEHBA, and DEHiBA. Complexation had negligible effect on the reaction of RH˙⁺ with TBP, for which a second-order rate coefficient (k) of (1.3 ± 0.1) × 10¹⁰ M^-1 s^-1 was measured. For DEHBA and DEHiBA, UO₂²⁺ complexation afforded an increase in their respective rate coefficients: k(RH˙⁺ + [UO₂(NO₃)₂(DEHBA)₂]) = (2.5 ± 0.1) × 10¹⁰ M^-1 s^-1 and k(RH˙⁺ + [UO₂(NO₃)₂(DEHiBA)₂]) = (1.6 ± 0.1) × 10¹⁰ M^-1 s^-1. This enhancement with complexation is indicative of an alternative RH˙⁺ reaction pathway, which is more readily accessible for [UO₂(NO₃)₂(DEHBA)₂] as it exhibited a much larger kinetic enhancement than [UO₂(NO₃)₂(DEHiBA)₂], 2.6× vs. 1.4×, respectively. Complementary quantum mechanical calculations suggests that the difference in reaction kinetic enhancement between TBP and DEHBA/DEHiBA is attributed to a combination of reaction pathway (electron/hole transfer vs. proton transfer) energetics and electron density distribution, wherein attendant nitrate counter anions effectively 'shield' TBP from RH˙⁺ electron transfer processes.

Speciation analysis the complexation of uranyl nitrate with tri-n-butyl phosphate in supercritical CO2

The complexation of solid uranyl nitrate with tri-n-butyl phosphate (TBP) in supercritical CO₂ is quite different from that of a liquid–liquid extraction system because fewer water molecules are involved. Here, the complexation mechanism was investigated by molecular dynamics simulation, emphasising on speciation distribution analysis. In the anhydrous uranyl nitrate system, poly-core uranyl-TBP species [UO₂(NO₃)₂]₂·3TBP and [UO₂(NO₃)₂]₃·3TBP were formed in addition to the predominant [UO₂(NO₃)₂]·1TBP and [UO₂(NO₃)₂]·2TBP species. The poly-core species was mainly constructed via the linkage of U Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O⋯U contributed by pre-developed [UO₂(NO₃)₂]·1TBP species. However, in the hydrated uranyl nitrate system, TBP·[UO₂(NO₃)₂]·H₂O species form, preventing the formation of the poly-core species. The complexation developed differently depending on the TBP to the uranyl nitrate ratio, the solute densities and the participation of water. It suggested that the kinetically favoring species would gradually convert into the thermodynamically stable species [UO₂(NO₃)₂]·2TBP by ligand exchange.

Uranyl nitrate complex with TBP in supercritical CO₂ could form 1 : 1 and 1 : 2 species, and further generate 2 : 3 and 3 : 3 poly-core uranyl species.

Key Factors Determining Efficiency of Liquid-Liquid Extraction: Implications from Molecular Dynamics Simulations of Biphasic Behaviors of CyMe4-BTPhen and Its Am(III) Complexes

CyMe₄-BTPhen (2,9-bis(5,5,8,8-tetramethyl-5,6,7,8-tetrahydro-1,2,4-benzotriazin-3-yl)-1,10-phenanthroline, denoted as L) has been considered as a promising extractant in lanthanide(III)/actinide(III) separation. Vast endeavors in its application put forward a compelling need on the understanding of the underlying mechanism in the liquid-liquid extraction. To address the issue of its dynamics in biphasic systems, we carried out molecular dynamics (MD) simulations of L and its complexes with a heavy f-block metal ion, americium(III) (Am³⁺) in "oil"/water binary solvents. Two types of organic phases have been considered, differing in the presence of octanol in the bulk n-dodecane or not, and the distribution of the solutes and their interfacial behaviors have been investigated. Two of the key factors that determine the efficiency of a liquid-liquid extraction protocol were delineated and discussed, that is, the appropriate ligand to enhance the lipophilicity of AmL complexes and appropriate way to form ion pairs to minimize the attraction between the complexes and aqueous phase. The simulations showed that the charge states of both ligand and AmL complexes were strongly correlated with their phase behavior, and the migration of neutral species was driven by van der Waals interactions while that of charged species by electrostatic interactions, indicating stronger lipophilicity of the former than the latter. The presence of octanol facilitated the migration of the ligand from the interface to the organic phase via hydrogen bond between its polar head and the ligand or the AmL complexes and constituted a polar core in the organic phase. This work bridged the widely used liquid-liquid extraction technique in chemistry to a fundamental chemical concept, that is, minimization of hydrophilicity and maximization of lipophilicity to facilitate phase transfer from the aqueous phase to the organic phase, and is expected to improve the understanding of dynamics of ligands and their complexes with metal ions and to contribute to the development of efficient protocols for phase transfer of target species.

Interfacial structure in the liquid-liquid extraction of rare earth elements by phosphoric acid ligands: a molecular dynamics study

Solvent extraction (SX), wherein two immiscible liquids, one containing the extractant molecules and the other containing the solute to be extracted are brought in contact to effect the phase transfer of the solute, underpins metal extraction and recovery processes. The interfacial region is of utmost importance in the SX process, since besides thermodynamics, the physical and chemical heterogeneity at the interface governs the kinetics of the process. Yet, a fundamental understanding of this heterogeneity and its implications for the extraction mechanism are currently lacking. We use molecular dynamics (MD) simulations to study the liquid-liquid interface under conditions relevant to the SX of Rare Earth Elements (REEs) by a phosphoric acid ligand. Simulations revealed that the extractant molecules and varying amounts of acid and metal ions partitioned to the interface. The presence of these species had a significant effect on the interfacial thickness, hydrogen bond life times and orientations of the water molecules at the interface. Deprotonation of the ligands was essential for the adsorption of the metal ions at the interface, with these ions forming a number of different complexes at the interface involving one to three extractant molecules and four to eight water molecules. Although the interface itself was rough, no obvious 'finger-like' water protrusions penetrating the organic phase were seen in our simulations. While the results of our work help us gain fundamental insights into the sequence of events leading to the formation of a variety of interfacial complexes, they also emphasize the need to carry out a more detailed atomic level study to understand the full mechanism of extraction of REEs from the aqueous to organic phases by phosphoric acid ligands.

Strong influence of weak hydrogen bonding on actinide-phosphonate complexation: accurate predictions from DFT followed by experimental validation

Among the varied classes of weak hydrogen bond, the CHO type is one of immense interest as it governs the finer structures of biological and chemical molecules, hence determining their functionalities. In the present work, this weak hydrogen bond has been shown to strongly influence the complexation behaviour of uranyl nitrate [UO2(NO3)2] with diamyl-H-phosphonate (DAHP) and its branched isomer disecamyl-H-phosphonate (DsAHP). The structures of the bare ligands and complexes have been optimized by density functional theory (DFT) calculations. Surprisingly, despite having the same chemical composition the branched UO2(NO3)2·2DsAHP complex shows a remarkably higher stability (by ∼14 kcal mol-1) compared to the UO2(NO3)2·2DAHP complex. Careful inspection of the optimized structures reveals the existence of multiple CHO hydrogen-bonding interactions between the nitrate oxygens or U[double bond, length as m-dash]O oxygens and the α-hydrogens in the alkyl chains of the ligands. Comparatively stronger such bonds are found in the UO2(NO3)2·2DsAHP complex. The binding free energies associated with the complexes are computed and favoured superior binding energetics for the more stable UO2(NO3)2·2DsAHP complex. Calculations involving diisoamyl-H-phosphonate (DiAHP) and its complexes have also been performed. Theoretical predictions are experimentally tested by carrying out the extraction of U(vi) from nitric acid media using these ligands. DAHP, DsAHP and DiAHP are synthesised, characterised by NMR and evaluated for their physicochemical properties viz. viscosity, density and aqueous solubility. It was experimentally discovered that indeed DsAHP complexation with uranyl nitrate is more favoured. H-phosphonates are generically classified as acidic extractants owing to the formation of an enol tautomer at lower acidities, hence complexing the metal ion by proton exchange. Our experiments interestingly reveal a neutral ligand characteristic for DsAHP alone which is generically an acidic extractant. Furthermore, the enol tautomer of H-phosphonates that governs their extraction profiles at low acidities is also explored by DFT and the anomalous pH dependent complexation trend of DsAHP could be successfully explained. The extractions of Pu(iv) and Th(iv) have also been carried out in addition to U(vi). Solvent extraction behaviour of Am(iii) was also studied with all three ligands; the positive binding energies computed for the Am(iii) complexation corroborate with our experimental results on the poor extraction of Am(iii).

Microscopic Behaviors of Tri- n-Butyl Phosphate, n-Dodecane, and Their Mixtures at Air/Liquid and Liquid/Liquid Interfaces: An AMBER Polarizable Force Field Study

In solvent extraction processes for recovering metal ions from used nuclear fuel, as well as other industrial applications, a better understanding of the metal complex phase transfer phenomenon would greatly aid ligand design and process optimization. We have approached this challenge by utilizing the classical molecular dynamics simulations technique to gain visual appreciation of the vapor/liquid and liquid/liquid interface between tri- n-butyl phosphate (TBP) and n-dodecane with air and water. In this study, we successfully reparameterized polarizable force fields for TBP and n-dodecane that accurately reproduced several of their thermophysical properties such as density, heat of vaporization, and dipole moment. Our models were able to predict the surface and interfacial tension of different systems when compared to experimental results that were also performed by us. Through this study, we gained atomistic understanding of the behaviors of TBP and n-dodecane at the interface against air and water, useful in further computational studies of such systems. Finally, our studies indicate that the initial configuration of a simulation may have a large effect on the final result.

UO22+ structure in solvent extraction phases resolved at molecular and supramolecular scales: a combined molecular dynamics, EXAFS and SWAXS approach

We developed a polarizable force field for unraveling the UO₂²⁺ structure in both aqueous and solvent extraction phases.

Structure, Dynamics, and Adsorption of Charged Guest within the Nanocavity of Polymer-Functionalized Neutral Macrocyclic Host

Host-guest encapsulation has been widely applied for purification and seizing of the metal ions. Macrocyclic crown ethers are one of the most popular hosts in the field of host-guest chemistry, which on functionalization with polymers are employed as an effective adsorbent. In spite of their vast applications, the microscopic information about their sensing mechanism toward cations/molecules is very scarce. Therefore, the present study is focused on the molecular insights of ion-exchange mechanism within the cavity of crown ether-functionalized polymers using molecular dynamics (MD) simulations. This present study investigates the molecular-level events of chloromethylated polystyrene (CMPS) bearing dibenzo-18-crown-6 (DB18C6) in the aqueous and acidic environment, which has been found to be particularly successful in sensing of various alkali and alkali earth metal ions. A strategy has been envisaged to design a crown ether-based functionalized polymeric resin, which exhibits good match of properties with the in-house-synthesized resin. The MD studies well capture the experimentally observed Langmuir-type adsorption isotherms of Li⁺ ions on crown ether-grafted polymer resins. The presence of acid reduces the adsorption of Li⁺ ions due to the competition with H₃O⁺ ions. In addition, the results revealed that the "adsorption in crown cavity" follows a dual residence time function. To the best of our knowledge, this is the first report on the adsorption isotherm of functionalized crown ether using MD simulations. The structure and dynamics of binding sites were explored using radial distribution functions and diffusion coefficients. All of these effects have been studied for different Li⁺-ion concentrations, acid concentrations, and counterions as well as different lengths of polymer chains and degrees of polymerization. Overall, the present study provides insights into and quantitative information about adsorption on the CMPS-DB18C6 resin, which might be useful in myriads of host-guest-based adsorption experiments.

Molecular Facts on the Structure and Dynamics of Electrolyte Species in Cu-Cl Cycle for Hydrogen Generation: An Insight from Molecular Dynamic Simulations

The Cu complex, which is the key chemical species in well-known Cu-Cl hybrid thermochemical cycles and also in numerous metal hydrometallurgical and sedimentary deposit processes, displays a wide variety of structural and dynamical characteristics that are further complicated by the presence of multiple oxidation states of Cu ions with different coordination chemistries, therefore they are difficult to explore from experiments alone. In this article, an attempt has been made to understand the coordination behavior of the Cu complex using MD simulations. The study provides compelling evidence of the experimentally observed multiple stoichiometries of Cu ions, i.e., 1:6:0, 1:5:1, and 1:4:2 for Cu⁺:H₂O:Cl^- and 1:6:0 for Cu²⁺:H₂O:Cl^-. The presence of the anionic Cu complex, [Cu⁺Cl₂]^-·2H₂O, [Cu⁺Cl₂]^-·3H₂O, [Cu²⁺Cl₃]^-·H₂O, and [Cu²⁺Cl₃]^-·2H₂O, was captured in the presence of excess chloride ions. Furthermore, the probability distribution profiles have been estimated to determine the most possible complex in the considered systems. The results establish structural and dynamical reformation of the Cu complex with change in the salt concentration or variation in the solvent medium in which they are dissolved. Moreover, the structure and kinetics of the Cu ions in the Cu-Cl electrolyzer have been explored over a large range of the electric field by extending the simulated systems for varied strengths of the electric fields. It has been observed that with an increase in the strength of the electric field, the water molecules lose their coordination strength with central Cu ions, which, on the other hand, results in a significant change in the structure of the captured complex. The diffusion dynamics of the ions is altered while applying the electric field, which is furthermore modified while increasing the strength of electric field beyond a critical limit. In fact, the diffusion mechanism of the ions was seen to be transformed from Brownian-like to linear motion and then to hopping diffusion with the increasing strength of the electric field. To the best of our knowledge, this is the first time when the multiple oxidation states of the Cu ion are explored using MD simulations, and the coexisting pictures of the multiple coordinations and the solvent effects have been clearly revealed. Also to date, the present article is the first one to report the insights of the structure and the dynamics of the ions in the Cu-Cl electrolyzer over a wide range of the electric field. The present studies will be very helpful in understanding the mechanism involved in numerous metal hydrometallurgical and sedimentary deposit processes and to comprehend the analogies involved in the electrode reactions of the Cu-Cl cycle for hydrogen generation.

Molecular dynamics simulation for the test of calibrated OPLS-AA force field for binary liquid mixture of tri-iso-amyl phosphate and n-dodecane

Tri-isoamyl phosphate (TiAP) has been proposed to be an alternative for tri-butyl phosphate (TBP) in the Plutonium Uranium Extraction (PUREX) process. Recently, we have successfully calibrated and tested all-atom optimized potentials for liquid simulations using Mulliken partial charges for pure TiAP, TBP, and dodecane by performing molecular dynamics (MD) simulation. It is of immense importance to extend this potential for the various molecular properties of TiAP and TiAP/n-dodecane binary mixtures using MD simulation. Earlier, efforts were devoted to find out a suitable force field which can explain both structural and dynamical properties by empirical parameterization. Therefore, the present MD study reports the structural, dynamical, and thermodynamical properties with different mole fractions of TiAP-dodecane mixtures at the entire range of mole fraction of 0-1 employing our calibrated Mulliken embedded optimized potentials for liquid simulation (OPLS) force field. The calculated electric dipole moment of TiAP was seen to be almost unaffected by the TiAP concentration in the dodecane diluent. The calculated liquid densities of the TiAP-dodecane mixture are in good agreement with the experimental data. The mixture densities at different temperatures are also studied which was found to be reduced with temperature as expected. The plot of diffusivities for TiAP and dodecane against mole fraction in the binary mixture intersects at a composition in the range of 25%-30% of TiAP in dodecane, which is very much closer to the TBP/n-dodecane composition used in the PUREX process. The excess volume of mixing was found to be positive for the entire range of mole fraction and the excess enthalpy of mixing was shown to be endothermic for the TBP/n-dodecane mixture as well as TiAP/n-dodecane mixture as reported experimentally. The spatial pair correlation functions are evaluated between TiAP-TiAP and TiAP-dodecane molecules. Further, shear viscosity has been computed by performing the non-equilibrium molecular dynamics employing the periodic perturbation method. The calculated shear viscosity of the binary mixture is found to be in excellent agreement with the experimental values. The use of the newly calibrated OPLS force field embedding Mulliken charges is shown to be equally reliable in predicting the structural and dynamical properties for the mixture without incorporating any arbitrary scaling in the force field or Lennard-Jones parameters. Further, the present MD simulation results demonstrate that the Stokes-Einstein relation breaks down at the molecular level. The present methodology might be adopted to evaluate the liquid state properties of an aqueous-organic biphasic system, which is of great significance in the interfacial science and technology.

Mechanism of phase transfer of uranyl ions: a vibrational sum frequency generation spectroscopy study on solvent extraction in nuclear reprocessing

VSFG study on interfaces of TBP/uranyl aqueous solutions shows that uranyl ion does not form complexes with TBP at the interface, and proposes a new extraction mechanism.

TBP Assisted Uranyl Extraction in Water-Dodecane Biphasic System: Insights from Molecular Dynamics Simulation

In the PUREX (Plutonium Uranium Recovery by Extraction Process) process, the extraction of uranyl ion from dissolver solution to the organic phase is influenced by co extraction of the other species, such as water and nitric acid and it is assumed that the presence of water or acid droplets in the organic phase intensifies the coordination mechanism of TBP. The present study illustrates the uranyl extraction from the aqueous phase to the organic phase using molecular dynamics (MD) simulation. Here, we consider the biphasic systems to gain insights into the characteristics of the interface and humidity of the organic phase under different acidic and neutral conditions. MD being a force field method, can’t satisfactorily model the bond making and breaking process therefore a priori choice has been made concerning the different status of proton for the acidic phase. Further, the importance of charge species transferability during uranyl-TBP complexation have been investigated considering two different models of uranyl nitrate; united UO2(NO3)2 complex and separate UO2 2+ and NO3 – ions. From the results, it is recommended to use the ionic uranyl model with separate UO2 2+ and NO3 – to study the structural and dynamical properties of extracted uranyl ions in the organic phase. Also, it was noticed that extracted uranyl ions in the organic phase are not completely dehydrated but are surrounded by water molecules. In other words the results show co extraction of other species such as water and acid molecules to the organic phase. Most remarkably, the present study evident that the neutral HNO3 effectively represents the acidity effect for the receiving phase in terms of acid/water extraction and their aggregation to form water droplet, especially when ionic model of uranyl nitrate is considered.

Ionic hydration-induced evolution of decane–water interfacial tension

We show that ionic hydration is responsible for the non-monotonic variation of the interfacial tension with increasing ionic concentration.

Network analysis and percolation transition in hydrogen bonded clusters: nitric acid and water extracted by tributyl phosphate

Network analysis of hydrogen bonded clusters formed in simulation by extraction of nitric acid and water by TBP interprets cluster topologies and identifies the mechanism for third phase formation.