Hyperion: A New Computational Tool for Relativistic Ab Initio Hyperfine Coupling

Tris-Silanide f-Block Complexes: Insights into Paramagnetic Influence on NMR Chemical Shifts

The paramagnetism of f-block ions has been exploited in chiral shift reagents and magnetic resonance imaging, but these applications tend to focus on 1H NMR shifts as paramagnetic broadening makes less sensitive nuclei more difficult to study. Here we report a solution and solid-state (ss) 29Si NMR study of an isostructural series of locally D 3h -symmetric early f-block metal(III) tris-hypersilanide complexes, [M{Si(SiMe3)3}3(THF)2] (1-M; M = La, Ce, Pr, Nd, U); 1-M were also characterized by single crystal and powder X-ray diffraction, EPR, ATR-IR, and UV-vis-NIR spectroscopies, SQUID magnetometry, and elemental analysis. Only one SiMe3 signal was observed in the 29Si ssNMR spectra of 1-M, while two SiMe3 signals were seen in solution 29Si NMR spectra of 1-La and 1-Ce. This is attributed to dynamic averaging of the SiMe3 groups in 1-M in the solid state due to free rotation of the M-Si bonds and dissociation of THF from 1-M in solution to give the locally C 3v -symmetric complexes [M{Si(SiMe3)3}3(THF) n ] (n = 0 or 1), which show restricted rotation of M-Si bonds on the NMR time scale. Density functional theory and complete active space self-consistent field spin-orbit calculations were performed on 1-M and desolvated solution species to model paramagnetic NMR shifts. We find excellent agreement of experimental 29Si NMR data for diamagnetic 1-La, suggesting n = 1 in solution and reasonable agreement of calculated paramagnetic shifts of SiMe3 groups for 1-M (M = Pr and Nd); the NMR shifts for metal-bound 29Si nuclei could only be reproduced for diamagnetic 1-La, showing the current limitations of pNMR calculations for larger nuclei.

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.

Relativistic Douglas-Kroll-Hess calculations of hyperfine interactions within first-principles multireference methods

A relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas-Kroll-Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods, including spin-orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field formalism with restricted active space state interaction for spin-orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativistic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs), such as Tb(III)Pc2 and Tb(II)(CpiPr5)2 without ligand truncation using well-converged basis sets. In particular, for the divalent SMM, which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

Relativistic Quantum Chemical Investigation of Actinide Covalency Measured by Electron Paramagnetic Resonance Spectroscopy

We investigate actinide covalency effects in two [AnCptt3] (An = Th, U) complexes recently studied with pulsed electron paramagnetic resonance spectroscopy, using the Hyperion package to obtain relativistic hyperfine coupling constants from relativistic multiconfigurational wave functions. 1H and 13C HYSCORE simulations using the computed parameters show excellent agreement with the experimental data, highlighting the accuracy of modern relativistic ab initio methods. The extent of covalency indicated from the calculations on [ThCptt3] is in agreement with the original report based on traditional spectral fitting methods, while the covalency in [UCptt3] is found to be previously overestimated. The latter is due to the paramagnetic spin-orbit effect that arises naturally in a relativistic theory of hyperfine coupling and yet was not accounted for in the original study, thus highlighting the necessity of relativistic approaches for the interpretation of magnetic resonance data pertaining to actinides.

Unified construction of relativistic Hamiltonians

It is shown that the four-component (4C), quasi-four-component (Q4C), and exact two-component (X2C) relativistic Hartree-Fock equations can be implemented in a unified manner by making use of the atomic nature of the small components of molecular 4-spinors. A model density matrix approximation can first be invoked for the small-component charge/current density functions, which gives rise to a static, pre-molecular mean field to be combined with the one-electron term. As a result, only the nonrelativistic-like two-electron term of the 4C/Q4C/X2C Fock matrix needs to be updated during the iterations. A "one-center small-component" approximation can then be invoked in the evaluation of relativistic integrals, that is, all atom-centered small-component basis functions are regarded as extremely localized near the position of the atom to which they belong such that they have vanishing overlaps with all small- or large-component functions centered at other nuclei. Under these approximations, the 4C, Q4C, and X2C mean-field and many-electron Hamiltonians share precisely the same structure and accuracy. Beyond these is the effective quantum electrodynamics Hamiltonian that can be constructed in the same way. Such approximations lead to errors that are orders of magnitude smaller than other sources of errors (e.g., truncation errors in the one- and many-particle bases as well as uncertainties of experimental measurements) and are, hence, safe to use for whatever purposes. The quaternion forms of the 4C, Q4C, and X2C equations are also presented in the most general way, based on which the corresponding Kramers-restricted open-shell variants are formulated for "high-spin" open-shell systems.

Metal-carbon bonding in early lanthanide substituted cyclopentadienyl complexes probed by pulsed EPR spectroscopy

We examine lanthanide (Ln)-ligand bonding in a family of early Ln3+ complexes [Ln(Cptt)3] (1-Ln, Ln = La, Ce, Nd, Sm; Cptt = C5H3tBu2-1,3) by pulsed electron paramagnetic resonance (EPR) methods, and provide the first characterization of 1-La and 1-Nd by single crystal XRD, multinuclear NMR, IR and UV/Vis/NIR spectroscopy. We measure electron spin T1 and Tm relaxation times of 12 and 0.2 μs (1-Nd), 89 and 1 μs (1-Ce) and 150 and 1.7 μs (1-Sm), respectively, at 5 K: the T1 relaxation of 1-Nd is more than 102 times faster than its valence isoelectronic uranium analogue. 13C and 1H hyperfine sublevel correlation (HYSCORE) spectroscopy reveals that the extent of covalency is negligible in these Ln compounds, with much smaller hyperfine interactions than observed for equivalent actinide (Th and U) complexes. This is corroborated by ab initio calculations, confirming the predominant electrostatic nature of the metal-ligand bonding in these complexes.

Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods

An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. 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.

Dissecting the Molecular Origin of g-Tensor Heterogeneity and Strain in Nitroxide Radicals in Water: Electron Paramagnetic Resonance Experiment versus Theory

Nitroxides are common EPR sensors of microenvironmental properties such as polarity, numbers of H-bonds, pH, and so forth. Their solvation in an aqueous environment is facilitated by their high propensity to form H-bonds with the surrounding water molecules. Their g- and A-tensor elements are key parameters to extracting the properties of their microenvironment. In particular, the gxx value of nitroxides is rich in information. It is known to be characterized by discrete values representing nitroxide populations previously assigned to have different H-bonds with the surrounding waters. Additionally, there is a large g-strain, that is, a broadening of g-values associated with it, which is generally correlated with environmental and structural micro-heterogeneities. The g-strain is responsible for the frequency dependence of the apparent line width of the EPR spectra, which becomes evident at high field/frequency. Here, we address the molecular origin of the gxx heterogeneity and of the g-strain of a nitroxide moiety (HMI: 2,2,3,4,5,5-hexamethylimidazolidin-1-oxyl, C9H19N2O) in water. To treat the solvation effect on the g-strain, we combined a multi-frequency experimental approach with ab initio molecular dynamics simulations for structural sampling and quantum chemical EPR property calculations at the highest realistically affordable level, including an explicitly micro-solvated HMI ensemble and the embedded cluster reference interaction site model. We could clearly identify the distinct populations of the H-bonded nitroxides responsible for the gxx heterogeneity experimentally observed, and we dissected the role of the solvation shell, H-bond formation, and structural deformation of the nitroxide in the creation of the g-strain associated with each nitroxide subensemble. Two contributions to the g-strain were identified in this study. The first contribution depends on the number of hydrogen bonds formed between the nitroxide and the solvent because this has a large and well-understood effect on the gxx-shift. This contribution can only be resolved at high resonance frequencies, where it leads to distinct peaks in the gxx region. The second contribution arises from configurational fluctuations of the nitroxide that necessarily lead to g-shift heterogeneity. These contributions cannot be resolved experimentally as distinct resonances but add to the line broadening. They can be quantitatively analyzed by studying the apparent line width as a function of microwave frequency. Interestingly, both theory and experiment confirm that this contribution is independent of the number of H-bonds. Perhaps even more surprisingly, the theoretical analysis suggests that the configurational fluctuation broadening is not induced by the solvent but is inherently present even in the gas phase. Moreover, the calculations predict that this broadening decreases upon solvation of the nitroxide.