Magnetic Resonance Properties of Actinyl Carbonate Complexes and Plutonyl(VI)-tris-nitrate

Quantifying Covalency and Environmental Effects in RASSCF-Simulated O K-Edge XANES of Uranyl

A RASSCF approach to simulate the O K-edge XANES spectra of uranyl is employed, utilizing three models that progressively improve the representation of the local crystal environment. Simulations successfully reproduce the observed three-peak profile of the experimental spectrum and confirm peak assignments made by Denning. The [UO2Cl4]2- model offers the best agreement with experiment, with peak positions (to within 1 eV) and relative peak separations accurately reproduced. Establishing a direct link between a specific electronic transition and peak intensity is complicated, as a large number of possible transitions can contribute to the overall peak profile. Furthermore, a relationship between oxygen character in the antibonding orbital and the strength of the transition breaks down when using a variety of orbital composition approaches at larger excitation energy. Covalency analysis of the U-O bond in both the ground- and excited-state reveals a dependence on the crystal environment. Orbital composition analysis reveals an underestimation of the uranium contribution to ground-state bonding orbitals when probing O K-edge core-excited states, regardless of the uranyl model employed. However, improving the environmental model provides core-excited state electronic structures that are better representative of that of the ground-state, validating their use in the determination of covalency and bonding.

Paramagnetic Properties of [AnIV(NO3)6]2- Complexes (An = U, Np, Pu) Probed by NMR Spectroscopy and Quantum Chemical Calculations

Actinide +IV complexes with six nitrates [AnIV(NO3)6]2- (An = Th, U, Np, and Pu) have been studied by 15N and 17O NMR spectroscopy in solution and first-principles calculations. Magnetic susceptibilities were evaluated experimentally using the Evans method and are in good agreement with the ab initio values. The evolution in the series of the crystal field parameters deduced from ab initio calculations is discussed. The NMR paramagnetic shifts are analyzed based on ab initio calculations. Because the cubic symmetry of the complex quenches the dipolar contribution, they are only of Fermi contact origin. They are evaluated from first-principles based on a complete active space/density functional theory (DFT) strategy, in good accordance with the experimental one. The ligand hyperfine coupling constants are deduced from paramagnetic shifts and calculated using unrestricted DFT. The latter are decomposed in terms of the contribution of molecular orbitals. It highlights two pathways for the delocalization of the spin density from the metallic open-shell 5f orbitals to the NMR active nuclei, either through the valence 5f hybridized with 6d to the valence 2p molecular orbitals of the ligands, or by spin polarization of the metallic 6p orbitals which interact with the 2s-based molecular orbitals of the ligands.

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.

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.

Bounding [AnO2]2+ (An = U, Np) covalency by simulated O K-edge and An M-edge X-ray absorption near-edge spectroscopy

Restricted active space simulations are shown to accurately reproduce and characterise both O K-edge and U M4,5-edge spectra of uranyl in excellent agreement with experimental peak positions and are extended to the Np analogue. Analysis of bonding orbital composition in the ground and O K-edge core-excited states demonstrates that metal contribution is underestimated in the latter. In contrast, An M4/5-edge core-excited states produce bonding orbital compositions significantly more representative of those in the ground state. Quantum Theory of Atoms in Molecules analysis is employed to explain the discrepancy between K- and M-edge data and demonstrates that the location of the core-hole impacts the pattern of electron localisation in core-excited states. An apparent contradiction to this behaviour in neptunyl is rationalised in terms interelectronic repulsion between the unpaired 5f electron and the excited core-electron.

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.

Quadratic Spin-Orbit Mechanism of the Electronic g-Tensor

Understanding how the electronic g-tensor is linked to the electronic structure is desirable for the correct interpretation of electron paramagnetic resonance spectra. For heavy-element compounds with large spin-orbit (SO) effects, this is still not completely clear. We report our investigation of quadratic SO contributions to the g-shift in heavy transition metal complexes. We implemented third-order perturbation theory in order to analyze the contributions arising from frontier molecular spin orbitals (MSOs). We show that the dominant quadratic SO term─spin-Zeeman (SO²/SZ)─generally makes a negative contribution to the g-shift, irrespective of the particular electronic configuration or molecular symmetry. We further analyze how the SO²/SZ contribution adds to or subtracts from the linear orbital-Zeeman (SO/OZ) contribution to the individual principal components of the g-tensor. Our study suggests that the SO²/SZ mechanism decreases the anisotropy of the g-tensor in early transition metal complexes and increases it in late transition metal complexes. Finally, we apply MSO analysis to the investigation of g-tensor trends in a set of closely related Ir and Rh pincer complexes and evaluate the influence of different chemical factors (the nuclear charge of the central atom and the terminal ligand) on the magnitudes of the g-shifts. We expect our conclusions to aid the understanding of spectra in magnetic resonance investigations of heavy transition metal compounds.

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.

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.

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.

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.

Synthesis, Structural, and Electronic Properties of K4PuVIO2(CO3)3(cr): An Environmentally Relevant Plutonium Carbonate Complex

The chemical properties of actinide materials are often predefined and described based on the data available for isostructural species. This is the case for potassium plutonyl (Pu^VI) carbonate, K₄Pu^VIO₂(CO₃)_3(cr), a complex relevant for nuclear technology and the environment, of which the crystallographic and thermodynamic properties of which are still lacking. We report here the synthesis and characterization of Pu^VI achieved by single-crystal X-ray diffraction analysis and high-energy-resolution fluorescence-detected X-ray absorption near-edge structure at the Pu M₄-edge coupled with electronic structure calculations. The crystallographic properties of Pu^VI are compared with isostructural uranium (U) and neptunium (Np) compounds. Actinyl (An^VI) axial bond lengths, [O-An^VI-O]²⁺, are correlated between solid, K₄An^VIO₂(CO₃)_3(cr), and aqueous, [An^VIO₂(CO₃)₃]^4-_(aq) species for the U^VI-Np^VI-Pu^VI series. The spectroscopic data are compared to KPu^VO₂CO_3(cr) and Pu^IVO_2(cr) to tackle the trend in the electronic structure of Pu^VI regarding the oxidation state changes and local structural modifications around the Pu atom.

31P MAS NMR and DFT study of crystalline phosphate matrices

We present the study of the phosphorus local environment by using ³¹P MAS NMR in a series of seven double monophosphates M^IIM^IV(PO₄)₂ (M^II and M^IV being divalent and tetravalent cations, respectively) of yavapaiite and low-yavapaiite type crystal structures. Solid-state and cluster DFT calculations were found to be efficient for predicting the ³¹P isotropic chemical shift and chemical shift anisotropy. To achieve this performance, however, a proper computational optimisation of the experimental structural data was required. From the three optimisation methods tested, the full optimisation provided the best reference structure for the calculation of the NMR parameters of the studied phosphates. Also, a better prediction of the chemical shifts was possible by using a correction to the GIPAW calculated shielding.

Electronic Structure and Chemical Bonding of [AmO2(H2O) n ]2+/1

Systematic americyl-hydration cations were investigated theoretically to understand the electronic structures and bonding in [(AmO2)(H2O) n ]2+/1+ (n = 1-6). We obtained the binding energy using density functional theory methods with scalar relativistic and spin-orbit coupling effects. The geometric structures of these species have been investigated in aqueous solution via an implicit solvation model. Computational results reveal that the complexes of five equatorial water molecules coordinated to americyl ions are the most stable due to the enhanced ionic interactions between the AmO2 2+/1+ cation and multiple oxygen atoms as electron donors. As expected, Am-Owater bonds in such series are electrostatic in nature and contain a generally decreasing covalent character when hydration number increases.

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.

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.

NMR analysis of an Fe(i)–carbene complex with strong magnetic anisotropy

A paramagnetic, easy-plane anisotropic Fe^I complex, bearing cyclic-alkyl(amino) carbene (cAAC) ligands, is studied by means of NMR and DFT.

Modeling the electronic states and magnetic properties derived from the f1 configuration in lanthanocene and actinocene compounds

The electronic structure and magnetic properties of a series of Kramers ions with f¹ configuration in axial symmetry have been analyzed with a combination of theoretical methods: ab initio relativistic wavefunction methods as well as a crystal-field (CF) model with parameters extracted from the ab initio calculations.

A Nonorthogonal State-Interaction Approach for Matrix Product State Wave Functions

We present a state-interaction approach for matrix product state (MPS) wave functions in a nonorthogonal molecular orbital basis. Our approach allows us to calculate, for example, transition and spin-orbit coupling matrix elements between arbitrary electronic states, provided that they share the same one-electron basis functions and size of the active orbital space, respectively. The key element is the transformation of the MPS wave functions of different states from a nonorthogonal to a biorthonormal molecular orbital basis representation, by exploiting a sequence of nonunitary transformations, following a proposal by Malmqvist [Int. J. Quantum Chem. 1986, 30, 479]. This is well-known for traditional wave function parametrizations but has not yet been exploited for MPS wave functions.

Kohn–Sham calculations of NMR shifts for paramagnetic 3d metal complexes: protocols, delocalization error, and the curious amide proton shifts of a high-spin iron(ii) macrocycle complex

Ligand chemical shifts (pNMR shifts) are analyzed using DFT. A large difference in the amide proton shifts of a high-spin Fe(ii) complex arises from O → Fe dative bonding which only transfers β spin density to the metal.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Theoretical study of the electronic spectra of neutral and cationic NpO and NpO2

The electronic spectra of neutral NpO and NpO2 as well as of their mono- (NpO(+), NpO2(+)) and dications (NpO(2+), NpO2(2+)) were studied using multiconfigurational relativistic quantum chemical calculations at the complete active space self-consistent field/CASPT2 level of theory taking into account spin-orbit coupling. The active space included 16 orbitals: all the 7s, 6d, and 5f orbitals of neptunium together with selected orbitals of oxygen. The vertical excitation energies on the ground state geometries have been computed up to ca. 35,000 cm(-1). The gas-phase electronic spectra were evaluated on the basis of the computed Einstein coefficients at 298 K and 3000 K. The computed vertical transition energies show good agreement with previous condensed-phase results on NpO2(+) and NpO2(2+).

Temperature dependence of contact and dipolar NMR chemical shifts in paramagnetic molecules

Using a recently proposed equation for NMR nuclear magnetic shielding for molecules with unpaired electrons [A. Soncini and W. Van den Heuvel, J. Chem. Phys. 138, 021103 (2013)], equations for the temperature (T) dependent isotropic shielding for multiplets with an effective spin S equal to 1/2, 1, 3/2, 2, and 5/2 in terms of electron paramagnetic resonance spin Hamiltonian parameters are derived and then expanded in powers of 1/T. One simplifying assumption used is that a matrix derived from the zero-field splitting (ZFS) tensor and the Zeeman coupling matrix (g-tensor) share the same principal axis system. The influence of the rhombic ZFS parameter E is only investigated for S = 1. Expressions for paramagnetic contact shielding (from the isotropic part of the hyperfine coupling matrix) and pseudo-contact or dipolar shielding (from the anisotropic part of the hyperfine coupling matrix) are considered separately. The leading order is always 1/T. A temperature dependence of the contact shielding as 1/T and of the dipolar shielding as 1/T(2), which is sometimes assumed in the assignment of paramagnetic chemical shifts, is shown to arise only if S ≥ 1 and zero-field splitting is appreciable, and only if the Zeeman coupling matrix is nearly isotropic (Δg = 0). In such situations, an assignment of contact versus dipolar shifts may be possible based only on linear and quadratic fits of measured variable-temperature chemical shifts versus 1/T. Numerical data are provided for nickelocene (S = 1). Even under the assumption of Δg = 0, a different leading order of contact and dipolar shifts in powers of 1/T is not obtained for S = 3/2. When Δg is not very small, dipolar and contact shifts both depend in leading order in 1/T in all cases, with sizable contributions in order 1/T(n) with n = 2 and higher.

Calculating NMR Chemical Shifts for Paramagnetic Metal Complexes from First-Principles

Paramagnetic effects on NMR shifts (pNMR) for paramagnetic metal complexes are calculated from first-principles, without recourse to spin Hamiltonian parameters. A newly developed code based on complete active space (CAS) and restricted active space (RAS) techniques in conjunction with treating spin-orbit (SO) coupling via state interaction is applied to (13)C NMR shifts of actinyl tris-carbonate complexes, specifically [UO2(CO3)3](5-) and [NpO2(CO3)3](4-). The experimental pNMR shifts as well as the sizable difference of the (13)C NMR shift for these iso-electronic species are well reproduced by the calculations. Approximations to the pNMR shift equations using spin Hamiltonian parameters or the magnetic susceptibility are calculated for the same systems at the same level of theory, and it is shown how the approximations relate to the ab initio data.

Single-ion 4f element magnetism: an ab-initio look at Ln(COT)2−

The orbital and spin magnetizations of the ground state of the Ln(COT)₂⁻ series are investigated at the ab-initio level.