Protactinium and the intersection of actinide and transition metal chemistry

Extension of the D3 and D4 London dispersion corrections to the full actinides series

Efficient dispersion corrections are an indispensable component of modern density functional theory, semi-empirical quantum mechanical, and even force field methods. In this work, we extend the well established D3 and D4 London dispersion corrections to the full actinides series, francium, and radium. To keep consistency with the existing versions, the original parameterization strategy of the D4 model was only slightly modified. This includes improved reference Hirshfeld atomic partial charges at the ωB97M-V/ma-def-TZVP level to fit the required electronegativity equilibration charge (EEQ) model. In this context, we developed a new actinide data set called AcQM, which covers the most common molecular actinide compound space. Furthermore, the efficient calculation of dynamic polarizabilities that are needed to construct CAB6 dispersion coefficients was implemented into the ORCA program package. The extended models are assessed for the computation of dissociation curves of actinide atoms and ions, geometry optimizations of crystal structure cutouts, gas-phase structures of small uranium compounds, and an example extracted from a small actinide complex protein assembly. We found that the novel parameterizations perform on par with the computationally more demanding density-dependent VV10 dispersion correction. With the presented extension, the excellent cost-accuracy ratio of the D3 and D4 models can now be utilized in various fields of computational actinide chemistry and, e.g., in efficient composite DFT methods such as r2SCAN-3c. They are implemented in our freely available standalone codes (dftd4, s-dftd3) and the D4 version will be also available in the upcoming ORCA 6.0 program package.

Proton-Coupled Electron Transfer at the Pu5+/4+ Couple

The synthesis and solution and solid-state characterization of [Pu4+(NPC)4], 1-Pu, (NPC = [NPtBu(pyrr)2]-; tBu = C(CH3)3; pyrr = pyrrolidinyl) and [Pu3+(NPC)4][K(2.2.2.-cryptand)], 2-Pu, is described. Cyclic voltammetry studies of 1-Pu reveal a quasi-reversible Pu4+/3+ couple, an irreversible Pu5+/4+ couple, and a third couple evincing a rapid proton-coupled electron transfer (PCET) reaction occurring after the electrochemical formation of Pu5+. The chemical identity of the product of the PCET reaction was confirmed by independent chemical synthesis to be [Pu4+(NPC)3(HNPC)][B(ArF5)4], 3-Pu, (B(ArF5)4 = tetrakis(2,3,4,5,6-pentafluourophenyl)borate) via two mechanistically distinct transformations of 1-Pu: protonation and oxidation. The kinetics and thermodynamics of this PCET reaction are determined via electrochemical analysis, simulation, and density functional theory. The computational studies demonstrate a direct correlation between the changing nature of 5f and 6d orbital participation in metal-ligand bonding and the electron density on the Nim atom with the thermodynamics of the PCET reaction from Np to Pu, and an indirect correlation with the roughly 5-orders of magnitude faster Pu PCET compared to Np for the An5+ species.

Chronicles of plutonium peroxides: spectroscopic characterization of a new peroxo compound of Pu(IV)

Although hydrogen peroxide (H2O2) has been highly used in nuclear chemistry for more than 75 years, the preparation and literature description of tetravalent actinide peroxides remain surprisingly scarce. A new insight is given in this topic through the synthesis and thorough structural characterization of a new peroxo compound of Pu(IV).

Chemdex: quantification and distributions of valence numbers, oxidation numbers, coordination numbers, electron numbers, and covalent bond classes for the elements

The paper introduces Chemdex, a freely accessible web-based database of over 70 000 compounds characterised by crystallography from across the periodic table. Its software calculates for an atom of interest within each compound classifications including valence number, oxidation number, coordination number, electron number, several covalent bond classifications, and the attached atom set. Users may explore distributions of these classifications by percentages and heat map displays for individual elements or sets of elements, or in several cases for one classification plotted against a second. These properties often display clear periodicity. Based upon distributions across the periodic table of valence numbers, electron numbers, coordination numbers, and attached atom data suggestions are made regarding the placement of hydrogen in the periodic table, membership of group 3 in the periodic table, locations of the early actinoids in the periodic table, and assignments of certain elements as metalloids.

Periodic trends in trivalent actinide halides, phosphates, and arsenates

Due to the limited abundance of the actinide elements, computational methods, for now, remain an exclusive avenue to investigate the periodic trends across the actinide series. As every actinide element can exhibit a +3-oxidation state, we have explored model systems of gas-phase actinide trihalides, phosphates, and arsenates across the series to capture the periodic trends. By doing so, we were able to capture the periodic trends down the halogen series as well, and for the first time we are reporting a study on actinide astatides. Using scalar and spin-orbit relativistic Density Functional Theory (DFT) calculations, we have explored the variations in bond lengths, bond angles, and the charges on actinides (An). Despite the use of different sets of ligands, the trends remain similar. The properties of trivalent Pa, U, Np, and Pu are nearly identical; similar ionic radii could be the reason. The actinide elements show a tendency to exhibit a pre-Pu and a post-Cm behaviour, with Am acting as a switch. This could be due to the change in the behaviour from d-f-type to f-filling/d-type at around Pu-Cm in the actinides as already proposed in the previous literature. Bond lengths in the AnX3 increase down the halide series, and the atomic charges decrease on the actinide elements.

Planar Four-Membered Diboron Actinide Compound with Double Möbius Aromaticity

The Möbius rule predicts that a planar four-membered metallacycle can be aromatic with four mobile electrons, but such a simple ring has escaped recognition because it usually favors Hückel anti-aromaticity. Here, we report that a quasi-square four-membered actinide compound (Pa₂B₂) is doubly Möbius aromatic. Chemical bonding analyses reveal that this diboron protactinium molecule has four delocalized π electrons in addition to four delocalized σ electrons, satisfying the 4n Möbius rule for both σ and π components. Energetically, the block-localized wavefunction method, which is the simplest variant of ab initio valence bond theory, shows that the delocalization energy for the π and σ electrons reaches up to 65.0 and 72.3 kcal/mol, respectively, while the extra cyclic resonance energy (ECRE) amounts to 45 kcal/mol. The large positive ECRE values strongly confirm the unprecedented double Möbius aromaticity in Pa₂B₂. We anticipate that this new type of aromatic molecule can enrich the concept of Möbius aromaticity and open a new avenue for actinide compounds.

Prediction of the structures and heats of formation of MO2, MO3, and M2O5 for M = V, Nb, Ta, Pa

Structures for the mono-, di-, and tri-bridge isomers of M₂O₅ as well as those for the MO₂ and MO₃ fragments for M = V, Nb, Ta, and Pa were optimized at the density functional theory (DFT) level. Single point CCSD(T) calculations extrapolated to the complete basis set (CBS) limit at the DFT geometries were used to predict the energetics. The lowest energy dimer isomer was the di-bridge for M = V and Nb and the tri-bridge for M = Ta and Pa. The di-bridge isomers were predicted to be composed of MO₂⁺ and MO₃^- fragments, whereas the mono- and tri-bridge are two MO₂⁺ fragments linked by an O^2-. The heats of formation of M₂O₅ dimers, as well as MO₂ and MO₃ neutral and ionic species were predicted using the Feller-Peterson-Dixon (FPD) approach. The heats of formation of the MF₅ species were calculated to provide additional benchmarks. Dimerization energies to form the M₂O₅ dimers are predicted to become more negative going down group 5 and range from -29 to -45 kcal mol^-1. The ionization energies (IEs) for VO₂ and TaO₂ are essentially the same at 8.75 eV whereas the IEs for NbO₂ and PaO₂ are 8.10 and 6.25 eV, respectively. The predicted adiabatic electron affinities (AEAs) range from 3.75 eV to 4.45 eV for the MO₃ species and vertical detachment energies from 4.21 to 4.59 eV for MO₃^-. The calculated MO bond dissociation energies increase from 143 kcal mol^-1 for M = V to ∼170 kcal mol^-1 for M = Nb and Ta to ∼200 kcal mol^-1 for M = Pa. The M-O bond dissociation energies are all similar ranging from 97 to 107 kcal mol^-1. Natural bond analysis provided insights into the types of chemical bonds in terms of their ionic character. Pa₂O₅ is predicted to behave like an actinyl species dominated by the interactions of approximately linear PaO₂⁺ groups.

Protactinium and Actinium Monohydrides: A Theoretical Study on Their Spectroscopic and Thermodynamic Properties

Spectroscopic and thermodynamics properties including bond dissociation energies (BDEs), adiabatic electron affinities (AEAs), and ionization energies (IEs) have been predicted for AcH and PaH using the Feller-Peterson-Dixon composite approach. Comparisons with previous studies on ThH and UH were performed to identify possible trends in the actinide series. Multireference CASPT2 calculations were used to predict the spin-orbit effects and obtain potential energy curves for the low-lying Ω states around the equilibrium distance as well as the vertical detachment energies (VDEs) from AcH^- and PaH^- to excited states of the neutral species. The calculated AEA for AnH (An = Ac, Th, Pa, U) showed that the AEA increases from AcH (0.425 eV) to ThH (0.820 eV) and decreases to PaH (0.781 eV) and to UH (0.457 eV), whereas the IE values are 5.887 eV (AcH), 6.181 eV (ThH), 6.204 eV (PaH), and 6.182 eV (UH). The ground state of AcH, AcH^-, PaH, and PaH^- are predicted to be¹Σ⁺₀,²Π_3/2, ³H₄, and ⁴I_9/2, respectively. The BDEs for AcH and PaH are 276.4 and 237.2 kJ/mol, and those for AcH^- and PaH^- are 242.8 and 239.8 kJ/mol, respectively. The natural bond analysis shows a significant ionic character, An⁺H^-, in the bonding of the neutral hydrides.

On the volatility of protactinium in chlorinating and brominating gas media

A multi-target recoil chamber technique was applied to study online chemical properties of protactinium in chlorinating and brominating gas media using 226Pa (T 1/2 = 1.8 min) decaying by alpha emission (74%) and β+/EC decay (26%). A 58 MeV proton beam passing 15 × 50 μg/cm2 thick 232Th targets enabled production of 226Pa formed in the reaction 232Th(p,7n). Isothermal gas chromatography in quartz columns allowed for the determination of adsorption enthalpies of oxohalides and pure halides of Pa5+ compounds. On the basis of empirical correlations, these adsorption enthalpies (ΔH0 ads) could be converted to sublimation enthalpies (ΔH0 subl). Resulting values for the assumed compounds PaCl5, PaOCl3, PaBr5, and PaOBr3 were 113 ± 15, 329 ± 16, 165 ± 5 and 235 ± 17 kJ/mol, respectively. These values are rather similar to known ΔH0 subl data for group-5 elements Nb, Ta and Db in support of the assumption that Pa is a pseudo-group 5 element.

Understanding Periodic and Non-periodic Chemistry in Periodic Tables

The chemical elements are the "conserved principles" or "kernels" of chemistry that are retained when substances are altered. Comprehensive overviews of the chemistry of the elements and their compounds are needed in chemical science. To this end, a graphical display of the chemical properties of the elements, in the form of a Periodic Table, is the helpful tool. Such tables have been designed with the aim of either classifying real chemical substances or emphasizing formal and aesthetic concepts. Simplified, artistic, or economic tables are relevant to educational and cultural fields, while practicing chemists profit more from "chemical tables of chemical elements." Such tables should incorporate four aspects: (i) typical valence electron configurations of bonded atoms in chemical compounds (instead of the common but chemically atypical ground states of free atoms in physical vacuum); (ii) at least three basic chemical properties (valence number, size, and energy of the valence shells), their joint variation across the elements showing principal and secondary periodicity; (iii) elements in which the (sp)8, (d)10, and (f)14 valence shells become closed and inert under ambient chemical conditions, thereby determining the "fix-points" of chemical periodicity; (iv) peculiar elements at the top and at the bottom of the Periodic Table. While it is essential that Periodic Tables display important trends in element chemistry we need to keep our eyes open for unexpected chemical behavior in ambient, near ambient, or unusual conditions. The combination of experimental data and theoretical insight supports a more nuanced understanding of complex periodic trends and non-periodic phenomena.

Bridging the Transuranics with Uranium(IV) Sulfate Aqueous Species and Solid Phases

Isolating isomorphic compounds of tetravalent actinides (i.e., Th^IV, U^IV, Np^IV, and Pu^IV) improve our understanding of the bonding behavior across the series, in addition to their relationship with tetravalent transition metals (Zr and Hf) and lanthanides (Ce). Similarities between these tetravalent metals are particularly illuminated in their hydrolysis and condensation behavior in aqueous systems, leading to polynuclear clusters typified by the hexamer [M^IV₆O₄(OH)₄]¹²⁺ building block. Prior studies have shown the predominance and coexistence of smaller species for Th^IV (monomers, dimers, and hexamers) and larger species for U^IV, Np^IV, and Pu^IV (including 38-mers and 70-mers). We show here that aqueous uranium(IV) sulfate also displays behavior similar to that of Th^IV (and Zr^IV) in its isolated solid-phase and solution speciation. Two single-crystal X-ray structures are described: a dihydroxide-bridged dimer (U₂) formulated as U₂(OH)₂(SO₄)₃(H₂O)₄ and a monomer-linked hexamer framework (U-U₆) as (U(H₂O)_3.5)₂U₆O₄(OH)₄(SO₄)₁₀(H₂O)₉. These structures are similar to those previously described for Th^IV. Moreover, cocrystallization of monomer and dimer and of dimer and monomer-hexamer phases for both Th^IV (prior) and U^IV (current) indicates the coexistence of these species in solution. Because it was not possible to effectively study the sulfate-rich solutions via X-ray scattering from which U₂ and U-U₆ crystallized, we provide a parallel solution speciation study in low sulfate conditions, as a function of the pH. Raman spectroscopy, UV-vis spectroscopy, and small-angle X-ray scattering of these show decreasing sulfate binding, increased hydrolysis, increased species size, and increased complexity, with increasing pH. This study describes a bridge across the first half the actinide series, highlighting U^IV similarities to Th^IV, in addition to the previously known similarities to the transuranic elements.

Relativistic Effects Stabilize the Planar Wheel-like Structure of Actinide-Doped Gold Clusters: An@Au7 (An = Th to Cm)

Despite the chemistry of actinide-ligand bonding is continuing and of burgeoning interest, investigations of the chemical bonding of bimetallic complexes involving transuranics remain relatively less, and there are rarely studies on the bonding features between actinide and coinage metals (CM). We present a systematic research on the series of An@Au₇ (An = Th to Cm), UCM₇ (CM = Cu, Ag, Au), and WAu₇ clusters to investigate the unique geometries, electronic structures, and chemical bonding between An 5f6d orbitals and CM ns orbitals, and to find their periodicity across the actinides and within the group of transition metals. A unique planar wheel-like structure for An@Au₇ clusters with the help of actinide metals encapsulation via spin-orbit coupling, resulting in An(III). Instead, the transition-metal (TM) element W retains its usual six-gold-coordination structure in WAu₇, thus forcing the seventh Au out of plane. The An-CM interactions, depending on the ion radii, become stronger with the increase of the atomic number of the actinide metals, as well as the CM. These results show that the presence of actinides in clusters can lead to unique electronic and geometrical structures.

Density functional theory investigations on the coordination of Pa(v) with N,N-dialkylamide

The density functional theory (DFT) method was used to study the coordination of a series of N,N-dialkylamides with Pa(v) to shed light on the inherent principles for screening amide extractants of Pa(v) from aqueous solution.

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

Carbon monoxide activation by atomic thorium: ground and excited state reaction pathways

Multi-reference configuration interaction (MRCI) and single reference coupled cluster calculations are performed for the ThCO and OThC isomers. Scalar and spin-orbit relativistic effects are considered through a relativistic pseudopotential and the coupling of MRCI wavefunctions via the Breit-Pauli spin-orbit Hamiltonian. Optimized geometries, excitation energies, and vibrational frequencies are reported for both isomers. Full potential energy profiles are constructed for the Th+CO reaction and the conversion of the produced ThCO to OThC. Linear ThCO was found to be more stable than the highly ionic bent OThC system by about 4 kcal mol^-1. The interconversion barrier is estimated to be around 30 kcal mol^-1. Our results are in agreement with earlier experimental data for the two isomers. The lowest lying states of Th do not populate f-orbitals and resemble the electronic structure of Ti. Therefore, the ability of the two atoms to activate the C[triple bond, length as m-dash]O bond is compared. OTiC is found to be about 40 kcal mol^-1 less stable than TiCO revealing the efficiency of Th and possibly other f-block elements to activate multiple chemical bonds as opposed to d-block metals.

Protactinium(V) in aqueous solution: a light actinide without actinyl moiety

This review highlights recent data on the complexation of Pa(V) with inorganic (fluoride and sulphate) and organic (oxalate, nitrilotriacetate, diethylenetriaminepentaacetate) ligands in solution. New thermodynamic parameters relative to the complexation of Pa(V) with sulphate are presented. The review also includes gas phase and theoretical studies focused on the interaction of Pa(V) in the dioxo and oxo forms with water.

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.

Theoretical prediction of some layered Pa2O5 phases: structure and properties

Density functional theory (DFT) was used to predict and study protactinium pentoxide (Pa₂O₅), which presents a fluorite and layered protactinium oxide-type structure. Although the layered structure has been observed with the isostructural transition Nb and Ta metal pentoxides experimentally, the detailed structure and properties of the layered Pa₂O₅ are not clear and understandable. Our theoretical prediction explored some possible stable structures of the Pa₂O₅ stoichiometry according to the existing M₂O₅ structures (where M is an actinide Np or transition Nb, Ta, and V metal) and replacing the M ions with protactinium ions. The structural, mechanical, thermodynamic and electronic properties including lattice parameters, bulk moduli, elastic constants, entropy and band gaps were predicted for all the simulated structures. Pa₂O₅ in the β-V₂O₅ structure was found to be a competitive structure in terms of stability, whereas Pa₂O₅ in the ζ-Nb₂O₅ structure was found to be the most stable overall. This is consistent with Sellers's experimental observations. In particular, Pa₂O₅ in the ζ-Nb₂O₅ structure is predicted to be charge-transfer insulators. Furthermore, we predict that ζ-Nb₂O₅-structured Pa₂O₅ is the most thermodynamically stable under ambient conditions and pressure.

Density functional theory (DFT) was used to predict and study protactinium pentoxide (Pa₂O₅), which presents a fluorite and layered protactinium oxide-type structure.

Theoretical studies on the oxidation states and electronic structures of actinide-borides: AnB12 (An = Th-Cm) clusters

As B12 clusters exhibit significant structural stability due to double aromaticity, metal doped-B12 clusters often prefer a half sandwich structure. Herein, we report a systematic theoretical study on the geometric and electronic structures, and chemical bonding of the half sandwich AnB12 (An = Th to Cm) clusters to explore the stability and extent of covalency of the An-B bonds of these actinide borides. We have shown that in the gas-phase clusters, the significant stability of AnB12 is determined by electrostatic and orbital interactions between the An 5f6d7s orbitals and π-type molecular orbitals from B 2p orbitals of the B12 unit. A change-over of An-B bond length from An = Th to Cm is found at An = Pa as a result of actinide contraction combined with weakening An-B bonding due to an energy decrease and orbital localization of the 5f orbitals. Consistently, the oxidation states of the An atoms at first increase from Th(f0)IV to Pa(f0)V, and then due to the 5f-AO contraction, they smoothly decline to U(f2)IV, Np(f4)III and Pu(f5)III, and then eventually to Am(f7)II but Cm(f7)III, both with a half-filled 5f shell.