Magnetic Properties and Electronic Structure of Neptunyl(VI) Complexes: Wavefunctions, Orbitals, and Crystal-Field Models

A Comprehensive Multireference Study of Excited-State Ni-Br Bond Homolysis in (dtbbpy)NiII(aryl)(Br)

The mechanism of visible light-driven Ni-C(aryl) bond homolysis in (2,2'-bipyridine)NiII(aryl)(halide) complexes, which play a crucial role in metallaphotoredox catalysis for cross-coupling reactions, has been well studied. Differently, the theoretical understanding of Ni-halide bond homolysis remains limited. In this study, we introduce a novel electronic structural framework to elucidate the mechanisms underlying photoinduced Ni-Br bond rupture in the (dtbbpy)NiII(aryl)(Br) complex. Using multireference ab initio calculations, we characterized the excited state potential energy surfaces corresponding to metal-to-ligand charge transfer (MLCT) and ligand-to-metal charge transfer (LMCT). Our calculations reveal that the Ni-Br dissociation, triggered by an external photocatalyst, begins with the promotion of Ni(II) to a 1MLCT excited state. This state undergoes intersystem crossing with repulsive triplet surfaces corresponding to the 3MLCT and Br-to-Ni 3LMCT states, resulting in Ni-Br bond breaking via the Dexter energy transfer mechanism. In the absence of a photocatalyst, the photoexcited Ni(II) favors Ni-C(aryl) homolysis, whereas the presence of a photocatalyst promotes Ni-Br dissociation. The Ni(III) species, resulting from the oxidation of Ni(II) by the photocatalyst, was found to be unproductive toward Ni-Br or Ni-C(aryl) activation.

Actinyl Electronic Structure Probed by XAS: The Role of Many-Body Effects

A detailed analysis of the wave functions for the M5 to 5f excitations in the linear actinyls, UO22+, NpO22+, and PuO22+, and the theoretical X-ray absorption spectra obtained with these wave functions in comparison with experimental M5-edge high-resolution X-ray absorption near-edge structure (HR-XANES) spectra is presented. The wave functions include full treatment of scalar and spin-orbit relativistic effects through the use of a Dirac-Coulomb Hamiltonian; many-body effects are included in determining the wave functions. The character of the excited states and of the active spaces to describe the wave functions for these states are investigated and analyzed. It is shown that the excited states cannot, in general, be described with a single configuration but have an essential multiconfiguration character. The characterization of the properties of the excited states and the X-ray absorption spectra was achieved through the use of novel methods.

Radiation-Induced Defects in Uranyl Trinitrate Solids

Actinides are inherently radioactive; thus, ionizing radiation is emitted by these elements can have profound effects on its surrounding chemical environment through the formation of free radical species. While previous work has noted that the presence of free radicals in the system impacts the redox state of the actinides, there is little atomistic understanding of how these metal cations interact with free radicals. Herein, we explore the effects of radiation (UV and γ) on three U(VI) trinitrate complexes, M[UO2(NO3)3] (where M=K+, Rb+, Cs+), and their respective nitrate salts in the solid state via electron paramagnetic resonance (EPR) and Raman spectroscopy paired with Density Functional Theory (DFT) methods. We find that the alkali salts form nitrate radicals under UV and γ irradiation, but also note the presence of additional degradation products. M[UO2(NO3)3] solids also form nitrate radicals and additional DFT calculations indicate the species corresponds to a change from the bidentate bound nitrate anion into a monodentate NO3 • radical. Computational studies also highlight the need to include the second sphere coordination environment around the [UO2(NO3)3]0,1 species to gain agreement between the experimental and predicted EPR signatures.

Consistent Second-Order Treatment of Spin-Orbit Coupling and Dynamic Correlation in Quasidegenerate N-Electron Valence Perturbation Theory

We present a formulation and implementation of second-order quasidegenerate N-electron valence perturbation theory (QDNEVPT2) that provides a balanced and accurate description of spin-orbit coupling and dynamic correlation effects in multiconfigurational electronic states. In our approach, the energies and wave functions of electronic states are computed by treating electron repulsion and spin-orbit coupling operators as equal perturbations to the nonrelativistic complete active-space wave functions, and their contributions are incorporated fully up to the second order. The spin-orbit effects are described using the Breit-Pauli (BP) or exact two-component Douglas-Kroll-Hess (DKH) Hamiltonians within spin-orbit mean-field approximation. The resulting second-order methods (BP2- and DKH2-QDNEVPT2) are capable of treating spin-orbit coupling effects in nearly degenerate electronic states by diagonalizing an effective Hamiltonian expanded in a compact non-relativistic basis. For a variety of atoms and small molecules across the entire periodic table, we demonstrate that DKH2-QDNEVPT2 is competitive in accuracy with variational two-component relativistic theories. BP2-QDNEVPT2 shows high accuracy for the second- and third-period elements, but its performance deteriorates for heavier atoms and molecules. We also consider the first-order spin-orbit QDNEVPT2 approximations (BP1- and DKH1-QDNEVPT2), among which DKH1-QDNEVPT2 is reliable but less accurate than DKH2-QDNEVPT2. Both DKH1- and DKH2-QDNEVPT2 hold promise as efficient and accurate electronic structure methods for treating electron correlation and spin-orbit coupling in a variety of applications.

Electronic Structure of Actinyls: Orbital Properties

A detailed analysis is presented for the covalent character of the orbitals in the actinyls: UO22+, NpO22+, and PuO22+. Both the initial, or ground state, GS, configuration and the excited configurations where a 3d electron is excited into the open valence, nominally the 5f shell, are considered. The orbitals are determined as fully relativistic, four component Dirac-Coulomb Hartree-Fock solutions. Several measures, which go beyond the commonly used population analyses, are used to characterize the covalent character of an orbital in order to obtain reliable estimates of the covalency. Although there are differences in the covalent character of the orbitals for the initial and excited configurations of the different actinyls, there is a surprising similarity in the covalent character for all of the states considered. This is true both between the initial and excited configurations as well as between the different actinyls. The analysis emphasizes the 5f covalent character in the closed shell bonding orbitals and the open shell antibonding orbitals since the focus is on characterizing orbitals needed in a many-body treatment of the actinyl wave functions. However, estimates are also made of the participation of the actinide 6d in the covalent bonding.

Multiconfiguration Pair-Density Functional Theory for Vertical Excitation Energies in Actinide Molecules

Modeling actinides with electronic structure theories is challenging because these systems present a strong ligand field and metal-ligand covalency. We systematically investigate the effectiveness of pair-density functional theory (PDFT) for the calculation of vertical excitation energies in An(III), [AnIIICl6]3-, and [AnVIO2]2+ (An = U, Np, Pu, and Am). We compare the performance of PDFT, hybrid PDFT, and multistate PDFT with traditional active-space methods followed by perturbation theory, like multistate CASPT2, and with experimental data. Overall, multistate PDFT gives quantitative agreement with multistate CASPT2 at a significantly reduced computational cost.

The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

Magnetic Properties of Tetravalent Pu in the Perovskites BaPuO3 and SrPuO3

BaPuO3 and SrPuO3 were synthesized, and their structures were refined in the orthorhombic space group Pbnm, a common distortion from the classic Pm3̅m cubic perovskite. Magnetic-susceptibility measurements, obtained as a function of temperature over the range of 1.8-320 K, exhibit temperature-dependent behavior, with evidence of long-range magnetic order at temperatures higher than their lanthanide and actinide analogues: BaPuO3 below 164(1) K and SrPuO3 below 76(1) K. Effective moments of 1.66(10)μB for BaPuO3 and 1.84(8)μB for SrPuO3 were obtained by fitting their paramagnetic susceptibilities using the Curie-Weiss law. Both are below the free-ion value of 2.68 μB expected for a Pu4+ 5I4 ground level. Ab initio wave function calculations, performed at the relativistic complete active space level including spin-orbit coupling and with an embedded cluster approach that neglects interactions between Pu centers, were used to generate embedded-cluster Pu4+ magnetic susceptibilities. The calculations agree well with experimental data at higher temperatures, providing evidence that a single-ion representation is sufficient to account for the observed paramagnetic behavior without the need to invoke charge transfer, disproportionation, strong covalent bonding, or other more complex electronic behavior.

Are Actinyl Cations Good Probes for Structure Determination in Solution by NMR?

We report on NMR spectroscopy, CAS-based method calculations, and X-ray diffraction of An^V and An^VI complexes with a neutral and slightly flexible TEDGA ligand. After checking that pNMR shifts mainly arise from pseudocontact interactions, we analyze pNMR shifts considering the axial and rhombic anisotropy of the actinyl magnetic susceptibilities. The results are compared to those of a previous study performed on [An^VIO₂]²⁺ complexes with dipicolinic acid. It is shown that 5f² cations (Pu^VI and Np^V) make very good candidates for determining the structure of actinyl complexes in solution by ¹H NMR spectroscopy as shown by the invariance of the magnetic properties to the equatorial ligands, as opposed to the Np^VI complexes with a 5f¹ configuration.

Relativistic Semistochastic Heat-Bath Configuration Interaction

In this work we present the extension of semistochastic heat-bath configuration interaction (SHCI) to work with any two-component and four-component Hamiltonian. The vertical detachment energy (VDE) of AuH₂^- and zero-field splitting (ZFS) of NpO₂²⁺ is calculated by correlating more than 100 spinors in both cases. This work demonstrates the capability of SHCI to treat problems where both relativistic effect and electron correlation are important.

Theoretical Study of the Excited States and Luminescent Properties of (H2O)nUO2Cl2 (n = 1-3)

The low-lying excited-state properties of the water-solvated UO₂Cl₂ complexes, i.e., (H₂O)_nUO₂Cl₂ (n = 1-3), below 33,000 cm^-1, are investigated based on the ab initio NEVPT2 and CCSD(T) with inclusion of scalar relativistic and spin-orbit coupling effects. The simulated luminescence spectral curves agree well with the experimental spectrum in aqueous solution at -120 °C. Water coordination is found to significantly affect the character of luminescent state, which is changed from the ³Φ_g state in UO₂Cl₂ to the ³Δ_g state in (H₂O)_2,3UO₂Cl₂. This is distinctly different from the observed unchanged nature of luminescent state in the cases of Ar coordination to UO₂Cl₂ and H₂O coordination to UO₂F₂ in the previous work. Furthermore, by combining with the theoretical results for the solvated UO₂F₂ system, the reason why water coordination does not remarkably change the spectral shape of UO₂Cl₂, as opposed to UO₂F₂, was explained based on the analysis of two key spectral parameters, O-U-O symmetrical vibrational frequency and U-O bond length elongation. The roles of ligand field and spin-orbit coupling in the determination of luminescent state character and spectral shape in uranyl dihalide complexes are deeply discussed and summarized. These results deepen our understanding of the luminescent properties of uranyl complexes in aqueous solution.

Simulating Spin-Orbit Coupling with Quasidegenerate N-Electron Valence Perturbation Theory

We present the first implementation of spin-orbit coupling effects in fully internally contracted second-order quasidegenerate N-electron valence perturbation theory (SO-QDNEVPT2). The SO-QDNEVPT2 approach enables the computations of ground- and excited-state energies and oscillator strengths combining the description of static electron correlation with an efficient treatment of dynamic correlation and spin-orbit coupling. In addition to SO-QDNEVPT2 with the full description of one- and two-body spin-orbit interactions at the level of two-component Breit-Pauli Hamiltonian, our implementation also features a simplified approach that takes advantage of spin-orbit mean-field approximation (SOMF-QDNEVPT2). The accuracy of these methods is tested for the group 14 and 16 hydrides, 3d and 4d transition metal ions, and two actinide dioxides (neptunyl and plutonyl dications). The zero-field splittings of group 14 and 16 molecules computed using SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the available experimental data. For the 3d transition metal ions, the SO-QDNEVPT2 method is significantly more accurate than SOMF-QDNEVPT2, while no substantial difference in the performance of two methods is observed for the 4d ions. Finally, we demonstrate that for the actinide dioxides the results of SO-QDNEVPT2 and SOMF-QDNEVPT2 are in good agreement with the data from previous theoretical studies of these systems. Overall, our results demonstrate that SO-QDNEVPT2 and SOMF-QDNEVPT2 are promising multireference methods for treating spin-orbit coupling with a relatively low computational cost.

A Coupled EXAFS-Molecular Dynamics Study on PuO2+ and NpO2+ Hydration: The Importance of Electron Correlation in Force-Field Building

The physicochemical properties of the monovalent actinyl cations, PuO₂⁺ and NpO₂⁺, in water have been studied by means of classical molecular dynamic simulations. A specific set of cation-water intermolecular potentials based on ab initio potential energy surfaces has been built on the basis of the hydrated ion concept. The TIP4P water model was adopted. Given the paramagnetic character of these actinyls, the cation-water interaction energies were computed from highly correlated wave functions using the NEVPT2 method. It is shown that the multideterminantal character of the wave function has a relevant effect on the main distances of the hydrated molecular cations. Several structural, dynamical, and energetic properties of the aqueous solutions have been obtained and analyzed. Structural RDF analysis gives An-O_yl distances of 1.82 and 1.84 Å and An-O(water) distances of 2.51 and 2.53 Å for PuO₂⁺ and NpO₂⁺ in water, respectively. Experimental EXAFS spectra from dilute aqueous solutions of PuO₂⁺ and NpO₂⁺ are revisited and analyzed, assuming tetra- and pentahydration of the actinyl cations. Simulated EXAFS spectra have been computed from the snapshots of the MD simulations. Good agreement with the experimental information available is found. The global analysis leads us to conclude that both PuO₂⁺ and NpO₂⁺ cations in water are stable pentahydrated aqua ions.

Near-infrared circular dichroism of the ytterbium DOTMA complex: an ab initio investigation

The electronic structure and circular dichroism spectra of the ytterbium(III) complex [Yb(DOTMA)]^- are calculated using complete and restricted active space self-consistent field wavefunction methods with the spin-orbit coupling treated by the state interaction approach. The influence of the dynamical correlation effect is then included via the 2nd order perturbation method. The experimental circular dichroism spectrum is well reproduced by calculations, both in terms of relative energy excitations and in terms of rotatory strength intensities. The results allow highlighting the mechanism that drives the chiroptical properties in Yb(III) complexes and reveal the importance of taking into account the 4f¹²5d¹ electronic configurations in the calculated wavefunctions to properly describe the chiroptical properties of the 4f-4f transitions.

X-ray absorption spectra of f-element complexes: insight from relativistic multiconfigurational wavefunction theory

X-ray absorption near edge structure (XANES) spectroscopy, coupled with ab initio calculations, has emerged as the state-of-the-art tool for elucidating the metal-ligand bonding in f-element complexes. This highlight presents recent efforts in calculating XANES spectra of lanthanide and actinide compounds with relativistic multiconfiguration wavefunction approaches that account for differences in donation bonding in the ground state (GS) versus a core-excited state (ES), multiplet effects, and spin-orbit-coupling. With the GS and ES wavefunctions available, including spin-orbit effects, an arsenal of chemical bonding tools that are popular among chemists can be applied to rationalize the observed intensities in terms of covalent bonding.

Spin-Orbit Natural Transition Orbitals and Spin-Forbidden Transitions

Natural transition orbitals (NTOs) are in widespread use for visualizing and analyzing electronic transitions. The present work introduces the analysis of formally spin-forbidden transitions with the help of complex-valued spin-orbit (SO) NTOs. The analysis specifically focuses on the components in such transitions that cause their intensity to be nonzero because of SO coupling. Transition properties such as transition dipole moments are partitioned into SO-NTO hole-particle pairs, such that contributions to the intensity from specific occupied and unoccupied orbitals are obtained. The method has been implemented within the restricted active space (RAS) self-consistent field wave function theory framework, with SO coupling treated by RAS state interaction. SO-NTOs have a broad range of potential applications, which is illustrated by the T₂-S₁ state mixing in pyrazine, spin-forbidden versus spin-allowed 4f-5d transitions in the Tb³⁺ ion, and the phosphorescence of tris(2-phenylpyridine) iridium [Ir(ppy)₃].

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.

Combined Experimental/Theoretical Study on the Luminescent Properties of Homoleptic/Heteroleptic Erbium(III) Anilate-Based 2D Coordination Polymers

The synthesis, structural and photophysical characterization, and theoretical studies on homo/heteroleptic neutral 2D-layered coordination polymers (CPs), obtained by combining the Er^III ion with chlorocyananilate (ClCNAn) and/or tetrafluoroterephthalate (F₄BDC) linkers, are herein reported. The structure of the heteroleptic Er^III-based CP, formulated as [Er₂(ClCNAn)₂(F₄BDC)(DMSO)₆]_n (1) is also reported. 1 crystallizes in the triclinic P1̅ space group, and the structure consists of neutral 2D layers formed by Er^III ions linked through the two linkers oriented in such a way that the neighboring 2D layers are eclipsed along the a axis, leading to parallelogram-like cavities. Photophysical measurements highlight the prominent role of chlorocyananilate linkers as optical antennas toward lanthanide ions, while wave-function-theory analysis supports the experimental findings, providing evidence for the effect of ligand substitution on the luminescence properties of homo/heteroleptic 2D CPs.

Computational and Spectroscopic Tools for the Detection of Bond Covalency in Pu(IV) Materials

Plutonium is used as a major component of new-generation nuclear fuels and of radioisotope batteries for Mars rovers, but it is also an environmental pollutant. Plutonium clearly has high technological and environmental importance, but it has an extremely complex, not well-understood electronic structure. The level of covalency of the Pu 5f valence orbitals and their role in chemical bonding are still an enigma and thus at the frontier of research in actinide science. We performed fully relativistic quantum chemical computations of the electronic structure of the Pu⁴⁺ ion and the PuO₂ compound. Using four different theoretical tools, it is shown that the 5f orbitals have very little covalent character although the 5f(_7/2) a_2u orbital with the highest orbital energy has the greatest extent of covalency in PuO₂. It is illustrated that the Pu M_4,5 edge high-energy resolution X-ray absorption near-edge structure (Pu M_4,5 HR-XANES) spectra cannot be interpreted in terms of dipole selection rules applied between individual 3d and 5f orbitals, but the selection rules must be applied between the total wavefunctions for the initial and excited states. This is because the states cannot be represented by single determinants. They are shown to involve major redistributions on the 5f electrons over the different 5f orbitals. These redistributions could be viewed as shake-up-like excitations in the 5f shell from the lowest orbital energy from J = 5f(_5/2) into higher orbital energy J = 5f(_7/2). We show that the second peak in the Pu M₄ edge and the high-energy shoulder of the Pu M₅ edge HR-XANES spectra probe the 5f(_7/2) a_2u orbital; thus, these spectral features are expected to change upon bond variations. We describe theoretical and spectroscopy tools, which can be applied for all actinide elements in materials with cubic structure.

Enhanced 5f-δ bonding in [U(C7H7)2]-: C K-edge XAS, magnetism, and ab initio calculations

5f covalency in [U(C₇H₇)₂]^- was probed with carbon K-edge X-ray absorption spectroscopy (XAS) and electronic structure theory. The results revealed U 5f orbital participation in δ-bonding in both the ground- and core-excited states; additional 5f ϕ-mixing is observed in the core-excited states. Comparisons with U(C₈H₈)₂ show greater δ-covalency for [U(C₇H₇)₂]^-.

Magnetic properties of organolanthanide(II) complexes, from the electronic structure and the crystal field effect

The magnetic properties of a series of organometallic complexes [LnCp3]- and Ln(CNT)2, where Cp = cyclopentadienyl and CNT = cyclononatetraenyl, of the lanthanide ions in the 2+ oxidation state, are theoretically studied in terms of the electronic structure obtained via multiconfigurational wave function-based methods. Calculations are performed for two groups of ion complexes selected based on their preferred electronic configuration 4fn+1 or 4fn5d1 (n is the number of f electrons in the 3+ ion). All the properties are discussed in terms of the electron density distribution of the ground state and ligand field effects. This analysis allows giving some molecular design strategies relevant to exploit the magnetic properties in applications like Single-Molecule Magnets (SMMs) for lanthanide ions in the 2+ oxidation state.

Spectroscopic characterization of neptunium(VI), plutonium(VI), americium(VI) and neptunium(V) encapsulated in uranyl nitrate hexahydrate

The coordination of crystalline products resulting from the co-crystallization of Np(vi), Pu(vi), Am(vi), and Np(v) with uranyl nitrate hexahydrate (UNH) has been revealed through solid-state spectroscopic characterization via diffuse reflectance UV-Vis-NIR spectroscopy, SEM-EDS, and extended X-ray absorption fine structure (EXAFS) spectroscopy. Density functional and multireference wavefunction calculations were performed to analyze the An(vi/v)O2(NO3)2·2H2O electronic structures and to help assign the observed transitions in the absorption spectra. EXAFS show a similar coordination between the U(VI) in UNH and Np(vi) and Pu(vi); while Am resulted in a similar coordination to Am(iii), as reduction of Am(vi) occurred prior to EXAFS data being obtained. The co-crystallization of the oxidized transuranic species-penta- and hexavalent-with UNH, represents a significant advance from not only a practical standpoint in providing an elegant solution for used nuclear fuel recycle, but also as an avenue to expand the fundamental understanding of the 5f electronic behavior in the solid-state.

Temperature Dependence of 1 H Paramagnetic Chemical Shifts in Actinide Complexes, Beyond Bleaney's Theory: The AnVI O2 2+ -Dipicolinic Acid Complexes (An=Np, Pu) as an Example

Actinide +VI complexes ( A n V I = U V I , N p V I and P u V I ) with dipicolinic acid derivatives were synthesized and characterized by powder XRD, SQUID magnetometry and NMR spectroscopy. In addition, N p V I and P u V I complexes were described by first principles CAS based and two-component spin-restricted DFT methods. The analysis of the ¹ H paramagnetic NMR chemical shifts for all protons of the ligands according to the X-rays structures shows that the Fermi contact contribution is negligible in agreement with spin density determined by unrestricted DFT. The magnetic susceptibility tensor is determined by combining SQUID, pNMR shifts and Evans' method. The SO-RASPT2 results fit well the experimental magnetic susceptibility and pNMR chemical shifts. The role of the counterions in the solid phase is pointed out; their presence impacts the magnetic properties of the N p V I complex. The temperature dependence of the pNMR chemical shifts has a strong 1 / T contribution, contrarily to Bleaney's theory for lanthanide complexes. The fitting of the temperature dependence of the pNMR chemical shifts and SQUID magnetic susceptibility by a two-Kramers-doublet model for the N p V I complex and a non-Kramers-doublet model for the P u V I complex allows for the experimental evaluation of energy gaps and magnetic moments of the paramagnetic center.

Luminescence, chiroptical, magnetic and ab initio crystal-field characterizations of an enantiopure helicoidal Yb(iii) complex

The electronic structure of a chiral Yb(iii)-based complex is fully determined by taking advantage of experimental magnetic, luminescence, and chiroptical (NIR-ECD and CPL) characterizations in combination with ab initio wavefunction calculations.

Exploring the Oxidation States of Neptunium with Schiff Base Coordination Complexes

A pair of neptunium Schiff base coordination complexes, Np^VIO₂L(MeOH) and Np^IVL₂ {H₂L = N,N'-bis[(4,4'-diethylamino)salicylidene]-1,2-phenylenediamine}, have been synthesized and analyzed by several characterization methods including single-crystal X-ray diffraction, electronic absorption, ¹H NMR, cyclic voltammetry, and theoretical interpretation. Structural analysis reveals that Np^VIO₂L(MeOH) and Np^IVL₂ are isomorphous with the previously reported U^VIO₂L(MeOH) and M^IVL₂ (M = Pu, Ce, U, Th) complexes, respectively, allowing for a direct comparison across the series. The reduction of Np^VIO₂L(MeOH) in situ or direct synthesis from a (Np^VO₂)⁺ source shows evidence of a pentavalent neptunyl (Np^VO₂L)_x^n- species as determined by UV/vis/NIR and ¹H NMR spectroscopy. The synthesis of (Np^VO₂L)_x^n- directly from a (Np^VO₂)⁺ starting material gives a similar spectrum. Theoretical analysis offers insight into the electronic structure for a better understanding of the bonding patterns and relative stability of the different oxidation states. Computational results show that the Np-L covalent interactions in Np^IVL₂ are similar to those in the Np^VIO₂L(MeOH) complex, indicating that neither the presence of the axial oxo ligands nor the oxidation state significantly modify the nature of the Np-L bonds.

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.

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.

Prediction of Binding Stability of Pu(IV) and PuO2(VI) by Nitrogen Tridentate Ligands in Aqueous Solution

Plutonium has potential applications in energy production in well-controlled nuclear reactors. Since nuclear power plants have great merit as environmentally friendly energy sources with a recyclable system, a recycling system for extracting Pu from spent fuels using suitable extractants has been proposed. Pu leakage is a potential environmental hazard, hence the need for chemical sensor development. Both extractants and chemical sensors involve metal–ligand interactions and to develop efficient extractants and chemical sensors, structural information about Pu ligands must be obtained by quantum calculations. Herein, six representative nitrogen tridentate ligands were introduced, and their binding stabilities were evaluated. The tridentate L6, which contains tri-pyridine chelate with benzene connectors, showed the highest binding energies for Pu(IV) and PuO₂(VI) in water. Analysis based on the quantum theory of atoms in molecular analysis, including natural population analysis and electron density studies, provided insight into the bonding characteristics for each structure. We propose that differences in ionic bonding characteristics account for the Pu-ligand stability differences. These results form a basis for designing novel extractants and organic Pu sensors.

Energy decomposition analysis based on broken symmetry unrestricted density functional theory

In this paper, the generalized Kohn-Sham energy decomposition analysis (GKS-EDA) scheme is extended to molecular interactions in open shell singlet states, which is a challenge for many popular EDA methods due to the multireference character. Based on broken symmetry (BS) unrestricted density functional theory with a spin projection approximation, the extension scheme, named GKS-EDA(BS) in this paper, divides the total interaction energy into electrostatic, exchange-repulsion, polarization, correlation, and dispersion terms. Test examples include the pancake bond in the phenalenyl dimer, the ligand interactions in the Fe(ii)-porphyrin complexes, and the radical interactions in dehydrogenated guanine-cytosine base pairs and show that GKS-EDA(BS) is a practical EDA tool for open shell singlet systems.

Ligand effects on electronic structure and bonding in U(iii) coordination complexes: a combined MCD, EPR and computational study

Spectroscopy and theory enable broader insight into electronic structure and bonding in U(iii) coordination complexes, focusing on systems with Tp* ligands.

Carbon–sulfur bond strength in methanesulfinate and benzenesulfinate ligands directs decomposition of Np(v) and Pu(v) coordination complexes

Adjusting intra-ligand bond strengths in actinide sulfinate complexes directs towards alternative cleavage of carbon–sulfur or actinide–sulfinate bonds.

Magnetic Coupling in the Ce(III) Dimer Ce2(COT)3

The monomer [Ce(COT)₂]^- and the dimer [Ce₂(COT)₃], with Ce(III) and COT = 1,3,5,7-cyclooctatetraenide, are studied by quantum chemistry calculations. Due to the large spin-orbit coupling, the ground state of the monomer is a strong mixing of σ and π states. The experimental isotropic coupling in the dimer was evaluated by Walter et al. to be J = -7 cm^-1 (with a Heisenberg Hamiltonian [Formula: see text]) with a small anisotropic coupling of 0.02 cm^-1. The coupling between the two Ce(III) in the dimer is calculated using CI methods. The low energy part of the spectra are modeled by spin Hamiltonians. All spin Hamiltonians parameters are deduced from ab initio calculations. g factors are calculated for both the pseudodoublet of the monomer and the pseudotriplet of the dimer and their sign have been determined. The magnetic coupling in the dimer is rationalized by a model based on crystal field theory. The kinetic and exchange contributions arising from the different configurations to the isotropic and anisotropic couplings are evaluated. It is shown that the main contribution to isotropic coupling is kinetic and originates from the f_σ-f_σ interaction due to the large transfer integral between those orbitals. However, the f_π-f_π interaction plays a non-negligible role. The anisotropic coupling originates from the difference of exchange energy of states arising from the f_σf_π configuration and is, in no matter, related to the anisotropy of the local magnetic moments as already pointed by van Vleck for a fictitious s-p system. The analysis of the natural orbitals evidences a superexchange mechanism through a σ_CH^* orbital of the bridging cycle favored by a local 4f_σ/5d_σ hybridization and that the δ type orbitals, both the HOMOs of the ligands and the virtual f_δ orbitals of the cerium atoms play an important polarization role, and to a less extend the π type orbitals, the HOMOs-1 of the ligands, and the metal f_π orbitals.

Elucidating cation–cation interactions in neptunyl dications using multi-reference ab initio theory

Understanding the binding mechanism in neptunyl clusters formed due to cation–cation interactions is of crucial importance in nuclear waste reprocessing and related areas of research.

Ab Initio Study of Covalency in the Ground versus Core-Excited States and X-ray Absorption Spectra of Actinide Complexes

Relativistic multireference ab initio wave function calculations within the restricted active space (RAS) framework were performed to calculate metal and ligand X-ray absorption (XAS) near-edge spectroscopy (XANES) intensities for the metal M_4,5 edges of [PuO₂(H₂O)₅]²⁺, [An^VIO₂]²⁺ (An = U, Np, Pu), and [AmCl₆]^3- and the Cl K edge of the Am complex. The extent of An(5f)-ligand bonding was determined via natural localized molecular orbital analyses of the relevant spin-orbit coupled multireference states. The calculated spectra are in good agreement with experiments and allow a detailed assignment of the observed spectral features. The XANES M_4,5-edge spectra are representative of the actinide orbital covalency in the probed core-excited states, which may be different from the ground-state covalency. An assignment of ground-state An orbital covalency based on XAS spectra should therefore be made with caution.

Structure and properties of [(4,6-tBu2C6H2O)2Se]2An(THF)2, An = U, Np, and their reaction with p-benzoquinone

The synthesis and characterization of U(iv) and Np(iv) selenium bis(phenolate) complexes are reported. The reaction of two equivalents of the U(iv) complex with p-benzoquinone results in the formation of a U(v)-U(v) species with a bridging reduced quinone. This represents a rare example of high-valent uranium chemistry as well as a rare example of a neptunium aryloxide complex.

Similar ligand-metal bonding for transition metals and actinides? 5f1 U(C7H7)2-versus 3d n metallocenes

U(C7H7)2- is a fascinating 5f1 complex whose metal-ligand bonding was assigned in the literature as being very similar to 3d7 cobaltocene, based on a crystal-field theoretical interpretation of the experimental magnetic resonance data. The present work provides an in-depth theoretical study of the electronic structure, bonding, and magnetic properties of the 5f1 U(C7H7)2-vs. 3d metallocenes with V, Co, and Ni, performed with relativistic wavefunction and density functional methods. The ligand to metal donation bonding in U(C7H7)2- is strong and in fact similar to that in vanadocene, in the sense that the highest occupied arene orbitals donate electron density into empty metal orbitals of the same symmetry with respect to the rotational axis (3dπ for V, 5fδ for U), but selectively with α spin (↑). For Co and Ni, the dative bonding from the ligands is β spin (↓) selective into partially filled 3dπ orbitals. In all systems, this spin delocalization triggers spin polarization in the arene σ bonding framework, causing proton spin densities opposite to those of the carbons. As a consequence, the proton spin densities and hyperfine coupling constants are negative for the Co and Ni complex, but positive for vanadocene. The of U(C7H7)2- is negative and similar to that of cobaltocene, but only because of the strong spin-orbit coupling in the actinocene, which causes to be opposite to the sign of the proton spin density. The study contributes to a better understanding of actinide 5f vs. transition metal 3d covalency, and highlights potential pitfalls when interpreting experimental magnetic resonance data in terms of covalent bonding for actinide complexes.

Ab initio calculations of heavy-actinide hexahalide compounds: do these heavy actinides behave like their isoelectronic lanthanide analogues?

Research on heavy actinides has experienced an increased interest in the last few years due to new synthetic techniques and recent technological advances that have allowed for obtaining important information even from very small samples. This area presents challenges not only from the experimental point of view but also from the theoretical perspective. This work deals with a multiconfigurational CASSCF and NEVPT2 benchmark study based on a two-step methodology that considers first correlation effects and then the spin-orbit coupling applied to berkelium (Bk), californium (Cf), einsteinium (Es) and fermium (Fm) hexahalides. Optical properties, such as f → d transitions and crystal-field parameters, have been calculated and rationalized. The results for these trivalent actinides indicate that the electronic structure of the low-lying states is reproduced accurately with small basis sets. The ground-state multiplets are isolated, in the same manner as their isoelectronic lanthanide counterparts. In the case of tetravalent berkelium, the picture is different regarding the electronic structure where crystal-field theory fails due to considerable ligand-to-metal charge transfer contributions to the ground state.

Electronic Structure and Properties of Berkelium Iodates

The reaction of ²⁴⁹Bk(OH)₄ with iodate under hydrothermal conditions results in the formation of Bk(IO₃)₃ as the major product with trace amounts of Bk(IO₃)₄ also crystallizing from the reaction mixture. The structure of Bk(IO₃)₃ consists of nine-coordinate Bk^III cations that are bridged by iodate anions to yield layers that are isomorphous with those found for Am^III, Cf^III, and with lanthanides that possess similar ionic radii. Bk(IO₃)₄ was expected to adopt the same structure as M(IO₃)₄ (M = Ce, Np, Pu), but instead parallels the structural chemistry of the smaller Zr^IV cation. Bk^III-O and Bk^IV-O bond lengths are shorter than anticipated and provide further support for a postcurium break in the actinide series. Photoluminescence and absorption spectra collected from single crystals of Bk(IO₃)₄ show evidence for doping with Bk^III in these crystals. In addition to luminescence from Bk^III in the Bk(IO₃)₄ crystals, a broad-band absorption feature is initially present that is similar to features observed in systems with intervalence charge transfer. However, the high-specific activity of ²⁴⁹Bk (t_1/2 = 320 d) causes oxidation of Bk^III and only Bk^IV is present after a few days with concomitant loss of both the Bk^III luminescence and the broadband feature. The electronic structure of Bk(IO₃)₃ and Bk(IO₃)₄ were examined using a range of computational methods that include density functional theory both on clusters and on periodic structures, relativistic ab initio wave function calculations that incorporate spin-orbit coupling (CASSCF), and by a full-model Hamiltonian with spin-orbit coupling and Slater-Condon parameters (CONDON). Some of these methods provide evidence for an asymmetric ground state present in Bk^IV that does not strictly adhere to Russel-Saunders coupling and Hund's Rule even though it possesses a half-filled 5f ⁷ shell. Multiple factors contribute to the asymmetry that include 5f electrons being present in microstates that are not solely spin up, spin-orbit coupling induced mixing of low-lying excited states with the ground state, and covalency in the Bk^IV-O bonds that distributes the 5f electrons onto the ligands. These factors are absent or diminished in other f⁷ ions such as Gd^III or Cm^III.

Uranyl Carbonate Complexes in Aqueous Solution and Their Ligand NMR Chemical Shifts and 17O Quadrupolar Relaxation Studied by ab Initio Molecular Dynamics

Dynamic structural effects, NMR ligand chemical shifts, and ¹⁷O NMR quadrupolar relaxation rates are investigated in the series of complexes UO₂²⁺, UO₂(CO₃)₃^4-, and (UO₂)₃(CO₃)₆^6-. Car-Parrinello molecular dynamics (CPMD) is used to simulate the dynamics of the complexes in water. NMR properties are computed on clusters extracted from the CPMD trajectories. In the UO₂²⁺ complex, coordination at the uranium center by water molecules causes a decrease of around 300 ppm for the uranyl ¹⁷O chemical shift. The final value of this chemical shift is within 40 ppm of the experimental range. The UO₂(CO₃)₃^4- and (UO₂)₃(CO₃)₆^6- complexes show a solvent dependence of the terminal carbonate ¹⁷O and ¹³C chemical shifts that is less pronounced than that for the uranyl oxygen atom. Corrections to the chemical shift from hybrid functionals and spin-orbit coupling improve the accuracy of chemical shifts if the sensitivity of the uranyl chemical shift to the uranyl bond length (estimated at 140 ppm per 0.1 Å from trajectory data) is taken into consideration. The experimentally reported trend in the two unique ¹³C chemical shifts is correctly reproduced for (UO₂)₃(CO₃)₆^6-. NMR relaxation rate data support large ¹⁷O peak widths, but remain below those noted in the experimental literature. Comparison of relaxation data for solvent-including versus solvent-free models suggest that carbonate ligand motion overshadows explicit solvent effects.

Electronic g Tensors in UV Complexes-A Computational Study

The theory and computation of EPR parameters from first principles has seen a great deal of development over the past two decades. In particular, various techniques for the computation of the electronic g tensor have been implemented in many quantum chemistry packages. These methods have been successfully applied to paramagnetic organic species and transition metal systems. The situation is less well-understood and established in the case of actinide-containing molecules and there is a dearth of experimental data available for validation of computations. In this study quantum chemical techniques have been used to evaluate the g tensor for U^V complexes, for which experimental data are available for comparison. The g tensors were calculated using relatively simple, state-averaged complete active space self-consistent field (SA-CASSCF) calculations. This approach is shown to be capable of providing useful accuracy. Aspects of the computations that should be refined to provide a more quantitative approach are discussed. The key features of the underlying electronic structure that influence the computed g values are delineated, providing a simple physical picture of these subtle molecular properties.