Next Article in Journal
A Clock Transition in the Cr7Mn Molecular Nanomagnet
Next Article in Special Issue
DFT Investigations of the Magnetic Properties of Actinide Complexes
Previous Article in Journal
Interfacial Spin Manipulation of Nickel-Quinonoid Complex Adsorbed on Co(001) Substrate
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling the Effect of Zero-Field Splitting on the 1H, 13C and 29Si Chemical Shifts of Lanthanide and Actinide Compounds

by
Austin W. Lloyd
*,
Helen M. Moylan
and
Joseph J. W. McDouall
School of Chemistry, The University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
Magnetochemistry 2019, 5(1), 3; https://doi.org/10.3390/magnetochemistry5010003
Submission received: 30 September 2018 / Revised: 17 December 2018 / Accepted: 20 December 2018 / Published: 11 January 2019
(This article belongs to the Special Issue Magnetic Properties of Complexes of Actinide Elements)

Abstract

:
The prediction of paramagnetic NMR (pNMR) chemical shifts in molecules containing heavy atoms presents a significant challenge to computational quantum chemistry. The importance of meeting this challenge lies in the central role that NMR plays in the structural characterisation of chemical systems. Hence there is a need for reliable assignment and prediction of chemical shifts. In a previous study [Trends in Physical Chemistry, 17, 25–57, (2017)] we looked at the computation of pNMR chemical shifts in lanthanide and actinide complexes using a spin Hamiltonian approach. In that study we were principally concerned with molecules with S = 1/2 ground states. In the present work we extend that study by looking at the effect of zero field splitting (ZFS) for six complexes with S = 3/2 ground states. It is shown that the inclusion of ZFS can produce substantial shifts in the predicted chemical shifts. The computations presented are typically sufficient to enable assignment of experimental spectra. However for one case, in which the peaks are closely clustered, the inclusion of ZFS re-orders the chemical shifts making assignment quite difficult. We also observe, and echo, the previously reported importance of including the paramagnetic spin-orbit hyperfine interaction for 13 C and 29 Si atoms, when these are directly bound to a heavy element and thus subject to heavy-atom-light-atom effects. The necessary computations are very demanding, and more work is needed to find theoretical and computational approaches that simplify the evaluation of this term. We discuss the computation of each term required in the spin Hamiltonian. The systems we study in this work are restricted to a single heavy atom ion (one Nd(III) and five U(III) complexes), but typify some of the computational complexity encountered in lanthanide and actinide containing molecules.

1. Introduction

Interest in lanthanide and actinide chemistry continues to expand as new applications are found and existing areas of research are developed. Examples include novel synthetic techniques, such as the development of reusable reagents in the synthesis of cyclic oxycarbon compounds [1], and the investigation into the nature of actinide containing single molecule magnets [2].
Radioactivity is the most exploited property of the actinides, as they potentially offer an efficient and low carbon alternative to fossil fuels for large scale energy generation. At the opposite end of the scale, nuclear powered batteries are used for powering small devices such as heart pacemakers [3]. Because of their industrial use, the reprocessing of nuclear waste to extract fissile material and reduce waste volumes is of paramount importance. This leads directly to a demand for better extractants, which in turn needs a better understanding of metal-ligand selectivity for actinide over lanthanide ions [4]. This is an area that benefits from improved and reliable computational modelling techniques.
Nuclear magnetic resonance spectroscopy (NMR) is a well-established and much used analytical tool for characterising molecules and providing insight into their constituent chemical bonds. Many experimental techniques exist for obtaining NMR spectra in diverse chemical and physical situations. The analysis and interpretation of complex spectra for most diamagnetic molecules, that lack heavy nuclei, are routine. Complementary to this experimental work, the development of computational tools for predicting NMR chemical shifts has also become common [5].
NMR studies of actinide species provide significant challenges to computational techniques. Firstly, for light atoms the effects of relativity are small, and the representation of electron correlation is significantly more important. However the relativistic correction to the energy of an atomic system approximately scales as ( Z c 1 ) 2 [6]. For heavy elements relativistic effects have a profound influence on the ordering of energy levels within their atoms and molecules, and a concomitant effect on the molecular properties of such systems. The computational treatment of relativistic effects begins with scalar relativistic corrections and these are often sufficient for the prediction of molecular geometries. For paramagnetic molecules the presence of the unpaired electronic spin requires that the coupling between orbital and spin angular momenta are treated reliably. This can also have a significant effect on the splitting of electronic energy levels. Equally, spin-dependent properties such as the Zeeman interaction, the hyperfine interaction and the spin-spin interaction will depend on the theoretical method used for their evaluation and require calibration against reliable experimental data. These considerations provide significant challenges for the theoretical formulation of paramagnetic NMR (pNMR) chemical shifts, as well as for their computational realisation. Calculating pNMR shifts of paramagnetic molecules remains a challenging problem that is still under much investigation [7,8].
Experimental spectra are normally reported as δ shifts. The δ shift is obtained with respect to a reference chemical shift, σ ref . It is conventional to separate the experimental paramagnetic chemical shift σ exp into three terms: the orbital contribution σ orb , which is approximated as the shift that would be observed if the species was diamagnetic, and the contact σ c and pseudocontact σ pc shifts arising from the unpaired electronic spin:
δ = σ ref σ exp
σ exp = σ orb + σ c + σ pc .
For simple systems, the contact term relates to the interaction of the unpaired electron spin density and the nuclear magnetic moment at the nucleus. The contact interaction is isotropic. The pseudocontact shift is the residual long-range interaction of the spin-dipoles and is anisotropic [9]. Qualitatively the contact shift is traditionally interpreted as arising from the transfer of spin density through bonds and may be associated with covalent interactions. The pseudocontact shift is interpreted as acting through space and is a dipole-dipole type interaction [10]. The theoretical justification for the decomposition into these terms weakens when relativistic treatments are included [11].
The theoretical framework of pNMR shifts can be formulated in a number of ways. The full pNMR shielding expression in a sum over states (SOS) formalism was derived in references [12,13]. In reference [14] multireference wave function calculations of pNMR shifts were performed, for the first time, directly employing the SOS formalism without reference to the spin Hamilitonian formalism. This study was extended to carbonate complexes of U, Np and Pu [15]. In references [14,15] the wave function code used was implemented as part of the OpenMolcas code but has not yet been released. Accordingly, in this work we shall adopt the spin Hamiltonian approach, which represents each of the three paramagnetic operators with a coupling matrix and is detailed in the next section. This approach has the advantage that the necessary components are readily available in many quantum chemistry packages; however, the main disadvantage is that the computational framework that is used to calculate each component has a significant effect on the results obtained. To achieve reliable results it is necessary to assess each component, and calibrate the most effective approach against an often-limited set of experimental values for similar chemical species. Detailed overviews of the evaluation of spin Hamiltonian parameters in the context of modern computational chemistry are available from Vaara [16], Bolvin [17] and Autschbach [18]. We shall mostly be concerned with studying the effect of zero field splitting (ZFS) on the predicted 1 H, 13 C, 29 Si shifts of six compounds. We will illustrate situations where the inclusion of ZFS causes a reordering in the predicted 1 H shifts, and where the inclusion of spin-orbit effects on the hyperfine term for the 13 C shifts of nuclei involved in a η -coordinated C-U bond dominate the hyperfine interaction. With the exception of the first example, our results are accurate enough to enable qualitative assignment of observed shifts.

2. The Spin Hamiltonian Approach to Paramagnetic NMR

Applying a magnetic field B to an isolated atom with nuclear spin I nuc results in each energy level splitting by 2 γ nuc | B | , where γ nuc is the magnetogyric ratio of the nucleus under investigation. In any electronic system, the circulation of electrons induced by the applied magnetic field generates an opposing magnetic field. This effect is represented by σ orb , the chemical shift tensor. For a diamagnetic molecule, the Hamiltonian corresponding to the energy shift, in atomic units, is given by
H NMR = γ nuc B · ( 1 σ orb ) · I n u c .
However for a system with unpaired electrons, there is a permanent magnetic moment associated with the electronic spin. It is energetically favourable for this to align with an applied magnetic field and the resulting (thermally averaged) dipole must be accounted for. The corresponding paramagnetic NMR Hamiltonian is given by [19,20,21,22]:
H pNMR = γ nuc B · ( 1 σ orb ) · I n u c + μ B B · g · S ˜ + S ˜ · D · S ˜ + S ˜ · A · I n u c .
The Zeeman coupling matrix, g , represents the interaction of the electronic spin-dipole with the applied magnetic field. An effective parameter, referred to as the pseudospin S ˜ , is introduced. In the absence of spin-orbit coupling S ˜ = S , which is an often-used approximation. For molecules where the spin-orbit interaction is strong, the effective degeneracy of the molecule is reduced so a smaller effective spin S ˜ is chosen based on the splitting pattern of the electronic states [17]. The interaction between the nuclear and electronic spin-dipoles is represented by the hyperfine coupling matrix, A . In a system with more than one unpaired electron, the interaction of the electron spin-dipoles is represented by the ZFS or D tensor.

2.1. The Orbital Shielding Tensor, σ o r b

For a diamagnetic molecule the orbital, or Ramsey, shielding matrix [23] is the only contribution to the observed chemical shift and as such, its isotropic value is readily extractable from NMR spectra. For a paramagnetic molecule it is not possible to directly extract this term and as a result, the reliability of any computational method must be assessed on diamagnetic analogues of the species under investigation. Because magnetic fields are defined by the action of an infinitesimal rotation of a vector potential and this operation is independent of the origin of the co-ordinate system used, the calculation of magnetic properties has an artificial dependence on the choice of origin, known as gauge dependence. The gauge dependence vanishes in the limit of the exact wave function. It has been demonstrated that for calculations targeting a single nucleus, the dependence is minimised by setting the origin at the nucleus. For molecules in general, choosing the centre of electronic charge as the origin minimises gauge dependence [24], but the most reliable approach is to employ gauge-including atomic orbitals (GIAO) [25,26], which remove the gauge dependence.
For actinide compounds, modelling relativistic effects is essential and many frameworks can be used. The relativistic Dirac equation can be satisfied by a four component wave function, partitioned into pairs of electronic and positronic spinors. The presence of the positronic solution hinders a variational approach, since simply minimising the energy is no longer possible due to an unbounded and infinite set of negative positronic states. For four component methods, this issue must be addressed, e.g., by kinetic balance [27]. However, since the positronic component is much smaller, efficient two component techniques have been developed to eliminate it. These techniques produce a spin-free (scalar) and a spin-dependent component. Observables in the Dirac framework are implicitly dependent on all four components, and the corresponding two component operator erroneously does not involve the full solution. This effect is called the picture change error [18] and can be addressed using a transformation of the operator equivalent to that used to transform the four component wave function to a two component wave function.
Scalar relativistic approaches such as the zero-order regular approximation (ZORA) [28], and the Douglas-Kroll Hamiltonian (DKH) [29] are implemented in a number of quantum chemistry packages. ZORA only holds for basis functions that are independent of the perturbation, which theoretically precludes using GIAOs, although in practice the error introduced is small for large basis sets [30].

2.2. The Zeeman Coupling Matrix, g

The Zeeman coupling matrix, g , more commonly referred to as the “g tensor”, (since it relates two quantities defined in space and spin coordinate frames respectively it is not a true tensor.) represents the effect of the magnetic field on the average electronic spin-dipole. With the addition of spin-orbit coupling, the electronic spin of a molecule is no longer well defined and the pseudospin parameter is introduced to describe the observed splitting of spin multiplets [19]. The theory and computation of spin-orbit coupling has a detailed review by Marian [31]. Because spin-orbit coupling has such a large effect on the effective spin-dipole, the reliable evaluation of the g tensor is directly dependent on how well this effect is modelled.
For cases of a weak spin-orbit interaction, density functional theory (DFT) and linear response theory (DFT-LRT) can successfully be used [32]. However DFT-LRT is unable to produce reliable results where degeneracies or low-lying excited states exist, or where the ground state is multiconfigurational in nature. Instead we have used a multiconfigurational approach that can describe the ground and excited states equally well, i.e., state-averaged complete active-space self-consistent field theory (SA-CASSCF). It is possible to include spin-orbit effects at a variational level. This type of approach is available in the SPOCK-CI program [33], which is a spin-orbit configuration interaction (CI) package. While this would be the preferred approach it is computationally demanding and perturbation techniques can be used to represent the interaction more efficiently.
One multiconfigurational approach is to use double perturbation theory and a SOS formalism to generate the g tensor directly. The double perturbation theory approach can be modified for degenerate systems, but for nearly degenerate systems, individual contributions to the SOS depend inversely on the energy difference between states, if this energy difference is small the resulting perturbation theory can be divergent. The spin-orbit term can also be applied using a quasi-degenerate perturbation theory (QDPT) approach which constructs a matrix over the manifold of spin states. The diagonal elements are taken as the energies of the spin-free states ( e.g., from a SA-CASSCF calculation), and the off-diagonal elements couple the spin-free states under the action of the spin-orbit operator. Diagonalisation gives the spin-orbit energy levels and wave functions.
In the case of a molecule with spin S = 1 2 , the lowest energy pair of spin-orbit wave functions must be degenerate in the absence of an applied magnetic field due to Kramers theorem [34], which states that the degeneracy of a state with half-integer spin must be at least two in the absence of a magnetic field. The energy splitting Δ E Zeeman upon application of a magnetic field is given by the difference in expectation values of the Kramers pair, represented by | Φ ± 1 / 2 , the subscript indicating the sign of the associated pseudospin, that is,
Δ E Zeeman = μ B Φ 1 / 2 | B · μ | Φ 1 / 2 μ B Φ 1 / 2 | B · μ | Φ 1 / 2 .
= μ B Φ 1 / 2 | B · g · S ˜ | Φ 1 / 2 μ B Φ 1 / 2 | B · g · S ˜ | Φ 1 / 2 .
In the above equation, the magnetic moment operator μ is related to the magnetogyric ratio for the electron, g e = 2.0023193 and the spin, S , and angular momentum, L , operators by μ = μ B ( L + g e S ) . Squaring Equations (5) and (6) gives a relationship between the related tensor G = g g T , called the Abragam-Bleaney tensor [35], and the matrix elements of μ in the Kramers pair basis. The principal values of the Zeeman matrix are the square roots of the eigenvalues of the Abragam-Bleaney tensor [36]. While this relationship is only valid for S = 1 2 , similar ones have been derived for S > 1 2 [37]. This process has two drawbacks, the resulting g matrix is necessarily symmetric, and because of the use of a square root precursor matrix, there is a sign indeterminacy for the principal values of g . For weak spin coupling, the principal values must be close to the free electron value, and therefore positive. Otherwise, the signs can be inferred by exploiting relationships between the sign of combinations of the principal values and the sign of combinations of the eigenvalues of the μ κ matrices (where κ represents a cartesian direction) [37].
The final consideration is the treatment of the spin-orbit operator. It is possible to use the full two-electron form of the operator, but as with two-electron repulsion integrals, this results in a considerable computational cost. It is often preferable to use a one-electron operator. Three operators are in common use, the first an empirical operator based on an effective nuclear charge calibrated against calculations on small molecules as developed by Kosekii [38,39,40]. The second is the spin-orbit mean field (SOMF) operator which forms an effective one-electron operator, similar to the Fock operator in Hartree–Fock theory. Finally the atomic mean-field approximation (AMFI) to the spin–orbit operator forms a similar one-electron operator, but discards the multicentre two-electron spin-orbit integrals, exploiting the short range of the two-electron spin-orbit operator while retaining the simpler one-centre, two-electron integrals [41,42].

2.3. The ZFS Tensor, D

The ZFS interaction primarily represents the interelectronic dipolar interaction, and can be split into a first order term, directly corresponding to spin-spin interaction and a second-order term representing corrections due to spin-orbit coupling [43]. When represented over a manifold of spin states S , M S , the first order term is diagonal. It is evaluated over the ground state eigenfunction. The second-order term includes off-diagonal elements with selection rule 1 Δ S , Δ M S 1 . This term must be evaluated over ground and excited states. Available computational approaches are similar to those used for the Zeeman matrix. Single reference wave function methods include a modified DFT coupled-perturbed approach to handle the spin-orbit perturbation [44]. In multireference approaches, the spin-spin contribution can be obtained from first-order perturbation theory, while the spin-orbit contribution requires second-order perturbation theory. However, as for the Zeeman matrix, near degeneracy often produces spurious results. A QDPT effective Hamiltonian approach gives more reliable ZFS tensors in such cases [45], and it is the latter approach that we have adopted. We shall outline below for each molecule the choice of active space used and states averaged in producing the ZFS tensor.

2.4. The Hyperfine Coupling Matrix, A

The hyperfine coupling constant (HFCC) matrix represents the interaction between the nuclear and the electronic spin-dipoles, and like the Zeeman matrix, the matrix for nucleus K can be related to an operator with matrix F K , κ = S ˜ κ · A K . For a Kramers doublet it is given by [17]
A K T A K κ υ = 2 g K 2 M p 2 ψ , ψ doublet ψ | F K , κ | ψ ψ | F K , υ | ψ
where g K is the g-factor for the target nucleus, and M p the mass of the proton. A wide range of approaches are available, such as DFT linear response and SOS double perturbation methods. However, the hyperfine interaction has strong inverse dependence on electron-nucleus separation and is strongly influenced by relativistic effects [6]. Relativistic corrections cannot be ignored, at least scalar relativistic treatments are needed for a reliable approach. However additional corrections for spin-orbit coupling are computationally demanding, since three sets of coupled-perturbed equations must be solved for each nucleus. In principle, it is possible to avoid the solution of the CP-SCF equations by adopting a 4-component approach [46,47], although this has an associated computational overhead in performing the 4-component calculation. For a light target nucleus the corrections are small and often can be safely ignored, except when the target nucleus is bonded to a heavy nucleus; the latter can have a substantial effect on the electronic structure of its neighbour, a phenomenon referred to as the ”heavy atom-light atom” (HALA) effect [48], and the spin-orbit correction to the hyperfine interaction must be included for reliable results.
A remaining issue in predicting isotropic chemical shifts is due to our calculations being based on a single geometry, ignoring vibrational and rotational motion. This is a valid approach for electronic transitions which operate on a much faster timescale than significant nuclear motion, but nuclear transitions are slow enough for the motion to be significant. The pNMR measurement is an average over nuclear motion. It is necessary to calculate the orbital shielding and hyperfine matrices for all nuclei and then average them in appropriate groupings to generate the averaged hyperfine matrices for each group of equivalent atoms.

2.5. Zero-Field Spitting—Theory

Ramsey [12,23,49] posited that because nuclear motion is much slower than electronic motion, the changes in the nuclear magnetic moment μ nuc changes much more slowly than the electronic magnetic moment μ , and can be treated as time-independent, effectively as a parameter. Assuming equilibrium is reached within experimental time scales, the observed NMR shielding tensor can be shown to be given by the second derivative of the electronic energy with respect to the components of the magnetic field and the nuclear magnetic moment at the limit of zero field and magnetic moment.
σ κ υ = 2 E B κ μ nuc , ν B 0 μ nuc 0
The resulting energy levels are often close enough to have significant occupation above the ground state, so Boltzmann averaging of observables is required. The Boltzmann average A 0 of a set of observables { A n } with n energy levels, { W n } is given by:
A 0 = n A n e W n / k B T n e W n / k B T
Following Moon and Patchkovskii’s summary [50] the approach we use is based on Kurland and McGarvey’s work [51] which allows the pNMR Hamiltonian to be split into Zeeman and hyperfine terms. The field-independent form of the shielding tensor [50] can be decomposed into diamagnetic and paramagnetic shielding tensors, σ orb and σ p respectively:
σ κ υ orb = 2 E B κ μ nuc , ν 0 σ κ υ p = 1 k B T E B κ · E μ nuc , υ 0 ,
Taking g and A as the following partial derivatives
g κ υ = 2 E B κ S ν B 0 A κ υ = γ nuc 2 E S κ μ nuc , ν μ nuc 0
it can be shown that, in the low spin-orbit coupling limit [10,11,13]
1 1 σ a b p = 1 k B T E B a · E μ nuc , b 0 = 1 γ nuc k B T ( g · S ) a ( S · A ) b 0 = 1 γ nuc k B T κ υ { x , y , z } g a κ Z κ υ A υ b ,
where Z is a 3 × 3 coupling matrix. In the absence of zero-field splitting, Z κ υ = 1 3 S ( S + 1 ) δ κ υ , where δ κ υ is the Kronecker delta and the resulting matrix is diagonal.
Pennanen and Vaara [22] formulated a practical approach to modelling ZFS in the context of modern quantum chemistry, the ”no coupling approximation.” This requires the formation of the ( 2 S + 1 ) × ( 2 S + 1 ) matrix, S · D · S , which is then diagonalised to obtain eigenvalues E n and eigenstates | n , the latter being linear combinations of pure spin states | S M S . Soncini and van den Heuvel [13] presented a fully coupled approach, which can be expressed as:
Z κ υ = ω 1 n m eigenstates q n m n | S κ | m m | S υ | n where ω = n e E n / k B T
with q n m being a Boltzmann weighting term
q n m = e E n / k B T E m = E n k B T E n E m ( e E n / k B T e E m / k B T ) otherwise
Pierre Curie observed that the magnetisation of a paramagnetic material is approximately inversely proportional to the temperature. If the occupation of the ground state ( i.e., its Boltzmann weight) is nearly one, this is observed for a doublet and in the no coupling approximation. In the coupling described in Equations (12) and (14), the collection of terms corresponding to E m = E n have approximate 1 / T dependence and the combined contribution is referred to as the Curie term. The remaining terms have an approximate 1 / T 2 dependence.
The derivation used for both approaches assumes weak spin-orbit coupling. Van den Heuvel and Soncini presented an extended approach [52] which is valid for strong spin-orbit coupling. Because of its potential as an improved treatment, an outline of the approach is presented in Appendix A. The approach depends on using Equation A1 to form irreducible tensor operator (ITO) representations of the μ and F K matrices introduced in Section 2.2 and Section 2.4 respectively.

2.6. Contact and Pseudocontact Shifts

The separation of the paramagnetic chemical shift into contact and pseudocontact shifts provides insight into to the chemical environment of the target nucleus; however, as additional effects are included in the model, the separation becomes increasingly difficult and debatable. The simplest approach is to decompose the Zeeman and hyperfine coupling matrices into isotropic and anisotropic components, and evaluate their nature and relative sizes in terms of their scaling with respect to the fine structure constant α [22,53].
In the decomposition of the Zeeman matrix, the only zeroth-order term corresponds to the magnetogyric ratio of the free electron, the remaining terms relate the g-shift, Δ g = g g e I 3 which is decomposed into an isotropic scalar g-shift, Δ g iso = 1 3 trace Δ g , and the remaining anisotropic matrix, Δ g ˜ = Δ g Δ g iso I 3 . The hyperfine matrix has no first order terms, but is entirely second-order in the absence of relativistic corrections and we will denote the scalar relativistic (strictly the scalar relativistic corrections implicitly include O ( α 4 ) terms, but to simplify the analysis we will consider A SR to be a O ( α 2 ) term.) matrix as A SR . This is decomposed into its isotropic scalar A iso = 1 3 trace A SR and an anisotropic spin-dipolar matrix A SD = A SR A iso I 3 , which is symmetric. The relativistic corrections are entirely fourth order, and are separated into an isotropic scalar A iso Δ SO and an anisotropic A aniso Δ SO matrix. The anisotropic component is further split into a symmetric matrix A Δ SO , SYM = 1 2 A aniso Δ SO + A aniso Δ SO T and an asymmetric matrix A Δ SO , AS = A aniso Δ SO A Δ SO , SYM . These decompositions are outlined in Table 1.
The product of the matrices, g · A , has fifteen terms, six of which are sixth-order in α . In this work we use two methods to qualitatively divide the paramagnetic chemical shift into contact and pseudocontact contributions [50]. The first method is to relate the terms derived from isotropic scalars of the hyperfine matrix to the contact shift and those from the anisotropic matrices to the pseudocontact shift [10,54]. The second is to truncate the expansion to fourth order, resulting in nine remaining terms. The remaining difference is that all contributions with an anisotropic component, most notably the term A iso Δ g ˜ (a contact term in method 1), are considered to relate to the pseudospin [11,55,56]. We will refer to these decompositions as methods 1 and 2 respectively, as outlined in Table 2. Method 2 is elaborated on by Kaupp [7], who presents a separation of terms for a doublet. In EPR terminology, A iso Δ SO is sometimes referred to as the pseudocontact term.
The anisotropic matrices are traceless, so without ZFS, Z is diagonal, so only the isotropic hyperfine/isotropic Zeeman and the anisotropic hyperfine/anisotropic Zeeman terms produce an interaction with a non-zero trace. In the absence of spin-orbit coupling these interactions can unambiguously be related to the contact and pseudocontact terms, respectively. However the nature of the fourth order spin-orbit corrections to the hyperfine is unclear and the application of zero-field splitting results in non-diagonal elements in Z which couple isotropic and anisotropic terms. As a result, the separation of magnetic chemical shifts into contact, bond-mediated effects and pseudocontact, dipolar-interaction effects is not strictly consistent. With this in mind, we will use both methods to analyse our results.
The truncation to fourth order terms used in the second method is based on the assumption that the order in α determines the size of the contributions. This implies that the spin-orbit corrections to the hyperfine are small compared to the scalar relativistic hyperfine. If these, computationally demanding, corrections are not included, no sixth-order terms are present, so the final shift is the same for both decompositions. However in cases of hyperfine calculations on heavy nuclei or on nuclei that are bonded to heavy nuclei, this assumption is not valid: spin-orbit corrections typically scale in terms of the atomic number to the fourth power (Z 4 ) making the spin-orbit correction considerably larger. In such cases, the isotropic spin-orbit correction may be larger than the isotropic value of the uncorrected hyperfine matrix [17].

2.7. An Example of a Spin-Orbit Description of a Molecular Ground State

Our calculations allow for mixing of electronic states under the action of the spin-orbit operator. The number of electronic states included in each calculation was motivated by the quality of the Zeeman matrix produced, as discussed in our previous paper [57]. As an example, for our first actinide molecule (Section 4.2) we performed state-averaged CASSCF calculations over 14 quartet and 25 doublet states (i.e., being a subset of the 35 quartet and 112 doublet states implicitly included in the CASSCF(7,3) configuration state function (CSF) expansion.) to generate the Zeeman and ZFS matrices. The resulting ground state wave function, corrected for spin-orbit coupling is detailed in Table 3.
Chibotaru states that the spin Hamiltonian approach is applicable when a QDPT-SOC calculation is convergent for a manifold of spin states [37], and later he discusses the differences between the fictitious pseudospin and the spin. In the limiting cases of weak SOC the pseudospin can be set equal to the spin. In the limiting case of strong SOC, the degeneracy of the states is broken in low symmetry molecules, and the pseudospin mirrors a doublet, S ˜ = 1/2 [17]. Chibotoaru considers actinide ionic complexes to fall in the intermediate range [37], where the model Hamiltonian cannot be defined by the total angular momentum because of the mixture of states that we describe above; however, he does still associate the pseudospin with a manifold of electronic states, and so this is the approach we employ. Furthermore, in systems for which S 3 / 2 the calculation ideally should include third order tensorial spin, Zeeman and Hyperfine matrices [52], as we summarise in the Appendix A. However, it must be noted that to implement this scheme would require the generalised hyperfine tensors in Equation (7), that are not currently generally available. We are not aware of the use of the generalised hyperfine operators by any other researchers.
In Table 4 we list the relative energies of the first five Kramers pairs for all molecules studied in this work.

3. Previous Work and the ZFS Tensor

In a recent study we looked at the pNMR shifts of 14 f-element compounds [57]. The present work studies the effects of zero field splitting on the six f 3 compounds from that study. One compound contains neodymium as a sample lanthanide, detailed in Section 4.1 and shown in Figure 1. The remaining five are uranium compounds detailed in Section 4.2, Section 4.3, Section 4.4, Section 4.5 and Section 4.6. Equation (4) shows that the pNMR shift depends explicitly on the four magnetic interactions ( σ orb , g , D and A ) and implicitly on the geometry at which the various terms are evaluated. We summarise the findings of our previous study in justifying the approach we have adopted in evaluating each of the necessary terms. More details can be found in the aforementioned paper.

3.1. Geometry Optimisation, R

Geometry optimisation was performed using GAUSSIAN 09 [58]. Since geometries largely depend on the ground state electronic structure, it was considered reasonable to use DFT with the PBE0 [59] functional. Due to the presence of uranium, relativistic effects are important, but it is sufficient to only consider scalar relativistic effects. We prefer an all-electron approach and previously created a program in which implements the ZORA [27,60], and applies the correction to the one-electron Hamiltonian in atomic blocks, removing the gauge dependence issue. This program was interfaced to the GAUSSIAN programs and used in all geometry optimisations. We used the segmented all-electron relativistically corrected (SARC) [61,62] basis set for the heavy atoms with Def2-SVP [63] for the remaining atoms. This approach showed reliable agreement against a test set of crystal structure data for a set of 11 molecules, producing correlation coefficients of R 2 = 1.00 for bond length and R 2 = 0 . 99 for bond angles [57].

3.2. Orbital Shielding, σ o r b

The orbital components of the chemical shifts were obtained using GAUSSIAN 09, DFT linear response theory (PBE0), the ZORA relativistic treatment outlined above, and GIAO [25,64]. A Def2-TZVP basis was used for 1 H chemical shifts, and an uncontracted Def2-SVP basis for the 13 C, otherwise the basis used was as in the previous section. This approach was calibrated against 20 1 H and 32 13 C chemical shifts of 4 diamagnetic molecules, producing correlation coefficients for 1 H and C 13 of R 2 = 0.99 and R 2 = 0.98 respectively [57]. Deviations were noted when a light atom under investigation was directly bonded to a heavy atom. These cases are subject to the HALA effect working through the chemical bond and as such, need a more elaborate relativistic treatment [48].

3.3. Hyperfine Coupling Matrix, A

The reported A tensors were calculated using the ORCA 3.0.3 program [65], DFT (PBE0), the Douglas–Kroll–Hess transformation to second-order (DKH2) [66,67]. Analysis of the effect that picture change has on these hyperfine coupling coefficients (without spin-orbit correction) was performed in our previous paper [57], and except for 1 H calculations, picture change effects were found to be significant. Therefore they were included as implemented in ORCA and documented by Sandhoefer et al. [68]. The second-order DKH transformation is applied to the Fermi contact, spin-dipole and Zeeman operators. The transformation of the paramagnetic spin-orbit term was treated to first order only. Relevant experimental results are extremely scarce, but the approach was compared against [U(C 7 H 7 ) 2 ] , showing reasonable accuracy.
For the nuclei of interest, a Def2-TZVP basis was used with three additional uncontracted s orbitals with large exponents added to allow improved description of the electronic environment close to the nucleus, otherwise the basis used was as in the previous section. The exponents were derived using a geometric progression starting from the largest s function exponents of the existing basis, following the requirements of basis set design for hyperfine coupling calculations outlined by Chipman [69]. Because of the high computational cost, the calculation generally did not include spin-orbit corrections to the hyperfine calculation. We indicate in the text the two molecules for which this term was evaluated. In these cases the calculation used the same methods as detailed above for the other hyperfine terms but included solution of the relevant coupled-perturbed Kohn–Sham equations to obtain the derivative of the spin density, followed by contraction with the nuclear (spin)-electron (orbit) integrals.

3.4. Zeeman Coupling Matrix, g

It is possible to calculate g tensors with DFT and use linear response theory to include spin-orbit coupling as a first-order perturbation [17,70]; however, this approach is only effective for weak spin-orbit coupling and in the absence of low-lying excited states. For the molecules under investigation here neither condition is true, and the ground states are all clearly degenerate. As a result a multiconfigurational approach is necessary. ORCA 3.0.3 was used to perform a series of State-Averaged Complete Active-Space Self-Consistent Field (SA-CASSCF) calculations, using an active space of 3 electrons and 7 orbitals, i.e., the f electrons of the heavy element distributed in the relevant manifold of f orbitals. Relativistic effects were added using DKH2. Spin-orbit coupling was included by utilising the SOMF operator [71], using the centre of electronic charge as the origin of the molecule to minimise gauge dependence [72]. A SARC basis set was used for uranium, with a Def2-SVP for hydrogen and the remainder using a Def2-TZVP. The g tensor was evaluated within the effective Hamiltonian approach. The number of states included in the SA-CASSCF calculations were explored empirically, by increasing the number of states, until some convergence in the principal values of the g tensor was achieved [17,60].

3.5. The ZFS Tensor, D

For efficiency, the D tensor was calculated at the same level of theory as detailed in the previous section to obtain g .
We have written a Matlab script that uses the approaches detailed in Section 2.5 to model the effect of zero-field splitting on pNMR shifts. We have separated the results into contact and pseudocontact components as outlined in Section 2.6. This script was validated with test data (http://www.ens-lyon.fr/crmn/pnmr/wp-content/uploads/2013/12/2014-02-Mariapfarr-presentation-Juha-Vaara.pdf) and a generally available program (http://ja01.chem.buffalo.edu/ jochena/downloads/downloads.html) provided by Autschbach [10,54].

4. Results

4.1. 1 H pNMR Shifts of Neodymium Tetrapyridyl Appended Cyclen

The first molecule under investigation is a neodymium 1,4,7,10-tetrakis(2-pyridylmethyl)- 1,4,7,10-tetraazacyclododecane complex [ NdL py ] 3 + with a triflate CF 3 SO 3 counterion. This was synthesised and studied by Natrajan et al. [73], who produced 500Mz 1 H NMR spectra of the species in a D 2 O solvent at room temperature. The molecule is shown in Figure 1. The observed shifts were 0.2 , 2.4 , 4.6 , 5.2 , 6, 8, 8.3 , 8.4 , 10.8 and 12.4 ppm and were not assigned in the paper. In our previous work [57], we found that inclusion of the triflate counterion produced a significant change in the geometry of the complex and so needs to be included in modelling this structure.
The calculation of the g and D matrices used 21 quartet and 40 doublet states from a SA-CASSCF(3,7) calculation. The matrices are listed in Table S1 of the ESI. The resulting isotropic g value is negative, a phenomenon observed for states with significant spin-orbit coupling as discussed in reference [37,74].
We find in the absence of ZFS corrections from our calculations a range of shifts spanning 2.56 16.51 ppm. Inclusion of ZFS produces a wider spread of values which are compared in Figure 2. Table 5 outlines the details of the components.
It is noticeable that adding ZFS has very little effect on the contact part of the shielding, but the size (and composition) of the pseudocontact term significantly changes. Examining the pseudocontact components, the term derived from the isotropic hyperfine changes only slightly and the interaction between the anisotropic Zeeman and the spin-dipole hyperfine term is negligible. The change derives from the interaction of the spin-dipole part of the hyperfine and the isotropic g term (decomposed into the free electron part and the shift). While the predicted chemical shift changes by only a few ppm, it is sufficient to shift the predicted order of the peaks and as a result assignment of all ten 1 H shifts are not possible. We have omitted the spin-orbit correction to the hyperfine coupling matrix as the calculations are extremely demanding. However it is unlikely that these 1 H shifts will be strongly influenced by this omission since the protons are well removed from the heavy atom.

4.2. 1 H pNMR Shifts of Tris(2,6-bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridine) Uranium(III), U(MeBTP) 3 3 +

2,6-Bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridine (BTP) is an effective chelating agent and has been investigated as a reagent for processing waste containing uranium(III) species [4,75], with the family of ligands showing preference for coordinating selectively with trivalent actinides over trivalent lanthanides due to the “much softer character of the N atoms in this weakly basic ligand” [4]. It is thought that actinide/lanthanide selectivity results from the difference in the nature of the f orbitals, the lanthanide 4f orbitals displaying “core-like” behaviour and participating weakly in bonding, whereas the more extended 5f orbitals of actinide species participate more in covalent bonding [76]. As such there is much interest in better understanding the nature of this U-N bond. The experimental results that we are using for comparison are 200 MHz 1 H NMR spectra taken at 30 C in deuterated pyridine [4]. The molecule has four chemically equivalent types of hydrogen atoms. The numbering used to refer to these sets of protons is shown in Figure 3 and the stereochemistry in Figure 4. All hydrogen atoms are relatively distant from the U-N bond, and as such we expect that spin-orbit component of the hyperfine interaction is not significant, hence they have not been included here.
The SA-CASSCF(3,7) calculation of the g and D matrices used 14 quartet and 25 doublet states. The matrices are listed in Table S2 of the ESI. Our results are shown in Table 6 and Table 7. For each distinct type of hydrogen we show the effect of the ZFS correction. We list the contribution to the paramagnetic shift from each component of g · Z · A and a qualitative partition of contact and pseudocontact terms. Overall the paramagnetic shift is calculated to be slightly larger than the orbital contribution, except for the H ( 2 ) atoms. In the absence of the ZFS correction H ( 1 ) and H ( 3 ) proton chemical shifts are predicted to be larger in magnitude that the observed shifts, although the signs of the shifts are predicted correctly. The ZFS correction reduces the magnitudes in the direction of the experimental values. In this system including the ZFS produces shifts that are closer to the observed values but are still too large in magnitude. The experimental shifts span a range of about 52 ppm, while the neglect of ZFS altogether gives a span of calculated values of 117 ppm. The ZFS correction reduces the range to 90 ppm. This is shown in Figure 5. The contact and pseudocontact terms vary according to the method of partitioning them, but the contact shift is consistently predicted to be larger than the pseudocontact shift in all cases. However for H ( 2 ) and H ( 3 ) the values are much closer in magnitude.

4.3. 1 H and 13 C pNMR Shifts of Tris(pentamethylcyclopentadienyl) Uranium (III), U(C 5 Me 5 ) 3

It was originally thought impossible to synthesise an organometallic complex consisting of a single atom and three η 5 -bonded pentamethylcyclopentadienyl rings due to the steric hindrance presented by the methyl groups. However a samarium complex was synthesised in 1991 followed by a uranium complex shortly after. The optimised geometry of the complex is shown in Figure 6, with the methyl carbons residing above the plane of the ring due to steric effects. The paper reporting the synthesis [77] also reported the chemical shifts observed in 1 H and 13 C spectra taken at 500MHz, 25 C in a deuterated benzene solvent and these act as our baseline for comparison.
The SA-CASSCF(3,7) calculation of the g and D matrices used 28 quartet and 55 doublet states. The matrices are listed in Table S3 of the ESI. Our results for the methyl hydrogen atoms are shown in Table 8, again reflecting the treatment of ZFS, the individual contributions to the paramagnetic shift and the qualitative partition of contact and pseudocontact terms. Overall the paramagnetic shift is calculated to be significantly smaller (14–17%) than the orbital contribution.
In the absence of the ZFS correction, the predicted chemical shift was 3.02 ppm, with the correction the predicted chemical shift was 3.36 ppm, hence the full ZFS correction only contributes 0.34 ppm to the final computed shift. This is unlikely to be significant compared to the errors associated with the approximations used to generate the pNMR matrices, and this result remains a good match for the experimental value. Both the contact and pseudocontact shielding terms are positive, although the two decompositions used disagree with respect to the relative magnitudes of the contact and pseudocontact shifts. It must be noted that the hyperfine couplings have been evaluated using a DFT method. The choice of exchange-correlation functional can introduce a variation in the predicted hyperfine operators. It is difficult currently to deal with this property using purely wave functional techniques, which would in principle would allow a degree of systematic refinement. Given the size and complexity of the systems studied here, the approach is a suitable starting point. However, as wave function methods become more available for these operators further refinements of the computations can and should be sought [78].
Table 9 presents the shifts for 13 C pNMR shifts for the ring and methyl carbon atoms.
Due to the computational expense of calculating the spin-orbit corrections to the hyperfine coupling matrix and the fact that 13 C is a light atom, we initially did not apply this correction to our calculations. However this complex contains a set of η -coordinated carbon–uranium bonds, and the HALA effects alluded to earlier are likely to be in operation. Additionally, in the work of Autschbach et al. on actinide carbonate molecules [14,15], it was noted that the spin-orbit correction dominated the hyperfine interaction. Accordingly, we computed the spin-orbit correction to the hyperfine tensor using the methodology detailed in Section 3.3, and report our findings in Table 10. The format of this table is similar to the previous ones, but the individual contributions have been expanded. The table includes a comparison with the values before the spin-orbit correction to the hyperfine matrix was added and the bottom section highlights the two decomposition methods. Table 8 shows that the predicted isotropic chemical shielding (and chemical shift) of the two decomposition methods is different. This is due to the neglect of the O ( α 6 ) terms that method 2 uses, which is not valid when spin-orbit effects are significant. With reference to the method 1 decomposition, we note that the effect of the inclusion of the spin-orbit hyperfine term is less marked on the methyl carbon atoms than for the ring carbons. There is now qualitative agreement with experiment for the latter. For the methyl carbons, the effect is smaller, perhaps due to their being further removed from the heavy nucleus.
Figure 7 shows the deviation from the averaged value for the isotropic paramagnetic shielding of the ring and methyl 13 C nuclei. Overall the ring 13 C range from 63 – 92ppm and the methyl 13 C range from 74 – 94ppm. Note that the experimental values of the shifts lie within both ranges.

4.4. 29 Si pNMR Shifts of U(N(SiMe t Bu 2 ) 2 ) 3

Our next molecule is a single molecule magnet which, unusually, has a trigonal planar coordination around the uranium atom when in the solid state, instead of the more common trigonal pyramidal geometry that is seen in other three-coordinate actinide compounds. This is due to the large steric hindrance introduced by the t-butyl groups. Steric hindrance is a commonly used tool to prevent the molecules oligomerising, and there is much interest in understanding how the magnetic properties of the molecule are affected by the ligand. For comparison with experiment, we use the 73 Hz 29 Si chemical shift of 296.04 ppm taken at 25 C in a solution with deuterated benzene [2]. Our calculation refers to an optimised gas phase structure [57], and is shown in Figure 8.
The SA-CASSCF(3,7) calculation of the g and D matrices used 21 quartet and 39 doublet states. The matrices are listed in Table S4 of the ESI. Table 11 shows the resulting ZFS calculations. The orbital and paramagnetic contributions are predicted to have a similar size, and the contact interaction is positive and much larger than the pseudocontact, so the paramagnetic shift is largely due to increased spin density at the nucleus. Without applying ZFS, the chemical shift predicted was 18.54 ppm less than the experimental value. Using the full coupling method, the predicted chemical shift reduced by 46.19 ppm, which is further from the experimental value. Given our findings for the 13 C shifts in U(C 5 Me 5 ) 3 , the spin-orbit correction to the hyperfine tensor cannot be neglected. Due to computational cost we have only obtained the corrected tensors for two silicon nuclei. The hyperfine matrices of the two silicon atoms considered result in very different isotropic shielding, specifically 559.96 and 371.12 ppm. We note that in our optimised structure, the uranium-silicon distance of the first atom is 3.31 Angstroms whereas the second is 3.66 Angstroms. The magnitude of the change in the isotropic shielding for the silicon atoms that we have computed clearly indicates that it is essential to carry out the average over the full set of silicon atoms to obtain a reliable estimate of the shifts. This remains to be done and we are exploring more computationally efficient ways to obtain the spin-orbit hyperfine interaction.

4.5. 29 Si pNMR Shifts of Tris(trimethylsilyltetramethylcyclopentadienyl)uranium(lll), U(C 5 Me 4 SiMe 3 ) 3 .

Tris(trimethylsilyltetramethylcyclopentadienyl)uranium(lll) demonstrates even more steric crowding than the pentamethylcyclopentadienyl ring studied in Section 4.3, and was part of a systematic study of 29 Si shifts of similar compounds, focused on discovering trends [79] in the chemical shifts of uranium-silicon compounds. The SI for reference [79] includes a graph that may indicate very weak correlation between U-Si distance and 29 Si, but cautions against treating their results as definitive, citing the complexity involved when studying paramagnetic species and effects that are in opposition. That study reported an experimental value of 155 ppm ( 99.2 MHz spectra, taken at 298 K in deuterated benzene with an iodide counter anion). The structure is shown in Figure 9. The SA-CASSCF(3,7) calculation of the g and D matrices used 21 quartet and 39 doublet states. The matrices are listed in Table S5 of the ESI. Table 12 shows the resulting ZFS calculations. Our calculated value before splitting was + 18.92 ppm from the observed value, applying ZFS using the full coupling approach resulted in a small shift of + 1.73 ppm, so in this case ZFS has very little effect on the total shift. For silicon nuclei, the implication of our results presented in Section 4.3 is that spin-orbit effects must also be important here. Further investigation by including spin-orbit corrections to the hyperfine matrices is needed to clarify if this is the case. Although the situation is slightly different in that the silicon atoms are not directly bound to the heavy atom and this may ameliorate the relative importance of including the spin-orbit part of the hyperfine interaction.

4.6. 29 Si pNMR Shifts of the Mixed-Sandwich uranium(III) Compound, U( η -C 5 Me 4 Et)( η -(1,4-C 8 H 6 (Si i Pr 3 ) 2 )

Mixed-sandwich organo-uranium complexes have been used for more environmentally friendly synthesis of oxocarbo rings from carbon dioxide and carbon monoxide with the ability to reform the complex after use. Similar uses for olefin processing exist, and there is interest in understanding this process better [1]. The structure of the compound studied is shown in Figure 10, and the cited study reported an experimental value of 136.70 ppm ( 99.2 MHz spectra, taken at 298 K in deuterated benzene with an iodide counter anion).
The SA-CASSCF(3,7) calculation of the g and D matrices used 13 quartet and 32 doublet states. The matrices are listed in Table S6 of the ESI. Inspection of our results in Table 13 shows that the contact term is substantially larger than the pseudocontact, and as for the example in Section 4.5, we interpret this as implying a strong interaction through the C–Si bond.

4.7. Summary

Considering our six examples as a whole, Figure 11 shows the calculated values of pNMR shifts of the five uranium molecules plotted against the experimental values. Since this is a small sample, firm conclusions are not possible but, perhaps not surprisingly, the correlation coefficients improve as the treatment of ZFS improves; from R 2 = 0.86 for the treatment without ZFS to R 2 = 0.90 with of ZFS.

5. Conclusions

In this paper we have demonstrated the importance of including ZFS for assigning and predicting the pNMR chemical shifts of the lanthanide and actinide f 3 species. With the exception of the first example, our results are accurate enough to assign the observed shifts, and the first example demonstrates that applying ZFS can reorder shifts that are tightly grouped. However we have also demonstrated the importance of including spin-orbit coupling corrections to the hyperfine coupling matrix, especially when studying the chemical shifts of nuclei directly bonded to heavy nuclei. These necessary corrections are computationally demanding, and more work is needed to find theoretical and computational approaches that simplify applying this correction without sacrificing accuracy. The systems we have used as examples in this work typify the complexity of lanthanide and actinide containing molecules and the difficulties they pose for the computational predictions of their magnetic spectra. The modelling remains challenging and incomplete and requires further refinement as theoretical and computational progress is made.

Supplementary Materials

Calculated EPR data for the six molecules (Zeeman, ZFS and hyperfine parameters) is available online at https://www.mdpi.com/2312-7481/5/1/3/s1, Table S1: EPR data for [Nd LPy]3+. Atom numbering of protons refers to Figure 1 of the main text. Hyperfine values averaged over chemically equivalent atoms, Table S2: EPR data for tris(2,6-bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridine) uranium(III), U(MeBTP)3+. Atom numbering of protons refers to Figure 3 of the main text. Hyperfine values averaged over chemically equivalent atoms, Table S3: EPR data for tris(pentamethylcyclopentadienyl) uranium (III), U(C5Me5)3. Hyperfine values averaged over chemically equivalent atoms, Table S4: EPR data for U(N(SiMetBu2)2)3. Hyperfine values averaged over chemically equivalent atoms, Table S5: EPR data for Tris(trimethylsilyltetramethylcyclopentadienyl)uranium(lll), U(C5Me4SiMe3)3. Hyperfine values averaged over chemically equivalent atoms, Table S6: EPR data for U( η -C5Me4Et)( η -(1,4-C8H6(SiiPr3)2). Hyperfine values averaged over chemically equivalent atoms.

Author Contributions

Original draft preparation, A.W.L.; reviewing and editing, A.W.L. and J.J.W.M.; geometry optimisation, H.M.M.; pNMR matrix calibration and calculation, H.M.M.; application of ZFS, A.W.L.; spin-orbit corrected hyperfine calculations, J.J.W.M. and A.W.L.

Funding

This research received no external funding.

Acknowledgments

H.M.M. thanks the Nuclear FiRST Doctoral Training Center and the EPRSC for the award of a studentship during which some of this work was carried out. The authors also thank IT services (University of Manchester) for the provision of time on the Computational Shared Facility (CSF).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
CASSCFComplete active space self-consistent field theory
CIConfiguration interaction
DFTDensity functional theory
DKHDouglas-Kroll-Hess
GIAOGauge including atomic orbitals
HALAHeavy atom-light atom
LRTLinear response theory
PBE0Hybrid Perdew–Burke-Ernzerhof (PBE) functional: Exchange energy 75% functional, 25% Hartree-Fock value.
QDPTQuasi degenerate perturbation theory
SARCScalar all-electron relativistic ( basis sets ).
SOMFSpin-orbit, mean field
SOSSum over states
ZFSZero Field Splitting
ZORAZeroth-order regular approximation

Appendix A. Relaxing the Weak Spin-Orbit Coupling Assumption

We outline Van den Heuvel and Soncini’s [52] improved approach to zero-field splitting, applicable for systems with large spin-orbit coupling. The set of spin tensors { S q k } form a complete basis for C ( 2 S + 1 ) × ( 2 S + 1 ) , so it follows that any 2 S + 1 by 2 S + 1 matrix X can be expressed as a linear combination of S q k matrices [12,81]
X = k = 0 2 S q = k k ( 1 ) q X q k S q k , where X q k = 2 k + 1 S | | S k | | S 2 Trace ( S q k X ) .
Ignoring zero-field splitting, the pNMR contribution to the shift is given by
1 1 σ κ υ P = μ B g nuc k B T · 1 2 S + 1 k = 0 2 S q = k k g q κ k A q υ k * S | | S k | | S 2 2 k + 1
ZFS is accounted for by considering a projection P of the ZFS Hamiltonian
P H Z F S = k k S q = 2 k 2 k ( 1 ) q D q 2 k S q 2 k where P = q = | S , S q = | S , S q q |
again the procedure is to diagonalise the ( 2 S + 1 ) × ( 2 S + 1 ) matrix S · D · S . The resulting 2 S + 1 eigenvectors grouped into n λ groups of degenerate | λ state with energy E λ can be used to calculate the paramagnetic shift as in the equation below
σ κ υ p = μ B ω g N μ N k , k = 0 2 S q = k k q = k k g q κ k λ , λ Q k q k q ( λ , λ ) a q κ k
where
1 1 Q k q k q ( λ , λ ) = e E n / k B T k T { λ } , { λ } λ | S q k | λ λ | S q k | λ E m = E n 2 e E n / k B T E n E m e { λ } , { λ } λ | S q k | λ λ | S q k | λ otherwise
The main barrier to using Equation (16) is that the A q k terms are derived from generalised hyperfine property tensors that are currently not generally available.

References

  1. Tsoureas, N.; Summerscales, O.T.; Cloke, F.G.N.; Roe, S.M. Steric effects in the reductive coupling of CO by mixed-sandwich uranium(III) complexes. Organometallics 2013, 32, 1352–1362. [Google Scholar] [CrossRef]
  2. Goodwin, C.A.; Tuna, F.; McInnes, E.J.; Liddle, S.T.; McMaster, J.; Vitorica-Yrezabal, I.J.; Mills, D.P. A structurally authenticated trigonal planar actinide complex. Chem. A Eur. J. 2014, 20, 14579–14583. [Google Scholar] [CrossRef] [PubMed]
  3. Parsonnet, V.; Driller, J.; Cook, D.; Rizvi, S.A. Thirty-one years of clinical experience with nuclear-powered pacemakers. PACE Pacing Clin. Electrophysiol. 2006, 29, 195–200. [Google Scholar] [CrossRef] [PubMed]
  4. Iveson, P.B.; Rivière, C.; Guillaneux, D.; Nierlich, M.; Thuéry, P.; Ephritikhine, M.; Madic, C. Selective complexation of uranium(III) over cerium(III) by 2,6-bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridines: 1H NMR and X-ray crystallography studies. Chem. Commun. 2001, 1512–1513. [Google Scholar] [CrossRef]
  5. Helgaker, T.; Jaszunski, M.; Ruud, K. Ab initio methods for the calculation of NMR shielding and indirect spin-spin coupling constants. Chem. Rev. 1999, 99, 293–352. [Google Scholar] [CrossRef] [PubMed]
  6. Autschbach, J. Relativistic calculations of magnetic resonance parameters: background and some recent developments. Philos. Trans. A. Math. Phys. Eng. Sci. 2014, 372, 20120489. [Google Scholar] [CrossRef] [PubMed]
  7. Kaupp, M.; Köhler, F.H. Combining NMR spectroscopy and quantum chemistry as tools to quantify spin density distributions in molecular magnetic compounds. Coord. Chem. Rev. 2009, 253, 2376–2386. [Google Scholar] [CrossRef]
  8. Patchkovskii, S.; Autschbach, J.; Ziegler, T. Curing difficult cases in magnetic properties prediction with self-interaction corrected density functional theory. J. Chem. Phys. 2001, 115, 26–42. [Google Scholar] [CrossRef]
  9. Bertini, I.; Luchinat, C.; Parigi, G. Magnetic susceptibility in paramagnetic NMR. Prog. Nucl. Magn. Reson. Spectrosc. 2002, 40, 249–273. [Google Scholar] [CrossRef]
  10. Martin, B.; Autschbach, J. Temperature dependence of contact and dipolar NMR chemical shifts in paramagnetic molecules. J. Chem. Phys. 2015, 142, 054108. [Google Scholar] [CrossRef] [Green Version]
  11. Vaara, J.; Rouf, S.A.; Mares, J. Magnetic couplings in the chemical shift of paramagnetic NMR. J. Chem. Theory Comput. 2015, 11, 4840–4849. [Google Scholar] [CrossRef] [PubMed]
  12. Van den Heuvel, W.; Soncini, A. NMR chemical shift as analytical derivative of the Helmholtz free energy. J. Chem. Phys. 2013, 138, 054113. [Google Scholar] [CrossRef] [PubMed]
  13. Soncini, A.; Van den Heuvel, W. Communication: Paramagnetic NMR chemical shift in a spin state subject to zero-field splitting. J. Chem. Phys. 2013, 138, 021103. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Gendron, F.; Sharkas, K.; Autschbach, J. Calculating NMR chemical shifts for paramagnetic metal complexes from first-principles. J. Phys. Chem. Lett. 2015, 6, 2183–2188. [Google Scholar] [CrossRef]
  15. Gendron, F.; Autschbach, J. Ligand NMR Chemical Shift Calculations for Paramagnetic Metal Complexes: 5f1 vs 5f2 Actinides. J. Chem. Theory Comput. 2016, 12, 5309–5321. [Google Scholar] [CrossRef]
  16. Vaara, J. Chemical shift in paramagnetic systems. Sci. Technol. At. Mol. Condens. Matter Biol. Syst. 2013, 3, 41–67. [Google Scholar] [CrossRef]
  17. Bolvin, H.; Autschbach, J. Relativistic methods for calculating electron paramagnetic resonance (EPR) parameters. In Handbook of Relativistic Quantum Chemistry; Liu, W., Ed.; Springer: Berlin/Heidelberg, Germany, 2017; pp. 725–763. [Google Scholar]
  18. Autschbach, J. NMR calculations for paramagnetic molecules and metal complexes. Annu. Rep. Comput. Chem. 2015. [Google Scholar] [CrossRef]
  19. Pryce, M.H.L. A modified perturbation procedure for a problem in paramagnetism. Proc. Phys. Soc. Sect. A 1950, 63, 25–29. [Google Scholar] [CrossRef]
  20. Abragam, A.; Pryce, M.H.L. Theory of the nuclear hyperfine structure of paramagnetic resonance spectra in crystals. Proc. R. Soc. A Math. Phys. Eng. Sci. 1951, 205, 135–153. [Google Scholar] [CrossRef]
  21. McConnell, H.M.; Chesnut, D.B. Theory of isotropic hyperfine interactions in π-electron radicals. J. Chem. Phys. 1958, 28, 107–117. [Google Scholar] [CrossRef]
  22. Pennanen, T.O.; Vaara, J. Nuclear magnetic resonance chemical shift in an arbitrary electronic spin state. Phys. Rev. Lett. 2008, 100, 133002. [Google Scholar] [CrossRef]
  23. Ramsey, N.F. Magnetic shielding of nuclei in molecules. Physica 1951, 17, 303–307. [Google Scholar] [CrossRef]
  24. Facelli, J.C. Chemical shift tensors: Theory and application to molecular structural problems. Prog. Nucl. Magn. Reson. Spectrosc. 2011, 58, 176–201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Ditchfield, R. Self-consistent perturbation theory of diamagnetism. Mol. Phys. 1974, 27, 789–807. [Google Scholar] [CrossRef]
  26. Schreckenbach, G.; Ziegler, T. Calculation of NMR shielding tensors using gauge-including atomic orbitals and modern density functional theory. J. Phys. Chem. 1995, 99, 606–611. [Google Scholar] [CrossRef]
  27. McDouall, J.J.W. Computational Quantum Chemistry; Theoretical and Computational Chemistry Series; Royal Society of Chemistry: Cambridge, UK, 2013. [Google Scholar]
  28. Filatov, M. On representation of the Hamiltonian matrix elements in relativistic regular approximation. Chem. Phys. Lett. 2002, 365, 222–231. [Google Scholar] [CrossRef] [Green Version]
  29. Reiher, M. Douglas-Kroll-Hess theory: A relativistic electrons-only theory for chemistry. Theor. Chem. Acc. 2006, 116, 241–252. [Google Scholar] [CrossRef]
  30. Autschbach, J.; Govind, N.; Atta-Fynn, R.; Bylaska, E.J.; Weare, J.W.; de Jong, W.A. Computational Tools for Predictive Modeling of Properties in Complex Actinide Systems. In Computational Methods in Lanthanide and Actinide Chemistry; John Wiley & Sons Ltd.: Chichester, UK, 2015; pp. 299–342. [Google Scholar]
  31. Marian, C.M. Spin-orbit coupling in molecules. In Reviews in Computational Chemistry; Wiley-Blackwell: Hoboken, NJ, USA, 2001; pp. 99–204. [Google Scholar]
  32. Neese, F. Prediction of electron paramagnetic resonance g values using coupled perturbed Hartree–Fock and Kohn–Sham theory. J. Chem. Phys. 2001, 115, 11080–11096. [Google Scholar] [CrossRef]
  33. Kleinschmidt, M.; Tatchen, J.; Marian, C.M. SPOCK.CI: A multireference spin-orbit configuration interaction method for large molecules. J. Chem. Phys. 2006, 124, 124101. [Google Scholar] [CrossRef]
  34. Klein, M.J. On a degeneracy theorem of Kramers. Am. J. Phys. 1952. [Google Scholar] [CrossRef]
  35. Gerloch, M.; McMeeking, R.F. Paramagnetic properties of unsymmetrical transition-metal complexes. J. Chem. Soc. Dalt. Trans. 1975, 0, 2443. [Google Scholar] [CrossRef]
  36. Bolvin, H. An alternative approach to the g-matrix: Theory and applications. ChemPhysChem 2006, 7, 1575–1589. [Google Scholar] [CrossRef] [PubMed]
  37. Chibotaru, L.F. Ab initio methodology for pseudospin Hamiltonians of anisotropic magnetic complexes. In Advances in Chemical Physics Volume 153; Wiley-Blackwell: Hoboken, NJ, USA, 2013; pp. 397–519. [Google Scholar]
  38. Koseki, S.; Schmidt, M.W.; Gordon, M.S. MCSCF/6-31G(d,p) calculations of one-electron spin-orbit coupling constants in diatomic molecules. J. Phys. Chem. 1992, 96, 10768–10772. [Google Scholar] [CrossRef] [Green Version]
  39. Koseki, S.; Gordon, M.S.; Schmidt, M.W.; Matsunaga, N. Main group effective nuclear charges for spin-orbit calculations. J. Phys. Chem. 1995, 99, 12764–12772. [Google Scholar] [CrossRef]
  40. Koseki, S.; Schmidt, M.W.; Gordon, M.S. Effective nuclear charges for the first- through third-row transition metal elements in spin-orbit calculations. J. Phys. Chem. A 1998, 102, 10430–10435. [Google Scholar] [CrossRef]
  41. Heß, B.A.; Marian, C.M.; Wahlgren, U.; Gropen, O. A mean-field spin-orbit method applicable to correlated wavefunctions. Chem. Phys. Lett. 1996, 251, 365–371. [Google Scholar] [CrossRef]
  42. Vahtras, O.; Engström, M.; Schimmelpfennig, B. Electronic g-tensors obtained with the mean-field spin-orbit Hamiltonian. Chem. Phys. Lett. 2002, 351, 424–430. [Google Scholar] [CrossRef]
  43. Neese, F. Calculation of the zero-field splitting tensor on the basis of hybrid density functional and Hartree-Fock theory. J. Chem. Phys. 2007, 127, 164112. [Google Scholar] [CrossRef]
  44. Neese, F. Analytic derivative calculation of electronic g-tensors based on multireference configuration interaction wavefunctions. Mol. Phys. 2007, 105, 2507–2514. [Google Scholar] [CrossRef]
  45. Neese, F. ORCA—An Ab Initio, DFT and Semiempirical SCF-MO Package (Manual), 3.0.3 ed.; Max-Planck-Institute for Chemical Energy Conversion: Mülheim an der Ruhr, Germany, 2013; p. 595. [Google Scholar]
  46. Malkin, E.; Repiský, M.; Komorovský, S.; Mach, P.; Malkina, O.L.; Malkin, V.G. Effects of finite size nuclei in relativistic four-component calculations of hyperfine structure. J. Chem. Phys. 2011, 134, 044111. [Google Scholar] [CrossRef]
  47. Rouf, S.A.; Mares, J.; Vaara, J. Relativistic Approximations to Paramagnetic NMR Chemical Shift and Shielding Anisotropy in Transition Metal Systems. J. Chem. Theory Comput. 2017, 13, 3731–3745. [Google Scholar] [CrossRef] [PubMed]
  48. Autschbach, J. Relativistic effects on NMR. Sci. Technol. At. Mol. Condens. Matter Biol. Syst. 2013, 3, 69–117. [Google Scholar] [CrossRef]
  49. Ramsey, N.F. The internal diamagnetic field correction in measurements of the proton magnetic moment. Phys. Rev. 1950, 77. [Google Scholar] [CrossRef]
  50. Moon, S.; Patchkovskii, S. First-principles calculations of paramagnetic NMR shifts. In Calculation of NMR and EPR Parameters: Theory and Applications; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2004; pp. 325–338. [Google Scholar]
  51. Kurland, R.J.; McGarvey, B.R. Isotropic NMR shifts in transition metal complexes: The calculation of the fermi contact and pseudocontact terms. J. Magn. Reson. 1970, 2, 286–301. [Google Scholar] [CrossRef]
  52. Van den Heuvel, W.; Soncini, A. NMR chemical shift in an electronic state with arbitrary degeneracy. Phys. Rev. Lett. 2012, 109, 073001. [Google Scholar] [CrossRef] [PubMed]
  53. Pennanen, T.O.; Vaara, J. Density-functional calculations of relativistic spin-orbit effects on nuclear magnetic shielding in paramagnetic molecules. J. Chem. Phys. 2005, 123, 174102. [Google Scholar] [CrossRef] [PubMed]
  54. Autschbach, J.; Patchkovskii, S.; Pritchard, B. Calculation of hyperfine tensors and paramagnetic NMR shifts using the relativistic zeroth-order regular approximation and density functional theory. J. Chem. Theory Comput. 2011, 7, 2175–2188. [Google Scholar] [CrossRef]
  55. Rouf, S.A.; Mares, J.; Vaara, J. 1H chemical shifts in paramagnetic Co(II) pyrazolylborate complexes: A first-principles study. J. Chem. Theory Comput. 2015, 11, 1683–1691. [Google Scholar] [CrossRef]
  56. Liimatainen, H.; Pennanen, T.O.; Vaara, J. 1H chemical shifts in nonaxial, paramagnetic chromium(III) complexes—Application of novel pNMR shift theory. Can. J. Chem. 2009, 87, 954–964. [Google Scholar] [CrossRef]
  57. Moylan, H.M.; McDouall, J.J.W. Computational study of chemical shifts in paramagnetic f-element compounds. Trends Phys. Chem. 2018, 17, 25–57. [Google Scholar]
  58. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 09, Revision B.01; Gaussian, Inc.: Wallingford, CT, USA, 2010. [Google Scholar]
  59. Adamo, C.; Barone, V. Toward reliable density functional methods without adjustable parameters: The PBE0 model. J. Chem. Phys. 1999, 110, 6158. [Google Scholar] [CrossRef]
  60. Moylan, H.M.; McDouall, J.J.W. Electronic g tensors in U(V) complexes—A computational study. Chem. A Eur. J. 2017, 23, 7798–7808. [Google Scholar] [CrossRef] [PubMed]
  61. Pantazis, D.A.; Neese, F. All-electron scalar relativistic basis sets for the lanthanides. J. Chem. Theory Comput. 2009, 5, 2229–2238. [Google Scholar] [CrossRef]
  62. Pantazis, D.A.; Neese, F. All-electron scalar relativistic basis sets for the actinides. J. Chem. Theory Comput. 2011, 7, 677–684. [Google Scholar] [CrossRef]
  63. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297. [Google Scholar] [CrossRef] [PubMed]
  64. Wolinski, K.; Hinton, J.F.; Pulay, P. Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations. J. Am. Chem. Soc. 1990, 112, 8251–8260. [Google Scholar] [CrossRef]
  65. Neese, F. The ORCA program system. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 73–78. [Google Scholar] [CrossRef]
  66. Hess, B.A. Applicability of the no-pair equation with free-particle projection operators to atomic and molecular structure calculations. Phys. Rev. A 1985, 32, 756–763. [Google Scholar] [CrossRef]
  67. Jansen, G.; Hess, B.A. Revision of the Douglas-Kroll transformation. Phys. Rev. A 1989. [Google Scholar] [CrossRef]
  68. Sandhoefer, B.; Kossmann, S.; Neese, F. 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. J. Chem. Phys. 2013, 138, 104102. [Google Scholar] [CrossRef]
  69. Chipman, D.M. The spin polarization model for hyperfine coupling constants. Theor. Chim. Acta 1992, 82, 93–115. [Google Scholar] [CrossRef]
  70. Schreckenbach, G.; Ziegler, T. Calculation of the g tensor of electron paramagnetic resonance spectroscopy using gauge-including atomic orbitals and density functional theory. J. Phys. Chem. A 1997, 101, 3388–3399. [Google Scholar] [CrossRef]
  71. Neese, F. Efficient and accurate approximations to the molecular spin-orbit coupling operator and their use in molecular g-tensor calculations. J. Chem. Phys. 2005, 122. [Google Scholar] [CrossRef] [PubMed]
  72. Luzanov, A.; Babich, E.; Ivanov, V. Gauge-invariant calculations of magnetic properties in semiempirical approaches. Application to full-CI π-electron scheme. J. Mol. Struct. 1994, 311, 211–220. [Google Scholar] [CrossRef]
  73. Natrajan, L.S.; Khoabane, N.M.; Dadds, B.L.; Muryn, C.A.; Pritchard, R.G.; Heath, S.L.; Kenwright, A.M.; Kuprov, I.; Faulkner, S. Probing the structure, conformation, and stereochemical exchange in a family of lanthanide complexes Derived from Tetrapyridyl-Appended Cyclen. Inorg. Chem. 2010, 49, 7700–7709. [Google Scholar] [CrossRef] [PubMed]
  74. Chibotaru, L.F.; Ungur, L. Negative g factors, Berry phases, and magnetic properties of complexes. Phys. Rev. Lett. 2012, 109, 246403. [Google Scholar] [CrossRef] [PubMed]
  75. Kolarik, Z.; Mullich, U.; Gassner, F. Extraction of Am(III) and Eu(III) nitrates by 2-6-di-(5,6-dipropyl-1,2,4-triazin-3-yl)pyridines. Solvent Extr. Ion Exch. 1999, 17, 1155–1170. [Google Scholar] [CrossRef]
  76. Cotton, S.A. Lanthanide and Actinide Chemistry, 2nd ed.; Wiley-Blackwell: Hoboken, NJ, USA, 2006; p. 263. [Google Scholar]
  77. Evans, W.J.; Forrestal, K.; Ziller, J. Activity of [Sm(C5Me5)3] in ethylene polymerization and synthesis of [U(C5Me5)3], the first tris(pentamethylcyclopentadienyl) 5f-element complex. Angew. Chem. Int. Ed. 1997, 36, 774–776. [Google Scholar] [CrossRef]
  78. Autschbach, J.; Ziegler, T. Double perturbation theory: a powerful tool in computational coordination chemistry. Coord. Chem. Rev. 2003, 238–239, 83–126. [Google Scholar] [CrossRef]
  79. Windorff, C.J.; Evans, W.J. 29Si NMR spectra of silicon-containing uranium complexes. Organometallics 2014, 33, 3786–3791. [Google Scholar] [CrossRef]
  80. Siladke, N.A.; Ziller, J.W.; Evans, W.J. Synthesis, structure, and reactivity of the sterically crowded complex (C5Me4SiMe3)3U. Z. Anorg. Allg. Chem. 2010. [Google Scholar] [CrossRef]
  81. Griffith, J. Some investigations in the theory of open-shell ions. Mol. Phys. 1960, 3, 477–483. [Google Scholar] [CrossRef]
Figure 1. Structure of the cyclen complex [ NdL py ] 3 + . Labels refer to the 10 chemically distinct H atoms, with 1,3,5 representing equatorial CH 2 protons and 2,4,6 representing axial protons. Colours indicate atom type (grey C, blue N, red O, yellow S, light blue F).
Figure 1. Structure of the cyclen complex [ NdL py ] 3 + . Labels refer to the 10 chemically distinct H atoms, with 1,3,5 representing equatorial CH 2 protons and 2,4,6 representing axial protons. Colours indicate atom type (grey C, blue N, red O, yellow S, light blue F).
Magnetochemistry 05 00003 g001
Figure 2. Experimental and computed pNMR spectra for [ NdL py ] 3 + with and without ZFS at 298.15K. For numbering of protons see Figure 1.
Figure 2. Experimental and computed pNMR spectra for [ NdL py ] 3 + with and without ZFS at 298.15K. For numbering of protons see Figure 1.
Magnetochemistry 05 00003 g002
Figure 3. (a) The MeBTP ligand, labels identify the four types of chemically equivalent 1 H, (b) The full structure of Tris(2,6-bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridine) uranium(III).
Figure 3. (a) The MeBTP ligand, labels identify the four types of chemically equivalent 1 H, (b) The full structure of Tris(2,6-bis(5,6-dialkyl-1,2,4-triazin-3-yl)pyridine) uranium(III).
Magnetochemistry 05 00003 g003
Figure 4. Stereochemistry of U(MeBTP) 3 3 + , showing the three planar ligands, staggered with approximate 3-fold rotational symmetry.
Figure 4. Stereochemistry of U(MeBTP) 3 3 + , showing the three planar ligands, staggered with approximate 3-fold rotational symmetry.
Magnetochemistry 05 00003 g004
Figure 5. Experimental and predicted pNMR shifts of U(MeBTP) 3 3 + with and without ZFS at 303K. For numbering of protons see Figure 3. The order is consistent in all cases.
Figure 5. Experimental and predicted pNMR shifts of U(MeBTP) 3 3 + with and without ZFS at 303K. For numbering of protons see Figure 3. The order is consistent in all cases.
Magnetochemistry 05 00003 g005
Figure 6. The structure and stereochemistry of U(C 5 Me 5 ) 3 .
Figure 6. The structure and stereochemistry of U(C 5 Me 5 ) 3 .
Magnetochemistry 05 00003 g006
Figure 7. Difference of the predicted isotropic shielding of individual 13 C nuclei from the averaged value at 298 K.
Figure 7. Difference of the predicted isotropic shielding of individual 13 C nuclei from the averaged value at 298 K.
Magnetochemistry 05 00003 g007
Figure 8. (a) Chemical structure of U(N(SiMe t Bu 2 ) 2 ) 3 , (b) Stereochemistry of U(N(SiMe t Bu 2 ) 2 ) 3 with uranium and silicon nuclei shown as spheres and the remaining non-hydrogen nuclei in wireframe. Stated distances are in Angstroms.
Figure 8. (a) Chemical structure of U(N(SiMe t Bu 2 ) 2 ) 3 , (b) Stereochemistry of U(N(SiMe t Bu 2 ) 2 ) 3 with uranium and silicon nuclei shown as spheres and the remaining non-hydrogen nuclei in wireframe. Stated distances are in Angstroms.
Magnetochemistry 05 00003 g008
Figure 9. Tris(trimethylsilylcyclopentadienyl)uranium(lll), U(C 5 Me 4 SiMe 3 ) 3 .
Figure 9. Tris(trimethylsilylcyclopentadienyl)uranium(lll), U(C 5 Me 4 SiMe 3 ) 3 .
Magnetochemistry 05 00003 g009
Figure 10. U( η -C 5 Me 4 Et)( η -(1,4-C 8 H 6 (Si i Pr 3 ) 2 ).
Figure 10. U( η -C 5 Me 4 Et)( η -(1,4-C 8 H 6 (Si i Pr 3 ) 2 ).
Magnetochemistry 05 00003 g010
Figure 11. Calculated results against experimental results.
Figure 11. Calculated results against experimental results.
Magnetochemistry 05 00003 g011
Table 1. Decomposition of the components of the Zeeman and hyperfine coupling matrices by isotropic and anisotropic components, I 3 indicates a 3 × 3 identity matrix.
Table 1. Decomposition of the components of the Zeeman and hyperfine coupling matrices by isotropic and anisotropic components, I 3 indicates a 3 × 3 identity matrix.
Matrix O(1) ComponentO( α 2 ) ComponentsO( α 4 ) Components
g g e I 3 Δ g iso I 3 , Δ g ˜
A A iso I 3 , A SD A iso Δ SO I 3 , A Δ SO , SYM , A Δ SO , AS
Table 2. Classification of components of the g · A matrix.
Table 2. Classification of components of the g · A matrix.
HyperfineZeemanOrder in α Decomposition
Method 1Method 2
Isotropic A iso Free electron g e 2ndContactContact 1
Scalar Iso. g-shift Δ g iso 4thContactContact 1
Relativistic Aniso. g-shift Δ g ˜ 4thContactPseudocontact 2
Isotropic A iso Δ SO Free electron g e 4thContactContact 1
SO Correction Iso. g-shift Δ g iso 6thContact-
Aniso. g-shift Δ g ˜ 6thContact-
Anisotropic A SD Free electron g e 2ndPseudocontactPseudocontact 2
Scalar Iso. g-shift Δ g iso 4thPseudocontactPseudocontact 2
Relativistic Aniso. g-shift Δ g ˜ 4thPseudocontactPseudocontact 3
Symmetric A Δ SO , SYM Free electron g e 4thPseudocontactPseudocontact 2
Anisotropic Iso. g-shift Δ g iso 6thPseudocontact-
SO Correction Aniso. g-shift Δ g ˜ 6thPseudocontact-
Asymmetric A Δ SO , AS Free electron g e 4thPseudocontactPseudocontact
Anisotropic Iso. g-shift Δ g iso 6thPseudocontact-
SO Correction Aniso. g-shift Δ g ˜ 6thPseudocontact-
1 For a doublet is a pure contact term, relates isotropic spin density at the atom [7]. 2 For a doublet is a pure spin-dipolar term. 3 For a doublet relates to spin density in p-type orbitals.
Table 3. One of the Kramers pair wave functions of U(MeBTP) 3 3 + , spin-orbit corrected, SA-CASSCF(3,7) wave function, with 14 quartet and 25 doublet states. The second state can be inferred, via Kramers degeneracy, exchanging the weights of each + M S state with that of the corresponding M S state. States with weight less than 0.01 are omitted.
Table 3. One of the Kramers pair wave functions of U(MeBTP) 3 3 + , spin-orbit corrected, SA-CASSCF(3,7) wave function, with 14 quartet and 25 doublet states. The second state can be inferred, via Kramers degeneracy, exchanging the weights of each + M S state with that of the corresponding M S state. States with weight less than 0.01 are omitted.
Quartet ( a ) Doublet ( b )
| 3 2 , 3 2 | 3 2 , 1 2 | 3 2 , 1 2 | 3 2 , 3 2 | 1 2 , 1 2 | 1 2 , 1 2
WeightRootWeightRootWeightRootWeightRootWeightRootWeightRoot
0.07610.06820.02110.04810.01250.0113
0.04820.02030.02420.08930.01770.0196
0.06430.04240.01930.02140.030100.0129
0.01850.01490.0515 0.01710
0.01760.012120.01111
( a ) Energies of first 12 quartet spin-free SA-CASSCF states, measured from lowest root are (1) 0 cm 1 , (2) 447 cm 1 , (3) 448.1 cm 1 , (4) 671.4 cm 1 , (5) 684.7 cm 1 , (6) 773.7 cm 1 , (7) 1457.9 cm 1 , (8) 1712.7 cm 1 , (9) 1784.5 cm 1 , (10) 1985.8 cm 1 , (11) 2121.2 cm 1 , (12) 2164.9 cm 1 . ( b ) Energies of first 10 doublet states, measured from lowest root in (a) are (1) 2971.7 cm 1 , (2) 3022.4 cm 1 , (3) 3079.2 cm 1 , (4) 3124.8 cm 1 , (5) 3766.2 cm 1 , (6) 3788.9 cm 1 , (7) 4100.9 cm 1 , (8) 4174.8 cm 1 , (9) 4421.0 cm 1 , (10) 4439.1 cm 1 .
Table 4. Relative energies (cm 1 ) for the 10 lowest spin-orbit states of the six compounds studied here, showing the increased gap between first and second excited pairs compared to that of the first pair and the ground state pair.
Table 4. Relative energies (cm 1 ) for the 10 lowest spin-orbit states of the six compounds studied here, showing the increased gap between first and second excited pairs compared to that of the first pair and the ground state pair.
States[ NdL py ] 3 + U(MeBTP) 3 3 + U(C 5 Me 5 ) 3 U(N(SiMe t Bu 2 ) 2 ) 3 U(C 5 Me 4 SiMe 3 ) 3 U( η -C 5 Me 4 Et)
* ( η -(1,4-C 8 H 6 (Si i Pr 3 ) 2 )
1–2 000  0  0     0
3–441  4458019061    374
5–6139 2601365432144    841
7–8218 387249511011034    1190
9–10265 464255212691183    1255
Table 5. Calculated 1 H pNMR shifts of [ NdL py ] 3 + organised in increasing δ . For numbering of protons see Figure 1. The reference chemical shift is 31 . 58 ppm (TMS) and the temperature is 298.15K.
Table 5. Calculated 1 H pNMR shifts of [ NdL py ] 3 + organised in increasing δ . For numbering of protons see Figure 1. The reference chemical shift is 31 . 58 ppm (TMS) and the temperature is 298.15K.
Component H ( 9 ) H ( 7 ) H ( 6 ) H ( 4 ) H ( 1 )
A g Order No ZFS ( a ) ZFS ( b ) No ZFSZFSNo ZFSZFSNo ZFSZFSNo ZFSZFS
A iso g e 2 −0.07−0.070.440.443.633.622.512.511.041.04
Δ g iso 4 0.120.12−0.70−0.70−5.77−5.76−3.99−3.99−1.66−1.66
Δ g ˜ 4 0.000.000.000.000.00−0.010.00-0.010.000.00
A SD g e 2 0.00−10.860.003.330.002.780.000.100.000.97
Δ g iso 4 0.0017.270.00−5.300.00−4.420.00−0.170.00−1.54
Δ g ˜ 4 6.416.60−1.95−2.01−1.68-1.720.010.00−0.56−0.57
Orbital Contribution σ orb 27.5928.4028.1123.6623.00
ContactMethod 1 ( c ) σ c 0.040.04−0.26−0.26−2.14−2.15−1.48−1.49−0.62−0.61
Method 2 0.040.04−0.26−0.26−2.14−2.15−1.48−1.48−0.62−0.61
PseudocontactMethod 1 σ pc 6.4113.01−1.95−3.98−1.68−3.370.01−0.07−0.56−1.14
Method 2 6.4113.00−1.95−3.98−1.68−3.370.01−0.08−0.56−1.14
Total isotropic shielding σ iso = σ orb + σ c + σ pc 34.0540.6426.1824.1624.2922.5922.1922.1221.8221.24
Delta σ ref σ iso −2.46−9.065.407.427.298.999.399.469.7610.35
Component H ( 3 ) H ( 10 ) H ( 2 ) H ( 5 ) H ( 8 )
A g Order No ZFSZFSNo ZFSZFSNo ZFSZFSNo ZFSZFSNo ZFSZFS
A iso g e 2 9.869.8420.5120.493.283.28−2.90−2.8925.9625.93
Δ g iso 4 −15.67−15.65−32.62−32.58−5.22−5.214.614.60−41.29−41.24
Δ g ˜ 4 0.00−0.020.00−0.050.00−0.010.000.010.00−0.06
A SD g e 2 0.00−2.270.00−2.410.002.260.008.470.00−2.02
Δ g iso 4 0 0.003.610.003.830.00−3.600.00−13.470.003.22
Δ g ˜ 4 1.811.831.421.46−1.36−1.39−5.16−5.301.201.24
Orbital Contribution σ orb 23.1429.0423.7727.3529.11
Contact Method 1 σ c −5.81−5.83−12.11−12.14−1.94−1.941.711.72−15.33−15.38
Method 2 −5.81−5.81−12.11−12.09−1.94−1.931.711.71−15.33−15.31
PseudocontactMethod 1 σ pc 1.813.171.422.88−1.36−2.74−5.16−10.301.202.41
Method 2 1.813.151.422.83−1.36−2.74−5.16−10.291.202.38
Total isotropic shielding σ iso = σ orb + σ c + σ pc 19.1320.4818.3519.7820.4819.1023.9018.7714.9816.18
Delta σ ref σ iso 12.4511.1013.2311.8011.1012.487.6812.8116.6015.41
( a ) The No ZFS approach uses Equation (11) with Z = 1 3 S ( S + 1 ) I , ( b ) The full coupling approach uses Equation (11) with Z defined in Equations 12 and (14), ( c ) Methods 1 and 2 are outlined in Table 2, ( d ) The observed shifts were 0.2 , 2.4 , 4.6 , 5.2 , 6, 8, 8.3 , 8.4 , 10.8 and 12.4 ppm and were not assigned in the paper.
Table 6. Calculated 1 H pNMR shifts for the methyl H atoms of U(MeBTP) 3 3 + in ppm at 303K). For numbering of protons see Figure 3. The final, ZFS corrected, computational values are highlighted.
Table 6. Calculated 1 H pNMR shifts for the methyl H atoms of U(MeBTP) 3 3 + in ppm at 303K). For numbering of protons see Figure 3. The final, ZFS corrected, computational values are highlighted.
Component Methyl H ( 1 ) Methyl H ( 2 )
A g Order No ZFSZFS No ZFSZFS
A iso g e 2 −113.43−69.47 28.7917.63
Δ g iso 4 166.55102.00 −42.28−25.89
Δ g ˜ 4 −0.008.85 0.00−2.25
A SD g e 2 0.002.26 0.006.82
Δ g iso 4 0.00−3.32 −0.00−10.01
Δ g ˜ 4 −2.12−0.85 −6.27−2.54
Ave. Calc. U-H distance, Angstrom7.16 6.15
Orbital Contribution28.96 29.88
Reference shift31.58 31.58
ContactMethod 153.1241.38 −13.48−10.51
Method 253.1232.53 -13.48-8.26
PseudocontactMethod 1−2.12−1.91 −6.27−5.73
Method 2−2.126.94 −6.27−7.98
Total isotropic shielding79.9668.43 10.1313.65
Delta −48.38−36.85 30.4617.93
Experimental δ shift [4]−18.04 4.59
Table 7. Calculated 1 H pNMR shifts for the ring H atoms of U(MeBTP) 3 in ppm at 303K. For numbering see Figure 5. The final, ZFS corrected, computational values are highlighted.
Table 7. Calculated 1 H pNMR shifts for the ring H atoms of U(MeBTP) 3 in ppm at 303K. For numbering see Figure 5. The final, ZFS corrected, computational values are highlighted.
Component Ortho-Pyridine H ( 3 ) Para-Pyridine H ( 4 )
A g Order No ZFSZFS No ZFSZFS
A iso g e 2 −49.74−30.46 155.8595.45
Δ g iso 4 73.0344.72 −228.84−140.15
Δ g ˜ 4 −0.003.88 0.00−12.16
A SD g e 2 0.00−10.72 0.00−15.21
Δ g iso 4 0.0015.73 0.0022.34
Δ g ˜ 4 10.014.03 14.485.83
Calculated U–H distance, Angstrom5.52 6.38
Orbital Contribution21.75 22.56
Reference shift31.58 31.58
ContactMethod 123.2918.14 −72.99−56.86
Method 223.2914.26 −72.99−44.70
PseudocontactMethod 110.019.04 14.4812.96
Method 210.0112.92 14.480.80
Total isotropic shielding55.0548.94 −35.95−21.35
Delta −23.47−17.36 67.5352.92
Experimental δ shift [4]5.64 33.95
Table 8. Calculated 1 H pNMR shifts for U(C 5 Me 5 ) 3 in ppm at 298K. The final, ZFS corrected, computational values are highlighted.
Table 8. Calculated 1 H pNMR shifts for U(C 5 Me 5 ) 3 in ppm at 298K. The final, ZFS corrected, computational values are highlighted.
AgOrder No ZFSZFS
A iso g e 2 −6.02−5.13
Δ g iso 4 7.836.68
Δ g ˜ 4 −0.001.12
A SD g e 2 −0.00−2.36
Δ g iso 4 0.003.07
Δ g ˜ 4 2.831.61
Orbital Contribution29.95
ContactMethod 11.822.67
Method 21.821.55
PseudocontactMethod 12.832.32
Method 22.833.44
Total isotropic shielding34.6034.94
Reference shift, 31.58 ppm (TMS)
Delta −3.02−3.36
Experimental δ shift [77]−0.90
Table 9. Calculated 13 C pNMR shifts for U(C 5 Me 5 ) 3 in ppm with no spin-orbit corrections to the hyperfine term at 298 K). The final, ZFS corrected, computational values are highlighted.
Table 9. Calculated 13 C pNMR shifts for U(C 5 Me 5 ) 3 in ppm with no spin-orbit corrections to the hyperfine term at 298 K). The final, ZFS corrected, computational values are highlighted.
Component CyclopentadienylMethyl
A g Order No ZFSZFS No ZFSZFS
A iso g e 2 256.34218.66 −343.69−293.18
Δ g iso 4 −333.83−284.76 447.59381.80
Δ g ˜ 4 −0.00−47.79 0.0064.07
A SD g e 2 −0.01−64.23 −0.01−11.63
Δ g iso 4 0.0183.65 0.0115.14
Δ g ˜ 4 77.5043.90 13.767.88
Orbital Contribution57.01 173.62
Reference shift189.03 189.03
ContactMethod 1−77.49−113.89 103.89152.69
Method 2−77.49−66.10 103.8988.62
PseudocontactMethod 177.5163.31 13.7611.39
Method 277.5115.53 13.7675.46
Total isotropic shielding57.026.44 291.27337.70
Delta 132.00182.59 -102.25−148.68
Experimental δ shift [77]324.80 −86.70
Table 10. Calculated 13 C pNMR shifts for U(C 5 Me 5 ) 3 in ppm with spin-orbit corrections to the hyperfine term at 298 K). The final, ZFS corrected, computational values are highlighted.
Table 10. Calculated 13 C pNMR shifts for U(C 5 Me 5 ) 3 in ppm with spin-orbit corrections to the hyperfine term at 298 K). The final, ZFS corrected, computational values are highlighted.
Component CyclopentadienylMethyl
A g Order No ZFSZFS No ZFSZFS
Isotropic Hyperfine Contribution
A iso g e 2 256.34218.66 −343.69−293.18
Δ g iso 4 −333.83−284.76 447.59381.80
Δ g ˜ 4 −0.00−47.79 0.0064.07
A iso Δ SO g e 4 74.4463.50 −5.38−4.59
Δ g iso 6 −96.94−82.69 7.005.97
Δ g ˜ 6 −0.00−13.88 0.001.00
Anisotropic Hyperfine Contribution
A SD g e 2 −0.01−64.23 −0.01−11.63
Δ g iso 4 0.0183.65 0.0115.14
Δ g ˜ 4 77.5043.90 13.767.88
A Δ SO , SYM g e 4 0.00127.76 0.0010.78
Δ g iso 6 −0.00−166.38 −0.00−14.04
Δ g ˜ 6 −152.63−86.94 −12.76−7.31
A Δ SO , AS g e 4 −0.000.00 −0.00−0.00
Δ g iso 6 −0.00−0.00 −0.000.00
Δ g ˜ 6 −0.00−0.01 0.000.00
Orbital Contribution57.01 173.62
Reference shift189.03 189.03
Method 1
Contact * A no SOC −77.49−113.89 103.89152.69
A with SOC −99.99−146.96 105.52155.08
Pseudocontact A no SOC 77.5163.31 13.7611.39
A with SOC −75.13−62.25 1.000.82
Isotropic shielding A no SOC 57.026.44 291.27337.70
A with SOC −118.11−133.18 280.14329.52
Delta A no SOC 132.00182.59 −102.25−148.68
A with SOC 307.14341.23 −91.11−140.50
A g Order No ZFSZFS No ZFSZFS
Method 2
Contact A no SOC −77.49−66.10 103.8988.62
A with SOC −3.05−2.61 137.4684.04
Pseudocontact A no SOC 77.5115.53 13.7675.46
A with SOC 77.51143.29 13.7686.24
Isotropic shielding A no SOC 57.026.44 291.27337.70
A with SOC 131.47197.69 324.84343.9
Delta A no SOC 132.00182.59 −102.25−148.68
A with SOC 57.57−8.66 −135.81−154.87
Experimental δ shift [77]324.80 −86.70
* Rows marked A—no SOC refer to the results before the hyperfine matrix has been corrected with spin-orbit effects, A—with SOC are results including this correction.
Table 11. Calculated 29 Si pNMR shifts of U(N(SiMe t Bu 2 ) 2 ) 3 in ppm at 298 K. The final, ZFS corrected, computational values are highlighted.
Table 11. Calculated 29 Si pNMR shifts of U(N(SiMe t Bu 2 ) 2 ) 3 in ppm at 298 K. The final, ZFS corrected, computational values are highlighted.
AgOrder No ZFSZFS
A iso g e 2 −624.79−608.68
Δ g iso 4 889.64866.71
Δ g ˜ 4 −0.0048.19
A SD g e 2 −0.00−6.14
Δ g iso 4 −0.008.74
Δ g ˜ 4 37.2639.49
ContactMethod 1264.85306.22
Method 2264.85258.03
PseudocontactMethod 137.2642.09
Method 237.2690.28
Orbital Contribution344.04
Reference shift331.60
Total isotropic shielding646.15692.34
Delta −314.54 360.74
Experimental δ shift [2]−296.00
Table 12. Calculated 29 Si pNMR shifts of U(C 5 Me 5 · SiMe 3 ) 3 in ppm at 298 K. The final, ZFS corrected, computational values are highlighted.
Table 12. Calculated 29 Si pNMR shifts of U(C 5 Me 5 · SiMe 3 ) 3 in ppm at 298 K. The final, ZFS corrected, computational values are highlighted.
AgOrder No ZFSZFS
A iso g e 2 −389.29−376.21
Δ g iso 4 532.93515.02
Δ g ˜ 4 0.004.55
A SD g e 2 0.000.99
Δ g iso 4 0.00−1.36
Δ g ˜ 4 −19.09−20.19
Ave. Calc. U–H distance, Angstrom4.16
ContactMethod 1143.64143.36
Method 2143.64138.81
PseudocontactMethod 1−19.09−20.56
Method 2−19.09−16.01
Orbital Contribution358.12
Reference shift346.58
Total isotropic shielding482.66480.93
Delta −136.08 134.35
Experimental δ shift [79,80]−155.00
Table 13. Calculated 29 Si pNMR shifts of U(C 5 Me 5 · Et)(1,4-C 8 H 6 (Si i Pr) in ppm at 298 K. The final, ZFS corrected, computational values are highlighted.
Table 13. Calculated 29 Si pNMR shifts of U(C 5 Me 5 · Et)(1,4-C 8 H 6 (Si i Pr) in ppm at 298 K. The final, ZFS corrected, computational values are highlighted.
AgOrder No ZFSZFS
A iso g e 2 −324.51−298.57
Δ g iso 4 512.84471.85
Δ g ˜ 4 0.0030.66
A SD g e 2 0.00−6.60
Δ g iso 4 −0.0010.43
Δ g ˜ 4 5.666.25
ContactMethod 1188.33203.95
Method 2188.33173.28
PseudocontactMethod 15.6610.08
Method 25.6640.74
Orbital Contribution326.64
Reference shift346.34
Total isotropic shielding520.64540.66
Delta −174.29 194.32
Experimental δ shift [1]−136.70

Share and Cite

MDPI and ACS Style

Lloyd, A.W.; Moylan, H.M.; McDouall, J.J.W. Modelling the Effect of Zero-Field Splitting on the 1H, 13C and 29Si Chemical Shifts of Lanthanide and Actinide Compounds. Magnetochemistry 2019, 5, 3. https://doi.org/10.3390/magnetochemistry5010003

AMA Style

Lloyd AW, Moylan HM, McDouall JJW. Modelling the Effect of Zero-Field Splitting on the 1H, 13C and 29Si Chemical Shifts of Lanthanide and Actinide Compounds. Magnetochemistry. 2019; 5(1):3. https://doi.org/10.3390/magnetochemistry5010003

Chicago/Turabian Style

Lloyd, Austin W., Helen M. Moylan, and Joseph J. W. McDouall. 2019. "Modelling the Effect of Zero-Field Splitting on the 1H, 13C and 29Si Chemical Shifts of Lanthanide and Actinide Compounds" Magnetochemistry 5, no. 1: 3. https://doi.org/10.3390/magnetochemistry5010003

APA Style

Lloyd, A. W., Moylan, H. M., & McDouall, J. J. W. (2019). Modelling the Effect of Zero-Field Splitting on the 1H, 13C and 29Si Chemical Shifts of Lanthanide and Actinide Compounds. Magnetochemistry, 5(1), 3. https://doi.org/10.3390/magnetochemistry5010003

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop