Calculation of Hyperfine Tensors and Paramagnetic NMR Shifts Using the Relativistic Zeroth-Order Regular Approximation and Density Functional Theory

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

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.

Zero-field splitting parameters within exact two-component theory and modern density functional theory using seminumerical integration

An efficient implementation of zero-field splitting parameters based on the work of Schmitt et al. [J. Chem. Phys. 134, 194113 (2011)] is presented. Seminumerical integration techniques are used for the two-electron spin-dipole contribution and the response equations of the spin-orbit perturbation. The original formulation is further generalized. First, it is extended to meta-generalized gradient approximations and local hybrid functionals. For these functional classes, the response of the paramagnetic current density is considered in the coupled-perturbed Kohn-Sham equations for the spin-orbit perturbation term. Second, the spin-orbit perturbation is formulated within relativistic exact two-component theory and the screened nuclear spin-orbit (SNSO) approximation. The accuracy of the implementation is demonstrated for transition-metal and diatomic main-group compounds. The efficiency is assessed for Mn and Mo complexes. Here, it is found that coarse integration grids for the seminumerical schemes lead to drastic speedups while introducing clearly negligible errors. In addition, the SNSO approximation substantially reduces the computational demands and leads to very similar results as the spin-orbit mean field Ansatz.

Lewis Acid Probe for Basicity of Sulfide Electrolytes Investigated by 11B Solid-State NMR

Sulfide-based solid-state lithium-ion batteries (SSLIB) have attracted a lot of interest globally in the past few years for their high safety and high energy density over the traditional lithium-ion batteries. However, sulfide electrolytes (SEs) are moisture-sensitive which pose significant challenges in the material preparation and cell manufacturing. To the best of our knowledge, there is no tool available to probe the types and the strength of the basic sites in sulfide electrolytes, which is crucial for understanding the moisture stability of sulfide electrolytes. Herein, we propose a new spectral probe with the Lewis base indicator BBr3 to probe the strength of Lewis basic sites on various sulfide electrolytes by 11B solid-state NMR spectroscopy (11B-NMR). The active sulfur sites and the corresponding strength of the sulfide electrolytes are successfully evaluated by the proposed Lewis base probe. The probed strength of the active sulfur sites of a sulfide electrolyte is consistent with the results of DFT (density functional theory) calculation and correlated with the H2S generation rate when the electrolyte was exposed in moisture atmosphere. This work paves a new way to investigate the basicity and moisture stability of the sulfide electrolytes.

Nature of NMR Shifts in Paramagnetic Octahedral Ru(III) Complexes with Axial Pyridine-Based Ligands

In recent decades, transition-metal coordination compounds have been extensively studied for their antitumor and antimetastatic activities. In this work, we synthesized a set of symmetric and asymmetric Ru(III) and Rh(III) coordination compounds of the general structure (Na⁺/K⁺/PPh₄⁺/LH⁺) [trans-M^IIIL(eq)_nL(ax)₂]^- (M = Ru^III or Rh^III; L(eq) = Cl, n = 4; L(eq) = ox, n = 2; L(ax) = 4-R-pyridine, R = CH₃, H, C₆H₅, COOH, CF₃, CN; L(ax) = DMSO-S) and systematically investigated their structure, stability, and NMR properties. ¹H and ¹³C NMR spectra measured at various temperatures were used to break down the total NMR shifts into the orbital (temperature-independent) and hyperfine (temperature-dependent) contributions. The hyperfine NMR shifts for paramagnetic Ru(III) compounds were analyzed in detail using relativistic density functional theory (DFT). The effects of (i) the 4-R substituent of pyridine, (ii) the axial trans ligand L(ax), and (iii) the equatorial ligands L(eq) on the distribution of spin density reflected in the "through-bond" (contact) and the "through-space" (pseudocontact) contributions to the hyperfine NMR shifts of the individual atoms of the pyridine ligands are rationalized. Further, we demonstrate the large effects of the solvent on the hyperfine NMR shifts and discuss our observations in the general context of the paramagnetic NMR spectroscopy of transition-metal complexes.

Mononuclear Oxidovanadium(IV) Complexes with BIAN Ligands: Synthesis and Catalytic Activity in the Oxidation of Hydrocarbons and Alcohols with Peroxides

Reactions of VCl3 with 1,2-Bis[(4-methylphenyl)imino]acenaphthene (4-Me-C6H4-bian) or 1,2-Bis[(2-methylphenyl)imino]acenaphthene (2-Me-C6H4-bian) in air lead to the formation of [VOCl2(R-bian)(H2O)] (R = 4-Me-C6H4 (1), 2-Me-C6H4 (2)). Thes complexes were characterized by IR and EPR spectroscopy as well as elemental analysis. Complexes 1 and 2 have high catalytic activity in the oxidation of hydrocarbons with hydrogen peroxide and alcohols with tert-butyl hydroperoxide in acetonitrile at 50 °С. The product yields are up to 40% for cyclohexane. Of particular importance is the addition of 2-pyrazinecarboxylic acid (PCA) as a co-catalyst. Oxidation proceeds mainly with the participation of free hydroxyl radicals, as evidenced by taking into account the regio- and bond-selectivity in the oxidation of n-heptane and methylcyclohexane, as well as the dependence of the reaction rate on the initial concentration of cyclohexane.

The evolution of paramagnetic NMR as a tool in structural biology

Paramagnetic NMR data contain extremely accurate long-range information on metalloprotein structures and, when used in the frame of integrative structural biology approaches, they allow for the retrieval of structural details to a resolution that is not achievable using other techniques. Paramagnetic data thus represent an extremely powerful tool to refine protein models in solution, especially when coupled to X-ray or cryoelectron microscopy data, to monitor the formation of complexes and determine the relative arrangements of their components, and to highlight the presence of conformational heterogeneity. More recently, theoretical and computational advancements in quantum chemical calculations of paramagnetic NMR observables are progressively opening new routes in structural biology, because they allow for the determination of the structure within the coordination sphere of the metal center, thus acting as a loupe on sites that are difficult to observe but very important for protein function.

Paramagnetic NMR Shielding Tensors Based on Scalar Exact Two-Component and Spin-Orbit Perturbation Theory

The temperature-dependent Fermi-contact and pseudocontact terms are important contributions to the paramagnetic NMR shielding tensor. Herein, we augment the scalar-relativistic (local) exact two-component (X2C) framework with spin-orbit perturbation theory including the screened nuclear spin-orbit correction for the EPR hyperfine coupling and g tensor to compute these temperature-dependent terms. The accuracy of this perturbative ansatz is assessed with the self-consistent spin-orbit two-component and four-component treatments serving as reference. This shows that the Fermi-contact and pseudocontact interaction is sufficiently described for paramagnetic NMR shifts; however, larger deviations are found for the EPR spectra and the principle components of the EPR properties of heavy elements. The impact of the perturbative treatment is further compared to that of the density functional approximation and the basis set. Large-scale calculations are routinely possible with the multipole-accelerated resolution of the identity approximation and the seminumerical exchange approximation, as shown for [CeTi₆O₃(OiPr)₉(salicylate)₆].

Electron correlation and vibrational effects in predictions of paramagnetic NMR shifts

Electronic structure calculations are fundamentally important for the interpretation of nuclear magnetic resonance (NMR) spectra from paramagnetic systems that include organometallic and inorganic compounds, catalysts, or metal-binding sites in proteins. Prediction of induced paramagnetic NMR shifts requires knowledge of electron paramagnetic resonance (EPR) parameters: the electronic g tensor, zero-field splitting D tensor, and hyperfine A tensor. The isotropic part of A, called the hyperfine coupling constant (HFCC), is one of the most troublesome properties for quantum chemistry calculations. Yet, even relatively small errors in calculations of HFCC tend to propagate into large errors in the predicted NMR shifts. The poor quality of A tensors that are currently calculated using density functional theory (DFT) constitutes a bottleneck in improving the reliability of interpretation of the NMR spectra from paramagnetic systems. In this work, electron correlation effects in calculations of HFCCs with a hierarchy of ab initio methods were assessed, and the applicability of different levels of DFT approximations and the coupled cluster singles and doubles (CCSD) method was tested. These assessments were performed for the set of selected test systems comprising an organic radical, and complexes with transition metal and rare-earth ions, for which experimental data are available. Severe deficiencies of DFT were revealed but the CCSD method was able to deliver good agreement with experimental data for all systems considered, however, at substantial computational costs. We proposed a more computationally tractable alternative, where the A was computed with the coupled cluster theory exploiting locality of electron correlation. This alternative is based on the domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) method. In this way the robustness and reliability of the coupled cluster theory were incorporated into the modern formalism for the prediction of induced paramagnetic NMR shifts, and became applicable to systems of chemical interest. This approach was verified for the bis(cyclopentadienyl)vanadium(II) complex (Cp₂V; vanadocene), and the metal-binding site of the Zn²⁺ → Co²⁺ substituted superoxide dismutase (SOD) metalloprotein. Excellent agreement with experimental NMR shifts was achieved, which represented a substantial improvement over previous theoretical attempts. The effects of vibrational corrections to orbital shielding and hyperfine tensor were evaluated and discussed within the second-order vibrational perturbation theory (VPT2) framework.

Discovery of a Dysprosium Metallocene Single-Molecule Magnet with Two High-Temperature Orbach Processes

Magnetic bistability in single-molecule magnets (SMMs) is a potential basis for new types of nanoscale information storage material. The standard model for thermally activated relaxation of the magnetization in SMMs is based on the occurrence of a single Orbach process. Here, we show that incorporating a phosphorus atom into the framework of the dysprosium metallocene [(Cp^iPr5)Dy(Cp^PEt4)]⁺[B(C₆F₅)₄]⁻ (Cp^iPr5 is penta-isopropylcyclopentadienyl, Cp^PEt4 is tetraethylphospholyl) leads to the occurrence of two distinct high-temperature Orbach processes, with energy barriers of 1410(10) and 747(7) cm^–1, respectively. These barriers provide experimental evidence for two different spin–phonon coupling regimes, which we explain with the aid of ab initio calculations. The strong and highly axial crystal field in this SMM also allows magnetic hysteresis to be observed up to 70 K, using a scan rate of 25 Oe s^–1. In characterizing this SMM, we show that a conventional Debye model and consideration of rotational contributions to the spin–phonon interaction are insufficient to explain the observed phenomena.

A phospholyl-ligated dysprosium metallocene single-molecule magnet shows two high-temperature Orbach relaxation processes with effective energy barriers of 1410(10) and 747(7) cm⁻¹, and magnetic hysteresis up to 70 K when using a scan rate of 25 Oe s⁻¹.

Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance

We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C₅Me₅)₂Y(μ-S)₂Mo(μ-S)₂Y(C₅Me₅)₂]^-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C₆Cl₅)₄]^- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.

p-Block metal complexes with bis(pyrazol-1-yl)acetato ligands

The hetereoscorpionate ligands bis(pyrazol-1-yl)acetic acid (Hbpza) and bis(3,5-dimethylpyrazol-1-yl)acetic acid (Hbdmpza) are reacted with [Sn(OAc)₂] or [Sn(acac)₂] to yield the corresponding Sn(II) complexes. A single crystal X-ray determination for the in solution monomeric complex [Sn(bpza)₂] (1a) revealed a dinuclear molecular structure [Sn₂(bpza)₄] (1b), with κ³-N,N,O-coordinated bpza ligands at each of the Sn(II) and two bpza ligands μ-bridging between the Sn(II) centres. The molecular structure of [Sn(bdmpza)₂] (2) exhibits a homoleptic bisligand complex with both ligands displaced by the free electron pair, which is proven by DFT calculations. Oxidation of complex 2 in an attempt to synthesize a homoleptic Sn(IV) complex leads to the formation of [Sn(bdmpza)F₃] (3). In addition, homoleptic bisligand gallium complexes [Ga(bdmpza)₂][ClO₄] (4) and [Ga(bpza)₂][GaCl₄] (5) were prepared and characterized by ⁷¹Ga NMR and IR spectroscopy as well as by X-ray crystallography.

Hyperfine Coupling Constants in Local Exact Two-Component Theory

We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within the one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, that is, the derivatives of the unitary decoupling transformation are considered. This requires the solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened nuclear spin-orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12-17 transition-metal compounds. The error introduced by DLU is negligible, and the DLU-X2C Hamiltonian accurately reproduces its four-component "fully" relativistic parent results. Functionals with a large amount of Hartree-Fock exchange such as CAM-QTP-02 and ωB97X-D are generally favorable. The pure density functional r²SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-ζ bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C₆Cl₅)₄]^- with more than 4700 primitive basis functions and four rare-earth single-molecule magnets: [La(OAr*)₃]^-, [Lu(NR₂)₃]^-, [Lu(OAr*)₃]^-, and [TbPc₂]^-. Here, the results with the spin-orbit DLU-X2C Hamiltonian are in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.

Paramagnetic NMR Shielding Tensors and Ring Currents: Efficient Implementation and Application to Heavy Element Compounds

We present an efficient implementation of paramagnetic NMR shielding tensors and shifts in a nonrelativistic and scalar-relativistic density functional theory framework. For the latter, we make use of the scalar exact two-component Hamiltonian in its local approximation, and generally we apply the well established (multipole-accelerated) resolution of the identity approximation and the seminumerical exchange approximation. The perturbed density matrix of a paramagnetic NMR shielding calculation is further used to study the magnetically induced current density and ring currents of open-shell systems as illustrated for [U@Bi₁₂]^3-. [U@Bi₁₂]^3- features delocalized highest occupied molecular orbitals and sustains a net diatropic ring current of ca. 18 nA/T through the Bi₁₂ torus similar to the all-metal aromatic heavy-element cluster [Th@Bi₁₂]^4-.

Crystal and Substituent Effects on Paramagnetic NMR Shifts in Transition-Metal Complexes

Nuclear magnetic resonance (NMR) spectroscopy of paramagnetic molecules provides detailed information about their molecular and electron-spin structure. The paramagnetic NMR spectrum is a very rich source of information about the hyperfine interaction between the atomic nuclei and the unpaired electron density. The Fermi-contact contribution to ligand hyperfine NMR shifts is particularly informative about the nature of the metal–ligand bonding and the structural arrangements of the ligands coordinated to the metal center. In this account, we provide a detailed experimental and theoretical NMR study of compounds of Cr(III) and Cu(II) coordinated with substituted acetylacetonate (acac) ligands in the solid state. For the first time, we report the experimental observation of extremely paramagnetically deshielded ¹³C NMR resonances for these compounds in the range of 900–1200 ppm. We demonstrate an excellent agreement between the experimental NMR shifts and those calculated using relativistic density-functional theory. Crystal packing is shown to significantly influence the NMR shifts in the solid state, as demonstrated by theoretical calculations of various supramolecular clusters. The resonances are assigned to individual atoms in octahedral Cr(acac)₃ and square-planar Cu(acac)₂ compounds and interpreted by different electron configurations and magnetizations at the central metal atoms resulting in different spin delocalizations and polarizations of the ligand atoms. Further, effects of substituents on the ¹³C NMR resonance of the ipso carbon atom reaching almost 700 ppm for Cr(acac)₃ compounds are interpreted based on the analysis of Fermi-contact hyperfine contributions.

The ligand NMR shifts in paramagnetic acetylacetonato Cr(III) and Cu(II) complexes have been predicted and measured in the solid state and interpreted by relativistic DFT calculations. The effects of the metal atom, ligand, and crystal packing on the spin delocalization and polarization reflected in the Fermi-contact contribution to the hyperfine interaction are rationalized.

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.

A Quantum Chemistry View on Two Archetypical Paramagnetic Pentacoordinate Nickel(II) Complexes Offers a Fresh Look on Their NMR Spectra

Quantum chemical methods for calculating paramagnetic NMR observables are becoming increasingly accessible and are being included in the inorganic chemistry practice. Here, we test the performance of these methods in the prediction of proton hyperfine shifts of two archetypical high-spin pentacoordinate nickel(II) complexes (NiSAL-MeDPT and NiSAL-HDPT), which, for a variety of reasons, turned out to be perfectly suited to challenge the predictions to the finest level of detail. For NiSAL-MeDPT, new NMR experiments yield an assignment that perfectly matches the calculations. The slightly different hyperfine shifts from the two “halves” of the molecules related by a pseudo-C₂ axis, which are experimentally divided into two well-defined spin systems, are also straightforwardly distinguished by the calculations. In the case of NiSAL-HDPT, for which no X-ray structure is available, the quality of the calculations allowed us to refine its structure using as a starting template the structure of NiSAL-MeDPT.

State-of-the-art quantum chemical methods and paramagnetism-tailored NMR experiments provide a deep insight on the relation between the spectra and the electronic structure for two paramagnetic pentacoordinate nickel(II) complexes.

Electron-Nucleus Hyperfine Coupling Calculated from Restricted Active Space Wavefunctions and an Exact Two-Component Hamiltonian

Exact two-component (X2C) relativistic nuclear hyperfine magnetic field operators were incorporated in X2C ab initio wavefunction calculations at the multireference restricted active space (RAS) level for calculations of nuclear hyperfine magnetic properties. Spin-orbit coupling was treated via RAS state interaction (SO-RASSI). The method was tested by calculations of electron-nucleus hyperfine coupling constants. The approach, implemented in the OpenMolcas program, overcomes a major limitation of a previous SO-RASSI implementation for hyperfine coupling that relied on nonrelativistic hyperfine operators [J. Chem. Theor. Comput. 2015, 11, 538-549] and therefore had limited applicability. Results from calculations on systems with light and heavy main group elements, transition metals, lanthanides, and one actinide complex demonstrate reasonably good agreement with experimental data, where available, as long as the active space can generate sufficient spin polarization.

A Sterically Hindered Derivative of 2,1,3-Benzotelluradiazole: A Way to the First Structurally Characterised Monomeric Tellurium-Nitrogen Radical Anion

Interaction of the tetradentate redox-active 6,6'-[1,2-phenylenebis(azanediyl)]bis(2,4-di-tert-butylphenol) (H₄ L) with TeCl₄ leads to neutral diamagnetic compound TeL (1) in high yield. The molecule of 1 has a nearly planar TeN₂ O₂ fragment, which suggests the formulation of 1 as Te^II L^2- , in agreement with the results of DFT calculations and QTAIM and NBO analyses. Reduction of 1 with one equivalent of [CoCp₂ ] leads to quantitative formation of the paramagnetic salt [CoCp₂ ]⁺ [1]^.- , which was characterised by single-crystal XRD. The solution EPR spectrum of [CoCp₂ ]⁺ [1]^.- at room temperature features a quintet due to splitting on two equivalent ¹⁴ N nuclei. Below 150 K it turns into a broad singlet line with two weak satellites due to the splitting on the ¹²⁵ Te nucleus. Two-component relativistic DFT calculations perfectly reproduce the a(¹⁴ N) HFI constants and A_∥ (¹²⁵ Te) value responsible for the low-temperature satellite splitting. Calculations predict that the additional electron in 1^.- is localised mainly on L, while the spin density is delocalised over the whole molecule with significant localisation on the Te atom (≥30 %). All these data suggest that 1^.- can be regarded as the first example of a structurally characterised monomeric tellurium-nitrogen radical anion.

Solution of a Puzzle: High-Level Quantum-Chemical Treatment of Pseudocontact Chemical Shifts Confirms Classic Semiempirical Theory

A recently popularized approach for the calculation of pseudocontact shifts (PCSs) based on first-principles quantum chemistry (QC) leads to different results than the classic "semiempirical" equation involving the susceptibility tensor. Studies that attempted a comparison of theory and experiment led to conflicting conclusions with respect to the preferred theoretical approach. In this Letter, we show that after inclusion of previously neglected terms in the full Hamiltonian, one can deduce the semiempirical equations from a rigorous QC-based treatment. It also turns out that in the long-distance limit, one can approximate the complete A tensor in terms of the g tensor. By means of Kohn-Sham density functional theory calculations, we numerically confirm the long-distance expression for the A tensor and the theoretically predicted scaling behavior of the different terms. Our derivation suggests a computational strategy in which one calculates the susceptibility tensor and inserts it into the classic equation for the PCS.

NWChem: Past, present, and future

Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.

Density Functional Calculations of Electron Paramagnetic Resonance g- and Hyperfine-Coupling Tensors Using the Exact Two-Component (X2C) Transformation and Efficient Approximations to the Two-Electron Spin-Orbit Terms

A two-component quasirelativistic density functional theory implementation of the computation of hyperfine and g-tensors at exact two-component (X2C) and Douglas-Kroll-Hess method (DKH) levels in the Turbomole code is reported and tested for a series of smaller 3d¹, 4d¹, and 5d¹ complexes, as well as for some larger 5d⁷ Ir and Pt systems in comparison with earlier four-component matrix-Dirac-Kohn-Sham results. A main emphasis is placed on efficient approximations to the two-electron spin-orbit contributions, comparing an existing implementation of two variants of Boettger's "scaled nuclear spin-orbit" (SNSO) approximation in the code with a newly implemented atomic mean-field spin-orbit (AMFSO) approximation. The different variants perform overall comparably well with the four-component data. The AMFSO approximation has the added advantage of being able to include the spin-other-orbit contributions arising from the Gaunt term of relativistic electron-electron interactions. These are of comparably larger importance for the 3d complexes than for their heavier homologues. The excellent agreement between X2C and four-component electron paramagnetic resonance parameter results provides the opportunity to treat large systems efficiently and accurately with the computationally more expedient two-component quasirelativistic methodology.

Pseudocontact shifts and paramagnetic susceptibility in semiempirical and quantum chemistry theories

Pseudocontact shifts are traditionally described as a function of the anisotropy of the paramagnetic susceptibility tensor, according to the semiempirical theory mainly developed by Kurland and McGarvey [J. Magn. Reson. 2, 286-301 (1970)]. The paramagnetic susceptibility tensor is required to be symmetric. Applying point-dipole approximation to the quantum chemistry theory of hyperfine shift, pseudocontact shifts are found to scale with a non-symmetric tensor that differs by a factor g^T/g_e from the paramagnetic susceptibility tensor derived within the semiempirical framework. We analyze the foundations of the Kurland-McGarvey pseudocontact shift expression and recall that it is inherently based on the Russell-Saunders (LS) coupling approximation for the spin-orbit coupling. We show that the difference between the semiempirical and quantum chemistry pseudocontact shift expressions arises directly from the different treatment of the orbital contribution to the hyperfine coupling.

Paramagnetic NMR in solution and the solid state

The field of paramagnetic NMR has expanded considerably in recent years. This review addresses both the theoretical description of paramagnetic NMR, and the way in which it is currently practised. We provide a review of the theory of the NMR parameters of systems in both solution and the solid state. Here we unify the different languages used by the NMR, EPR, quantum chemistry/DFT, and magnetism communities to provide a comprehensive and coherent theoretical description. We cover the theory of the paramagnetic shift and shift anisotropy in solution both in the traditional formalism in terms of the magnetic susceptibility tensor, and using a more modern formalism employing the relevant EPR parameters, such as are used in first-principles calculations. In addition we examine the theory first in the simple non-relativistic picture, and then in the presence of spin-orbit coupling. These ideas are then extended to a description of the paramagnetic shift in periodic solids, where it is necessary to include the bulk magnetic properties, such as magnetic ordering at low temperatures. The description of the paramagnetic shift is completed by describing the current understanding of such shifts due to lanthanide and actinide ions. We then examine the paramagnetic relaxation enhancement, using a simple model employing a phenomenological picture of the electronic relaxation, and again using a more complex state-of-the-art theory which incorporates electronic relaxation explicitly. An additional important consideration in the solid state is the impact of bulk magnetic susceptibility effects on the form of the spectrum, where we include some ideas from the field of classical electrodynamics. We then continue by describing in detail the solution and solid-state NMR methods that have been deployed in the study of paramagnetic systems in chemistry, biology, and the materials sciences. Finally we describe a number of case studies in paramagnetic NMR that have been specifically chosen to highlight how the theory in part one, and the methods in part two, can be used in practice. The systems chosen include small organometallic complexes in solution, solid battery electrode materials, metalloproteins in both solution and the solid state, systems containing lanthanide ions, and multi-component materials used in pharmaceutical controlled-release formulations that have been doped with paramagnetic species to measure the component domain sizes.

Interplay of Through-Bond Hyperfine and Substituent Effects on the NMR Chemical Shifts in Ru(III) Complexes

The links between the molecular structure and nuclear magnetic resonance (NMR) parameters of paramagnetic transition-metal complexes are still relatively unexplored. This applies particularly to the contact term of the hyperfine contribution to the NMR chemical shift. We report combining experimental NMR with relativistic density functional theory (DFT) to study a series of Ru(III) complexes with 2-substituted β-diketones. A series of complexes with systematically varied substituents was synthesized and analyzed using ¹H and ¹³C NMR spectroscopy. The NMR spectra recorded at several temperatures were used to construct Curie plots and estimate the temperature-independent (orbital) and temperature-dependent (hyperfine) contributions to the NMR shift. Relativistic DFT calculations of electron paramagnetic resonance and NMR parameters were performed to interpret the experimental observations. The effects of individual factors such as basis set, density functional, exact-exchange admixture, and relativity are analyzed and discussed. Based on the calibration study in this work, the fully relativistic Dirac-Kohn-Sham (DKS) method, the GIAO approach (orbital shift), the PBE0 functional with the triple-ζ valence basis sets, and the polarizable continuum model for describing solvent effects were selected to calculate the NMR parameters. The hyperfine contribution to the total paramagnetic NMR (pNMR) chemical shift is shown to be governed by the Fermi-contact (FC) term, and the substituent effect (H vs Br) on the through-bond FC shifts is analyzed, interpreted, and discussed in terms of spin-density distribution, atomic spin populations, and molecular-orbital theory. In contrast to the closed-shell systems of Rh(III), the presence of a single unpaired electron in the open-shell Ru(III) analogs significantly alters the NMR resonances of the ligand atoms distant from the metal center in synergy with the substituent effect.

Relativistic DFT Calculations of Hyperfine Coupling Constants in 5d Hexafluorido Complexes: [ReF6 ]2- and [IrF6 ]2

The performance of relativistic density functional theory (DFT) methods has been investigated for the calculation of the recently measured hyperfine coupling constants of hexafluorido complexes [ReF6 ]2- and [IrF6 ]2- . Three relativistic methods were employed at the DFT level of theory: the 2-component zeroth-order regular approximation (ZORA) method, in which the spin-orbit coupling was treated either variationally (EV ZORA) or as a perturbation (LR ZORA), and the 4-component Dirac-Kohn-Sham (DKS) method. The dependence of the results on the basis set and the choice of exchange-correlation functional was studied. Furthermore, the effect of varying the amount of Hartree-Fock exchange in the hybrid functionals was investigated. The LR ZORA and DKS methods combined with DFT led to very similar deviations (about 20 %) from the experimental values for the coupling constant of complex [ReF6 ]2- by using hybrid functionals. However, none of the methods were able to reproduce the large anisotropy of the hyperfine coupling tensor of complex [ReF6 ]2- . For [IrF6 ]2- , the EV ZORA and DKS methods reproduced the experimental tensor components with deviations of ≈10 and ≈5 % for the hybrid functionals, whereas the LR ZORA method predicted the coupling constant to be around one order of magnitude too large owing to the combination of large spin-orbit coupling and very low excitation energies.

Cyclic Distannene or Bis(stannylene) with a Ferrocenyl Backbone: Synthesis, Structure, and Coordination Chemistry

1,1'-Dilithioferrocene was reacted with 2 equiv of isopropyl (Ar*) or methyl (Ar') substituted terphenyl tin(II) chloride. Reaction product 1, carrying the bulkier terphenyl substituent Ar*, displays a bis(stannylene) structure in the solid state without formation of a tin-tin bond. Temperature-dependent solution ¹¹⁹Sn NMR spectroscopy, however, revealed a dynamic interplay between bis(stannylene) (100 °C) and cyclic distannene (-80 °C). In contrast to 1, the less bulky Ar' substituent results in a cyclic distannene 2. On the basis of temperature-dependent ¹¹⁹Sn NMR spectroscopy the Sn-Sn bond of compound 2 was preserved up to 100 °C. Both compounds were further characterized by solid-state ¹¹⁹Sn NMR spectroscopy as well as ¹¹⁹Sn and ⁵⁷Fe Mössbauer spectroscopy. 1 reacted as a chelating ligand with nickel and palladium complexes [Ni(cod)₂] and [Pd(nbe)₃] (nbe = norbornene). In the resulting coordination compounds the nonstabilized stannylene acts as a donor as well as an acceptor ligand and shows a dynamic switch from donor to acceptor behavior in the monopalladium complex.

Towards reliable references for electron paramagnetic resonance parameters based on quantum chemistry: the case of verdazyl radicals

We present an efficient and accurate computational procedure to calculate properties measurable by EPR spectroscopy. We simulate a molecular dynamics (MD) trajectory by employing the quantum mechanically derived force field (QMDFF) [S. Grimme, J. Chem. Theory Comput., 2014, 10, 4497] and sample the trajectory at different time steps. For each snapshot EPR properties are calculated with a hybrid density functional theory (DFT) method. EPR spectra are simulated based on the averaged results. We applied the strategy to a number of previously published and novel verdazyl radicals, for which we recorded EPR spectra. The resulting simulated spectra are compatible with experiment already before employing an additional fitting step, in contrast to those from single point electronic-structure calculations. After the refinement, the experimental data are excellently reproduced, and the fitted EPR parameters do not deviate much from the calculated ones. This provides confidence in ascribing a direct physical meaning to the refined data in terms of experimental EPR parameters rather than merely considering them as mathematical fit parameters. We also find that couplings to hydrogen nuclei have a significant influence on the spectra of verdazyl radicals.

Silver complexes stabilized by large silanethiolate ligands - crystal structures and luminescence properties

Bulky silanethiolate and disiladithiolate ligands were applied to synthesize one mononuclear and three trinuclear silver complexes including two cyclic "microclusters" and a linear tri-nuclear silanethiolate complex. All obtained compounds are characterized by X-ray diffraction and FT-IR. NMR and emission spectroscopies were used where possible. The first trinuclear anionic silver thiolate is structurally characterized. The influence of the different charge of cyclic silver complexes as well as the overall ligand environment on the structural properties is demonstrated. The impact of the different synthetic routes on the final structures of the obtained clusters - cyclic or linear - is discussed.

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.

Pseudo-Contact NMR Shifts over the Paramagnetic Metalloprotein CoMMP-12 from First Principles

Long-range pseudo-contact NMR shifts (PCSs) provide important restraints for the structure refinement of proteins when a paramagnetic metal center is present, either naturally or introduced artificially. Here we show that ab initio quantum-chemical methods and a modern version of the Kurland-McGarvey approach for paramagnetic NMR (pNMR) shifts in the presence of zero-field splitting (ZFS) together provide accurate predictions of all PCSs in a metalloprotein (high-spin cobalt-substituted MMP-12 as a test case). Computations of 314 ¹³ C PCSs using g- and ZFS tensors based on multi-reference methods provide a reliable bridge between EPR-parameter- and susceptibility-based pNMR formalisms. Due to the high sensitivity of PCSs to even small structural differences, local structures based either on X-ray diffraction or on various DFT optimizations could be evaluated critically by comparing computed and experimental PCSs. Many DFT functionals provide insufficiently accurate structures. We also found the available 1RMZ PDB X-ray structure to exhibit deficiencies related to binding of a hydroxamate inhibitor. This has led to a newly refined PDB structure for MMP-12 (5LAB) that provides a more accurate coordination arrangement and PCSs.

A computational workflow for designing silicon donor qubits

Developing devices that can reliably and accurately demonstrate the principles of superposition and entanglement is an on-going challenge for the quantum computing community. Modeling and simulation offer attractive means of testing early device designs and establishing expectations for operational performance. However, the complex integrated material systems required by quantum device designs are not captured by any single existing computational modeling method. We examine the development and analysis of a multi-staged computational workflow that can be used to design and characterize silicon donor qubit systems with modeling and simulation. Our approach integrates quantum chemistry calculations with electrostatic field solvers to perform detailed simulations of a phosphorus dopant in silicon. We show how atomistic details can be synthesized into an operational model for the logical gates that define quantum computation in this particular technology. The resulting computational workflow realizes a design tool for silicon donor qubits that can help verify and validate current and near-term experimental devices.