A Quasi-relativistic Density Functional Theory Study of the Actinyl(VI, V) (An = U, Np, Pu) Complexes with a Six-Membered Macrocycle Containing Pyrrole, Pyridine, and Furan Subunits

Computational study of core modified dipyriamethyrin for the competitive complexation of Am3+/Cm3+ from their trichlorides

In the process of handling and storage of radioactive actinides it is essential to selectively sequester the minor actinides, such as Am and Cm, through a competitive complexation process. Herein we computationally designed two core modified ligands (L21- and L3) through systematic oxygen substitution at the NH sites of dipyriamethyrin (L1_2H), a hexadentate expanded porphyrin, and studied their competitive complexation towards trivalent actinides (An = Am/Cm) from their trichlorides using density functional theory (DFT). We observed shorter An-N bonds and longer An-O bonds in complexes based on core modified ligands (L21- and L3). The An-Cl bond length increases with increasing axial coordination number (i.e., from L12- to L3) to accommodate the ligands. All the bonds were identified to be electrostatic in nature. L12- exhibits shorter bonds and larger bond orders on complexing with Am than with Cm. On moving from complexes of L21- to L3, the An-N bond lengths are shortened, while An-O bond lengths become larger. Between the complexes of Am and Cm, there is marginal difference in their bond distances with L21- and L3. Charge analysis shows ligand to metal charge transfer during coordination, with back-donation from An to N/O and Cl. The calculated spin-density analysis indicates that An remains in its trivalent oxidation state on complexation, while orbital occupation analysis shows that the 5f and 6d orbitals are involved in bonding; this was confirmed by molecular orbital (MO) analysis that shows the complexes of L21- and L3 to exhibit higher degeneracy in their overlapping MOs. Further, the energy decomposition analysis (EDA) confirms that all ionic bonds are primarily due to electrostatic contributions, where the orbital contributions increase from L12- to L3 complexes and maximum covalency was observed in Cm complexes due to the energy matching between the 5f orbitals of Cm and the 2p orbitals of N and Cl, compared to Am. To confirm the competitiveness in the complexation of the ligand towards Am vs. Cm, the thermodynamic parameters were analysed for the ligand and metal substitution reactions. L12- shows more affinity towards Am than Cm, while L21- and L3 prefer Cm.

Pyridine Diphosphonate Ligand for Stabilization of Tetravalent Uranium and Neptunium in Aqueous Medium under Aerobic Conditions

Though uranium is usually present in its +6 oxidation state (as uranyl ion) in aqueous solutions, its conversion to oxidation states such as +4 or +5 is a challenging task. Electrochemical reduction and axial oxo activation are the preferred methods to get stable unusual oxidation states of uranium in an aqueous medium. In previous studies, dicarboxylic acid has been used to stabilize UO2+ in aqueous alkaline solutions. In the present work, a diphosphonate ligand was chosen due to its higher complexing ability compared to that of the carboxylate ligands. Neptunium complexation studies with 2,6-pyridinediphosphonic acid (PyPOH) indicated the formation of different species at different pH values and the complexation facilitates disproportionation of NpO2+ to Np4+ and NpO22+ at pH 2. Hexavalent actinides form insoluble complexes in aqueous media at pH = 2, as confirmed by UO22+ complexation studies. The in situ complexation-driven precipitation resulted in conversion to pure Np4+ in aqueous media as the Np4+-PyPOH complex. A strong complexing ability of the PyPOH ligand toward the Np4+ ion is also seen for the stabilization of the electrochemically generated U4+ in aqueous medium under aerobic conditions. The U4+-PyPOH complex was found to be stable for 3 months. Raman, UV-vis, fluorescence, and cyclic voltametric studies along with density functional theory (DFT) calculations were done to get structural insights into the PyPOH complexes of actinides in different oxidation states.

Periodic Trends in the Stabilization of Actinyls in Their Higher Oxidation States Using Pyrrophen Ligands

Owing to the prominent existence and unique chemistry of actinyls, their complexation with suitable ligands is of significant interest. The complexation of high-valent actinyl moieties (An = U, Np, Pu and Am) with the acyclic sal-porphyrin analogue called "pyrrophen" (L₍₁₎) and its dimethyl derivative (L₍₂₎) with four nitrogen and two oxygen donor atoms was studied using relativistic density functional theory. Based on the periodic trends, the [U^VO₂-L₍₁₎/L₍₂₎]^1- complexes show shorter bond lengths and higher bond orders that increase across the series of pentavalent actinyl complexes mainly due to the localization of the 5f orbitals. Among the hexavalent complexes, the [U^VIO₂-L₍₁₎/L₍₂₎] complexes have the shortest bonds. Following the uranyl complex, due to the plutonium turn, the [Am^VIO₂-L₍₁₎/L₍₂₎] complexes exhibit comparable properties with those of the former. Charge analysis suggests the complexation to be facilitated through ligand-to-metal charge transfer (LMCT) mainly through σ donation. Thermodynamic feasibility of complexation was modeled using hydrated actinyl moieties in aqueous medium and was found to be spontaneous. The dimethylated pyrrophen (L₍₂₎) shows higher magnitudes of thermodynamic parameters indicating increased feasibility compared to the unsubstituted ligand (L₍₁₎). Energy decomposition analysis (EDA) along with extended transition-state-natural orbitals for chemical valence theory (ETS-NOCV) analysis shows that the dominant electrostatic contributions decrease across the series and are counteracted by Pauli repulsion. Slight but considerable covalency is provided to hexavalent actinyl complexes by orbital contributions; this was confirmed by molecular orbital (MO) analysis that suggests strong covalency in americyl (VI) complexes. In addition to the pentavalent and hexavalent actinyl moieties, heptavalent actinyl species of neptunyl, plutonyl, and americyl were studied. Beyond the influence of the charges, the geometric and electronic properties point to the stabilization of neptunyl (VII) in the pyrrophen ligand environment, while the others shift to a lower (+VI) and relatively stable OS on complexation.

Unusually Kinetically Inert Monocationic Neptunyl Complex with a Fluorescein-Modified 1,10-Phenanthroline-2,9-dicarboxylate Ligand: Specific Separation and Detection in Gel Electrophoresis

We found a singly charged Np(V)O₂⁺ complex with unprecedented kinetic inertness in aqueous solution, one million times slower than the widely accepted fast kinetics of neptunyl complexes. An inert NpO₂⁺ complex with a fluorescent 1,10-phenanthroline-2,9-dicarboxylate derivative was found by kinetic selection using polyacrylamide gel electrophoresis (PAGE) from a small chemical library. Autoreduction from Np(VI)O₂²⁺ to Np(V)O₂⁺ via complexation was observed. A remarkably small spontaneous dissociation rate constant of 8 × 10^-6 s^-1 (half-life of 23 h) was determined using PAGE. Selective detection of Np(V)O₂⁺ was achieved in PAGE with a detection limit of 68 pmol dm^-3 (17 fg). This system was successfully applied to simulated radioactive waste samples. Our finding that electron-rich NpO₂⁺ forms a uniquely inert complex with no strong electrostatic interaction reveals a new aspect of actinide chemistry for developing a novel separation system of real radioactive material samples.

Contribution of Molecular Dynamics in pNMR for the Structural Determination of AnV and AnVI Complexes in Solution

In this study, we propose to use classical molecular dynamics (MD) coupled with ¹H NMR spectroscopy to study the conformations of different actinyl An^VI (An = U, Np, and Pu) and An^V (An = Np) complexes with tetra-ethyl dyglicolamide (TEDGA) ligands in order to have a better representation of such complexes in solution. Molecular dynamics simulations showed its effectiveness in interpreting the experiments by the calculation of geometric factors needed for the determination of magnetic properties of these complexes. We demonstrated that different conformations of the An^V and An^VI complexes with TEDGA exist in solution with different coordination modes, which is experimentally confirmed by ¹H NMR and EXAFS spectroscopies. Furthermore, MD simulations provide additional insights into the structures of complexes in solution since conformations with fast exchanges, which are not accessible from NMR experiments, have been observed by MD simulations.

Chemical bonding in actinyl(V/VI) dipyriamethyrin complexes for the actinide series from americium to californium: a computational investigation

The separation of minor actinides in their dioxocation (i.e., actinyl) form in high-valence oxidation states requires efficient ligands for their complexation. In this work, we evaluate the complexation properties of actinyls including americyl, curyl, berkelyl, and californyl in their pentavalent and hexavalent oxidation states with the dipyriamethyrin ligand (L) using density functional theory calculations. The calculated bond parameters show shorter AnO_yl bonds with covalent character and longer An-N bonds with ionic character. The bonding between the actinyl cation and the ligand anion shows a flow of charges from the ligand to actinyl in all [An^V/VIO₂-L]^1-/0 complexes. However, across the series, backdonation of charges from the metal to the ligand becomes prominent and stabilizes the complexes. The thermodynamic parameters in the gas phase and solution suggest that the complex formation reaction is spontaneous for [Cf^V/VIO₂-L]^1-/0 complexes and spontaneous at elevated temperatures (>298.15 K) for all other complexes. Spin-orbit corrections have a quantitative impact while the overall trend remains the same. Energy decomposition analysis (EDA) reveals that the interaction between actinyl and the ligand is mainly due to electrostatic contributions that decrease from Am to Cf along with an increase in orbital contributions due to the backdonation of charges from the actinyl metal center to the ligand that greatly stabilizes the Cf complex. The repulsive Pauli energy contribution is observed to increase in the case of [An^VO₂-L]^1- complexes from Am to Cf while a decrease is observed among [An^VIO₂-L]⁰ complexes, showing minimum repulsion in [Cf^VIO₂-L]⁰ complex formation. Overall, the hexavalent actinyl complexes show greater stability (increasing from Am to Cf) than their pentavalent counterparts.

Theoretical Insights into the Reduction Mechanism of Np(VI) with Phenylhydrazine

Effectively adjusting and controlling the valence state of neptunium from the spent fuel reprocessing process is essential to separating neptunium. Hydrazine and its derivatives as free-salt reductants have been experimentally demonstrated to effectively reduce Np(VI) to Np(V). We have theoretically investigated the reduction mechanisms of Np(VI) with hydrazine and three derivatives (HOC₂H₄N₂H₃, CH₃N₂H₃, and CHON₂H₃) in previous works. Herein, we further explored the reduction reaction of Np(VI) with phenylhydrazine (C₆H₅N₂H₃) including the free radical ion mechanism and the free radical mechanism. Potential energy profiles (PEPs) indicate that the rate-determining step of both mechanisms is the first stage. Moreover, for the free radical ion mechanism, phenylhydrazine possesses better reduction ability to Np(VI) compared to HOC₂H₄N₂H₃, CH₃N₂H₃, and CHON₂H₃, which falls completely in line with the experimental results. Additionally, the analyses of the quantum theory of atoms in molecules (QTAIM), natural bond orbitals (NBOs), electron localization function (ELF), and localized molecular orbitals (LMOs) have been put forward to elucidate the bonding evolution for the structures of the reaction pathways. This work offers insights into the reduction mechanism of Np(VI) with phenylhydrazine from the theory point of view and contributes to design more high-efficiency reductants for the separation of U/Np and Np/Pu in spent fuel reprocessing.

Computational Study of Actinyl Ion Complexation with Dipyriamethyrin Macrocyclic Ligands

Relativistic density functional theory has been employed to characterize [AnO₂(L)]^0/-1 complexes, where An = U, Np, Pu, and Am, and L is the recently reported hexa-aza porphyrin analogue, termed dipyriamethyrin, which contains six nitrogen donor atoms (four pyrrolic and two pyridine rings). Shorter axial (An═O) and longer equatorial (An-N) bond lengths are observed when going from An^VI to An^V. The actinide to pyrrole nitrogen bonds are shorter as compared to the bonds to the pyridine nitrogens; the former also play a dominant role in the formation of the actinyl (VI and V) complexes. Natural population analysis shows that the pyrrole nitrogen atoms in all the complexes carry higher negative charges than the pyridine nitrogens. Upon binding actinyl ions with the ligand a significant ligand-to-metal charge transfer takes place in all the actinyl (VI and V) complexes. The formation energy of the actinyl(VI,V) complexes in the gas-phase is found to decrease in the order of UO₂L > PuO₂L > NpO₂L > AmO₂L. This trend is consistent with results for the formation of complexes in dichloromethane solution. The calculated ΔG and ΔH values are negative for all the complexes. Energy decomposition analysis (EDA) indicates that the interactions between actinyl(V/VI) and ligand are mainly controlled by electrostatic components over covalent orbital interactions, and the covalent character gradually decreases from U to Am for both pentavalent and hexavalent actinyl complexes.

Theoretical Study on the Reduction Mechanism of Np(VI) by Hydrazine Derivatives

The key to effective separation of neptunium from the spent fuel reprocessing process is to adjust and control its valence state. Hydrazine and its derivatives have been experimentally confirmed to be effective salt-free reductants for reducing Np(VI) to Np(V). We theoretically studied the reduction reactions of Np(VI) with three hydrazine derivatives (2-hydroxyethyl hydrazine (HOC₂H₄N₂H₃), methyl hydrazine (CH₃N₂H₃), and formyl hydrazide (CHON₂H₃)) and obtained the free radical ion mechanism and the free radical mechanism. Their potential energy profiles (PEPs) suggest that the free radical mechanism is the most probable reaction. Based on the energy barrier of the free radical ion mechanism, the trend of the reduction ability of the three hydrazine derivatives is HOC₂H₄N₂H₃ > CH₃N₂H₃ > CHON₂H₃, which is in excellent agreement with the experimental results. Lastly, the analyses of natural bond orbitals (NBOs), quantum theory of atoms-in-molecules (QTAIM), and electron localization function (ELF) have been carried out to explore the bonding evolution of the structures along the reaction pathways. This work provides an insight into the reduction mechanism of Np(VI) with hydrazine derivatives from the theoretical perspective and helps to design more effective reductants for the separation of U/Np and Np/Pu in spent fuel reprocessing.

Separation of actinides from lanthanides associated with spent nuclear fuel reprocessing in China: current status and future perspectives

Developing necessary reprocessing techniques to meet the remarkable increase of spent nuclear fuels (SNFs) is crucial for the sustainable development of nuclear energy. This review summarizes recent research progresses related to the SNF reprocessing in China, with an emphasis on actinides separation over lanthanides through three different techniques, hydrometallurgical reprocessing, pyrometallurgical processes, and selective crystallization based separation. Some future perspectives with respect to advanced actinide separation are also given.

Uranyl-Bound Tetra-Dentate Non-Innocent Ligands: Prediction of Structure and Redox Behaviour Using Density Functional Theory

Computational methods have been applied to understand the reduction potentials of [UO₂ -salmnt-L] complexes (L=pyridine, DMSO, DMF and TPPO), and their redox behavior is compared with previous experiments in dichloromethane solution. Since the experimental results were inconclusive regarding the influence of the uranyl-bound tetra-dentate 'salmnt' ligand, here we will show that salmnt acts as a redox-active ligand and exhibits non-innocent behavior to interfere with the otherwise expected one-electron metal (U) reduction. We have employed two approaches to determine the uranyl (VI/V) reduction potentials, using a direct study of one-electron reduction processes and an estimation of the overall reduction using isodesmic reactions. Hybrid density functional theory (DFT) methods were combined with the Conductor-like Polarizable Continuum Model (CPCM) to account for solvation effects. The computationally predicted one-electron reduction potentials for the range of [UO₂ -salmnt-L] complexes are in excellent agreement with shoulder peaks (∼1.4 eV) observed in the cyclic voltammetry experiments and clearly correlate with ligand reduction. Highly conjugated pi-bonds stabilize the ligand based delocalized orbital relative to the localized U f-orbitals, and as a consequence, the ligand traps the incoming electron. A second reduction step results in metal U(VI) to U(V) reduction, in good agreement with the experimentally assigned uranyl (VI/V) reduction potentials.

Novel diamide ligands with a central carbonyl group and their comparative evaluation with the diglycolamide ligand: synthesis, extraction, DFT and chromatographic studies

The extraction behaviour of U(VI), Th(IV) and Nd(III) was investigated as a function of nitric acid concentration for diamide based extractants, namely, N,N,N′,N′-tetraoctyl-3-carbonylpentanediamide (TOCPDA) and 4-carbonyl-heptanedioic acid bis-dioctylamide (CHADA). In addition, the distribution ratio was also measured for Pu(IV) and Sr(II) with 1.1 M CHADA in n-dodecane. These extractants were synthesized by adopting simple acid, amine coupling reaction with DCC (dicyclohexylcarbodiimide) and DMAP (N,N′-dimethylaminopyridine) as the coupling agent. The newly synthesized extractants were characterized by FT-IR, NMR, Mass, CHNS and HPLC. The extraction results indicated that CHADA shown has better extraction behavior for U(VI) compared to TOCPDA. In addition, CHADA coated HPLC column was examined for the retention behaviour of U(VI), Th(IV), and Nd(III). Computation studies based on density functional theory (DFT) were carried out to understand the complexing behaviour of U(VI), Pu(IV) and Sr(II) with CHADMA and TMCPDA.

Uncovering the impact of 'capsule' shaped amine-type ligands on Am(iii)/Eu(iii) separation

Separation of trivalent actinides (An(iii)) and lanthanides (Ln(iii)) in spent nuclear fuel reprocessing is extremely challenging mainly owing to their similar chemical properties. Two amine-type reagents, tetrakis(2-pyridyl-methyl)-1,2-ethylenediamine (TPEN) and its hydrophobic derivative N,N,N',N'-tetrakis((4-butoxypyridin-2-yl)methyl)-ethylenediamine (TBPEN), have been identified to possess a selectivity for Am(iii) over Eu(iii). In this work, the structures, bonding nature, and thermodynamic behaviors of the Am(iii) and Eu(iii) complexes with these two ligands have been systematically studied via scalar relativistic density functional theory (DFT) calculations. According to Mayer bond order and the quantum theory of atoms in molecules (QTAIM) analyses, interactions between the ligands and metal cations exhibit some degree of covalent character with relatively more covalency for Am(iii) complexes. In comparison with TPEN, TBPEN has better extractability but worse separation ability for Am(iii) and Eu(iii). Four nitrogen atoms in pyridine moieties may be responsible for the different extraction abilities of TPEN and TBPEN, while two nitrogen atoms in amine chains of these ligands appear to play more important roles in the separation of Am(iii)/Eu(iii). These different extraction behaviors may be attributed to the longer and thinner 'capsule' shaped TBPEN ligand compared to TPEN. Our study might provide new insights into understanding the selectivity of the amine-type ligands toward minor actinides, and pave the way for designing new TPEN derivatives for extraction and separation of An(iii)/Ln(iii).

Complexation of a macrocyclic ligand, 2,6-di (N-methyl)formamide-calix[4]pyridine, with Eu(III) and extraction of Eu(III) and Am(III)

Complexation of a new macrocyclic compound, 2,6-dimethylformamide-calix[4]pyridine (L1 ), with Eu(III) was studied by spectrophotometry. Stability constants of the Eu(III)/L1 complex in different solvents were determined. The results reveal that L1 forms moderately strong complexes with Eu(III) and other lanthanides in aprotic solvents and shows little binding ability with transition metals. Moreover, the binding strength of L1 weakens significantly in protic solvents. Using 2-bromodecanoic acid as the synergistic reagent, L1 extracts Am(III) and Eu(III) successfully with a separation factor of SFAm/Eu=1.3, and the distribution ratios of Am(III) and Eu(III) increases as the aqueous acidity is decreased. DFT computational studies were conducted to corroborate the solvent extraction data, and compare the coordination properties of Am(III)/Eu(III) complexes with L1 and a related, 2,6-diformamide-calix[4]pyridine (L2 ). The computational results suggest that L2 could form stable complexes [ML]3+ and ML(NO3)3 [where M represent Am(III) or Eu(III)] in aqueous phase, in sharp contrast to the case of L1 where such complexes in aqueous phase are not stable.

A combined DFT and molecular dynamics study of U(VI)/calcite interaction in aqueous solution

Here we present a combined DFT and molecular dynamics study of uranyl (U(VI)) interaction mechanisms with the calcite (104) surface in aqueous solution. The roles of three anion ligands (CO₃^2-, HCO₃^-, OH^-) and solvation effect in U(VI) interaction with calcite have been evaluated. According to our calculations, water adsorbed on the calcite (104) surface prefers to exist in molecular state rather than dissociative state. Energy analysis indicate that the positively charged uranyl species prefers to form surface complexes on the surface, while neutral uranyl species may bind with the surface via both surface complexing and ion exchange reactions of U(VI)→Ca(II). In contrast, the negatively charged uranyl species prefer to interact with the surface via ion exchange reactions of U(VI)→Ca(II), and the one with UO₂(CO₃)₂(H₂O)^2- as the reactant becomes the most favorable one in energy. We also found that uranyl adsorption increases the hydrophilicability of the (104) surface to different extents, where the UO₂(CO₃)₃Ca₂ species contributes to the largest degree of energy changes (-53kcal/mol). Our calculations proved that the (104) surface also has the ability to immobilize U(VI) via either surface complexing or ion exchange mechanisms under different pH values.

Ligand size dependence of U–N and U–O bond character in a series of uranyl hexaphyrin complexes: quantum chemical simulation and density based analysis

A quantum chemical and density based analysis of bonding in uranyl hexaphyrin complexes, looking for trends in stability and covalency.

Synthesis and characterization of a dipyriamethyrin–uranyl complex

A bench stable uranyl complex of a hexaazadipyriamethryin macrocyclic ligand has been prepared and characterized.