Relativistic spin-orbit effects on hyperfine coupling tensors by density-functional theory

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.

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

Herein we describe Hyperion, a new program for computing relativistic picture-change-corrected magnetic resonance parameters from scalar relativistic active space wave functions, with or without spin–orbit coupling (SOC) included a posteriori. Hyperion also includes a new orbital decomposition method for assisting active space selection for calculations of hyperfine coupling. For benchmarking purposes, we determine hyperfine coupling constants of selected alkali metal, transition metal, and lanthanide atoms, based on complete active space self-consistent field spin–orbit calculations in OpenMolcas. Our results are in excellent agreement with experimental data from atomic spectroscopy as well as theoretical predictions from four-component relativistic calculations.

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.

Theoretical analysis of the long-distance limit of NMR chemical shieldings

After some years of controversy, it was recently demonstrated how to obtain the correct long-distance limit [point-dipole approximation (PDA)] of pseudo-contact nuclear magnetic resonance chemical shifts from rigorous first-principles quantum mechanics [Lang et al., J. Phys. Chem. Lett. 11, 8735 (2020)]. This result confirmed the classical Kurland-McGarvey theory. In the present contribution, we elaborate on these results. In particular, we provide a detailed derivation of the PDA both from the Van den Heuvel-Soncini equation for the chemical shielding tensor and from a spin Hamiltonian approximation. Furthermore, we discuss in detail the PDA within the approximate density functional theory and Hartree-Fock theories. In our previous work, we assumed a relatively crude effective nuclear charge approximation for the spin-orbit coupling operator. Here, we overcome this assumption by demonstrating that the derivation is also possible within the fully relativistic Dirac equation and even without the assumption of a specific form for the Hamiltonian. Crucial ingredients for the general derivation are a Hamiltonian that respects gauge invariance, the multipolar gauge, and functional derivatives of the Hamiltonian, where it is possible to identify the first functional derivative with the electron number current density operator. The present work forms an important foundation for future extensions of the Kurland-McGarvey theory beyond the PDA, including induced magnetic quadrupole and higher moments to describe the magnetic hyperfine field.

A Method for the Variational Calculation of Hyperfine-Resolved Rovibronic Spectra of Diatomic Molecules

An algorithm for the calculation of hyperfine structure and spectra of diatomic molecules based on the variational nuclear motion is presented. The hyperfine coupling terms considered are Fermi-contact, nuclear spin-electron spin dipole-dipole, nuclear spin-orbit, nuclear spin-rotation, and nuclear electric quadrupole interactions. Initial hyperfine-unresolved wave functions are obtained for a given set of potential energy curves and associated couplings by a variation solution of the nuclear-motion Schrödinger equation. Fully hyperfine-resolved parity-conserved rovibronic Hamiltonian matrices for a given final angular momentum, F, are constructed and then diagonalized to give hyperfine-resolved energies and wave functions. Electric transition dipole moment curves can then be used to generate a hyperfine-resolved line list by applying rigorous selection rules. The algorithm is implemented in Duo, which is a general program for calculating spectra of diatomic molecules. This approach is tested for NO and MgH, and the results are compared to experiment and shown to be consistent with those given by the well-used effective Hamiltonian code PGOPHER.

A method to calculate the NMR spectra of paramagnetic species using thermalized electronic relaxation

For paramagnetic species, it has been long understood that the hyperfine interaction between the unpaired electrons and the nucleus results in a nuclear magnetic resonance (NMR) peak that is shifted by a paramagnetic shift, rather than split by the coupling, due to an averaging of the electronic magnetic moment caused by electronic relaxation that is fast in comparison to the hyperfine coupling constant. However, although this feature of paramagnetic NMR has formed the basis of all theories of the paramagnetic shift, the precise theory and mechanism of the electronic relaxation required to predict this result has never been discussed, nor has the assertion been tested. In this paper, we show that the standard semi-classical Redfield theory of relaxation fails to predict a paramagnetic shift, as does any attempt to correct for the semi-classical theory using modifications such as the inhomogeneous master equation or Levitt-di Bari thermalization. In fact, only the recently-introduced Lindbladian theory of relaxation in magnetic resonance [J.Magn.Reson., 310, 106645 (2019)] is able to correctly predict the paramagnetic shift tensor and relaxation-induced linewidth in pNMR. Furthermore, this new formalism is able to predict the NMR spectra of paramagnetic species outside the high-temperature and weak-order limits, and is therefore also applicable to dynamic nuclear polarization. The formalism is tested by simulations of five case studies, which include Fermi-contact and spin-dipolar hyperfine couplings, g-anisotropy, zero-field splitting, high and low temperatures, and fast and slow electronic relaxation.

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.

Hyperfine-Coupling Tensors from Projected Hartree-Fock Theory

We assess the calculation of hyperfine coupling (HFC) tensors by different variants of Projected Hartree-Fock (PHF) theory. For a set of small main-group S = 1/2 radicals (BO, CO⁺, CN, AlO, vinyl, methyl, ethynyl), spin-symmetry as well as complex-conjugation and point-group symmetry are first broken in a reference determinant, and then variationally restored, in the frame of the modern formulation of PHF theory. Historically, PHF theory was basically restricted to the restoration of spin symmetry from an unrestricted HF determinant (conserving S_z symmetry). This afforded unsatisfactory HFCs. We obtain far better results for isotropic (and anisotropic) HFCs when the variational energy is further lowered by working with generalized determinants that completely break spin symmetry, and when additional symmetries are used. Specifically, complex-conjugation projection recovers a substantial fraction of the dynamical correlation energy in small molecules, and the detailed equations for combined complex-conjugation, spin- and point-group projection in the density-matrix/diagonalization formulation of PHF theory are here reported for the first time. The compact representation of the PHF wave function allows for a straightforward evaluation of the spin-density matrix and of HFC tensors with little effort. The promising performance of PHF theory may motivate the application of post-PHF methods to the calculation of HFC tensors.

Paramagnetic Pyrazolylborate Complexes Tp2M and Tp*2M: 1H, 13C, 11B, and 14N NMR Spectra and First-Principles Studies of Chemical Shifts

The paramagnetic pyrazolylborates Tp₂M and Tp*₂M (M = Cu, Ni, Co, Fe, Mn, Cr, V) as well as [Tp₂M]⁺ and [Tp*₂M]⁺ (M = Fe, Cr, V) have been synthesized and their NMR spectra recorded. The ¹H signal shift ranges vary from ∼30 ppm (Cu(II) and V(III)) to ∼220 ppm (Co(II)), and the ¹³C signal shift ranges from ∼180 ppm (Fe(III)) to ∼1150 ppm (Cr(II)). The ¹¹B and ¹⁴N shifts are ∼360 and ∼730 ppm, respectively. Both negative and positive shifts have been observed for all nuclei. The narrow NMR signals of the Co(II), Fe(II), Fe(III), and V(III) derivatives provide resolved ¹³C,¹H couplings. All chemical shifts have been calculated from first-principles on a modern version of Kurland-McGarvey theory which includes optimized structures, zero-field splitting, and g tensors, as well as signal shift contributions. Temperature dependence in the Fe(II) spin-crossover complex results from the equilibrium of the ground singlet and the excited quintet. We illustrate both the assignment and analysis capabilities, as well as the shortcomings of the current computational methodology.

Probing basis set requirements for calculating hyperfine coupling constants

A series of basis sets optimized for the calculation of the hyperfine coupling constant is proposed. The pcH-n basis sets are defined in qualities from double-ζ to pentuple-ζ for the elements H to Ar. They are derived from the polarization consistent basis sets by addition of two tight s-functions and one tight p-, d-, and f-function and are shown to provide an exponential convergence toward the complete basis set limit, and they have significantly lower basis set errors than other commonly used basis sets for a given ζ quality. The pcH basis sets display very similar basis set convergence with a range of density functional theory methods and may also be suitable for wave function based methods.

Magnetic susceptibility and paramagnetism-based NMR

The magnetic interactions between the nuclear magnetic moment and the magnetic moment of unpaired electron(s) depend on the structure and dynamics of the molecules where the paramagnetic center is located and of their partners. The long-range nature of the magnetic interactions is thus a reporter of invaluable information for structural biology studies, when other techniques often do not provide enough data for the atomic-level characterization of the system. This precious information explains the flourishing of paramagnetism-assisted NMR studies in recent years. Many paramagnetic effects are related to the magnetic susceptibility of the paramagnetic metal. Although these effects have been known for more than half a century, different theoretical models and new approaches have been proposed in the last decade. In this review, we have summarized the consequences for NMR spectroscopy of magnetic interactions between nuclear and electron magnetic moments, and thus of the presence of a magnetic susceptibility due to metals, and we do so using a unified notation.

What are the methodological and theoretical prospects for paramagnetic NMR in structural biology? A glimpse into the crystal ball

NMR spectroscopy is very sensitive to the presence of unpaired electrons, which perturb the NMR chemical shifts, J splittings and nuclear relaxation rates. These paramagnetic effects have attracted increasing attention over the last decades, and their use is expected to increase further in the future because they can provide structural information not easily achievable with other techniques. In fact, paramagnetic data provide long range structural restraints that can be used to assess the accuracy of crystal structures in solution and to improve them by simultaneous refinements with the X-ray data. They are also precious for obtaining information on the conformational variability of biomolecular systems, possibly in conjunction with SAXS and/or DEER data. We foresee that new tools will be developed in the next years for the simultaneous analysis of the paramagnetic data with data obtained from different techniques, in order to take advantage synergistically of the information content of all of them. Of course, the use of the paramagnetic data for structural purposes requires the knowledge of the relationship between these data and the molecular coordinates. Recently, the equations commonly used, dating back to half a century ago, have been questioned by first principle quantum chemistry calculations. Our prediction is that further theoretical/computational improvements will essentially confirm the validity of the old semi-empirical equations for the analysis of the experimental paramagnetic data.

How Do Nuclei Couple to the Magnetic Moment of a Paramagnetic Center? A New Theory at the Gauntlet of the Experiments

The recent derivation, based on pure quantum chemistry (QC) first-principles, of the pseudocontact shifts (PCSs) caused by a paramagnetic metal center on far away nuclei has cast doubts on the validity of the semiempirical (SE) theory, predicting PCSs to arise from the metal magnetic susceptibility anisotropy. The SE theory has been used and applied countless times, especially in the last 2 decades, to obtain structural information on proteins containing paramagnetic metal ions. We show here that the QC and SE predictions can be directly tested against experiments, provided a suitable macromolecular system is used. The SE approach yields a good prediction of the experimental PCSs while the QC one does not. It appears that the classic theory is able to grasp satisfactorily the underlying physics.

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.

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.

Ab initio paramagnetic NMR shifts via point-dipole approximation in a large magnetic-anisotropy Co(ii) complex

Transition metal complexes can possess a large magnetic susceptibility anisotropy, facilitating applications such as paramagnetic tags or shift agents in nuclear magnetic resonance (NMR) spectroscopy. Due to its g-shift and zero-field splitting (ZFS) we demonstrate on a Co(ii) clathrochelate with an aliphatic 16-carbon chain, a modern approach for ab initio calculation of paramagnetic susceptibility. Due to its large anisotropy, large linear dimension but relatively low number of atoms, the chosen complex is especially well-suited for testing the long-range point-dipole approximation (PDA) for the pseudocontact shifts (PCSs) of paramagnetic NMR. A static structure of the complex is used to compare the limiting long-distance PDA with full first-principles quantum-mechanical calculation. A non-symmetric formula for the magnetic susceptibility tensor is necessary to be consistent with the latter. Comparison with experimental shifts is performed by conformational averaging over the chain dynamics using Monte Carlo simulation. We observe satisfactory accuracy from the rudimentary simulation and, more importantly, demonstrate the fast applicability of the ab initio PDA.

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.

Paramagnetic solid-state NMR of materials

Solid-state NMR (SSNMR) experiments are mainly performed on diamagnetic materials. However, the study of paramagnetic materials offers access to additional information in combination with the unpaired electron. This article discusses the experimental complications arising when considering paramagnetic materials such as large shift effects and large linewidths. A look at the theoretical description as well as experimental methods are shown and accompanied by techniques for signal assignment, also with the help of quantum-chemical calculations.

Solution and solid state structures and magnetism of a series of linear trinuclear compounds with a hexacoordinate LnIII and two terminal NiII centers

Reported are the syntheses, structures and magnetic properties, also by NMR spectroscopy in solution, of a series of 13 linear trinuclear 3d-4f compounds with a lanthanide(iii) surrounded by two Ni^II ions, NiLn^III, where the central Ln^III is hexacoordinate. For three of the crystal structures, an additional H₂O molecule is coordinated to the central Ln^III ion, leading to a monocapped trigonal prismatic structure. However, NMR spectroscopy indicates that in solution, these complexes also have a hexacoordinate Ln^III center. The solution magnetic anisotropies, determined by NMR spectroscopy, indicate that the axial components of the anisotropies are relatively small and that the Dy^III derivative might therefore not exhibit single molecule magnetism. The axial anisotropies determined by NMR spectroscopy are in good agreement with the expectations based on the distorted trigonal prismatic ligand field.

Coupled-Perturbed Scheme for the Calculation of Electronic g-Tensors with Local Hybrid Functionals

A coupled-perturbed Kohn-Sham (CPKS) scheme for calculating second-order magnetic properties has been developed for the case of general occupied-orbital-dependent (OOD) exchange-correlation functionals involving the exact-exchange energy density. The origin of the coupling terms in the functional derivatives of OOD functionals with respect to the orbitals has been thoroughly analyzed, and general expressions for the resulting coupling terms have been obtained. The generalized CPKS scheme thus obtained has been implemented within the MAG-ReSpect code and tested in calculations of electronic g-tensors with local hybrid functionals. Compared to previously tested global hybrids, like B3LYP, thermochemically optimized local hybrids provide only little to moderate improvement for test sets of main-group radicals and paramagnetic transition-metal complexes. Closer analyses point to possible areas in which the fundamentally more flexible local hybrids may be improved for the property at hand.

QM/MM investigations of organic chemistry oriented questions

About 35 years after its first suggestion, QM/MM became the standard theoretical approach to investigate enzymatic structures and processes. The success is due to the ability of QM/MM to provide an accurate atomistic picture of enzymes and related processes. This picture can even be turned into a movie if nuclei-dynamics is taken into account to describe enzymatic processes. In the field of organic chemistry, QM/MM methods are used to a much lesser extent although almost all relevant processes happen in condensed matter or are influenced by complicated interactions between substrate and catalyst. There is less importance for theoretical organic chemistry since the influence of nonpolar solvents is rather weak and the effect of polar solvents can often be accurately described by continuum approaches. Catalytic processes (homogeneous and heterogeneous) can often be reduced to truncated model systems, which are so small that pure quantum-mechanical approaches can be employed. However, since QM/MM becomes more and more efficient due to the success in software and hardware developments, it is more and more used in theoretical organic chemistry to study effects which result from the molecular nature of the environment. It is shown by many examples discussed in this review that the influence can be tremendous, even for nonpolar reactions. The importance of environmental effects in theoretical spectroscopy was already known. Due to its benefits, QM/MM can be expected to experience ongoing growth for the next decade.In the present chapter we give an overview of QM/MM developments and their importance in theoretical organic chemistry, and review applications which give impressions of the possibilities and the importance of the relevant effects. Since there is already a bunch of excellent reviews dealing with QM/MM, we will discuss fundamental ingredients and developments of QM/MM very briefly with a focus on very recent progress. For the applications we follow a similar strategy.

Curie-type paramagnetic NMR relaxation in the aqueous solution of Ni(ii)

The magnetic field of the Curie spin manifests itself as both the pNMR shielding tensor and Curie relaxation, in analogy with CSA relaxation theory.

Predicting the paramagnet-enhanced NMR relaxation of H₂ encapsulated in endofullerene nitroxides by density-functional theory calculations

We have investigated the structure and nuclear magnetic resonance (NMR) spectroscopic properties of some dihydrogen endofullerene nitroxides by means of density-functional theory (DFT) calculations. Quantum versus classical roto-translational dynamics of H₂ have been characterized and compared. Geometrical parameters and hyperfine couplings calculated by DFT have been input to the Solomon-Bloembergen equations to predict the enhancement of the NMR longitudinal relaxation of H₂ due to coupling with the unpaired electron. Estimating the rotational correlation time via computed molecular volumes leads to a fair agreement with experiment for the simplest derivative; the estimate is considerably improved by recourse to the calculation of the diffusion tensor. For the other more flexible congeners, the agreement is less good, which may be due to an insufficient sampling of the conformational space. In all cases, relaxation by Fermi contact and Curie mechanisms is predicted to be negligible.

Probing the C₆₀ triplet state coupling to nuclear spins inside and out

The photoexcitation of functionalized fullerenes to their paramagnetic triplet electronic state can be studied by pulsed electron paramagnetic resonance (EPR) spectroscopy, whereas the interactions of this state with the surrounding nuclear spins can be observed by a related technique: electron nuclear double resonance (ENDOR). In this study, we present EPR and ENDOR studies on a functionalized exohedral fullerene system, dimethyl[9-hydro (C60-Ih)[5,6]fulleren-1(9H)-yl]phosphonate (DMHFP), where the triplet electron spin has been used to hyperpolarize, couple and measure two nuclear spins. We go on to discuss the extension of these methods to study a new class of endohedral fullerenes filled with small molecules, such as H₂@C₆₀, and we relate the results to density functional calculations.

Derivation and assessment of relativistic hyperfine-coupling tensors on the basis of orbital-optimized second-order Møller-Plesset perturbation theory and the second-order Douglas-Kroll-Hess transformation

The accurate calculation of hyperfine-coupling tensors requires a good description of the electronic spin density, especially close to and at the nucleus. Thus, dynamic correlation as well as relativistic effects have to be included in the quantum-chemical calculation of this quantity. In this paper, orbital-optimized second-order Møller-Plesset perturbation theory (MP2) is combined with the second-order Douglas-Kroll-Hess (DKH) transformation to yield an efficient and accurate ab initio method for the calculation of hyperfine couplings for larger molecules including heavy elements. Particular attention is paid to the derivation of the hyperfine-coupling tensor in the DKH framework. In the presence of a magnetic field, the DKH-transformation is not unique. Two different versions can be found in the literature. In this paper, a detailed derivation of one-electron contributions to the hyperfine-coupling tensor as they arise in linear-response theory is given for both DKH-transformations. It turns out that one of the two variants produces divergent hyperfine-coupling constants. The possibility to remove this divergence through a physically motivated finite-nucleus model taking into account the different extent of charge and magnetization distribution is discussed. Hyperfine-coupling values obtained at the orbital-optimized MP2 level with second-order DKH corrections for the non-divergent variant are presented. The influence of a Gaussian nucleus model is studied. The method is compared to four-component, high-accuracy calculations for a number of cations and atoms. Comparison to B3LYP and B2PLYP is made for a set of transition-metal complexes of moderate size.

First-principles calculation of parameters of electron paramagnetic resonance spectroscopy in solids

The hyperfine A-tensor and Zeeman g-tensor parameterize the interaction of an 'effective' electron spin with the magnetic field due to the nuclear spin and the homogeneous external magnetic field, respectively. The A- and g-tensors are the quantities of primary interest in electron paramagnetic resonance (EPR) spectroscopy. In this paper, we review our work [E.S. Kadantsev, T. Ziegler, J. Phys. Chem. A 2008, 112, 4521; E. S. Kadantsev, T. Ziegler, J. Phys. Chem. A 2009, 113, 1327] on the calculation of these EPR parameters under periodic boundary conditions (PBC) from first-principles. Our methodology is based on the Kohn-Sham DFT (KS DFT), explicit usage of Bloch basis set made up of numerical and Slater-type atomic orbitals (NAOs/STOs), and is implemented in the 'full potential' program BAND. Our implementation does not rely on the frozen core approximation. The NAOs/STOs basis is well suited for the accurate representation of the electron density near the nuclei, a prerequisite for the calculation of highly accurate hyperfine parameters. In the case of g-tensor, our implementation is based on the method of Van Lenthe et al. [E. van Lenthe, P. E. S. Wormer, A. van der Avoird, J. Chem. Phys. 1997, 107, 2488] in which the spin-orbital coupling is taken into account variationally. We demonstrate the viability of our scheme by calculating EPR parameters of paramagnetic defects in solids. We consider the A-tensor of 'normal' and 'anomalous' muonium defect in IIIA-VA semiconductors as well as the S2 anion radical in KCl host crystal lattice.

Zero-point vibrational corrections to isotropic hyperfine coupling constants in polyatomic molecules

The present work addresses isotropic hyperfine coupling constants in polyatomic systems with a particular emphasis on a largely neglected, but a posteriori significant, effect, namely zero-point vibrational corrections. Using the density functional restricted-unrestricted approach, the zero-point vibrational corrections are evaluated for the allyl radical and four of its derivatives. In addition for establishing the numerical size of the zero-point vibrational corrections to the isotropic hyperfine coupling constants, we present simple guidelines useful for identifying hydrogens for which such corrections are significant. Based on our findings, we critically re-examine the computational procedures used for the determination of hyperfine coupling constants in general as well as the practice of using experimental hyperfine coupling constants as reference data when benchmarking and optimizing exchange-correlation functionals and basis sets for such calculations.

Validation of broken-symmetry density functional methods for the calculation of electron paramagnetic resonance parameters of dinuclear mixed-valence MnIVMnIII complexes

. The EPR parameters of a series of dinuclear manganese(III,IV) complexes with mono(μ-oxo), bis(μ-oxo), (μ-oxo)(μ-carboxylato), bis(μ-oxo)(μ-carboxylato), and (μ-oxo)bis(μ-carboxylato) bridges were studied by broken-symmetry density functional (DFT) methods. The influence of the exchange-correlation functional on the agreement with experiment has been evaluated systematically for g tensors; 55Mn, 14N, and 1H hyperfine coupling tensors; and Heisenberg exchange couplings. 14N and 1H hyperfine couplings, 55Mn hyperfine anisotropies, g tensors, and exchange couplings are well described by hybrid functionals with moderate exact-exchange admixtures such as B3LYP. The isotropic 55Mn hyperfine couplings require larger exact-exchange admixtures. However, the errors of the B3LYP calculations are systematic and may be corrected by a constant scaling factor, providing good predictive power for a wide range of EPR parameters with broken-symmetry DFT and standard functionals. The influence of terminal and bridging ligands on structure, spin-density distributions, and EPR parameters are evaluated systematically. Computed hyperfine and g tensors are not covariant to each other. This may have consequences for spectra simulations. The nature of the broken-symmetry state and the origin of its spin contamination were analyzed by an expansion into restricted determinants, based on paired orbitals.

1H chemical shifts in nonaxial, paramagnetic chromium(III) complexes — Application of novel pNMR shift theory

We present the first chemical application of the recent, general theory of the nuclear magnetic resonance shielding and chemical shift in paramagnetic compounds, to a set of nonaxial high-spin metallo-organic complexes. The theory is for the first time rigorous for systems of arbitrary spatial and spin symmetry, and introduces new structure to the isotropic, anisotropic but symmetric, and anisotropic and antisymmetric parts of the shielding tensor. We apply the theory using density functional calculations of the proton chemical shift in a family of nonaxial chromium(III) complexes possessing a quartet ground electronic spin state. We discuss the various contributions to the isotropic chemical shift, and compare the full theory to approximate forms appropriate to the doublet case on the one hand, and to the doublet case at the nonrelativistic limit, on the other hand. The performance of various exchange-correlation functionals in reproducing the recently measured experimental chemical shifts is evaluated.