Density Functional Studies of Actinyl Aquo Complexes Studied Using Small-Core Effective Core Potentials and a Scalar Four-Component Relativistic Method

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).

Coordination and hydrolysis of plutonium ions in aqueous solution using Car-Parrinello molecular dynamics free energy simulations

Car-Parrinello molecular dynamics (CPMD) simulations have been used to examine the hydration structures, coordination energetics, and the first hydrolysis constants of Pu(3+), Pu(4+), PuO2(+), and PuO2(2+) ions in aqueous solution at 300 K. The coordination numbers and structural properties of the first shell of these ions are in good agreement with available experimental estimates. The hexavalent PuO2(2+) species is coordinated to five aquo ligands while the pentavalent PuO2(+) complex is coordinated to four aquo ligands. The Pu(3+) and Pu(4+) ions are both coordinated to eight water molecules. The first hydrolysis constants obtained for Pu(3+) and PuO2(2+) are 6.65 and 5.70, respectively, all within 0.3 pH unit of the experimental values (6.90 and 5.50, respectively). The hydrolysis constant of Pu(4+), 0.17, disagrees with the value of -0.60 in the most recent update of the Nuclear Energy Agency Thermochemical Database (NEA-TDB) but supports recent experimental findings. The hydrolysis constant of PuO2(+), 9.51, supports the experimental results of Bennett et al. [Radiochim. Acta 1992, 56, 15]. A correlation between the pKa of the first hydrolysis reaction and the effective charge of the plutonium center was found.

Density functional theory calculations of the redox potentials of actinide(VI)/actinide(V) couple in water

The measured redox potential of an actinide at an electrode surface involves the transfer of a single electron from the electrode surface on to the actinide center. Before electron transfer takes place, the complexing ligands and molecules of solvation need to become structurally arranged such that the electron transfer is at its most favorable. Following the electron transfer, there is further rearrangement to obtain the minimum energy structure for the reduced state. As such, there are three parts to the total energy cycle required to take the complex from its ground state oxidized form to its ground state reduced form. The first part of the energy comes from the structural rearrangement and solvation energies of the actinide species before the electron transfer or charge transfer process; the second part, the energy of the electron transfer; the third part, the energy required to reorganize the ligands and molecules of solvation around the reduced species. The time resolution of electrochemical techniques such as cyclic voltammetry is inadequate to determine to what extent bond and solvation rearrangement occurs before or after electron transfer; only for a couple to be classed as reversible is it fast in terms of the experimental time. Consequently, the partitioning of the energy theoretically is of importance to obtain good experimental agreement. Here we investigate the magnitude of the instantaneous charge transfer through calculating the fast one electron reduction energies of AnO2(H2O)n(2+), where An = U, Np, and Pu, for n = 4-6, in solution without inclusion of the structural optimization energy of the reduced form. These calculations have been performed using a number of DFT functionals, including the recently developed functionals of Zhao and Truhlar. The results obtained for calculated electron affinities in the aqueous phase for the AnO2(H2O)5(2+/+) couples are within 0.04 V of accepted experimental redox potentials, nearly an order of magnitude improvement on previous calculated standard potentials E(0) values, obtained using both DFT and high level multireference approaches.

Oxo-functionalization and reduction of the uranyl ion through lanthanide-element bond homolysis: synthetic, structural, and bonding analysis of a series of singly reduced uranyl-rare earth 5f1-4f(n) complexes

The heterobimetallic complexes [{UO2Ln(py)2(L)}2], combining a singly reduced uranyl cation and a rare-earth trication in a binucleating polypyrrole Schiff-base macrocycle (Pacman) and bridged through a uranyl oxo-group, have been prepared for Ln = Sc, Y, Ce, Sm, Eu, Gd, Dy, Er, Yb, and Lu. These compounds are formed by the single-electron reduction of the Pacman uranyl complex [UO2(py)(H2L)] by the rare-earth complexes Ln(III)(A)3 (A = N(SiMe3)2, OC6H3Bu(t)2-2,6) via homolysis of a Ln-A bond. The complexes are dimeric through mutual uranyl exo-oxo coordination but can be cleaved to form the trimetallic, monouranyl "ate" complexes [(py)3LiOUO(μ-X)Ln(py)(L)] by the addition of lithium halides. X-ray crystallographic structural characterization of many examples reveals very similar features for monomeric and dimeric series, the dimers containing an asymmetric U2O2 diamond core with shorter uranyl U═O distances than in the monomeric complexes. The synthesis by Ln(III)-A homolysis allows [5f(1)-4f(n)]2 and Li[5f(1)-4f(n)] complexes with oxo-bridged metal cations to be made for all possible 4f(n) configurations. Variable-temperature SQUID magnetometry and IR, NIR, and EPR spectroscopies on the complexes are utilized to provide a basis for the better understanding of the electronic structure of f-block complexes and their f-electron exchange interactions. Furthermore, the structures, calculated by restricted-core or all-electron methods, are compared along with the proposed mechanism of formation of the complexes. A strong antiferromagnetic coupling between the metal centers, mediated by the oxo groups, exists in the U(V)Sm(III) monomer, whereas the dimeric U(V)Dy(III) complex was found to show magnetic bistability at 3 K, a property required for the development of single-molecule magnets.

Force field development for actinyl ions via quantum mechanical calculations: an approach to account for many body solvation effects

Advances in computational algorithms and methodologies make it possible to use highly accurate quantum mechanical calculations to develop force fields (pair-wise additive intermolecular potentials) for condensed phase simulations. Despite these advances, this approach faces numerous hurdles for the case of actinyl ions, AcO2(n+) (high-oxidation-state actinide dioxo cations), mainly due to the complex electronic structure resulting from an interplay of s, p, d, and f valence orbitals. Traditional methods use a pair of molecules (“dimer”) to generate a potential energy surface (PES) for force field parametrization based on the assumption that many body polarization effects are negligible. We show that this is a poor approximation for aqueous phase uranyl ions and present an alternative approach for the development of actinyl ion force fields that includes important many body solvation effects. Force fields are developed for the UO2(2+) ion with the SPC/Fw, TIP3P, TIP4P, and TIP5P water models and are validated by carrying out detailed molecular simulations on the uranyl aqua ion, one of the most characterized actinide systems. It is shown that the force fields faithfully reproduce available experimental structural data and hydration free energies. Failure to account for solvation effects when generating PES leads to overbinding between UO2(2+) and water, resulting in incorrect hydration free energies and coordination numbers. A detailed analysis of arrangement of water molecules in the first and second solvation shell of UO2(2+) is presented. The use of a simple functional form involving the sum of Lennard-Jones + Coulomb potentials makes the new force field compatible with a large number of available molecular simulation engines and common force fields.

Density functional theory investigation of the redox properties of tricyclopentadienyl- and phospholyluranium(IV) chloride complexes

The redox behavior of tricyclopentadienyl- and phospholyluranium(IV) chloride complexes L(3)UCl with L = C(5)H(5) (Cp), C(5)H(4)Me (MeCp), C(5)H(4)SiMe(3) (TMSCp), C(5)H(4)(t)Bu ((t)BuCp), C(5)Me(5) (Cp*), and C(4)Me(4)P (tmp), has been investigated using relativistic density functional theory calculations, with the solvent being taken into account using the conductor-like screening model. A very good linear correlation (r(2) = 0.99) has been obtained between the computed electron affinities of the L(3)UCl complexes and the experimental half-wave reduction potentials E(1/2) related to the U(IV)/U(III) redox systems. From a computational point of view, our study confirms the crucial importance of spin-orbit coupling and solvent corrections and the use of an extended basis set in order to achieve the best experiment-theory agreement. Considering oxidation of the uranium(IV) complexes, the instability of the uranium(V) derivatives [L(3)UCl](+) is revealed, in agreement with experimental electrochemical findings. The driving roles of both the electron-donating ability of the L ligand and the U 5f orbitals on the redox properties of the complexes are brought to light. Interestingly, we found and explained the excellent correlation between variations of the uranium Hirschfeld charges following U(IV)/U(III) electron capture and E(1/2). In addition, this work allowed one to estimate theoretically the half-wave reduction potential of [Cp*(3)UCl].

Production and properties of transuranium elements

We summarize historical perspective of the transuranium elements, neptunium (Np) through lawrencium (Lr), and recent progress on production, and nuclear and chemical properties of these elements. Exotic decay properties of heavy nuclei are also introduced. Chemical properties of transuranium elements in aqueous and solid states are summarized based on the actinide concept.

On the "yl" bond weakening in uranyl(VI) coordination complexes

The U-O(yl) triple bonds in the UO(2)(2+) aquo ion are known to be weakened by replacing the first shell water with organic or inorganic ligands. Weakening of the U-O(yl) bond may enhance the reactivity of "yl" oxygens and uranyl(VI) cation-cation interactions. Density functional theory calculations as well as previously published vibrational spectroscopic data have been used to study the origin of the U-O(yl) bond weakening in uranyl(VI) coordination complexes. Natural population analyses (NPA) revealed that the electron localization on the O(yl) 2p orbital is a direct measure of the U-O(yl) bond weakening, indicating that the bond weakening is correlated to the weakening of the U-O(yl) covalent bond and not that of the ionic bond. The Mulliken analysis gives poor results for uranium to ligand electron partitioning and is thus unreliable. Further analyses of molecular orbitals near the highest occupied molecular orbital (HOMO) show that both the σ and π donating abilities of the ligands may account for the U-O(yl) bond weakening. The mechanism of the bond weakening varies with coordinating ligand so that each case needs to be examined independently.