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

Periodic Bonding Behaviors of Actinide Atoms (An = Th-Cm) with Negatively Curved Nanographene

The nanographene with negative curvature has been extensively studied due to its interesting properties and potential applications. In the present work, we have performed all-electron scalar relativistic density functional theory (DFT) calculations to understand the periodic interaction mechanisms of actinide atoms (An = Th, Cm) with the TB8C nanographene. The encapsulated complexes (An@TB8C) were formed due to the octagonal vacancy in the TB8C nanographene. TB8C shows fairly high affinity toward An atoms, especially for Th and Pa. AIMD simulations further confirmed the effective trapping of An atom with TB8C. The partial covalent characters of An-C bonds in An@TB8C were revealed through various bond analysis methods. The 6d electrons of An play an important role in the participation of chemical bonds. The delocalization index (DI) is proposed as a useful descriptor in the study of bond strength involving the actinides. Electronic absorption spectra were simulated for further identification in the experiments. The current work has expanded the potential molecular properties and applications of nanographene.

High-valent actinyl (AnO2; An = U, Np and Pu) complexation with TEtraQuinoline (TEQ) ligand - a DFT study

A recently synthesised novel macrocycle composed of four quinoline units called TEtraquinoline (TEQ) (J. Am. Chem. Soc., 2023, 145(4), 2609-2618) was reported to exhibit transition metal complexation ability. Meanwhile, there is a growing interest in different binding motifs for radioactive and toxic actinides. In this study, we modelled high-valent actinyl (AnO2)n+, An = U, Np, Pu; n = 1, 2, 3 complexes of TEQ and studied their geometric and electronic properties using scalar relativistic density functional theory (SR-DFT). The calculated results showed that the equatorial An-N and axial AnO bonds were polar bonds with a high degree of covalence, the former being longer than the latter. Natural bond orbital (NBO) analysis of the An-N bond orders decreased from complexes of uranyl to plutonyl and from complexes of heptavalent to pentavalent actinyls. This was due to the localization of the 5f orbital in the heavier actinyl and the high charge on An. The charge analysis showed a ligand-to-metal charge transfer (LMCT) on complexation. It was interesting to observe metal-to-ligand spin delocalization only in the [AnVI/VIIO2L]2+/3+ complexes, where the spin density on An was observed to increase on complexation. Based on the assigned oxidation states, it was observed that the heptavalent neptunyl cation retained its formal high oxidation state on complexation with TEQ. The energetics associated with the formation reaction of all the actinyl-TEQ complexes suggest spontaneity at lower temperatures (i.e., lower than 298.15 K). The energy decomposition analysis (EDA) indicates that the electrostatic energy contributions were predominant in the [AnVO2L]+ complexes, while covalent (orbital) energy contributions were higher in the [AnVIO2L]2+ and [AnVIIO2L]3+ complexes. The extended transition state-natural orbitals for chemical valence (ETS-NOCV) analysis confirmed the prominent covalent character in [AnVIIO2L]3+ over [AnVIO2L]2+ and [AnVO2L]+and the back donation of charges from An to N that stabilizes TEQ.

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

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

High Coordination Numbers of Actinides (An) in AnC13+ Rings (An = Th and U)

An intriguingly high abundance of both ThC13+ and UC13+ cluster cations was observed in a previous mass spectrometry experiment; however, the structural identification of these cations has not yet been completed. In this study, we determined the lowest lying structures of ThC13+ and UC13+ clusters using an unbiased structural search method. The 13-coordinate planar ring configuration was the most stable for both ThC13+ and UC13+ cluster cations. The C-An bonds in ThC13+ and UC13+ show a small degree of covalency, originating from the overlap of the s, d, and f orbitals of the An atoms with C 2p orbitals of both π and σ characteristics. The infrared and electronic absorption spectra of the most favorable planar ring configurations were theoretically simulated to facilitate the identification of the molecular structures in future experiments. This study provides an in-depth understanding of the experimental mass spectra.

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

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

Trends in the Electronic Structure and Chemical Bonding of a Series of Porphyrinoid-Uranyl Complexes

In this paper, we have explored the relativistic density functional theory study on a series of deprotonated porphyrinoid (Lⁿ) complexes of uranyl to investigate the geometrical structures and chemical bonding. The ligands bound with uranyl in the 1:1 complexes [UO₂(Lⁿ)]^x (n = 4, 5, 6; x = 0, -1, -2), showing more thermodynamic stability for "in-cavity" structures of L⁵ and L⁶ than that of the "side-on" structure of L⁴ and an increase in stability with the increase of negative charges, L^2- < L^3- < L^4-. Among the six ligands, the cyclo[6]pyrrole presents the best selectivity toward uranyl. Based on chemical bonding analyses, the U-N_L bond in the in-cavity complexes adopts a typical dative N_L → U bond with mainly ionic bonding and significant covalency, which comes from the significant orbital interaction of U 5f_ϕ6d_δ7s hybrid AOs and N_L 2p-based MOs. This work provides a systematic understanding of the coordination chemistry in uranyl pyrrole-containing macrocycle complexes and the nature of chemical bonding in such systems, which may provide inspirations for the future design of synthetic targets that could be relevant to actinide separations or in the remediation of spent nuclear fuel.

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.

Effects of coordinating heteroatoms on molecular structure, thermodynamic stability and redox behavior of uranyl(vi) complexes with pentadentate Schiff-base ligands

Uranyl(vi) complexes with pentadentate N₃O₂-, N₂O₃- and N₂O₂S₁-donating Schiff base ligands, tBu,MeO–saldien–X²⁻ (X = NH, O and S), were synthesized and thoroughly characterized by ¹H NMR, IR, elemental analysis, and single crystal X-ray diffraction. The crystal structures of UO₂(tBu,MeO–saldien–X) showed that the U–X bond strength follows U–O ≈ U–NH > U–S. Conditional stability constants (β_X) of UO₂(tBu,MeO–saldien–X) in ethanol were investigated to understand the effect of X on thermodynamic stability. The log β_X decrease in the order of UO₂(tBu,MeO–saldien–NH) (log β_NH = 10) > UO₂(tBu,MeO–saldien–O) (log β_O = 7.24) > UO₂(tBu,MeO–saldien–S) (log β_S = 5.2). This trend cannot be explained only by Pearson's Hard and Soft Acids and Bases (HSAB) principle, but rather follows the order of basicity of X. Theoretical calculations of UO₂(tBu,MeO–saldien–X) suggested that the ionic character of U–X bonds decreases in the order of U–NH > U–O > U–S, while the covalency increases in the order U–O < U–NH < U–S. Redox potentials of all UO₂(tBu,MeO–saldien–X) in DMSO were similar to each other regardless of the difference in X. Spectroelectrochemical measurements and DFT calculations revealed that the center U⁶⁺ of each UO₂(tBu,MeO–saldien–X) undergoes one-electron reduction to afford the corresponding uranyl(v) complex. Consequently, the difference in X of UO₂(tBu,MeO–saldien–X) affects the coordination of tBu,MeO–saldien–X²⁻ with UO₂²⁺. However, the HSAB principle is not always prominent, but the Lewis basicity and balance between ionic and covalent characters of the U–X interactions are more relevant to determine the bond strengths.

The U–X bond strength and thermodynamic stability of uranyl(vi) complexes with pentadentate N₂O₂X₁-donating ligands (X = NH, O, S) are affected by the difference in X. In contrast, the X atom does not largely affect the redox behavior of the complexes.