Theoretical Study of Actinide(III)-DOTA Complexes

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.

Metal-Ligand Interactions in Scandium Complexes with Radiopharmaceutical Applications

The radioisotopes of scandium (43Sc, 44Sc, and 47Sc) are potential candidates for use in imaging and therapy both separately and as elementally matched pairs for radiotheranostics. In the present study the bonding interactions of Sc3+ with 18 hepta- to decadentate ligands are compared using density functional theory (DFT) calculations. The bonding analysis is based on the natural bond orbital (NBO) model. The main contributions to the bonding were assessed using natural energy decomposition analysis (NEDA). Most of the ligands have anionic character (charges from 2- to 8-); thus the electrical term determines the major differences in the interaction energies. However, interesting features were found in the covalent contributions manifested by the ligand → Sc3+ charge transfer (CT) interactions. Significant differences could be observed in the energetic contributions of the N and O donors to the total CT.

Accurate Evaluation of Dispersion Energies at Coupled Cluster Level to Understand the Substituent Effects in Am(III) and Eu(III) Complexes

The effect of cyclic and aromatic substituents on the complexation behavior of phosphine oxide ligands with Am(III) and Eu(III) was investigated at density functional theory (DFT) and domain-based local pair natural orbital coupled-cluster (DLPNO-CC) levels. Combining DFT with accurate coupled cluster methods, we have evaluated the dispersion energy contributions to the complexation energies for trivalent Am and Eu complexes for the first time. Irrespective of the nature of substituents on the P atom, the electronic structure of the P═O group remains identical in all of the ligands. The study reveals the importance of dispersion interactions during complexation and is estimated to be more significant for Am(III) than for Eu(III) complexes.

Understanding the Coordination Chemistry of AmIII/CmIII in the DOTA Cavity: Insights from Energetics and Electronic Structure Theory

The minor actinides Am/Cm show multiple possibilities for coordination, providing great opportunities for their extraction and adsorption separation. Herein, we report complexation in an aqueous medium of Am^III/Cm^III in the DOTA (H₄DOTA = 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid) cavity with axial ligands (OH^-, F^-, and H₂O), based on the energetics and electronic structure properties using density functional theory (DFT). The formation and substitution reactions of OH^--capped complexes are more likely to occur due to their enhanced hydration Gibbs free energies, followed by F^-, and then H₂O. Both the longer An-O_DOTA bond lengths and the larger bite angle (∠O-An-O) in the OH^--capped complexes reflect the enhanced coordination provided by the axial ligand, slightly less so for F^-. Energy decomposition analysis based on the electronic structure supports the preference for OH^--capped complexes with a near-perfect balance between attractive and repulsive contributions toward the interaction. Furthermore, molecular orbital analysis revealed that the frontier molecular orbitals of Am and Cm complexes are substantially different; that is, the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) compositions of the Am complexes are all contributed by 5f, while the HOMO and LUMO compositions of the Cm complexes are derived from 5f and 6d, respectively. Finally, the metal-exchange reactions demonstrate competitive complexation of DOTA toward Am^III over Cm^III for the OH^--capped system. These results imply the importance of coordination chemistry in actinide chemistry in general and specifically in Am^III/Cm^III solution chemistry.

Actinium coordination chemistry: A density functional theory study with monodentate and bidentate ligands

In the current study, the coordination chemistry of nine-coordinate Ac(III) complexes with 35 monodentate and bidentate ligands was investigated using density functional theory (DFT) in terms of their geometries, charges, reaction energies, and bonding interactions. The energy decomposition analysis with naturals orbitals for chemical valence (EDA-NOCV) and the quantum theory of atoms in molecules (QTAIM) were employed as analysis methods. Trivalent Ac exhibits the highest affinities toward hard acids (such as charged oxophilic donors, fluoride), so its classification as a hard acid is justified. Natural population analysis quantified the involvement of 5f orbitals on Ac to be about 30% of total valence electron natural configuration indicating that Ac is a member of the actinide series. Pearson correlation coefficients were used to study the pairwise correlations among the bond lengths, ΔG reaction energies, charges on Ac and donor atoms, and data from EDA-NOCV and QTAIM. Strong correlations and anticorrelations were found between Voronoi charges on donor atoms with ΔG, EDA-NOCV interaction energies and QTAIM bond critical point densities.

Theoretical Study of Complexes of Tetravalent Actinides with DOTA

1,4,7,10-Tetraazacyclododecane-N,N′,N″,N‴-tetraacetic acid (H4DOTA) is a prominent chelating ligand with potential applications in various fields, from radiotherapy to the separation of fission products. The present study explores the stability, structure, and bonding properties of its complexes with tetravalent actinides (An = Th, U, Np, Pu) using density functional theory and relativistic multireference calculations. Neutral complexes prefer to form symmetric (C4) structures with DOTA. The first coordination sphere of the actinide ions is readily saturated by a weakly bonded H2O ligand. The latter ligand reduces the molecular symmetry while exerting only marginal effects on the properties of the parent complex. An-ligand bonding is mainly electrostatic, but there are also significant charge-transfer contributions from DOTA to the An 6d/5f orbitals. The charge-transfer interactions and the covalent character of bonding increase gradually in the order of Th < U < Np < Pu, as indicated by analysis of the electron density distribution using the Quantum Theory of Atoms in Molecules.

Metal-ligand bonding in bispidine chelate complexes for radiopharmaceutical applications

The complexes of selected radionuclides relevant for nuclear medicine (InIII, BiIII, LuIII, AcIII and in addition LaIII for comparative purposes) with the octadentate (6,6′-((9-hydroxy-1,5-bis(methoxycarbonyl)-2,4-di(pyridin-2-yl)-3,7-diazabicyclo[3.3.1]nonane-3,7-diyl)bis(methylene))dipicolinic acid) ligand, H2bispa2, have been studied by density functional theory calculations modelling both isolated and aqueous solution conditions. The properties in focus are the encapsulation efficiency of the ligand for the different-size metals (M), the differences in bonding with the various MIII ions analysed using Bader’s atoms in molecules theory and the possibility and characteristics of nona- and decacoordination by H2O ligands. The computed results confirmed strong steric effects in the case of the In complex excluding higher than octacoordination. The studied properties depend strongly on the interplay of the sizes and electronic structures of the MIII ions. The computations support high stability of the complexes in aqueous solution, where also the solvation energies of the MIII ions (as dissociation products) play a significant role.

Chelating the Alpha Therapy Radionuclides 225Ac3+ and 213Bi3+ with 18-Membered Macrocyclic Ligands Macrodipa and Py-Macrodipa

The radionuclides 225Ac3+ and 213Bi3+ possess favorable physical properties for targeted alpha therapy (TAT), a therapeutic approach that leverages α radiation to treat cancers. A chelator that effectively binds and retains these radionuclides is required for this application. The development of ligands for this purpose, however, is challenging because the large ionic radii and charge-diffuse nature of these metal ions give rise to weaker metal-ligand interactions. In this study, we evaluated two 18-membered macrocyclic chelators, macrodipa and py-macrodipa, for their ability to complex 225Ac3+ and 213Bi3+. Their coordination chemistry with Ac3+ was probed computationally and with Bi3+ experimentally via NMR spectroscopy and X-ray crystallography. Furthermore, radiolabeling studies were conducted, revealing the efficient incorporation of both 225Ac3+ and 213Bi3+ by py-macrodipa that matches or surpasses the well-known chelators macropa and DOTA. Incubation in human serum at 37 °C showed that ∼90% of the 225Ac3+-py-macrodipa complex dissociates after 1 d. The Bi3+-py-macrodipa complex possesses remarkable kinetic inertness reflected by an EDTA transchelation challenge study, surpassing that of Bi3+-macropa. This work establishes py-macrodipa as a valuable candidate for 213Bi3+ TAT, providing further motivation for its implementation within new radiopharmaceutical agents.

Covalency of Trivalent Actinide Ions with Different Donor Ligands: Do Density Functional and Multiconfigurational Wavefunction Calculations Corroborate the Observed "Breaks"?

A comprehensive ab initio study of periodic actinide-ligand bonding trends for trivalent actinides is performed. Relativistic density functional theory (DFT) and complete active-space (CAS) self-consistent field wavefunction calculations are used to dissect the chemical bonding in the [AnCl₆]^3-, [An(CN)₆]^3-, [An(NCS)₆]^3-, [An(S₂PMe₂)₃], [An(DPA)₃]^3-, and [An(HOPO)]^- series of actinide (An = U-Es) complexes. Except for some differences for the early actinide complexes with DPA, bond orders and excess 5f-shell populations from donation bonding show qualitatively similar trends in 5f ⁿ active-space CAS vs DFT calculations. The influence of spin-orbit coupling on donation bonding is small for the tested systems. Along the actinide series, chemically soft vs chemically harder ligands exhibit clear differences in bonding trends. There are pronounced changes in the 5f populations when moving from Pu to Am or Cm, which correlate with previously noted "breaks" in chemical trends. Bonding involving 5f becomes very weak beyond Cm/Bk. We propose that Cm(III) is a borderline case among the trivalent actinides that can be meaningfully considered to be involved in ground-state 5f covalent bonding.

Metal–ligand interactions in complexes of cyclen-based ligands with Bi and Ac

The structural and bonding properties of Bi and Ac complexes with cyclen-based chelating ligands have been studied using relativistic DFT calculations in conjunction with TZ2P all-electron basis sets. Besides the parent cyclen ligand, the study has covered its extensions with pyridine-type (Lpy), carboxylate (DOTA, DOTPA), picolinate (MeDO2PA) and phosphonate (DOTMP) pendant arms. The effect of the cyclen ring size has been probed by increasing it from [12]aneN4 to [16]aneN4. Additional extensions in the DOTA complexes included the H2O ligand at the 9th coordination site as well as the p-SCN-Bn substituent (a popular linker to the targeting vector). The study focuses on the complex stability, the nature of bonding and the differences between Ac and Bi in the complexes. The metal–ligand interactions have been analysed by the Extended Transition State method combined with Natural Orbitals of Chemical Valence theory and Quantum Theory of Atoms in Molecules models.