Complexation of Trivalent Lanthanides with Hexaalkyl Nitrilotriacetamides into Methylimidazolium-Based Ionic Liquids: Spectroscopic, Electrochemical, Calorimetric, and Theoretical Insights

Exothermic, spontaneous inner-sphere complexation of trivalent lanthanides with N,N,N',N',N″,N″-hexaalkyl-substituted nitrilotriacetamides (HRNTAs) into 1-hexyl-3-methylimidazolium-bis(trifluoromethylsulfonyl)imide (C6mim NTf2) is reported with a predominant presence of ML and ML2 complexes having complexation constants, β1Eu 3.10 ± 0.02 ,β1Nd 2.67 ± 0.03 and β2Eu 5.22 ± 0.04, β2Nd 4.95 ± 0.02, respectively, for the n-butyl-substituted HRNTA (HBNTA); while those for the n-hexyl derivative (HHNTA) are β1Eu 4.27 ± 0.03, β1Nd 3.71 ± 0.03 and β2Eu 7.70 ± 0.03, β2Nd 7.18 ± 0.04, respectively. HHNTA shows better complexing ability; whereas the stronger complexation of the Nd3+ over Eu3+ is attributed to the lanthanide contraction. Furthermore, the nature of the ionic liquid also influences the extent of complexation with the trend: C4mim NTf2 > C6mim NTf2 > C8mim NTf2, which follows the order of their dielectric constants. Judd-Ofelt parameters were calculated from photoluminescence data to get an idea about the symmetry of the Eu3+ complexes. Electrochemical investigations give diffusion coefficient values of 1.17 × 10-7 and 8.26 × 10-8 cm2/s for the Eu3+ complexes of HBNTA and HHNTA, respectively. Changes in the spectral characteristics and peak positions are evidenced in the FTIR spectra on the complexation of Eu3+ with the HRNTA ligands in C6mim NTf2. Structure optimization for the complexes was performed by DFT computations.

Quantum tunnelling effect in the cis-trans isomerization of uranyl tetrahydroxide

The role of quantum tunnelling (QT) in the proton transfer kinetics of the uranyl tetrahydroxide (UTH, [UO2(OH)4]2-) cis to trans isomerization was computationally studied under three possible reaction pathways. The first pathway involved a direct proton transfer from the hydroxide ligand to the oxo atom. In the other two pathways, one or two water molecules were added to the second sphere. The first H2O, bound by hydrogen bonds to the ligands, acts as a bridge enabling a proton shuttling, a concerted hopping of a proton from the hydroxide to the oxo atom similar to the Grotthuss mechanism. In the third pathway, the second water molecule does not participate in the H-transfer chain, but works as an anchor for the first water molecule, limiting its movement and therefore enhancing the QT. Since experimentally the reaction occurs in water, the first two pathways (no water or one H2O) serve only as models of the gas phase behaviour, while the third pathway will always be thermodynamically and kinetically preferred. The effects were investigated in the gas phase as well as in a continuum aqueous model, including the H/D Kinetic Isotope Effect (KIE). The results indicate that at very low temperatures, QT is the only mechanism that permits the reaction kinetics, consistent with the large computed KIE. At higher temperatures, thermally activated tunnelling competes with the classical crossing over the potential barrier.

Understanding the Interaction of Uranyl Cation with Two C-Pivot Tripodal Amides: Synthesis, Complexation, Microcalorimetry, and DFT Studies

Complexation of uranyl ions with two structurally related C-pivotal tripodal amides with varying spacer lengths, synthesized for the first time, was studied by optical spectroscopy. In the tripodal amides, the coordination was through the carbonyl O atoms where the carbonyl groups were away from the central C-atom by three spacer atoms (LI) and four spacer atoms (LII), respectively. Increasing the spacer atoms going from LI to LII favors the complexation with the linear uranyl cations and results in stronger complex formation. The complexation heat between the uranyl cations and the two amide ligands was directly measured by microcalorimetric titrations. The complexation with both the ligands was driven by exothermic enthalpy and positive entropy changes. Formation of the complex proceeded by the replacement of water molecules from the primary coordination sphere of the uranyl cation. Both ligands formed bisolvated (ML2-type) complexes in which one unit of the ligand binds in a monodentate manner and the other in a bidentate mode. Density functional theory calculations further supported our experimental observations.

Synthesis, Structure, and Theoretical Calculations on NpO2Br42

The actinide-halogen complexes (AnO2X42-, X = Cl, Br, and I) are the simplest and most representative compounds for studying the bonding nature of actinides with ligands. In this work, we attempted to synthesize the crystals of NpO2X42- (X = Cl, Br, and I). The crystals of NpO2Cl42- and NpO2Br42- were successfully synthesized, in which the structure of NpO2Br42- was obtained for the first time. The crystal of NpO2I42- could not be obtained due to the rapid reduction of Np(VI) to Np(V) by I-. The molecular structures of NpO2Cl42- and NpO2Br42- were characterized by single-crystal X-ray diffraction and infrared, Raman, and UV-Vis-NIR absorption spectroscopy. The complexes of NpO2X42- (X = Cl, Br, and I) were also investigated by density functional theory calculations, and the calculated vibration frequencies and absorption features were comparable to the experimental results. Both the experimental results and theoretical calculations demonstrate the strengthened Np-O bonds and the weakened Np-X bonds across the NpO2X42- series; however, the population analysis on the frontier molecular orbitals (MOs) of NpO2X42- indicates a slight reduction in the Np-O bonding covalency and an enhancement in the Np-X bonding covalency from NpO2Cl42- to NpO2I42-. Results in this work have enriched the crystal database of the AnO2X42- family and provided insights into the bonding nature in the actinide complexes with soft- and hard-donor ligands.

Nature and coordination geometry of geologically relevant aqueous Uranium(VI) complexes up to 400 ºC: A review and new data

The structure of the uranyl aqua ion (UO₂²⁺) and a number of its inorganic complexes (specifically, UO₂Cl⁺, UO₂Cl₂⁰, UO₂SO₄⁰, [Formula: see text] , [Formula: see text] and UO₂OH₄^2-) have been characterised using X-Ray absorption spectroscopy/extended X-Ray absorption fine structure (XAS/EXAFS) at temperatures ranging from 25 to 326 ºC. Results of ab initio molecular dynamics (MD) calculations are also reported for uranyl in chloride and sulfate-bearing fluids from 25 to 400 ºC and 600 bar to 20 kilobar (kb). These results are reported alongside a comprehensive review of prior structural characterisation work with particular focus given to EXAFS works to provide a consistent and up-to-date view of the structure of these complexes under conditions relevant to U mobility in ore-forming systems and around high-grade nuclear waste repositories. Regarding reported EXAFS results, average equatorial coordination was found to decrease in uranyl and its sulfate and chloride complexes as temperature rose - the extent of this decrease differed between species and solution compositions but typically resulted in an equatorial coordination number of ∼3-4 at temperatures above 200 ºC. The [Formula: see text] complex was observed at temperatures from 25 to 247 ºC and exhibited no major structural changes over this temperature range. UO₂OH₄^2- exhibited only minor structural changes over a temperature range from 88 to 326 ºC and was suggested to manifest fivefold coordination with four hydroxyl molecules and one water molecule around its equator. Average coordination values derived from fits of the reported EXAFS data were compared to average coordination values calculated using the experimentally derived thermodynamic data for chloride complexes reported by Dargent et al. (2013) and Migdisov et al. (2018b), and for sulfate complexes reported by Alcorn et al. (2019) and Kalintsev et al. (2019). Sulfate EXAFS data were well described by available thermodynamic data, and chloride EXAFS data were described well by the thermodynamic data of Migdisov et al. (2018b), but not by the data of Dargent et al. (2013). The ab initio molecular dynamics calculations confirmed the trends in equatorial coordination observed with EXAFS and were also able to provide an insight into the effect of pressure in equatorial water coordination - for a given temperature, higher pressures appear to lead to a greater number of equatorially bound waters counteracting the temperature effect.

Redox Potentials of Uranyl Ions in Macrocyclic Complexes: Quantifying the Role of Counter-Ions

Several uranyl ions strapped with Schiff-base ligands in the presence of redox-innocent metal ions are synthesized, and their reduction potentials are recently estimated. The change in Lewis acidity of the redox-innocent metal ions contributes to ∼60 mV/pK_a unit quantified which is intriguing. Upon increasing the Lewis acidity of metal ions, the number of triflate molecules found near the metal ions also increases whose contributions toward the redox potentials remain poorly understood and not quantified until now. Most importantly, to ease the computational burden, triflate anions are often neglected in quantum chemical models due to their larger size and weak coordination to metal ions. Herein, we have quantified and dissected the individual contributions that arise alone from Lewis acid metal ions and from triflate anions with electronic structure calculations. The triflate anion contributions are large, in particular, for divalent and trivalent anions that cannot be neglected. It was presumed to be innocent, but we here show that they can contribute more than 50% to the predicted redox potentials, suggesting that their vital role in the overall reduction processes cannot be neglected.

Electrocatalytic and green system coupling strategy for simultaneous recovery and purification of uranium from uranium-containing wastewater

The recycling of uranium in wastewater is not only beneficial to the protection of ecological safety but also has great significance for the sustainable development of nuclear energy. However, there is no satisfactory method to recover and reuse uranium efficiently up to now. Here, we have developed an efficient and economical strategy that can achieve uranium recovery and direct reuse in wastewater. The feasibility analysis verified that the strategy still had good separation and recovery ability in acidic, alkaline, and high-salinity environments. The purity of uranium recovered from the separated liquid phase after electrochemical purification was up to about 99.95%. Ultrasonication could greatly increase the efficiency of this strategy, and 99.00% of high-purity uranium could be recovered within 2 h. We further improved the overall recovery rate by recovering the residual solid-phase uranium, and the overall recovery of uranium was increased to 99.40%. Moreover, the concentration of impurity ions in the recovered solution met the World Health Organization guidelines. In summary, the development of this strategy is of great importance for the sustainable use of uranium resources and environmental protection.

Influencing Bonding Interactions of the Neptunyl (V, VI) Cations with Electron-Donating and -Withdrawing Groups

Neptunium makes up the largest percentage of minor actinides found in spent nuclear fuel, yet separations of this element have proven difficult due to its rich redox chemistry. Developing new reprocessing techniques should rely on understanding how to control the Np oxidation state and its interactions with different ligands. Designing new ligands for separations requires understanding how to properly tune a system toward a desired trait through functionalization. Emerging technologies for minor actinide separations focus on ligands containing carboxylate or pyridine functional groups, which are desirable due to their high degree of functionalization. Here, we use DFT calculations to study the interactions of carboxylate and polypyridine ligands with the neptunyl cation [Np(V/VI)O2]+/2+. A systematic study is performed by varying the electronic properties of the carboxylate and polypyridine ligands through the inclusion of different electron-withdrawing and electron-donating R groups. We focus on how these groups can affect geometric properties, electronic structure, and bonding characterization as a function of the metal oxidation state and ligand character and discuss how these factors can play a role in neptunium ligand design principles.

Aggregation Behavior of Nitrilotriacetamide (NTAmide) Ligands in Thorium(IV) Extraction from Acidic Medium: Small-Angle Neutron Scattering, Fourier Transform Infrared, and Theoretical Studies

Two tripodal amides obtained from nitrilotriacetic acid with n-butyl and n-octyl alkyl chains (HBNTA(L_I) and HONTA(L_II), respectively) were studied for the extraction of Th(IV) ions from nitric acid medium. The effect of the diluent medium, i.e., n-dodecane alone and a mixture of n-dodecane and 1-decanol, onto aggregate formation were investigated using small angle neutron scattering (SANS) studies. In addition, the influence of the ligand structure, nitric acid, and Th(IV) loading onto ligand aggregation and third-phase formation tendency was discussed.The L_I/L_II exist as monomers (aggregarte radius for L_I: 6.0 Å; L_II:7.4 Å) in the presence of 1-decanol, whereas L_II forms dimers (aggregarte radius for L_II:9.3 Å; L_I does not dissolve in n-dodecane) in the absence of 1-decanol. The aggregation number increases for both the ligands after HNO₃ and Th(IV) loading. The maximum organic concentration (0.050 ± 0.004 M) of Th(IV) was reached without third-phase formation for 0.1 M L_I/L_II dissolved in 20% isodecanol +80% n-dodecane. The interaction of 1-decanol with L_II and HNO₃/Th(IV) with amidic oxygens of L_I/L_II results in shift of carbonyl stretching frequency, as shown by attenuated total reflectance-Fourier transform infrared (ATR-FTIR) studies. The structural and bonding information of the Th-L_I/L_II complex were derived from the density functional theoretical (DFT) studies. The molecular dynamics (MD) simulations suggested that the aggregation behavior of the ligand in the present system is governed by the population of hydrogen bonds by phase modifier around the ligand molecules. Although the theoretical studies suggested higher Gibbs free energy of complexation for Th⁴⁺ ions with L_I than L_II, the extraction was found to be higher with the latter, possibly due to the higher lipophilicity and solubility of the Th-L_II aggregate in the nonpolar media.

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.

Theoretically investigating the ability of phenanthroline derivatives to separate transuranic elements and their bonding properties

The bonding and separation properties of actinide Np3+, Pu3+, Am3+, and Cm3+ complexes formed with phenanthroline derivatives were studied using the DFT method.

Use of vibrational spectroscopy to identify the formation of neptunyl-neptunyl interactions: A paired Density Functional Theory and Raman spectroscopy study

Abstract: Actinyl-Actinyl interactions (AAIs) occur in pentavalent actinide systems, particularly for Np(V), and lead to complex vibrational signals that are challenging to analyze and interpret. Previous studies have focused on...

Disproportionation of the Uranyl(V) Coordination Complexes in Aqueous Solution through Outer-Sphere Electron Transfer

Among the linear actinyl(VI/V) cations, the uranyl(V) species are particularly intriguing because they are unstable and exhibit a unique behavior to undergo H⁺ promoted disproportionation in aqueous solution and form stable uranyl(VI) and U(IV) complexes. This study uses density functional theory (DFT) combined with the conductor-like polarizable continuum model approach to investigate [UO₂]^2+/+ to [U^IVO₂] reduction free energies (RFEs) and explores the stability of uranyl(V) complexes in aqueous solution through computing disproportionation free energies (DFEs) for an outer-sphere electron transfer process. In addition to the aqua complex (U1), another three commonly encountered ligands such as chloride (U2), acetate (U3), and carbonate (U4) in aqueous environmental conditions are taken into account. For the U1 complex, the computed 1e^- (V/IV) and 2e^- (VI/IV) RFEs are in good agreement with experiments. The computed DFEs reveal that the presence of H⁺ is imperative for the disproportionation to take place. Although the presence of the alkali cations favors the disproportionation to some extent, they cannot fully make the reaction thermodynamically feasible. For the anionic complexes, the high negative charge does not allow for the formation of a cation-cation encounter complex due to Coulombic repulsion. Furthermore, an additional factor is the ligand exchange reaction which is also an energy-demanding step. Therefore, the current study examined the Kern-Orlemann mechanism and our results validate the mechanism based on DFT computed DFEs and propose that for the anionic complexes, an outer-sphere electron transfer is highly probable and our computed protonation free energies further support this claim.

Role of O Substitution in Expanded Porphyrins on Uranyl Complexation: Orbital- and Density-Based Analyses

Search for new U(VI) sequestering macrocyclic ligands is an important area of research due to manifold applications. Besides hard- or soft-donor-based ligands, mixed-donor ligands are also gaining popularity in achieving optimized performances. However, how the combination of hard-soft-donor centers alters the bonding interactions with U(VI) is still not well-understood. Moreover, a consensus is yet to be reached on the nature and role of underlying covalent interactions in mixed N,O-donor ligands. In this work, using the relativistic density functional theory (DFT), we attempted to address these intriguing issues by investigating the subtle change in bonding characteristics of the uranyl ion upon binding with an expanded porphyrin, viz. sapphyrin, with subsequent O substitutions at the cavity. The results obtained from a range of modern analysis tools suggest that in the O-substituted sapphyrin variants, UO₂²⁺ prefers to bind with N over O, and an increase in the number of O-donor sites at the cavity prompts UO₂²⁺ to have a better interaction with the rest of the N-donor-centers. Although O donors are involved in more numbers of mixed molecular orbitals, the variation in the amplitude of overlap and the better σ-donation ability favor N to have stronger bonding interactions with uranyl. Molecular orbital (MO) and density of states (DOS) analyses show favorable participation of U(d), and the involvement of U(f) orbitals in bonding is of a low extent but non-negligible. Although electrostatic interaction dominates at U-O/N bonds in the equatorial plane, the quantum theory of atoms in molecules descriptors, MO analysis, and overlap-integral calculations confirm the presence of underlying near-degeneracy-driven covalent interactions.

Plutonium complexes in water: new approach to ab initio modeling

Geometry optimization and the electronic structure calculations of Pu Z+ complexes (Z = 3–6) in water solution have been performed, within the framework of the DMol3 and Relativistic Discrete-Variational (RDV) methods. For the simulation of Pu Z+ molecular environment in aqueous solution we used 22 and 32 water molecules randomly distributed around cation. To model the effect of bulk solvent environment we used COSMO (Conductor-like Screening Model) potential for water (ε = 78.54). The obtained results showed that this approach allows the modeling of water dissociation and the formation of hydrolysis products. Our previously suggested scheme for the calculation of interaction energies between selected fragments of multi-molecular systems provides the quantitative estimation of the interaction strengths between plutonium in various oxidation states and each ligand in the first and second coordination shells in water solution.

Highly efficient diglycolamide-functionalized dendrimers for the sequestration of tetravalent actinides: solvent extraction and theoretical studies

Two DGA-functionalized dendrimers were used for the extraction and bonding computations of tetravalent actinide ions Np⁴⁺ and Pu⁴⁺ from HNO₃ medium.

Evaluation of Chemical Bonding in Actinyl(VI/V) Oxo-Crown-Ether Complexes for Actinide Series from Uranium to Curium

The separation and management of nuclear waste is one of the problems that needs to be solved urgently, so finding a new radiation-proof and durable extractant to deal with nuclear waste is a difficult but desirable task. Since the successful isolation of the first pentavalent plutonium crown ether complex recently (Wang et al. CCS Chem. 2020, 2, 425-431), complexes with actinyl(V/VI) inserted into the cavity of 18-crown-6 ether (oxo-18C6), as well as their bonding character, need to be explored. Here we present a series of novel crown ether complexes containing actinyl(V/VI) and oxo-18C6 via computational prediction and analysis. On the basis of the calculations, actinyl(V/VI) are thermodynamically feasible and can be stabilized by oxo-18C6 ligand via six dative bonds between An ions and the oxo-18C6 O atoms in the "insertion" structure of [AnO₂(18C6)]^2+/+ complexes. The stability of actinyl(VI) species generally falls at minor actinides, ascribed to the reduced highest possible oxidation states of curium, which is mainly attributed to the mixing of bonding orbitals and non-bonding orbitals as well as the increase of occupation on partially 5f antibonding orbitals. It is found that the interactions between the actinyl(V/VI) and oxo-18C6 are mainly electronic interactions, with the well-known covalency contributions generally decreasing from uranium to curium due to energy degeneracy and spatial orbital contraction. This work would give a basic understanding of the coordination chemistry of actinyl(V/VI), which also provides inspirations on the design of new extractants for actinide separations.

Theoretical insights into selective separation of trivalent actinide and lanthanide by ester and amide ligands based on phenanthroline skeleton

Phenanthroline based ligands have shown potential performance for partitioning trivalent actinides from lanthanides. In this work, we have explored four ester and amide ligands based on the phenanthroline skeleton and elucidated the separation mechanism between Am(iii) and Eu(iii) ions. The molecular geometries and extraction reactions of the metal-ligand complexes were modeled by using scalar-relativistic density functional theory. The results show that the amide based ligands have stronger coordination ability with the metal ions than the corresponding ester based ligands. According to the thermodynamic results, ligands N,N'-diethyl-N,N'-ditolyl-2,9-diamide-1,10-phenanthroline (L2) and N,N'-(1,10-phenanthroline-2,9-diyl)bis(N-ethyl-P-methyl-N-(p-tolyl)phosphinic amide) (L4) appear to have the strongest complexing ability, which is supported by the result of electrostatic potential (ESP) and the M-OL bond orders. Moreover, ligand L2 has excellent selectivity for Am(iii)/Eu(iii) among the four ligands. Additionally, the bonding properties between the metal ions and the ligands reveal that the Am(iii)/Eu(iii) selectivity stems from the Am-N bonds with more covalent character, which is supported by the analysis of the hardness of the ligands and the bond orders. This work provides useful information for understanding the Am(iii)/Eu(iii) selectivity of phenanthroline derived ligands bearing ester and amide groups.

Investigation of Uranyl Sulfate Complexation under Hydrothermal Conditions by Quantitative Raman Spectroscopy and Density Functional Theory

Quantitative first and second formation constants of aqueous uranyl sulfate complexes were obtained from Raman spectra of solutions in fused silica capillary cells at 25 MPa, at temperatures ranging from 25 to 375 °C. Temperature-dependent values of the symmetric O-U-O vibrational frequencies of UO₂²⁺(aq), UO₂SO₄⁰(aq), and UO₂(SO₄)₂^2-(aq) were determined from the high-temperature spectra. Temperature-independent Raman scattering coefficients of UO₂²⁺(aq) were calculated directly from uranyl triflate spectra from 25 to 300 °C, while those of UO₂SO₄⁰(aq) and UO₂(SO₄)₂^2-(aq) were derived from spectroscopic data at 25 °C and concentrations calculated using the formation constants of Tian and Rao ( J. Chem. Thermodyn. 2009 , 41 , 569 - 574 ), together with the Specific Ion Interaction Theory (SIT) activity coefficient model. Chemical structures and vibrational frequencies predicted from Density Functional Theory (Gaussian 09) were employed to interpret the Raman spectra. Values of the cumulative formation constants ranged from log β₁ = 3.23 ± 0.08 and log β₂ = 4.22 ± 0.15 at 25 °C, to log β₁ = 12.35 ± 0.22 and log β₂ = 14.97 ± 0.02 at 350 °C. This is the first reported use of high-pressure fused silica capillary cells to determine formation constants of metal ligand complexes from their reduced isotropic Raman spectra under hydrothermal conditions.

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.

Ab initio molecular dynamics studies of hydroxide coordination of alkaline earth metals and uranyl

Ab initio molecular dynamics (AIMD) simulations of the Mg2+, Ca2+, Sr2+ and UO22+ ions in either a pure aqueous environment or an environment containing two hydroxide ions have been carried out at the density functional level of theory, employing the generalised gradient approximation via the PBE exchange-correlation functional. Calculated mean M-O bond lengths in the first solvation shell of the aquo systems compared very well to existing experimental and computational literature, with bond lengths well within values measured previously and coordination numbers in line with previously calculated values. When applied to systems containing additional hydroxide ions, the methodology revealed increased bond lengths in all systems. Proton transfer events (PTEs) were recorded and were found to be most prevalent in the strontium hydroxide systems, likely due to the low charge density of the ion and the consequent lack of hydroxide coordination. For all alkaline earths, intrashell PTEs which occurred outside of the first solvation shell were most prevalent. Only three PTEs were identified in the entire simulation data of the uranium dihydroxide system, indicating the clear impact of the increased charge density of the hexavalent uranium ion on the strength of metal-oxygen bonds in aqueous solution. Broadly, systems containing more charge dense ions were found to exhibit fewer PTEs than those containing ions of lower charge density.

Density functional theory (DFT) calculations of VI/V reduction potentials of uranyl coordination complexes in non-aqueous solutions

Of particular interest within the +6 uranium complexes is the linear uranyl(vi) cation and it forms numerous coordination complexes in solution and exhibits incongruent redox behavior depending on coordinating ligands. In this study, to determine the reduction potentials of uranyl complexes in non-aqueous solutions, a hybrid density functional theory (DFT) approach was used in which two different DFT functionals, B3LYP and M06, were applied. Bulk solvent effects were invoked through the conductor-like polarizable continuum model. The solute cavities were described with the united-atom Kohn-Sham (UAKS) cavity definition. Inside the cavity the dielectric constant matches the value of a vacuum and outside the cavity the dielectric constant value is the same as that of the solvent of interest, for example, dimethyl sulfoxide (DMSO), dimethylformamide (DMF), dichloromethane (DCM), acetonitrile and pyridine. With the help of the Nernst equation the calculated reduction potentials with respect to the ferrocene (Fc) reference electrode are converted into reduction free energies (RFEs). Uranyl complexes of organic ligands which range from mono- to hexa-dentate coordination modes were investigated in non-aqueous solutions of DMSO, DMF, DCM, acetonitrile and pyridine solutions. The effect of the spin-orbit correction and the reference electrode correction on the RFEs and various methods such as the direct method and the isodesmic reaction model were explored. Overall, our computational determination of RFEs of uranyl complexes in various non-aqueous solutions demonstrates that the RFEs can be obtained within ∼0.2 eV of experimental values.

Ion association with tetra-n-alkylammonium cations stabilizes higher-oxidation-state neptunium dioxocations

Extended-coordination sphere interactions between dissolved metals and other ions, including electrolyte cations, are not known to perturb the electrochemical behavior of metal cations in water. Herein, we report the stabilization of higher-oxidation-state Np dioxocations in aqueous chloride solutions by hydrophobic tetra-n-alkylammonium (TAA⁺) cations—an effect not exerted by fully hydrated Li⁺ cations under similar conditions. Experimental and molecular dynamics simulation results indicate that TAA⁺ cations not only drive enhanced coordination of anionic Cl^– ligands to Np^V/VI but also associate with the resulting Np complexes via non-covalent interactions, which together decrease the electrode potential of the Np^VI/Np^V couple by up to 220 mV (ΔΔG = −22.2 kJ mol⁻¹). Understanding the solvation-dependent interplay between electrolyte cations and metal–oxo species opens an avenue for controlling the formation and redox properties of metal complexes in solution. It also provides valuable mechanistic insights into actinide separation processes that widely use quaternary ammonium cations as extractants or in room temperature ionic liquids.

The electrochemical behaviour of redox-active metal cations foremost depends on the metal centre’s inner-sphere coordination environment. Here the authors show that electrolyte cations unexpectedly stabilize higher-oxidation-state neptunium dioxocations in water through extended-coordination sphere interactions.

Electronic Structure and Chemical Bonding of [AmO2(H2O) n ]2+/1

Systematic americyl-hydration cations were investigated theoretically to understand the electronic structures and bonding in [(AmO2)(H2O) n ]2+/1+ (n = 1-6). We obtained the binding energy using density functional theory methods with scalar relativistic and spin-orbit coupling effects. The geometric structures of these species have been investigated in aqueous solution via an implicit solvation model. Computational results reveal that the complexes of five equatorial water molecules coordinated to americyl ions are the most stable due to the enhanced ionic interactions between the AmO2 2+/1+ cation and multiple oxygen atoms as electron donors. As expected, Am-Owater bonds in such series are electrostatic in nature and contain a generally decreasing covalent character when hydration number increases.

Bonding trends across the series of tricarbonato-actinyl anions [(AnO2)(CO3)3]4- (An = U-Cm): the plutonium turn

Actinyl-tricarbonato anions [(AnO₂)(CO₃)₃]^4- (An = U-Cm) in various environments were investigated using theoretical approaches of quantum-mechanics, molecular-mechanics and cluster-models. Cations and solvent molecules in the 2^nd coordination sphere affect the equatorial An←O_eq bonds more than the axial An[triple bond, length as m-dash]O_ax bonds. Common actinide contraction is found for calculated and experimental axial bond lengths of ₉₂U to ₉₄Pu, though no longer for ₉₄Pu to ₉₆Cm. The tendency of U to Pu forming actinyl(vi) species dwindles away toward Cm, which already features the preferred An^III/Ln^III oxidation state of the later actinides and all lanthanides. The well known change from d-type to typical U-Pu-Cm type and then to f-type behavior is labeled as the plutonium turn, a phenomenon that is caused by f-orbital energy-decrease and f-orbital localization with increase of both nuclear charge and oxidation state, and a non-linear variation of effective f-electron population across the actinide series. Both orbital and configuration mixing and occupation of antibonding 5f type orbitals increase, weakening the AnO_ax bonds and reducing the highest possible oxidation states of the later actinides.

Systematic analysis of structural and topological properties: new insights into PuO2(H2O)n2+ (n = 1–6) complexes in the gas phase

In this study, the equilibrium, electronic structures, bonding and topological properties of PuO₂(H₂O)_n²⁺ (n = 1–6) complexes were systematically investigated.

Unusual intramolecular CH⋯O hydrogen bonding interaction between a sterically bulky amide and uranyl oxygen

An unusual intramolecular CH⋯O hydrogen bonding interaction between a sterically bulky amide and uranyl oxygen is found to selectively extract uranyl.

Exploring the Interaction Natures in Plutonyl (VI) Complexes with Topological Analyses of Electron Density

The interaction natures between Pu and different ligands in several plutonyl (VI) complexes are investigated by performing topological analyses of electron density. The geometrical structures in both gaseous and aqueous phases are obtained with B3LYP functional, and are generally in agreement with available theoretical and experimental results when combined with all-electron segmented all-electron relativistic contracted (SARC) basis set. The Pu- O y l bond orders show significant linear dependence on bond length and the charge of oxygen atoms in plutonyl moiety. The closed-shell interactions were identified for Pu-Ligand bonds in most complexes with quantum theory of atoms in molecules (QTAIM) analyses. Meanwhile, we found that some Pu-Ligand bonds, like Pu-OH(-), show weak covalent. The interactive nature of Pu-ligand bonds were revealed based on the interaction quantum atom (IQA) energy decomposition approach, and our results indicate that all Pu-Ligand interactions is dominated by the electrostatic attraction interaction as expected. Meanwhile it is also important to note that the quantum mechanical exchange-correlation contributions can not be ignored. By means of the non-covalent interaction (NCI) approach it has been found that some weak and repulsion interactions existed in plutonyl(VI) complexes, which can not be distinguished by QTAIM, can be successfully identified.

Paving the way for the synthesis of a series of divalent actinide complexes: a theoretical perspective

Recently, the +2 formal oxidation state in soluble molecular complexes for lanthanides (La-Nd, Sm-Lu) and actinides (Th and U) has been discovered [W. J. Evans, et al., J. Am. Chem. Soc., 2011, 133, 15914; J. Am. Chem. Soc., 2012, 134, 8420; J. Am. Chem. Soc., 2013, 135, 13310; Chem. Sci., 2015, 6, 517]. To explore the nature of the bonding and stabilities of the low-valent actinide complexes, a series of divalent actinide species, [AnCp'3](-) (An[double bond, length as m-dash]Th-Am, Cp' = [η(5)-C5H4(SiMe3)](-)) have been investigated in THF solution using scalar relativistic density functional theory. The electronic structures and electron affinity properties were systematically studied to identify the interactions between the +2 actinide ions and Cp' ligands. The ground state electron configurations for the [AnCp'3](-) species are [ThCp'3](-) 6d(2), [PaCp'3](-) 5f(2)6d(1), [UCp'3](-) 5f(3)6d(1), [NpCp'3](-) 5f(5), [PuCp'3](-) 5f(6), and [AmCp'3](-) 5f(7), respectively, according to the MO analysis. The total bonding energy decreases from the Th- to the Am-complex and the electrostatic interactions mainly dominate the bonding between the actinide atom and ligands. The electron affinity analysis suggests that the reduction reaction of AnCp'3→ [AnCp'3](-) should become increasingly facile across the actinide series from Th to Am, in accord with the known An(iii/ii) reduction potentials. This work expands the knowledge on the low oxidation state chemistry of actinides, and further motivates and guides the synthesis of related low oxidation state compounds of 5f elements.

New insights into the selectivity of four 1,10-phenanthroline-derived ligands toward the separation of trivalent actinides and lanthanides: a DFT based comparison study

To understand the correlation between the chemical structures of N-ligands and Ln/An separation capability, we selected four phenanthroline-derived ligands as examples to probe why some ligands work whereas others fail.

Acidity constants and redox potentials of uranyl ions in hydrothermal solutions

We report a first principles molecular dynamics (FPMD) study of the structures, acidity constants (pK_a) and redox potentials (E⁰) of uranyl (UO₂²⁺) from ambient conditions to 573 K.

Terminal U≡E (E = N, P, As, Sb, and Bi) bonds in uranium complexes: a theoretical perspective

The compound L-U-N [L = [N(CH2CH2NSiPr(i)3)3](3-), Pr(i) = CH(CH3)2] containing a terminal U-N triple bond has been synthesized and isolated successfully in experiments. To investigate the trend in the bonding nature of its pnictogen analogues, we have studied the L-U-E (E = N, P, As, Sb, and Bi) complexes using the scalar relativistic density functional theory. The terminal U-E multiple bond length increases in the order of U-N ≪ U-P < U-As < U-Sb < U-Bi, which can be supported by the hard and soft acids and bases (HSAB) theory. The U-E bond length, molecular orbital (MO), and natural bond orbital (NBO) reveal that the terminal U-E bonds should be genuine triple bonds containing one σ- and two π-bonding orbitals. Quantum theory of atoms in molecules (QTAIM) topological analysis and the electron localization function (ELF) suggest that the terminal U-E bond possesses covalent character and the covalency of U-E bonds decrease sharply when the terminal atom becomes heavier. This work presents a comparison about the bonding characteristic between the terminal U≡N bond and its heavier pnictogen (P, As, Sb, and Bi) analogues. It is expected that this work would shed light on the evaluation of the amount of 5f orbital participation in multiple bonds and further facilitate our deeper understanding of f-block elements.

Elucidating the structures and cooperative binding mechanism of cesium salts to the multitopic ion-pair receptor through density functional theory calculations

Quantum chemical studies predict the binding of Cesium salts to multitopic ion-pair receptor is through cooperative mechanism.

Effect of successive alkylation of N,N-dialkyl amides on the complexation behavior of uranium and thorium: solvent extraction, small angle neutron scattering, and computational studies

The effect of successive alkylation of the Cα atom adjacent to the carbonyl group in N,N-dialkyl amides (i.e., di(2-ethylhexyl)acetamide (D2EHAA), di(2-ethylhexyl)propionamide (D2EHPRA), di(2-ethylhexyl)isobutyramide (D2EHIBA), and di(2-ethylhexyl)pivalamide (D2EHPVA)) on the extraction behavior of hexavalent uranium (U(VI)) and tetravalent thorium (Th(IV)) ions has been investigated. These studies show that the extraction of Th(IV) is significantly suppressed compared to that of U(VI) with increased branching at the Cα atom adjacent to the carbonyl group. Small angle neutron scattering (SANS) studies showed an increased aggregation tendency in the presence of nitric acid and metal ions. D2EHAA showed more aggregation compared to its branched homologues, which explains its capacity for higher extraction of metal ions. These experimental observations were further supported by density function theory calculations, which provided structural evidence of differential binding affinities of these extractants for uranyl cations. The complexation process is primarily controlled by steric and electronic effects. Quantum chemical calculations showed that local hardness and polarizability can be extremely useful inputs for designing novel extractants relevant to a nuclear fuel cycle.

Theoretical investigation on multiple bonds in terminal actinide nitride complexes

A series of actinide (An) species of L-An-N compounds [An = Pa-Pu, L = [N(CH2CH2NSiPr(i)3)3](3-), Pr(i) = CH(CH3)2] have been investigated using scalar relativistic density functional theory (DFT) without considering spin-orbit coupling effects. The ground state geometric and electronic structures and natural bond orbital (NBO) analysis of actinide compounds were studied systematically in neutral and anionic forms. It was found that with increasing actinide atomic number, the bond length of terminal multiple An-N1 bond decreases, in accordance with the actinide contraction. The Mayer bond order of An-N1 decreases gradually from An = Pa to Pu, which indicates a decrease in bond strength. The terminal multiple bond for L-An-N compounds contains one σ and two π molecular orbitals, and the contributions of the 6d orbital to covalency are larger in magnitude than the 5f orbital based on NBO analysis and topological analysis of electron density. This work may help in understanding of the bonding nature of An-N multiple bonds and elucidating the trends and electronic structure changes across the actinide series. It can also shed light on the construction of novel An-N multiple bonds.

Investigations on preferential Pu(IV) extraction over U(VI) by N,N-dihexyloctanamide versus tri-n-butyl phosphate: evidence through small angle neutron scattering and DFT studies

Straight chain amide N,N-dihexyloctanamide (DHOA) has been found to be a promising alternative extractant to tri-n-butyl phosphate (TBP) for the reprocessing of irradiated uranium- and thorium-based fuels. Unlike TBP, DHOA displays preferential extraction of Pu(IV) over U(VI) at higher acidities (≥3 M HNO3) and poor extraction at lower acidities. Density functional theory (DFT) based calculations have been carried out on the structures and relative binding energies of U(VI) and Pu(IV) with the extractant molecules. These calculations suggest that the differential hardness of the two extractants is responsible for the preferential binding/complexation of TBP to uranyl, whereas the softer DHOA and the bulky nature of the extractant lead to stronger binding/complexation of DHOA to Pu(IV). In conjunction with quantum chemical calculations, small angle neutron scattering (SANS) measurements have also been performed for understanding the stoichiometry of the complex formed that leads to relatively lower extraction of Th(IV) (a model for Pu(IV)) as compared to U(VI) using DHOA and TBP as the extractants. The combined experimental and theoretical studies helped us to understand the superior complexation/extraction behavior of Pu(IV) over U(VI) with DHOA.

Promising density functional theory methods for predicting the structures of uranyl complexes

By examining the overall accuracy of different theoretical methods in predicting the U–X bond distances (of a series uranyl complexes), we found that both the global-hybrid meta-GGA functional of BB1K and the range-seperated LC-BLYP functional are fairly good (even better than the popular B3LYP method).