Zero-field splitting in nickel(II) complexes: A comparison of DFT and multi-configurational wavefunction calculations

Probing the dynamic behaviour and magnetic identification of seven coordinated Mn(II) complexes: a combined AIMD and multi-reference approach

We present an in-depth solution phase dynamics of rare seven coordinated pentagonal bipyramidal Mn(II) complexes, together with their binding affinity anticipated using ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT). Moreover, the simulations at different temperatures (25 °C and 90 °C) interpret the rigidity and stability of the ligands with Mn(II) ions. An intuitive approach for modulating the easy plane magnetic anisotropy of the mononuclear Mn(II) complex has been revealed by this work. In this regard, we have performed an extensive theoretical study based on the ab initio CASSCF/NEVPT2 method, exhibiting the presence of an easy plane magnetic anisotropy with a positive value of axial zero-field splitting (ZFS) parameter D. The complex's magnetic properties and electronic relaxation reveal that the rhombic ZFS term (E) can be modulated as the symmetry around the Mn(II) ion varies. The magnitude of the D-value increased with a more symmetrical equatorial ligand as found in the order of [Mn(pydpa)(H2O)] > [Mn(cbda)(H2O)]- > [Mn(dpaaa)(H2O)]- > [Mn(dpasam)(H2O)]-. Furthermore, we found that substituting the equatorial oxygen atom with heavier S and Se-donor atoms switches the sign of magnetic anisotropy for the Mn(II) complexes. The magnitude of the D-value increased when the energy levels of the ground state (GS) and the first excited state (ES) decreased. The observed magneto-structural correlation reveals that shortening the distance of the axial water molecule (Mn-O(w)) increases the D-value by an order of magnitude for the symmetrical [Mn(pydpa)(H2O)] complex. Overall, the combined analysis of solution phase dynamics of Mn(II) complexes and their magnetic characterization opens up new avenues in coordination chemistry, molecular magnetism, spin-crossover materials, and catalysis.

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.

Rapid and Accurate Prediction of the Axial Magnetic Anisotropy in Cobalt(II) Complexes Using a Machine-Learning Approach

Estimating the magnetic anisotropy for single-ion magnets is complex due to its multireference nature. This study demonstrates that deep neural networks (DNNs) can provide accurate axial magnetic anisotropy (D) values, closely matching the complete-active-space self-consistent-field (CASSCF) quality using density functional theory (DFT) data. We curated an 86-parameter database (UFF1) with electronic data from over 33000 cobalt(II) compounds. The DNN achieved an R2 of 0.906 and a mean absolute error of 18.1 cm-1 in comparison to reference CASSCF D values. Remarkably, it is 11 times more accurate than DFT methods and 7700 times faster. This approach hints at DNNs predicting the anisotropy in larger molecules, even when trained on smaller ligands.

Limitations on the D-Parameter in Ni(II) Complexes

A number of hexacoordinate, pentacoordinate, and tetracoordinate Ni(II) complexes have been investigated by applying ab initio CASSCF + NEVPT2 + SOC calculations and Generalized Crystal Field Theory. The geometry of the coordination polyhedron covers D4h, D3h, D2h, D2d, C4v, C3v, and C2v symmetry. The calculated spin-Hamiltonian parameters D and E were compared to the available experimental data. The limiting values of the D-parameter in the class of Ni(II) complexes are identified. Magnetic anisotropy in Ni(II) complexes, expressed by the axial zero-field splitting parameter D, seriously depends upon the ground and first excited electronic states. In hexacoordinate complexes, the ground electronic term is nondegenerate 3B1g for the D4h symmetry; D is slightly positive or negative. In tetracoordinate systems, D is only positive when the electronic ground state is nondegenerate 3A or 3B; this diverges on the τ4 path when oblate bisphenoid approaches the prolate geometry and a level crossing with 3E occurs. In pentacoordinate systems, D could be extremely negative when approaching a trigonal bipyramid (Addison index τ5 ∼ 1, ground state 3E″). In pentacoordinate Ni(II) complexes with the D3h and C3v symmetry of the coordination polyhedron, the ground electronic term is orbitally doubly degenerate which causes the D-parameter stays undefined. It is emphasized that one has to inspect compositions of the spin-orbit multiplets from the spin states |MS⟩ and check whether the weights confirm the expected spin-Hamiltonian picture: with D > 0, the ground state contains a dominant part of |0⟩ (close to 100%) whereas with D < 0 the spin-orbit doublet is formed of |±1⟩ with high weights (approaching 50 + 50%). The calculations show that the situations are not black and white, and the mixing of the states might be more complex especially when the rhombic zero-field splitting parameter E is in the play. In the case of the 3E ground term, six spin-orbit multiplets are formed by mixing six |MS⟩ states from the ground and quasi-degenerate excited states.

Magnetic Behavior of Trigonal (Bi-)pyramidal 3d8 Mononuclear Nanomagnets: The Case of [Ni(MDABCO)2Cl3]ClO4

This paper attempts to shed light on the origin of the magnetic behavior specific to trigonal bi- and pyramidal 3d8 mono- and polynuclear nanomagnets. The focus lies on entirely unraveling the system's intrinsic microscopic mechanisms and fundamental quantum mechanical relations governing the underlying electron dynamics. To this end, we develop a self-consistent approach to characterize, in great detail, all electron correlations and the ensuing fine structure of the energy spectra of a broad class of 3d8 systems. The mathematical framework is based on the multiconfigurational self-consistent field method and is devised to account for prospective quantum mechanical constraints that may confine the electron orbital dynamics while preserving the properties of all measurable quantities. We successfully characterize the experimentally observed magnetic anisotropy properties of a slightly distorted trigonal bipyramidal Ni2+ coordination complex, demonstrating that such compounds do not exhibit intrinsic huge zero-field splitting and inherent giant magnetic anisotropy. We reproduce qualitatively and quantitatively the behavior of the low-field magnetic susceptibility, magnetization, low-, and high-field electron paramagnetic resonance spectroscopy measurements and provide an in-depth analysis of the obtained results.

Fine Structure and the Huge Zero-Field Splitting in Ni2+ Complexes

We perform a thorough study of the ground state magnetic properties of nickel-based 3d8 complexes. This includes an in-depth analysis of the contribution of the crystal field, spin exchange and spin-orbit interactions to the ground state magnetic properties. Of particular interest to the current investigation are the presence and occurrence of non-trivial zero-field splitting. The study focuses on the cases of Ni2+ ideal octahedral, trigonal bipyramidal, square planar and tetrahedral geometries. We provide results for the complete energy spectrum, the fine structure related to the ground state and the second set of excited states, low-field magnetic susceptibility and magnetization. In addition, we examine the zero-field fine structure in square pyramidal, trigonal pyramidal and trigonal planar complexes. The obtained results unequivocally show that a moderate or highly coordinated 3d8 complex can neither exhibit spin-orbit-driven large and giant magnetic anisotropy nor a huge zero-field splitting. Moreover, in the trigonal bipyramidal coordination, a fine structure associated to the ground state cannot result from the spin-orbit coupling alone.

Single-Ion Magnets with Giant Magnetic Anisotropy and Zero-Field Splitting

The design of mononuclear molecular nanomagnets exhibiting a huge energy barrier to the reversal of magnetization have seen a surge of interest during the last few decades due to their potential technological applications. More specifically, single-ion magnets are peculiarly attractive by virtue of their rich quantum behavior and distinct fine structure. These are viable candidates for implementation as single-molecule high-density information storage devices and other applications in future quantum technologies. The present review presents the comprehensive state of the art in the topic of single-ion magnets possessing an eminent magnetization-reversal barrier, very slow magnetic relaxation and high blocking temperature. We turn our attention to the achievements in the synthesis of 3d and 4f single-ion magnets during the last two decades and discuss the observed magnetostructural properties underlying the anisotropy behavior and the ensuing remanence. Furthermore, we highlight the fundamental theoretical aspects to shed light on the complex behavior of these nanosized magnetic entities. In particular, we focus on key notions, such as zero-field splitting, anisotropy energy and quantum tunneling of the magnetization and their interdependence.

Computational design of magnetic molecules and their environment using quantum chemistry, machine learning and multiscale simulations

Having served as a playground for fundamental studies on the physics of d and f electrons for almost a century, magnetic molecules are now becoming increasingly important for technological applications, such as magnetic resonance, data storage, spintronics and quantum information. All of these applications require the preservation and control of spins in time, an ability hampered by the interaction with the environment, namely with other spins, conduction electrons, molecular vibrations and electromagnetic fields. Thus, the design of a novel magnetic molecule with tailored properties is a formidable task, which does not only concern its electronic structures but also calls for a deep understanding of the interaction among all the degrees of freedom at play. This Review describes how state-of-the-art ab initio computational methods, combined with data-driven approaches to materials modelling, can be integrated into a fully multiscale strategy capable of defining design rules for magnetic molecules.

Multiconfiguration Pair-Density Functional Theory for Chromium(IV) Molecular Qubits

Pseudotetrahedral organometallic complexes containing chromium(IV) and aryl ligands have been experimentally identified as promising molecular qubit candidates. Here we present a computational protocol based on multiconfiguration pair-density functional theory for computing singlet-triplet gaps and zero-field splitting (ZFS) parameters in Cr(IV) aryl complexes. We find that two multireference methods, multistate complete active space second-order perturbation theory (MS-CASPT2) and hybrid multistate pair-density functional theory (HMS-PDFT), perform better than Kohn-Sham density functional theory for singlet-triplet gaps. Despite the very small values of the ZFS parameters, both multireference methods performed qualitatively well. MS-CASPT2 and HMS-PDFT performed particularly well for predicting the trend in the ratio of the rhombic and axial ZFS parameters, |E/D|. We have also investigated the dependence and sensitivity of the calculated ZFS parameters on the active space and the molecular geometry. The methodologies outlined here can guide future prediction of ZFS parameters in molecular qubit candidates.

Unravelling the robustness of magnetic anisotropy of a nickelocene molecule in different environments: a first-principles-based study

Recent scanning tunneling spectroscopy with single metallocene molecule-functionalized tips have proved to be a powerful tool to probe and control individual spins and spin-spin exchange interactions due to the robustness of the magnetic properties of the metallocene molecule in different surroundings. However, accurate prediction of such robustness at a first-principles-based level by the conventional density functional theory (DFT) has remained challenging. In this paper, we have performed a detailed investigation of the evolution of electronic and magnetic properties of a nickelocene molecule (NiCp₂) in different environments, i.e., free-standing, adsorbed on Cu(100) and as a functionalized tip apex. Using an embedding method, which combines DFT and the complete active space self-consistent field (CASSCF) method recently developed, we demonstrate that the nickelocene molecule almost preserves its spin and magnetic anisotropy upon adsorption on Cu(100), and also in the position of the tip apex. In particular, the cyclic π* orbital of the Cp rings could hybridize with the singly occupied d_π orbitals of the Ni center of the molecule, protecting these orbitals from external states. Hence the molecular spin maintains S = 1, the same as in the free-standing case, and its magnetic anisotropy is also robust with energies of 3.56, 3.34, and 3.51 meV in free-standing, adsorbed on Cu(100), and functionalized tip apex states, respectively, in good agreement with previous theoretical and experimental results. This work thus provides a first-principles-based understanding of the relevant experiments. Such agreement between theoretical simulations and experimental measurements highlights the potential usefulness of the method for investigating the local electronic and spin states of organometallic molecule-surface composite systems.

Antisymmetric Exchange in a Real Copper Triangular Complex

The antisymmetric exchange, also known as the Dzyaloshinskii-Moriya interaction (DMI), is an effective interaction that may be at play in isolated complexes (with transition metals or lanthanides, for instance), nanoparticles, and highly correlated materials with adequate symmetry properties. While many theoretical works have been devoted to the analysis of single-ion zero-field splitting and to a lesser extent to symmetric exchange, only a few ab initio studies deal with the DMI. Actually, it originates from a subtle interplay between weak electronic interactions and spin-orbit couplings. This article aims to highlight the origin of this interaction from theoretical grounds in a real tri-copper(II) complex, capitalizing on previous methodological studies on bi-copper(II) model complexes. By tackling this three-magnetic-center system, we will first show that the multispin model Hamiltonian is appropriate for trinuclear (and likely for higher nuclearity) complexes, then that the correct application of the permutation relationship is necessary to explain the outcomes of the ab initio calculations, and finally, that the model parameters extracted from a binuclear model transfer well to the trinuclear complex. For a more theory-oriented purpose, we will show that the use of a simplified structural model allows one to perform more demanding electronic structure calculations. On this simpler system, we will first check that the previous transferability is still valid, prior to performing more advanced calculations on the derived two-magnetic-center model system. To this end, we will explain in detail the physics of the DMI in the copper triangle of interest, before advocating further theory/experiment efforts.

Interpretation of EPR and optical spectra of Ni(II) ions in crystalline lattices at ambient temperature

Many biologically important paramagnetic metal ions are characterized by electron paramagnetic resonance (EPR) spectroscopy to use as spin probes to investigate the structure and function of biomolecules. Though nickel(II) ions are an essential trace element and part of many biomolecules, the EPR properties are least understood. Herein, the EPR and optical absorption spectra measured at 300 K for Ni(II) ions diluted in two different diamagnetic hosts are investigated and reported. The EPR spectrum of a polycrystalline Ni/Mg(3-methylpyrazole)₆ (ClO₄ )₂ [Ni/MMPC] shows two transitions at X-band frequency (~9.5 GHz), suggesting the zero-field splitting parameter (D) is larger than the resonance field of the free electron (H_o ). This incomplete and complex spectrum is successfully analyzed to obtain EPR parameters. The EPR spectrum of the polycrystalline Ni/Zn(pyrazole)₆ (NO₃ )₂ [Ni/ZPN] shows a triplet spectrum indicating D < H_o . A detailed analysis of single-crystal EPR data yielded the spin Hamiltonian parameters. The optical absorption spectra are deconvoluted to understand the symmetry of the coordination environment in the complex.

Modern multireference methods and their application in transition metal chemistry

Transition metal chemistry is a challenging playground for quantum chemical methods owing to the simultaneous presence of static and dynamic electron correlation effects in many systems. Wavefunction based multireference (MR) methods constitute a physically sound and systematically improvable Ansatz to deal with this complexity but suffer from some conceptual difficulties and high computational costs. The latter problem partially arises from the unfavorable scaling of the Full Configuration Interaction (Full-CI) problem which in the majority of MR methods is solved for a subset of the molecular orbital space, the so-called active space. In the last years multiple methods such as modern variants of selected CI, Full-CI Quantum Monte Carlo (FCIQMC) and the density matrix renormalization group (DMRG) have been developed that solve the Full-CI problem approximately for a fraction of the computational cost required by conventional techniques thereby significantly extending the range of applicability of modern MR methods. This perspective review outlines recent advancements in the field of MR electronic structure methods together with the resulting chances and challenges for theoretical studies in the field of transition metal chemistry. In light of its emerging importance a special focus is put on the selection of adequate active spaces and the concomitant development of numerous selection aides in recent years.

Magnetic Anisotropy in a Cubane-like Ni4 Complex: An Ab Initio Perspective

Magnetic anisotropy, in the absence of an external magnetic field, relates to the degeneracy lift of energy levels. In the standard case of transition metal complexes, this property is usually modeled by an anisotropic spin Hamiltonian and one speaks of "zero-field splitting" (ZFS) of spin states. While the case of mononuclear complexes has been extensively described by means of ab initio quantum mechanical calculations, the literature on polynuclear complexes studied with these methodologies is rather scarce. In this work, advanced multiconfigurational wave function theory methods are applied to compute the ZFS of the ground S = 4 state of an actual tetranickel(II) complex, displaying a magnet behavior below 0.5 K. First, the isotropic couplings are computed in the absence of the spin-orbit coupling operator, in the full complex and also in clusters with only two active nickel(II) centers, confirming the occurrence of weak ferromagnetic couplings in this system. Second, the single-site magnetic anisotropies are computed on a cluster bearing only one active nickel(II) site, showing that the single-site anisotropy axes are not oriented in an optimal fashion for generating a large uniaxial molecular anisotropy. Furthermore, the possibility for involving only a few local orbital excited states in the calculation is assessed, actually opening the way for a consistent and manageable treatment of the ZFS of the ground S = 4 state. Third, multiconfigurational calculations are performed on the full complex, confirming the weak uniaxial anisotropy occurring for this state and also, interestingly, revealing a significant contribution of the lowest-lying orbitally excited S = 3 states. Overall, by comparison with the experiment, the reported results question the common habit of using only one structure, in particular derived from a crystallography experiment, to compute magnetic anisotropy parameters.

Structure and magnetic properties of a novel heteroheptanuclear metal string complex [Ni3Ru2Ni2(μ7-teptra)4(NCS)2](PF6)

We report the synthesis of a novel heteroheptanuclear metal string complex (HMSC) [Ni₃Ru₂Ni₂(μ₇-teptra)₄(NCS)₂](PF₆) 1 supported by tetra-pyridyl-tri-amine (H₃teptra) ligands. We employed X-ray diffraction and other spectroscopic techniques to characterize the complex. The observed remarkably short Ru-Ru distance of 2.2499(3) Å for 1 is indicative of a unique metal-metal interaction in the mixed-valence [Ru₂]⁵⁺ (S = 3/2) unit. The complex exhibits a relatively high magnetic moment value of 4.55 B.M. at 4 K, which increases rapidly to 6.00 B.M. at 30 K and remains at 6.11 B.M. from 50 to 300 K as shown by SQUID measurements, indicating a high spin (S≥ 3/2) system which is further supported by the analyses of EPR spectra at low temperatures. These magnetic behaviors can be ascribed to the result of spin-exchange interactions among multi-spin centers.

Field-dependent paramagnetic relaxation enhancement in solutions of Ni(II): What happens above the NMR proton frequency of 1 GHz?

An extended set of paramagnetic relaxation enhancement (PRE) data, up to the field of 32.9 Tesla, is reported for protons in an acidified aqueous solution of a Ni(II) salt in the presence and in the absence of added glycerol. For the 55% w/w glycerol sample, a distinct maximum in the PRE vs magnetic field curve is observed for the first time. The data are analysed using the Swedish slow-motion theory, including both the intramolecular (inner-sphere) and intermolecular (outer-sphere) contributions. The results indicate that estimating the outer-sphere part in the presence of the more efficient inner-sphere term is a difficult task.

On the magnetic properties of nanodiamonds: Electronic g-tensor calculations

The electronic g-tensor calculations are carried out for various paramagnetic defects introduced into hydrogenated diamond nanocrystal C₃₅H₃₆, showing that such a system can be successfully used to model magnetic properties of nanodiamonds (NDs) with paramagnetic centers containing no vacancies. In addition, it is revealed that, depending on the geometric positions in ND, paramagnetic centers of the same type produce noticeable variations of the g-tensor values. A side-by-side comparison of the performance of effective nuclear charge and spin-orbit mean field (SOMF) approaches indicates that the latter is more sensitive to the quality of basis sets, especially concerning diffuse functions, the inclusion of which is found to be nonbeneficial. What is more, the SOMF method also exhibits a much more pronounced gauge-origin dependence. Compared to electronic charge centroid, spin centers (SCs) demonstrate a superior suitability as gauge origins, providing a better agreement with diamagnetic and paramagnetic contributions of g-tensor obtained employing gauge-including atomic orbitals (GIAOs). Therefore, SCs can be recommended for the g-tensor calculations of NDs whenever GIAOs are not available.

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.

Modelling spin Hamiltonian parameters of molecular nanomagnets

Molecular nanomagnets encompass a wide range of coordination complexes possessing several potential applications. A formidable challenge in realizing these potential applications lies in controlling the magnetic properties of these clusters. Microscopic spin Hamiltonian (SH) parameters describe the magnetic properties of these clusters, and viable ways to control these SH parameters are highly desirable. Computational tools play a proactive role in this area, where SH parameters such as isotropic exchange interaction (J), anisotropic exchange interaction (Jx, Jy, Jz), double exchange interaction (B), zero-field splitting parameters (D, E) and g-tensors can be computed reliably using X-ray structures. In this feature article, we have attempted to provide a holistic view of the modelling of these SH parameters of molecular magnets. The determination of J includes various class of molecules, from di- and polynuclear Mn complexes to the {3d-Gd}, {Gd-Gd} and {Gd-2p} class of complexes. The estimation of anisotropic exchange coupling includes the exchange between an isotropic metal ion and an orbitally degenerate 3d/4d/5d metal ion. The double-exchange section contains some illustrative examples of mixed valance systems, and the section on the estimation of zfs parameters covers some mononuclear transition metal complexes possessing very large axial zfs parameters. The section on the computation of g-anisotropy exclusively covers studies on mononuclear Dy(III) and Er(III) single-ion magnets. The examples depicted in this article clearly illustrate that computational tools not only aid in interpreting and rationalizing the observed magnetic properties but possess the potential to predict new generation MNMs.

Recent progress in the electron paramagnetic resonance study of polymers

This review article provides an overview of the contemporary research based on a tailor-made technique to understand the paramagnetic behavior of different polymer classes.

Rapid and precise determination of zero-field splittings by terahertz time-domain electron paramagnetic resonance spectroscopy

Zero-field splitting (ZFS) parameters are fundamentally tied to the geometries of metal ion complexes. Despite their critical importance for understanding the magnetism and spectroscopy of metal complexes, they are not routinely available through general laboratory-based techniques, and are often inferred from magnetism data. Here we demonstrate a simple tabletop experimental approach that enables direct and reliable determination of ZFS parameters in the terahertz (THz) regime. We report time-domain measurements of electron paramagnetic resonance (EPR) signals associated with THz-frequency ZFSs in molecular complexes containing high-spin transition-metal ions. We measure the temporal profiles of the free-induction decays of spin resonances in the complexes at zero and nonzero external magnetic fields, and we derive the EPR spectra via numerical Fourier transformation of the time-domain signals. In most cases, absolute values of the ZFS parameters are extracted from the measured zero-field EPR frequencies, and the signs can be determined by zero-field measurements at two different temperatures. Field-dependent EPR measurements further allow refined determination of the ZFS parameters and access to the g-factor. The results show good agreement with those obtained by other methods. The simplicity of the method portends wide applicability in chemistry, biology and material science.

Assignment of solid-state 13C and 1H NMR spectra of paramagnetic Ni(II) acetylacetonate complexes aided by first-principles computations

Recent advances in computational methodology allowed for first-principles calculations of the nuclear shielding tensor for a series of paramagnetic nickel(II) acetylacetonate complexes, [Ni(acac)₂L₂] with L = H₂O, D₂O, NH₃, ND₃, and PMe₂Ph have provided detailed insight into the origin of the paramagnetic contributions to the total shift tensor. This was employed for the assignment of the solid-state ^1,2H and ¹³C MAS NMR spectra of these compounds. The two major contributions to the isotropic shifts are by orbital (diamagnetic-like) and contact mechanism. The orbital shielding, contact, as well as dipolar terms all contribute to the anisotropic component. The calculations suggest reassignment of the ¹³C methyl and carbonyl resonances in the acac ligand [Inorg. Chem.53, 2014, 399] leading to isotropic paramagnetic shifts of δ(¹³C) ≈ 800-1100 ppm and ≈180-300 ppm for ¹³C for the methyl and carbonyl carbons located three and two bonds away from the paramagnetic Ni(II) ion, respectively. Assignment using three different empirical correlations, i.e., paramagnetic shifts, shift anisotropy, and relaxation (T₁) were ambiguous, however the latter two support the computational results. Thus, solid-state NMR spectroscopy in combination with modern quantum-chemical calculations of paramagnetic shifts constitutes a promising tool for structural investigations of metal complexes and materials.

Interplay of spin-dependent delocalization and magnetic anisotropy in the ground and excited states of [Gd2@C78]- and [Gd2@C80]. J Chem Phys 2017;147:124305. [PMID: 28964020 DOI: 10.1063/1.5004183]

The magnetic properties and electronic structure of the ground and excited states of two recently characterized endohedral metallo-fullerenes, [Gd₂@C₇₈]^- (1) and [Gd₂@C₈₀]^- (2), have been studied by theoretical methods. The systems can be considered as [Gd₂]⁵⁺ dimers encapsulated in a fullerene cage with the fifteen unpaired electrons ferromagnetically coupled into an S = 15/2 high-spin configuration in the ground state. The microscopic mechanisms governing the Gd-Gd interactions leading to the ferromagnetic ground state are examined by a combination of density functional and ab initio calculations and the full energy spectrum of the ground and lowest excited states is constructed by means of ab initio model Hamiltonians. The ground state is characterized by strong electron delocalization bordering on a σ type one-electron covalent bond and minor zero-field splitting (ZFS) that is successfully described as a second order spin-orbit coupling effect. We have shown that the observed ferromagnetic interaction originates from Hund's rule coupling and not from the conventional double exchange mechanism. The calculated ZFS parameters of 1 and 2 in their optimized geometries are in qualitative agreement with experimental EPR results. The higher excited states display less electron delocalization, but at the same time they possess unquenched first-order angular momentum. This leads to strong spin-orbit coupling and highly anisotropic energy spectrum. The analysis of the excited states presented here constitutes the first detailed study of the effects of spin-dependent delocalization in the presence of first order orbital angular momentum and the obtained results can be applied to other mixed valence lanthanide systems.

Perspectives on paramagnetic NMR from a life sciences infrastructure

The effects arising in NMR spectroscopy because of the presence of unpaired electrons, collectively referred to as "paramagnetic NMR" have attracted increasing attention over the last decades. From the standpoint of the structural and mechanistic biology, paramagnetic NMR provides long range restraints that can be used to assess the accuracy of crystal structures in solution and to improve them by simultaneous refinements through NMR and X-ray data. These restraints also provide information on structure rearrangements and conformational variability in biomolecular systems. Theoretical improvements in quantum chemistry calculations can nowadays allow for accurate calculations of the paramagnetic data from a molecular structural model, thus providing a tool to refine the metal coordination environment by matching the paramagnetic effects observed far away from the metal. Furthermore, the availability of an improved technology (higher fields and faster magic angle spinning) has promoted paramagnetic NMR applications in the fast-growing area of biomolecular solid-state NMR. Major improvements in dynamic nuclear polarization have been recently achieved, especially through the exploitation of the Overhauser effect occurring through the contact-driven relaxation mechanism: the very large enhancement of the ¹³C signal observed in a variety of liquid organic compounds at high fields is expected to open up new perspectives for applications of solution NMR.

Behaviour of DFT-based approaches to the spin–orbit term of zero-field splitting tensors: a case study of metallocomplexes, MIII(acac)3 (M = V, Cr, Mn, Fe and Mo)

Zero-field splitting tensors of M^III(acac)₃ complexes are calculated using ab initio and DFT methods.

Ab initio modelling of magnetic anisotropy in Sr3NiPtO6

First principles calculations in the framework of Density Functional Theory (DFT) and wavefunction-based correlated methods have been performed to investigate in detail the magnetic anisotropy in Sr3NiPtO6. This material is known for the easy-plane anisotropy with a large anisotropy constant of about 7.5-9.3 meV. We find that by properly choosing the onsite Coulomb repulsion and exchange parameters, DFT can correctly explain the easy-plane magnetocrystalline anisotropy of the material, but the magnitude of the anisotropy constant is underestimated. On the other hand, a quantitative agreement with respect to experiments, both in the magnitude and direction of the magnetic anisotropy, can be recovered by using the wavefunction-based approach which is able to fully describe the multiplet physics. We also show that the presence of structural distortions of the local NiO6 coordination is crucial for stabilizing the magnetic anisotropy in this compound.

Origin of the Zero-Field Splitting in Mononuclear Octahedral Mn(IV) Complexes: A Combined Experimental and Theoretical Investigation

The aim of this work was to determine and understand the origin of the electronic properties of Mn(IV) complexes, especially the zero-field splitting (ZFS), through a combined experimental and theoretical investigation on five well-characterized mononuclear octahedral Mn(IV) compounds, with various coordination spheres (N6, N3O3, N2O4 in both trans (trans-N2O4) and cis configurations (cis-N2O4) and O4S2). High-frequency and -field EPR (HFEPR) spectroscopy has been applied to determine the ZFS parameters of two of these compounds, MnL(trans-N2O4) and MnL(O4S2). While at X-band EPR, the axial-component of the ZFS tensor, D, was estimated to be +0.47 cm(-1) for MnL(O4S2), and a D-value of +2.289(5) cm(-1) was determined by HFEPR, which is the largest D-magnitude ever measured for a Mn(IV) complex. A moderate D value of -0.997(6) cm(-1) has been found for MnL(trans-N2O4). Quantum chemical calculations based on two theoretical frameworks (the Density Functional Theory based on a coupled perturbed approach (CP-DFT) and the hybrid Ligand-Field DFT (LF-DFT)) have been performed to define appropriate methodologies to calculate the ZFS tensor for Mn(IV) centers, to predict the orientation of the magnetic axes with respect to the molecular ones, and to define and quantify the physical origin of the different contributions to the ZFS. Except in the case of MnL(trans-N2O4), the experimental and calculated D values are in good agreement, and the sign of D is well predicted, LF-DFT being more satisfactory than CP-DFT. The calculations performed on MnL(cis-N2O4) are consistent with the orientation of the principal anisotropic axis determined by single-crystal EPR, validating the calculated ZFS tensor orientation. The different contributions to D were analyzed demonstrating that the d-d transitions mainly govern D in Mn(IV) ion. However, a deep analysis evidences that many factors enter into the game, explaining why no obvious magnetostructural correlations can be drawn in this series of Mn(IV) complexes.

Magnetic couplings in the chemical shift of paramagnetic NMR

We apply the Kurland-McGarvey (J. Magn. Reson. 1970, 2, 286) theory for the NMR shielding of paramagnetic molecules, particularly its special case limited to the ground-state multiplet characterized by zero-field splitting (ZFS) interaction of the form S·D·S. The correct formulation for this problem was recently presented by Soncini and Van den Heuvel (J. Chem. Phys. 2013, 138, 054113). With the effective electron spin quantum number S, the theory involves 2S+1 states, of which all but one are low-lying excited states, between which magnetic couplings take place by Zeeman and hyperfine interactions. We investigate these couplings as a function of temperature, focusing on both the high- and low-temperature behaviors. As has been seen in work by others, the full treatment of magnetic couplings is crucial for a realistic description of the temperature behavior of NMR shielding up to normal measurement temperatures. At high temperatures, depending on the magnitude of ZFS, the effect of magnetic couplings diminishes, and the Zeeman and hyperfine interactions become effectively averaged in the thermally occupied states of the multiplet. At still higher temperatures, the ZFS may be omitted altogether, and the shielding properties may be evaluated using a doublet-like formula, with all the 2S+1 states becoming effectively degenerate at the limit of vanishing magnetic field. We demonstrate these features using first-principles calculations of Ni(II), Co(II), Cr(II), and Cr(III) complexes, which have ZFS of different sizes and signs. A non-monotonic inverse temperature dependence of the hyperfine shift is predicted for axially symmetric integer-spin systems with a positive D parameter of ZFS. This is due to the magnetic coupling terms that are proportional to kT at low temperatures, canceling the Curie-type 1/kT prefactor of the hyperfine shielding in this case.

Magnetic anisotropy of mononuclear Ni(II) complexes: on the importance of structural diversity and the structural distortions

Mononuclear Ni(II) complexes are particularly attractive in the area of single-molecule magnets as the axial zero-field splitting (D) for the Ni(II) complexes is in the range of -200 to +200 cm(-1) . Despite this advantage, very little is known on the origin of anisotropy across various coordination ligands, coordination numbers, and particularly what factors influence the D parameter in these complexes. To answer some of these questions, herein we have undertaken a detailed study of a series of mononuclear Ni(II) complexes with ab initio calculations. Our results demonstrate that three prominent spin-conserved low-lying d-d transitions contribute significantly to the D value. Variation in the sign and the magnitude of D values are found to correlate to the specific structural distortions. Apart from the metal-ligand bond lengths, two different parameters, namely, Δα and Δβ, which are correlated to the cis angles present in the coordination environment, are found to significantly influence the axial D values. Developed magneto-structural D correlations suggest that the D values can be enhanced significantly by fine tuning the structural distortion in the coordination environment. Calculations performed on a series of Ni(II) models with coordination numbers two to six unfold an interesting observation-the D parameter increases significantly upon a reduction in coordination number compared with a reference octahedral coordination. Besides, if high symmetry is maintained, even larger coordination numbers yield large D values.

An ab initio MO study of heavy atom effects on the zero-field splitting tensors of high-spin nitrenes: how the spin–orbit contributions are affected

The CASSCF and hybrid CASSCF–MRMP2 methods reproduce the ZFS tensors of spin-septet 2,4,6-trinitrenopyridines, focusing on the heavy atom effects on the spin–orbit terms of the tensors.

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.