Computational study of the interactions of tetravalent actinides (An = Th-Pu) with the α-Fe13 Keggin cluster

In recent years, evidence has emerged that actinide (An) uptake at the enhanced actinide removal plant (EARP) at the UK's Sellafield nuclear site occurs via An interactions with an α-Fe13 Keggin molecular cluster during ferrihydrite formation. We here study theoretically the substitution of aquo complexes of the actinides Th-Pu onto a Na-decorated α-Fe13 Keggin cluster using DFT at the PBE0 level. The optimised Pu-O and Pu-Fe distances are in good agreement with experiment and Na/An substitutions are significantly favourable energetically, becoming more so across the early 5f series, with the smallest and largest ΔrG° being for Th and Pu at -335.7 kJ mol-1 and -396.0 kJ mol-1 respectively. There is strong correlation between the substitution reaction energy and the ionic radii of the actinides (Δrε0R2 = 0.97 and ΔrG° R2 = 0.91), suggesting that the principal An-Keggin binding mode is ionic. Notwithstanding this result, Mulliken and natural population analyses reveal that covalency increases from Th-Pu in these systems, supported by analysis of the occupied Kohn-Sham molecular orbitals where enhanced An(5f)-O(2p) overlap is observed in the Np and Pu systems. By contrast, quantum theory of atoms in molecules analysis shows that U-Keggin binding is the most covalent among the five actinides, in keeping with previous studies.

Understanding the Complexation Behavior of Carbamoylphosphine Oxide Ligands with Representative f-Block Elements

The complexation behavior of carbamoylmethylphosphine oxide ligands (CMPO), a bifunctional phosphine oxide, and their substituted derivatives with Ce(III), Eu(III), Th(IV), U(VI), and Am(III) was probed at the density functional theory (DFT) level. The enhanced extraction of trivalent rare earth elements by the 2-diphenylphosphinylethyl derivative over the conventional CMPO ligand is identified due to the availability of an additional P═O donor group in the former. In addition, the orbital and dispersive interactions play a vital role in the preference of Th(IV) over U(VI) during extraction using CMPO ligands. The better complexing ability of ligands having long alkyl chain substituents at the P atom is justified due to the observed enhanced dispersion interactions in these systems.

Pressure-Induced Coordination Number Transition in Lanthanide Mellitate Coordination Polymers: Structure and Spectroscopy

High external pressure is found to induce a non-coordinated water molecule to bond to cerium in a previously studied mellitate coordination polymer, as determined by high-pressure single-crystal X-ray diffraction, resulting in a coordination number transition at 3.85 GPa from 9 to 9.5 where half the cerium ions are 10-coordinate. Also, bond length changes due to increased pressure are experimentally measured, whereas the cerium-carboxylate bond lengths overall change by -0.004(9) Å/GPa, the cerium-water bonds by -0.016(3) Å/GPa, and cerium-oxygen bonds overall by -0.010(6) Å/GPa, which corresponds well with theoretical bond length decreases determined for similar compounds. The high-pressure absorbance spectra of the analogous neodymium mellitate are examined and compared with the structural changes observed.

Tris(carbene)borates; alternatives to cyclopentadienyls in organolanthanide chemistry

The chemistry of the tris-carbene anion phenyltris(3-alkyl-imidazoline-2-yliden-1-yl)borate, [C3Me]- ligand, is initiated for f-block metal cations. Neutral, molecular complexes of the form Ln(C3)2I are formed for cerium(III), while a separated ion pair [Ln(C3)2]I forms for ytterbium(III). DFT/QTAIM computational analyses of the complexes and related tridentate tris(pyrazolyl)borate (Tp) - supported analogs demonstrates the anticipated strength of the σ donation and confirms greater covalency in the metal-carbon bonds of the [C3Me]- complexes in comparison with those in the TpMe,Me complexes. The DFT calculations demonstrate the crucial role of THF solvent in accurately reproducing the contrasting molecular and ion-pair geometries observed experimentally for the Ce and Yb complexes.

Covalent bond shortening and distortion induced by pressurization of thorium, uranium, and neptunium tetrakis aryloxides

Covalency involving the 5f orbitals is regularly invoked to explain the reactivity, structure and spectroscopic properties of the actinides, but the ionic versus covalent nature of metal-ligand bonding in actinide complexes remains controversial. The tetrakis 2,6-di-tert-butylphenoxide complexes of Th, U and Np form an isostructural series of crystal structures containing approximately tetrahedral MO₄ cores. We show that up to 3 GPa the Th and U crystal structures show negative linear compressibility as the OMO angles distort. At 3 GPa the angles snap back to their original values, reverting to a tetrahedral geometry with an abrupt shortening of the M-O distances by up to 0.1 Å. The Np complex shows similar but smaller effects, transforming above 2.4 GPa. Electronic structure calculations associate the M-O bond shortening with a change in covalency resulting from increased contributions to the M-O bonding by the metal 6d and 5f orbitals, the combination promoting MO₄ flexibility at little cost in energy.

Contrasting behaviour under pressure reveals the reasons for pyramidalization in tris(amido)uranium(III) and tris(arylthiolate) uranium(III) molecules

A range of reasons has been suggested for why many low-coordinate complexes across the periodic table exhibit a geometry that is bent, rather a higher symmetry that would best separate the ligands. The dominating reason or reasons are still debated. Here we show that two pyramidal UX₃ molecules, in which X is a bulky anionic ligand, show opposite behaviour upon pressurisation in the solid state. UN″₃ (UN3, N″ = N(SiMe₃)₂) increases in pyramidalization between ambient pressure and 4.08 GPa, while U(SAr)₃ (US3, SAr = S-C₆H₂-^tBu₃−2,4,6) undergoes pressure-induced planarization. This capacity for planarization enables the use of X-ray structural and computational analyses to explore the four hypotheses normally put forward for this pyramidalization. The pyramidality of UN3, which increases with pressure, is favoured by increased dipole and reduction in molecular volume, the two factors outweighing the slight increase in metal-ligand agostic interactions that would be formed if it was planar. The ambient pressure pyramidal geometry of US3 is favoured by the induced dipole moment and agostic bond formation but these are weaker drivers than in UN3; the pressure-induced planarization of US3 is promoted by the lower molecular volume of US3 when it is planar compared to when it is pyramidal.

The reasons for which many low-coordinate complexes exhibit bent geometry, rather than a higher symmetry, are still under debate. Here, the authors use high-pressure crystallography to examine whether low-coordinate f-block molecules become more planar or pyramidal under pressure; which happens is dictated by the dipole moment of the complex and the volume of the planar form.

Harnessing Greenhouse Gases Absorption by Doped Fullerenes with Externally Oriented Electric Field

In this work, a theoretical investigation of the effects caused by the doping of C₂₀ with silicon (Si) atom as well as the adsorption of CO, CO₂ and N₂ gases to C₂₀ and C₁₉Si fullerenes was carried out. In concordance with previous studies, it was found that the choice of the doping site can control the structural, electronic, and energetic characteristics of the C₁₉Si system. The ability of C₂₀ and C₁₉Si to adsorb CO, CO₂ and N₂ gas molecules was evaluated. In order to modulate the process of adsorption of these chemical species to C₁₉Si, an externally oriented electric field was included in the theoretical calculations. It was observed that C₁₉Si is highly selective with respect to CO adsorption. Upon the increase of the electric field intensity the adsorption energy was magnified correspondingly and that the interaction between CO and C₁₉Si changes in nature from a physical adsorption to a partial covalent character interaction.

Lanthanide Complexes Containing a Terminal Ln═O Oxo Bond: Revealing Higher Stability of Tetravalent Praseodymium versus Terbium

We report on the reactivity of gas-phase lanthanide-oxide nitrate complexes, [Ln(O)(NO₃)₃]^- (denoted LnO²⁺), produced via elimination of NO₂^• from trivalent [Ln^III(NO₃)₄]^- (Ln = Ce, Pr, Nd, Sm, Tb, Dy). These complexes feature a Ln^III-O^• oxyl, a Ln^IV═O oxo, or an intermediate Ln^III/IV oxyl/oxo bond, depending on the accessibility of the tetravalent Ln^IV state. Hydrogen atom abstraction reactivity of the LnO²⁺ complexes to form unambiguously trivalent [Ln^III(OH)(NO₃)₃]^- reveals the nature of the oxide bond. The result of slower reactivity of PrO²⁺ versus TbO²⁺ is considered to indicate higher stability of the tetravalent praseodymium-oxo, Pr^IV═O, versus Tb^IV═O. This is the first report of Pr^IV as more stable than Tb^IV, which is discussed with respect to ionization potentials, standard electrode potentials, atomic promotion energies, and oxo bond covalency via 4f- and/or 5d-orbital participation.

Stable Chelation of the Uranyl Ion by Acyclic Hexadentate Ligands: Potential Applications for 230U Targeted α-Therapy

Uranium-230 is an α-emitting radionuclide with favorable properties for use in targeted α-therapy (TAT), a type of nuclear medicine that harnesses α particles to eradicate cancer cells. To successfully implement this radionuclide for TAT, a bifunctional chelator that can stably bind uranium in vivo is required. To address this need, we investigated the acyclic ligands H₂dedpa, H₂CHXdedpa, H₂hox, and H₂CHXhox as uranium chelators. The stability constants of these ligands with UO₂²⁺ were measured via spectrophotometric titrations, revealing log β_ML values that are greater than 18 and 26 for the "pa" and "hox" chelators, respectively, signifying that the resulting complexes are exceedingly stable. In addition, the UO₂²⁺ complexes were structurally characterized by NMR spectroscopy and X-ray crystallography. Crystallographic studies reveal that all six donor atoms of the four ligands span the equatorial plane of the UO₂²⁺ ion, giving rise to coordinatively saturated complexes that exclude solvent molecules. To further understand the enhanced thermodynamic stabilities of the "hox" chelators over the "pa" chelators, density functional theory (DFT) calculations were employed. The use of the quantum theory of atoms in molecules revealed that the extent of covalency between all four ligands and UO₂²⁺ was similar. Analysis of the DFT-computed ligand strain energy suggested that this factor was the major driving force for the higher thermodynamic stability of the "hox" ligands. To assess the suitability of these ligands for use with ²³⁰U TAT in vivo, their kinetic stabilities were probed by challenging the UO₂²⁺ complexes with the bone model hydroxyapatite (HAP) and human plasma. All four complexes were >95% stable in human plasma for 14 days, whereas in the presence of HAP, only the complexes of H₂CHXdedpa and H₂hox remained >80% intact over the same period. As a final validation of the suitability of these ligands for radiotherapy applications, the in vivo biodistribution of their UO₂²⁺ complexes was determined in mice in comparison to unchelated [UO₂(NO₃)₂]. In contrast to [UO₂(NO₃)₂], which displays significant bone uptake, all four ligand complexes do not accumulate in the skeletal system, indicating that they remain stable in vivo. Collectively, these studies suggest that the equatorial-spanning ligands H₂dedpa, H₂CHXdedpa, H₂hox, and H₂CHXhox are highly promising candidates for use in ²³⁰U TAT.

Actinide Coordination Chemistry on Surfaces: Synthesis, Manipulation, and Properties of Thorium Bis(porphyrinato) Complexes

Actinide-based metal-organic complexes and coordination architectures encompass intriguing properties and functionalities but are still largely unexplored on surfaces. We introduce the in situ synthesis of actinide tetrapyrrole complexes under ultrahigh-vacuum conditions, on both a metallic support and a 2D material. Specifically, exposure of a tetraphenylporphyrin (TPP) multilayer to an elemental beam of thorium followed by a temperature-programmed reaction and desorption of surplus molecules yields bis(porphyrinato)thorium (Th(TPP)₂) assemblies on Ag(111) and hexagonal boron nitride/Cu(111). A multimethod characterization including X-ray photoelectron spectroscopy, scanning tunneling microscopy, temperature-programmed desorption, and complementary density functional theory modeling provides insights into conformational and electronic properties. Supramolecular assemblies of Th(TPP)₂ as well as individual double-deckers are addressed with submolecular precision, e.g., demonstrating the reversible rotation of the top porphyrin in Th(TPP)₂ by molecular manipulation. Our findings thus demonstrate prospects for actinide-based functional nanoarchitectures.

Covalency in AnCl3 (An = Th-No)

The geometric and electronic structures of AnCl3 are studied computationally using scalar relativistic, hybrid density functional theory (PBE0). The An-Cl bond lengths generally decrease across the 5f series, although there is a slight lengthening from Fm-Cl to No-Cl as the metal ions display increasing M(ii) character. Covalency in the An-Cl bond is studied using a wide range of metrics drawn from the Natural Bond Orbital, Natural Resonance Theory and Quantum Theory of Atoms-in-Molecules (QTAIM) methods, including bond order, orbital composition, orbital overlap and electron density topology data. Most metrics agree that the later An-Cl bonds are less ionic than might be anticipated on the basis of trends in the first half of the series, due to energy degeneracy-driven covalency in the β spin manifold; for example, the An-Cl QTAIM delocalisation index (bond order) for MdCl3 (0.88) is almost exactly the same as for NpCl3 (0.89). By contrast, the ratio of the kinetic to potential energy densities at the An-Cl bond critical points indicates that ionicity increases across the series, suggesting that the delocalisation index measures both orbital overlap and energy degeneracy-based covalency, while the bond critical point metric gauges only the former. Recalculation of all the data using the generalised gradient approximation PBE functional finds larger energy degeneracy-driven covalency in the later actinides than using hybrid DFT. Hence, we find that conclusions concerning the covalency of the An-Cl bond are dependent not only on the metric used to evaluate it, but also on the underlying electronic structure method.

Series of Tetravalent Actinide Amidinates: Structure Determination and Bonding Analysis

Two series of isostructural tetravalent actinide amidinates [AnX((S)-PEBA)₃] (An = Th, U, Np; X = Cl, N₃) bearing the chiral (S,S)-N,N'-bis(1-phenylethyl)benzamidinate ((S)-PEBA) ligand have been synthesized and thoroughly characterized in solid and in solution. This study expands the already reported tetravalent neptunium complexes to the lighter actinides thorium and uranium. Furthermore, a rare Ce(IV) amidinate [CeCl((S)-PEBA)₃] was synthesized to compare its properties to those of the analogous tetravalent actinide complexes. All compounds were characterized in the solid state using single-crystal XRD and infrared spectroscopy and in solution using NMR spectroscopy. Quantum chemical bonding analysis including also the isostructural Pa and Pu complexes was used to characterize the covalent contributions to any bond involving the metal cation. Th shows the least covalent character throughout the series, even substantially smaller than for the Ce complex. For U, Np, and Pu, similar covalent bonding contributions are found, but a natural population analysis reveals different origins. The 6d participation is the highest for U and decreases afterward, whereas the 5f participation increases continuously from Pa to Pu.

Correlating Electronic Structure and Magnetic Anisotropy in Actinide Complexes [An(COT)2], AnIII/IV = U, Np, and Pu

The electronic structures and magnetic anisotropies for compounds [An(COT)₂] (An = U^III/U^IV, Np^III/Np^IV and Pu^III/Pu^IV, COT = cyclooctatetraene) are characterized using scalar relativistic density functional theory calculations and second-order perturbation theory based on a complete active space self-consistent field reference including spin-orbit coupling. The degree of participation of 5f orbitals in actinide-ligand bonding and the associated metal-ligand covalency is found to trend as U > Np ≥ Pu for both the tetra-positive and tripositive An complexes. A spin-Hamiltonian analysis indicates only weak single-molecule magnet (SMM) characteristics for [U(COT)₂]^- and [Np(COT)₂] complexes and no significant SMM behavior for the other complexes. The weak SMM behavior in [U(COT)₂]^- and [Np(COT)₂] is attributed to a subtle interplay between local symmetry and ligand-field splitting. Such a result suggests that magnetic anisotropy in 5f³ ions can be modulated in general by electrostatic ligand field design. In particular, σ-donor ligands oriented 180 degrees relative to one another will have a maximal influence on the 5f-orbital ligand field splitting, while π donors like cyclopentadiene and COT generate ligand field influences that have more acute angles associated with corresponding atoms on the individual ligands. These observations rationalize the differences in SMM characteristics for [U(Bc^Me)₃] (Bc^Me- = dihydrobis(methylimidazolyl)borate) and [U(Bp^Me)₃] (Bp^Me- = dihydrobis(methylpyrazolyl)borate) and indicate strategies to design new actinide-based SMMs with high magnetic relaxation barriers.

Quantum chemical topology and natural bond orbital analysis of M–O covalency in M(OC6H5)4 (M = Ti, Zr, Hf, Ce, Th, Pa, U, Np)

V_XC(M,O): the exchange–correlation metric quantifies covalency between M and O atomic basins in M(OC₆H₅)₄ (M = Ti, Zr, Hf, Ce, Th, Pa, U, Np).

Systematic comparison of the structure of homoleptic tetradentate N2O2-type Schiff base complexes of tetravalent f-elements (M(IV) = Ce, Th, U, Np, and Pu) in solid state and in solution

A series of tetradentate N2O2-type Schiff base complexes with tetravalent 4f- and 5f-block metals, [M(salpn)2] (H2salpn = N,N'-disalicylidene-1,3-diaminopropane; M = Ce, Th, U, Np, and Pu), were prepared to systematically investigate their solid state structure, and their complexation behaviour in solution with the goal to investigate the subtle differences between 4f- and 5f-elements. X-ray diffraction revealed that all investigated metal cations form [M(salpn)2] complexes. All the complexes show the same ligand arrangement with meridional conformation, amongst which only Ce(iv) exhibits unique behaviour upon crystallisation. [Ce(salpn)2] crystallises in two less symmetric systems (P1[combining macron] or P21/n), whilst all the other [M(salpn)2] crystallise in a more symmetric orthorhombic system (Pban). Quantum chemical calculations suggest that the observed structural peculiarity of Ce(iv) stems from the geometrical flexibility due to the more "ionic" nature of bonds to the 4f element. 1H NMR measurements revealed that [M(salpn)2] forms two different species in solution with and without an additional solvent molecule, where the relative distribution of the two species depends mainly on the ionic radius of the metal centre. Again, Ce(iv) behaves differently from the tetravalent actinides with a higher ratio of the solvent-molecule-coordinated species than the ratio expected from its ionic radius. Hence, this study is successful in observing subtle differences between 4f- (i.e. Ce) and 5f-elements (actinides; Th, U, Np, and Pu) both in the solid state and in solution on an analytically distinguishable level, and in relating the observed subtle differences to their electronic structure.

Multiple Bonding in Lanthanides and Actinides: Direct Comparison of Covalency in Thorium(IV)- and Cerium(IV)-Imido Complexes

A series of thorium(IV)-imido complexes was synthesized and characterized. Extensive experimental and computational comparisons with the isostructural cerium(IV)-imido complexes revealed a notably more covalent bonding arrangement for the Ce═N bond compared with the more ionic Th═N bond. The thorium-imido moieties were observed to be 3 orders of magnitude more basic than their cerium congeners. More generally, these results provide unique experimental evidence for the larger covalent character of 4f⁰5d⁰ Ce(IV) multiple bonds compared to its 5f⁰6d⁰ Th(IV) actinide congener.

Reactions of cerium complexes with transition metal nitrides: synthesis and structure of heterometallic cerium complexes containing bridging catecholate ligands

Heterometallic cerium complexes containing bridging catecholate ligands have been synthesized from cerium complexes with Kläui's tripodal ligand and metal catecholates.