How accurate are electronic structure methods for actinoid chemistry?

Predissociation-based measurements of bond dissociation energies: US2, OUS, and USe

The uranium-containing molecules US2, OUS, and USe have been investigated using a pulsed laser ablation supersonic beam molecular source with time-of-flight mass spectrometric detection. Spectra have been recorded using the resonant two-photon ionization method over the spectroscopic range from 277 to 238 nm. These species have a myriad of excited electronic states in this spectroscopic region, leading to spectra that are highly congested and appear quasicontinuous. Sharp predissociation thresholds are observed, allowing precise bond dissociation energies to be measured. In the case of the triatomic molecules, it was necessary to use one laser for excitation and a delayed laser for ionization in order to observe a sharp predissociation threshold that allowed a precise bond dissociation energy to be measured. The resulting thermochemical values are D0(SU-S) = 4.910 ± 0.003 eV, D0(OU-S) = 5.035 ± 0.004 eV, and D0(USe) = 4.609 ± 0.009 eV. These results provide the first measurement of D0(USe) and reduce the error limits in the previous values of D0(SU-S) and D0(OU-S) by a factor of more than 70.

Gas-phase synthesis of [OU-X]+ (X = Cl, Br and I) from a UO22+ precursor using ion-molecule reactions and an [OUCH]+ intermediate

Difficulty in the preparation of gas-phase ions that include U in middle oxidation states(III,IV) have hampered efforts to investigate intrinsic structure, bonding and reactivity of model species. Our group has used preparative tandem mass spectrometry (PTMS) to synthesize a gas-phase U-methylidyne species, [OUCH]+, by elimination of CO from [UO2(CCH)]+ [M. J. van Stipdonk, I. J. Tatosian, A. C. Iacovino, A. R. Bubas, L. Metzler, M. C. Sherman and A. Somogyi, J. Am. Soc. Mass Spectrom., 2019, 30, 796-805], which has been used as an intermediate to create products such as [OUN]+ and [OUS]+ by ion-molecule reactions. Here, we investigated the reactions of [OUCH]+ with a range of alkyl halides to determine whether the methylidyne is a also a useful intermediate for production and study of the oxy-halide ions [OUX]+, where X = Cl, Br and I, formally U(IV) species for which intrinsic reactivity data is relatively scarce. Our experiments demonstrate that [OUX]+ is the dominant product ion generated by reaction [OUCH]+ with neutral regents such as CH3Cl, CH3CH2Br and CH2CHCH2I.

Multiconfiguration Pair-Density Functional Theory for Vertical Excitation Energies in Actinide Molecules

Modeling actinides with electronic structure theories is challenging because these systems present a strong ligand field and metal-ligand covalency. We systematically investigate the effectiveness of pair-density functional theory (PDFT) for the calculation of vertical excitation energies in An(III), [AnIIICl6]3-, and [AnVIO2]2+ (An = U, Np, Pu, and Am). We compare the performance of PDFT, hybrid PDFT, and multistate PDFT with traditional active-space methods followed by perturbation theory, like multistate CASPT2, and with experimental data. Overall, multistate PDFT gives quantitative agreement with multistate CASPT2 at a significantly reduced computational cost.

Route to Chemical Accuracy for Computational Uranium Thermochemistry

Benchmark spinor-based relativistic coupled-cluster calculations for the ionization energies of the uranium atom, the uranium monoxide molecule (UO), and the uranium dioxide molecule (UO₂) and for the bond dissociation energies of UO and UO₂ are reported. The accuracy of these calculations in the treatments of relativistic, electron-correlation, and basis-set effects is analyzed. The intrinsic convergence of the computed results and the favorable comparison with the experimental values demonstrate the unique applicability of the spinor representation of quantum-chemical methods to open-shell uranium-containing atomic and molecular species with uranium oxidation states ranging from U(0) to U(V).

Reactions of U+ with H2, D2, and HD Studied by Guided Ion Beam Tandem Mass Spectrometry and Theory

The kinetic energy-dependent reactions of the atomic actinide uranium cation (U⁺) with H₂, D₂, and HD were examined by guided ion beam tandem mass spectrometry. An average 0 K bond dissociation energy of D₀(U⁺ - H) = 2.48 ± 0.06 eV is obtained by analysis of the endothermic product ion cross sections. Quantum chemistry calculations were performed for comparison with experimental thermochemistry, including high-level CASSCF-CASPT2-RASSI calculations of the spin-orbit corrections. CCSD(T) and the CASSCF levels show excellent agreement with experiment, whereas B3LYP and PBE0 slightly overestimate and the M06 approach badly underestimates the bond energy for UH⁺. Theory was also used to investigate the electronic structures of the reaction intermediates and potential energy surfaces. The experimental product branching ratio for the reaction of U⁺ with HD indicates that these reactions occur primarily via a direct reaction mechanism, despite the presence of a deep-well for UH₂⁺ formation according to theory. The reactivity and hydride bond energy for U⁺ are compared with those for transition metal, lanthanide, and actinide cations, and periodic trends are discussed. These comparisons suggest that the 5f electrons on uranium are largely core and uninvolved in the reactive chemistry.

Guided Ion Beam and Quantum Chemical Investigation of the Thermochemistry of Thorium Dioxide Cations: Thermodynamic Evidence for Participation of f Orbitals in Bonding

Kinetic energy dependent reactions of ThO⁺ with O₂ are studied using a guided ion beam tandem mass spectrometer. The formation of ThO₂⁺ in the reaction of ThO⁺ with O₂ is observed to be slightly endothermic and also exhibits two obvious features in the cross section. These kinetic energy dependent cross sections were modeled to determine a 0 K bond dissociation energy of D₀(OTh⁺-O) = 4.94 ± 0.06 eV. This value is slightly larger but within experimental uncertainty of less precise previously reported experimental values. The higher energy feature in the ThO₂⁺ cross section was also analyzed and suggests formation of an excited state of the product ion lying 3.1 ± 0.2 eV above the ground state. Additionally, the thermochemistry of ThO₂⁺ was explored by quantum chemical calculations, including a full Feller-Peterson-Dixon (FPD) composite approach with correlation contributions up to CCSDT(Q) and four-component spin-orbit corrections, as well as more approximate CCSD(T) calculations including semiempirical estimates of spin-orbit energy contributions. The FPD approach predicts D₀(OTh⁺-O) = 4.87 ± 0.04 eV, in good agreement with the experimental value. Analogous FPD results for ThO⁺, ThO, and ThO₂ are also presented, including ionization energies for both ThO and ThO₂. The ThO₂⁺ bond energy is larger than those of its transition metal congeners, TiO₂⁺ and ZrO₂⁺, which can be attributed partially to an actinide contraction, but also to contributions from the participation of f orbitals on thorium that are unavailable to the transition metal systems.

Insight into the nature of M–C bonding in the lanthanide/actinide-biscarbene complexes: a theoretical perspective

The An/Ln–C bonding nature was explored using relativistic theory. Inclusion of Np and Pu extends understanding to later actinides bonding.

Actinide covalency measured by pulsed electron paramagnetic resonance spectroscopy

Our knowledge of actinide chemical bonds lags far behind our understanding of the bonding regimes of any other series of elements. This is a major issue given the technological as well as fundamental importance of f-block elements. Some key chemical differences between actinides and lanthanides-and between different actinides-can be ascribed to minor differences in covalency, that is, the degree to which electrons are shared between the f-block element and coordinated ligands. Yet there are almost no direct measures of such covalency for actinides. Here we report the first pulsed electron paramagnetic resonance spectra of actinide compounds. We apply the hyperfine sublevel correlation technique to quantify the electron-spin density at ligand nuclei (via the weak hyperfine interactions) in molecular thorium(III) and uranium(III) species and therefore the extent of covalency. Such information will be important in developing our understanding of the chemical bonding, and therefore the reactivity, of actinides.

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.

4-Component relativistic calculations of L3ionization and excitations for the isoelectronic species UO22+, OUN+and UN2

4-Component relativistic calculations explore uranium 2p_3/2ionization and excitation in the isoelectronic series UO₂²⁺, OUN⁺and UN₂.

Theoretical insight into the binding affinity enhancement of serine with the uranyl ion through phosphorylation

Various type complexes of uranyl ion with serine and phosphoserine were investigated which showed phosphorylation induced binding affinity enhancement of serine with the uranyl ion.

Dissecting the cation–cation interaction between two uranyl units

A theoretical study of the CCIs between two bare uranyl units and their spectroscopic characterization.

Theoretical study of the electronic spectra of neutral and cationic NpO and NpO2

The electronic spectra of neutral NpO and NpO2 as well as of their mono- (NpO(+), NpO2(+)) and dications (NpO(2+), NpO2(2+)) were studied using multiconfigurational relativistic quantum chemical calculations at the complete active space self-consistent field/CASPT2 level of theory taking into account spin-orbit coupling. The active space included 16 orbitals: all the 7s, 6d, and 5f orbitals of neptunium together with selected orbitals of oxygen. The vertical excitation energies on the ground state geometries have been computed up to ca. 35,000 cm(-1). The gas-phase electronic spectra were evaluated on the basis of the computed Einstein coefficients at 298 K and 3000 K. The computed vertical transition energies show good agreement with previous condensed-phase results on NpO2(+) and NpO2(2+).

U(III)-CN versus U(IV)-NC coordination in tris(silylamide) complexes

Treatment of the metallacycle [UN*2(N,C)] [N* = N(SiMe3)2; N,C = CH2SiMe2N(SiMe3)] with [HNEt3][BPh4], [HNEt3]Cl, and [pyH][OTf] (OTf = OSO2CF3) gave the cationic compound [UN*3][BPh4] (1) and the neutral complexes [UN*3X] [X = Cl (3), OTf (4)], respectively. The dinuclear complex [{UN*(μ-N,C)(μ-OTf)}2] (5) and its tetrahydrofuran (THF) adduct [{UN*(N,C)(THF)(μ-OTf)}2] (6) were obtained by thermal decomposition of 4. The successive addition of NEt4CN or KCN to 1 led to the formation of the cyanido-bridged dinuclear compound [(UN*3)2(μ-CN)][BPh4] (7) and the mononuclear mono- and bis(cyanide) complexes [UN*3(CN)] (2) and [M][UN*3(CN)2] [M = NEt4 (8), K(THF)4 (9)], while crystals of [K(18-crown-6)][UN*3(CN)2] (10) were obtained by the oxidation of [K(18-crown-6)][UN*3(CN)] with pyridine N-oxide. The THF adduct of 1, [UN*3(THF)][BPh4], and complexes 2-7, 9 and 10 were characterized by their X-ray crystal structure. In contrast to their U(III) analogues [NMe4][UN*3(CN)] and [K(18-crown-6)]2[UN*3(CN)2] in which the CN anions are coordinated to the metal center via the C atom, complexes 2 and 9 exhibit the isocyanide U-NC coordination mode of the cyanide ligand. This U(III)/U(IV) differentiation has been analyzed using density functional theory calculations. The observed preferential coordinations are well explained considering the electronic structures of the different species and metal-ligand bonding energies. A comparison of the different quantum descriptors, i.e., bond orders, NPA/QTAIM data, and energy decomposition analysis, has allowed highlighting of the subtle balance between covalent, ionic, and steric factors that govern the U-CN/NC bonding.

Plutonium and transplutonium element trioxides: molecular structures, chemical bonding, and isomers

The stability of the An·3O (An = Pu through Cf) lowest-energy isomers should decrease in the series Pu(vi)O₃ > Am(v)O₃ ≈ Bk(v)O₃ > Cm(iii)O(O₂) > Cf(iii)O(O₂).

Spectroscopy and structure of the simplest actinide bonds

Understanding the influence of electrons in partially filled f- and d-orbitals on bonding and reactivity is a key issue for actinide chemistry. This question can be investigated by using a combination of well-defined experimental measurements and theoretical calculations. Gas phase spectroscopic data are particularly valuable for the evaluation of theoretical models. Consequently, the primary objectives of our research have been to obtain gas phase spectra for small actinide molecules. To complement the experimental effort, we are investigating the potential for using relativistic ab initio calculations and semiempirical models to predict and interpret the electronic energy level patterns for f-element compounds. Multiple resonance spectroscopy and jet cooling techniques have been used to unravel the complex electronic spectra of Th and U compounds. Recent results for fluorides, sulfides, and nitrides are discussed.

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.

Synthesis and properties of uranium sulfide cations. An evaluation of the stability of thiouranyl, {S═U═S}2+

Atomic uranium cations, U(+) and U(2+), reacted with the facile sulfur-atom donor OCS to produce several monopositive and dipositive uranium sulfide species containing up to four sulfur atoms. Sequential abstraction of two sulfur atoms by U(2+) resulted in US2(2+); density functional theory computations indicate that the ground-state structure for this species is side-on η(2)-S2 triangular US2(2+), with the linear thiouranyl isomer, {S═U(VI)═S}(2+), some 171 kJ mol(-1) higher in energy. The result that the linear thiouranyl structure is a local minimum at a moderate energy suggests that it should be feasible to stabilize this moiety in molecular compounds.

Structure and other molecular properties of actinide trichlorides AnCl3 (An = Th-Cm)

The ground-state molecular properties of the trichlorides of light actinides (An = Th-Cm) have been predicted by state-of-the-art quantum chemical calculations. The ground electronic states have been determined by multireference calculations at the CASPT2 level including both scalar and spin-orbit relativistic effects. These studies supported the expected single-configuration character of ThCl3 and CmCl3 with their well-defined 6dσ/7s hybrid and 5f(7) configurations, respectively. In contrast, the intermediate actinides (PaCl3-AmCl3) with partly filled 5f shells have numerous very low-lying excited states and consequently a mixed character of the spin-orbit ground states. Apart from the planar ThCl3 the ground-state molecular geometries proved to be pyramidal with C(3v) symmetry. The gradually decreasing An-Cl bond distances reveal the actinide contraction known for the atomic and ionic radii of these actinide atoms. Other ground-state molecular properties as vibrational frequencies and natural charges have been obtained by density functional theory calculations using the B3LYP exchange-correlation functional in conjunction with small-core relativistic energy-consistent pseudopotentials for the actinides.

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.

Cation-π interaction of alkali metal ions with C24 fullerene: a DFT study

Using first principle calculations, we investigated cation-π interactions between alkali cations (Li(+), Na(+), and K(+)) and pristine C(24) or doped fullerenes of BC(23), and NC(23). The most suitable adsorption site is found to be atop the center of a six-membered ring of the exterior surface of C(24) molecule. Interaction energies of these cations decreased in the order: Li(+) > Na(+) > K(+), with values of -31.82, -22.36, and -15.68 kcal mol(-1), respectively. It was shown that the interaction energies are increased and decreased by impurity doping of B and N atoms in adjacent wall of adsorption site, depending on electron donating or receptivity of the doping atoms.