Compression of curium pyrrolidine-dithiocarbamate enhances covalency.

Curium is unique in the actinide series because its half-filled 5f 7 shell has lower energy than other 5f n configurations, rendering it both redox-inactive and resistant to forming chemical bonds that engage the 5f shell1-3. This is even more pronounced in gadolinium, curium's lanthanide analogue, owing to the contraction of the 4f orbitals with respect to the 5f orbitals4. However, at high pressures metallic curium undergoes a transition from localized to itinerant 5f electrons5. This transition is accompanied by a crystal structure dictated by the magnetic interactions between curium atoms5,6. Therefore, the question arises of whether the frontier metal orbitals in curium(III)-ligand interactions can also be modified by applying pressure, and thus be induced to form metal-ligand bonds with a degree of covalency. Here we report experimental and computational evidence for changes in the relative roles of the 5f/6d orbitals in curium-sulfur bonds in [Cm(pydtc)4]- (pydtc, pyrrolidinedithiocarbamate) at high pressures (up to 11 gigapascals). We compare these results to the spectra of [Nd(pydtc)4]- and of a Cm(III) mellitate that possesses only curium-oxygen bonds. Compared with the changes observed in the [Cm(pydtc)4]- spectra, we observe smaller changes in the f-f transitions in the [Nd(pydtc)4]- absorption spectrum and in the f-f emission spectrum of the Cm(III) mellitate upon pressurization, which are related to the smaller perturbation of the nature of their bonds. These results reveal that the metal orbital contributions to the curium-sulfur bonds are considerably enhanced at high pressures and that the 5f orbital involvement doubles between 0 and 11 gigapascal. Our work implies that covalency in actinides is complex even when dealing with the same ion, but it could guide the selection of ligands to study the effect of pressure on actinide compounds.

The Counterintuitive Relationship between Orbital Energy, Orbital Overlap, and Bond Covalency in CeF62- and CeCl62.

The 4f orbitals of Ce(IV) have shown appreciably enhanced covalent mixing with ligand orbitals relative to those of Ce(III). Here, X-ray spectroscopy, magnetic susceptibility measurements, and theoretical methods are used to investigate 4f covalency in CeF62- and CeCl62-. These techniques show covalent mixing between Ce 4f and F 2p orbitals to be about 25% less than mixing between Ce 4f and Cl 3p orbitals, placing CeF62- among the most ionic Ce(IV) compounds to-date. However, ligand field analysis using the experimental data shows significantly higher 4f orbital overlap with the F 2p orbitals compared to the Cl 3p. This result is counterintuitive since the Ce-F bonds display less 4f covalency despite their higher orbital overlap, and greater overlap is traditionally associated with enhanced bond covalency. The weaker covalency is attributed to the large energy gap between Ce 4f and F 2p orbitals strongly counteracting the higher orbital overlap. These results highlight that only a concerted consideration of both atomic orbital overlap and energy matching in f-element systems leads to an accurate picture of their bonding.

The resolution of the weak-exchange limit made rigorous, simple and general in binuclear complexes.

The correct interpretation of magnetic properties in the weak-exchange regime has remained a challenging task for several decades. In this regime, the effective exchange interaction between local spins is quite weak, of the same order of magnitude or smaller than the various anisotropic terms, which generates a complex set of levels characterized by spin mixing. Although the model multispin Hamiltonian in the absence of local orbital momentum, , is considered good enough to map the experimental energies at zero field and in the strong-exchange limit, theoretical works pointed out limitations of this simple model. This work revives the use of HMS from a new theoretical perspective, detailing point-by-point a strategy to correctly map the computational energies and wave functions onto HMS, thus validating it regardless of the exchange limit. We will distinguish two cases, based on experimentally characterized dicobalt(II) complexes from the literature. If centrosymmetry imposes alignment of the various rank-2 tensors constitutive of HMS in the first case, the absence of any symmetry element prevents such alignment in the second case. In such a context, the strategy provided herein becomes a powerful tool to rationalize the experimental magnetic data, since it is capable of fully and rigorously extracting the multispin model without any assumption on the orientation of its constitutive tensors. Furthermore, the strategy allows to question the use of the spin Hamiltonian approach by explicitly controlling the projection norms on the model space, which is showcased in the second complex where local orbital momentum could have occurred (distorted octahedra). Finally, previous theoretical data related to a known dinickel(II) complex is reinterpreted, clarifying initial wanderings regarding the weak exchange limit.

Equatorial Electronic Structure in the Uranyl Ion: Cs2UO2Cl4 and Cs2UO2Br4.

Electric field gradient (EFG) tensors in the equatorial plane of the linear UO22+ ion have been measured by nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) experiments and computed by relativistic Kohn-Sham methods with and without environment embedding for Cs2UO2Cl4 and Cs2UO2Br4. This approach expands the possibilities for probing the electronic structure in uranyl complexes beyond the strongly covalent U-O bonds. The combined analyses find that one of the two largest principal EFG tensor components at the halogen sites points along the U-X bond (X = Cl, Br), and the second is parallel to the UO22+ ion; in Cs2UO2Cl4, the components are nearly equal in magnitude, whereas in Cs2UO2Br4, due to short-range bromide-cesium interactions, the equatorial component is dominant for one pair of Br sites and the axial component is larger for the second pair. The directions and relative magnitudes of the field gradient principal axes are found to be sensitive to the σ and π electron donation by the ligands and the model of the environment. Chlorine-35 NQR spectra of 235U-depleted and 235U-enriched Cs2UO2Cl4 exhibited no uranium-isotope-dependent shift, but the resonance of the depleted sample displayed a 58% broader line width.

Probing Multiconfigurational States by Spectroscopy: The Cerium XAS L3 -edge Puzzle.

Invited for the cover of this issue are Prof. Jochen Autschbach and Dr. Dumitru-Claudiu Sergentu of State University of New York at Buffalo, and Dr. Corwin H. Booth of Lawrence Berkeley National Laboratory. The image depicts high-energy X-ray beams as lightnings probing Ce at the L3 edge in the iconic covalent-bonded Ce(C8 H8 )2 and in CeO2 . The mountain peaks in the background represent the double-peaked L3 edges. The peaks turn out to be intuitively interpreted in terms of localized orbitals and hence metal oxidation states. Read the full text of the article at 10.1002/chem.202100145.

Probing Multiconfigurational States by Spectroscopy: The Cerium XAS L3 -edge Puzzle.

The Ce L3 edge XAS spectra of CeO2 and cerocene [Ce(C8 H8 )2 ] were calculated with relativistic ab-initio multireference wavefunction approaches capable of reproducing the observed spectra accurately. The study aims to resolve the decades-long puzzle regarding the relationship between the number and relative intensities of the XAS peaks and the 4f electron occupation in the ground state (GS) versus the core-excited states (ESs). CeO2 and cerocene exemplify the different roles of covalent bonding and wavefunction configurational composition in the observed intensity patterns. Good agreement is found between the calculated GS 4f-shell occupations and the value derived from XAS measurements using peak areas (nf ). The identity of the two-peaked Ce L3 edge is fully rationalized from the perspective of the relaxed wavefunctions for the GS and core ESs. The states underlying the different peaks differ from each other in a surprisingly simple way that can be associated with 4f1 vs. 4f0 sub-configurations. Furthermore, part of one of the cerocene spectral peaks is associated with 4f2 sub-configurations. The pattern therefore reveals excited states that can be interpreted in terms of Ce IV and III oxidation numbers, as long assumed, with Ce II states additionally appearing in the cerocene spectrum. While this work demonstrates the rough accuracy of the conventional approach to determining nf from Ce L3 -edge XAS, limitations are highlighted in terms of the ultimate accuracy of this approach and the potential of observing new types of excited states. The need to determine the sources of nf by calculations, is stressed.

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 [AnCl6]3-, [An(CN)6]3-, [An(NCS)6]3-, [An(S2PMe2)3], [An(DPA)3]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 n 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.

Quantum chemical topology at the spin-orbit configuration interaction level: Application to astatine compounds.

We report a methodology that allows the investigation of the consequences of the spin-orbit coupling by means of the QTAIM and ELF topological analyses performed on top of relativistic and multiconfigurational wave functions. In practice, it relies on the "state-specific" natural orbitals (NOs; expressed in a Cartesian Gaussian-type orbital basis) and their occupation numbers (ONs) for the quantum state of interest, arising from a spin-orbit configuration interaction calculation. The ground states of astatine diatomic molecules (AtX with X = AtF) and trihalide anions (IAtI- , BrAtBr- , and IAtBr- ) are studied, at exact two-component relativistic coupled cluster geometries, revealing unusual topological properties as well as a significant role of the spin-orbit coupling on these. In essence, the presented methodology can also be applied to the ground and/or excited states of any compound, with controlled validity up to including elements with active 5d, 6p, and/or 5f shells, and potential limitations starting with active 6d, 7p, and/or 6f shells bearing strong spin-orbit couplings.

Ab Initio Analysis of Metal-Ligand Bonding in An(COT)2 with An=Th, U in Their Ground- and Core-Excited States.

Relativistic multireference ab initio wave function calculations with the restricted active space second-order perturbation theory (RASPT2) were performed on thorocene and uranocene to determine the actinide N4,5 -edge and carbon K-edge X-ray absorption near-edge structure (XANES) intensities and the metal-ligand orbital mixing in the ground state and core-excited states. Calculated spectral intensities show very good agreement with the experiments and therefore allow detailed and unambiguous assignment of the observed spectral features. φ-type covalent bonding or antibonding interactions are observed for thorocene in the core-excited states, though not in the ground state. This is because the molecular orbital of φ symmetry, which is the in-phase combination of the ligand Lφ and the Th 5fφ orbitals, can be populated with electrons in core-excited states, whereas it is essentially unoccupied in the ground state. For uranocene, the XANES spectra do not reveal much information beyond multiplet broadening, despite the presence of distinct peaks in the spectra. Every core-excited peak is best characterized by its own set of bond orbitals, as the excited state covalency is clearly different from the ground state covalency.

Probing the Electronic Structure of a Thorium Nitride Complex by Solid-State 15N NMR Spectroscopy.

The solid-state 15N NMR powder spectra of the thorium nitride complex, [K(18-crown-6)(THF)2][(R2N)3Th(µ-15N)Th(NR2)3] ([K][1-15N], R = SiMe3), and the thorium amide complex, [Th(NR2)3(15NH2)] (2-15N), were recorded. The spectrum for [K][1-15N] represents the first reported solid-state 15N NMR data for an actinide complex. The experimentally measured tensor spans are Ω = 847 ppm for [K][1-15N] and Ω = 237 ppm for 2-15N. Both shielding tensors exhibit axial symmetry, which for [K][1-15N] is consistent with a local rotational symmetry of its 15N-labeled nitride ligand. For 2-15N, the axial asymmetry can be rationalized by a quasi-free Th-NH2 bond rotation in the solid-state. Density functional theory calculations overestimate the tensor span somewhat for [K][1-15N], but provide isotropic shifts in good agreement with both the solid-state and solution values for both complexes. Natural localized molecular orbital analyses of the nuclear shielding reveal that the larger tensor span in [K][1-15N] vs 2-15N is primarily a consequence of more pronounced covalency of the σ(N-Th) bonds and large spin-orbit coupling due to significant Th 5f orbital contribution to those bonds, impacting the principal components of the shielding tensor perpendicular to the Th-N-Th axis. Overall, our analysis confirms the involvement of the 5f orbitals in Th-N multiple bonds and further demonstrates the value of solid-state NMR spectroscopy for interrogating actinide-ligand bonding.

Modern quantum chemistry with [Open]Molcas.

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

The Exceptional Diversity of Homoleptic Uranium-Methyl Complexes.

Homoleptic σ-bonded uranium-alkyl complexes have been a synthetic target since the Manhattan Project. The current study describes the synthesis and characterization of several unprecedented uranium-methyl complexes. Amongst these complexes, the first example of a homoleptic uranium-alkyl dimer, [Li(THF)4 ]2 [U2 (CH3 )10 ], as well as a seven-coordinate uranium-methyl monomer, {Li(OEt2 )Li(OEt2 )2 UMe7 Li}n were both crystallographically identified. The diversity of complexes reported herein provides critical insight into the structural diversity, electronic structure and bonding in uranium-alkyl chemistry.

Diuranium(IV) Carbide Cluster U2C2 Stabilized Inside Fullerene Cages.

Novel actinide cluster fullerenes, U2C2@Ih(7)-C80 and U2C2@D3h(5)-C78, were synthesized and fully characterized by mass spectrometry, single-crystal X-ray crystallography, UV-vis-NIR, nuclear magnetic resonance spectroscopy (NMR), X-ray absorption spectroscopy (XAS), Raman spectroscopy, IR spectroscopy, as well as density functional and multireference wave function calculations. The encapsulated U2C2 is the first example of a uranium carbide cluster featuring two U centers bridged by a C≡C unit. The U-C bond distances in these U2C2 clusters are in the range between 2.130 and 2.421 Å. While the U2C2 cluster in U2C2@C80 adopts a butterfly-shaped geometry with a U-C2-U dihedral angle of 112.7° and a U-U distance of 3.855 Å, the U-U distance in U2C2@C78 is 4.164 Å and the resulting U-C2-U dihedral angle is increased to 149.1°. The combined experimental and quantum-chemical results suggest that the formal U oxidation state is +4 in the U2C2 cluster, and each U center transfers three electrons to the C2n cage and one electron to C2. Different from the strong U═C covalent bonding reported for U2C@C80, the U-C bonds in U2C2 are less covalent and predominantly ionic. The C-C triple bond is somewhat weaker than in HCCH, and the C-C π bonds undergo donation bonding with the U centers. This work demonstrates that the combination of the unique encapsulation effect of fullerene cages and the variable oxidation states of actinide elements can lead to the stabilization of novel actinide clusters, which are not accessible by conventional synthetic methods.

Synthesis, Characterization, and Electrochemistry of the Homoleptic f Element Ketimide Complexes [Li]2[M(N═CtBuPh)6] (M = Ce, Th).

Reaction of [Ce(NO3)3(THF)4] with 6 equiv of Li(N═CtBuPh), followed by addition of 0.5 equiv of I2, affords the homoleptic Ce(IV) ketimide [Li]2[Ce(N═CtBuPh)6] (1), which can be isolated in 44% yield after workup. Similarly, reaction of [ThCl4(DME)2] (DME = 1,2-dimethoxyethane) with 6 equiv of Li(N═CtBuPh) in tetrahydrofuran affords the isostructural Th(IV) ketimide [Li]2[Th(N═CtBuPh)6] (2), which can be isolated in 53% yield after workup. Both 1 and 2 were fully characterized, including analysis by X-ray crystallography, allowing for a detailed structural and spectroscopic comparison. The electronic structures of 1 and 2 were also explored with density functional theory and multiconfigurational wave function calculations. Additionally, the redox chemistry of 1 was probed by cyclic voltammetry, which revealed a highly cathodic Ce(IV)/Ce(III) reduction potential, providing evidence for the ability of the ketimide ligand to stabilize high oxidation states of the lanthanides.

Magnetic circular dichroism spectra of transition metal complexes calculated from restricted active space wavefunctions.

First principles multiconfigurational restricted active space (RAS) self-consistent field (SCF) or configuration interaction (CI) approaches, augmented with a treatment of spin-orbit coupling by state interaction, were used to calculate the magnetic circular dichroism (MCD) , , and/or for closed- and open-shell transition metal complexes: PdCl42-, PdBr42-, AuCl4-, AuBr4-, MnO4-, CuCl42-, CuBr42-, and Fe(CN)63-. The were determined with a sum-over-states approach. It is shown that reasonably accurate MCD spectra can be obtained directly at a RAS level or at a RAS level augmented with corrections for the dynamic correlation. The sign and magnitude of the individual MCD terms can be unambiguously determined and assigned to particular electronic transitions.

Homoleptic Aryl Complexes of Uranium (IV).

The synthesis and characterization of sterically unencumbered homoleptic organouranium aryl complexes containing U-C σ-bonds has been of interest to the chemical community for over 70 years. Reported herein are the first structurally characterized, sterically unencumbered homoleptic uranium (IV) aryl-ate species of the form [U(Ar)6 ]2- (Ar=Ph, p-tolyl, p-Cl-Ph). Magnetic circular dichroism (MCD) spectroscopy and computational studies provide insight into electronic structure and bonding interactions in the U-C σ-bond across this series of complexes. Overall, these studies solve a decades-long challenge in synthetic uranium chemistry, enabling new insight into electronic structure and bonding in organouranium complexes.

Understanding and Controlling the Emission Brightness and Color of Molecular Cerium Luminophores.

Molecular cerium complexes are a new class of tunable and energy-efficient visible- and UV-luminophores. Understanding and controlling the emission brightness and color are important for tailoring them for new and specialized applications. Herein, we describe the experimental and computational analyses for series of tris(guanidinate) (1-8, Ce{(R2N)C(N iPr)2}3, R = alkyl, silyl, or phenyl groups), guanidinate-amide [GA, A = N(SiMe3)2, G = (Me3Si)2NC(N iPr)2], and guanidinate-aryloxide (GOAr, OAr = 2,6-di- tert-butylphenoxide) cerium(III) complexes to understand and develop predictive capabilities for their optical properties. Structural studies performed on complexes 1-8 revealed marked differences in the steric encumbrance around the cerium center induced by various guanidinate ligand backbone substituents, a property that was correlated to photoluminescent quantum yield. Computational studies revealed that consecutive replacements of the amide and aryloxide ligands by guanidinate ligand led to less nonradiative relaxation of bright excited states and smaller Stokes shifts. The results establish a comprehensive structure-luminescence model for molecular cerium(III) luminophores in terms of both quantum yields and colors. The results provide a clear basis for the design of tunable, molecular, cerium-based, luminescent materials.

The bonding picture in hypervalent XF3 (X = Cl, Br, I, At) fluorides revisited with quantum chemical topology.

Hypervalent XF3 (X = Cl, Br, I, At) fluorides exhibit T-shaped C2V equilibrium structures with the heavier of them, AtF3 , also revealing an almost isoenergetic planar D3h structure. Factors explaining this behavior based on simple "chemical intuition" are currently missing. In this work, we combine non-relativistic (ClF3 ), scalar-relativistic and two-component (X = Br - At) density functional theory calculations, and bonding analyses based on the electron localization function and the quantum theory of atoms in molecules. Typical signatures of charge-shift bonding have been identified at the bent T-shaped structures of ClF3 and BrF3 , while the bonds of the other structures exhibit a dominant ionic character. With the aim of explaining the D3h structure of AtF3 , we extend the multipole expansion analysis to the framework of two-component single-reference calculations. This methodological advance enables us to rationalize the relative stability of the T-shaped C2v and the planar D3h structures: the Coulomb repulsions between the two lone-pairs of the central atom and between each lone-pair and each fluorine ligand are found significantly larger at the D3h structures than at the C2v ones for X = Cl - I, but not with X = At. This comes with the increasing stabilization, along the XF3 series, of the planar D3h structure with respect to the global T-shaped C2v minima. Hence, we show that the careful use of principles that are at the heart of the valence shell electron pair repulsion model provides reasonable justifications for stable planar D3h structures in AX3 E2 systems. © 2017 Wiley Periodicals, Inc.

Scrutinizing "Invisible" astatine: A challenge for modern density functionals.

The main-group 6p elements did not receive much attention in the development of recent density functionals. In many cases it is still difficult to choose among the modern ones a relevant functional for various applications. Here, we illustrate the case of astatine species (At, Z = 85) and we report the first, and quite complete, benchmark study on several properties concerning such species. Insights on geometries, transition energies and thermodynamic properties of a set of 19 astatine species, for which reference experimental or theoretical data has been reported, are obtained with relativistic (two-component) density functional theory calculations. An extensive set of widely used functionals is employed. The hybrid meta-generalized gradient approximation (meta-GGA) PW6B95 functional is overall the best choice. It is worth noting that the range-separated HSE06 functional as well as the old and very popular B3LYP and PBE0 hybrid-GGAs appear to perform quite well too. Moreover, we found that astatine chemistry in solution can accurately be predicted using implicit solvent models, provided that specific parameters are used to build At cavities. © 2016 Wiley Periodicals, Inc.

Advances on the Determination of the Astatine Pourbaix Diagram: Predomination of AtO(OH)2 (-) over At(-) in Basic Conditions.

It is generally assumed that astatide (At(-) ) is the predominant astatine species in basic aqueous media. This assumption is questioned in non-complexing and non-reductive aqueous solutions by means of high-pressure anion-exchange chromatography. Contrary to what is usually believed, astatide is found to be a minor species at pH=11. A different species, which also bears a single negative charge, becomes predominant when the pH is increased beyond 7. Using competition experiments, an equilibrium constant value of 10(-6.9) has been determined for the formation of this species from AtO(OH) with the exchange of one proton. The identification of this species, AtO(OH)2 (-) , is achieved through relativistic quantum mechanical calculations, which rule out the significant formation of the AtO2 (-) species, while leading to a hydrolysis constant of AtO(OH) in excellent agreement with experiment when the AtO(OH)2 (-) species is considered. Beyond the completion of the Pourbaix diagram of astatine, this new information is of interest for the development of (211) At radiolabeling protocols.