The Absolute Shielding Constants of Heavy Nuclei: Resolving the Enigma of the (119)Sn Absolute Shielding

Relativistic DFT Calculations of Cadmium and Selenium Solid-State NMR Spectra of CdSe Nanocrystal Surfaces

Solid-state NMR spectra have been used to probe the structure of CdSe nanocrystals and propose detailed models of their surface structures. Density functional theory (DFT)-optimized cluster models that represent probable molecular structures of carboxylate-coordinated surface sites have been proposed. However, to the best of our knowledge, 113Cd and 77Se chemical shifts have not been calculated for these surface models. We performed relativistic DFT calculations of cadmium and selenium magnetic shielding tensors on model compounds with previously measured solid-state NMR spectra with (i) the four-component Dirac-Kohn-Sham (DKS) Hamiltonian and (ii) the scalar and (iii) spin-orbit levels within the ZORA Hamiltonian. Molecular clusters with Cd and Se sites in varying bonding environments were used to model CdSe (100) and CdSe(111) surfaces capped with carboxylic acid ligands. Our calculations identify the observed 113Cd isotropic chemical shifts δ(iso) of -465, -318, and -146 ppm arising from CdSeO3, CdSe2O2, and CdSe3O surface groups, respectively, with very good agreement with experimental measurements. The 113Cd chemical shifts linearly decrease with the number of O-neighbors. The calculated spans (δ11 - δ33) encompass the experimental values for CdSe3O and CdSe2O2 clusters but are slightly larger than the measured value for CdSeO3 clusters. Relativistic DFT calculations predicted a one-bond 113Cd-77Se scalar coupling of 258 Hz, which is in good agreement with the experimental values of 250 Hz. With a dense coverage of carboxylic acid ligands, the CdSe (100) surface shows a distribution of Cd-Se bond lengths and J-couplings. Relativistic DFT simulations thus aid in interpretation of NMR spectra of CdSe nanocrystals and related nanomaterials.

Relativistic study of parity-violating nuclear spin-rotation tensors

We present a four-component relativistic approach to describe the effects of the nuclear spin-dependent parity-violating (PV) weak nuclear forces on nuclear spin-rotation (NSR) tensors. The formalism is derived within the four-component polarization propagator theory based on the Dirac-Coulomb Hamiltonian. Such calculations are important for planning and interpretation of possible future experiments aimed at stringent tests of the standard model through the observation of PV effects in NSR spectroscopy. An exploratory application of this theory to the chiral molecules H₂X₂ (X = ¹⁷O, ³³S, ⁷⁷Se, ¹²⁵Te, and ²⁰⁹Po) illustrates the dramatic effect of relativity on these contributions. In particular, spin-free and spin-orbit effects are even of opposite signs for some dihedral angles, and the latter fully dominate for the heavier nuclei. Relativistic four-component calculations of isotropic nuclear spin-rotation constants, including parity-violating electroweak interactions, give frequency differences of up to 4.2 mHz between the H₂Po₂ enantiomers; on the nonrelativistic level of theory, this energy difference is 0.1 mHz only.

Performance of the LRESC Model on top of DFT Functionals for Relativistic NMR Shielding Calculations

The linear response within the elimination of the small component model (LRESC) is an insightful and computationally efficient method for including relativistic effects on molecular properties like the nuclear magnetic shielding constants, spin-rotation constant, g-tensor, and electric field gradient of heavy atom containing molecules with atoms belonging up to the sixth row of the periodic table. One of its main advantages is its capacity to analyze the electronic origin of the different relativistic correcting terms. Until now, it was always applied on top of Hartree-Fock ground-state wave functions (LRESC/HF) to calculate and analyze NMR shieldings. In this work, we show the performance of the LRESC formalism on top of some density functional theory (DFT) functionals to compute tin shielding constants in SnX₄ (X = H, F, Cl, Br, I) molecular systems. We analyze the performance of each LRESC/DFT scheme on reproducing the electronic mechanisms of the shieldings, taking as a benchmark the results of relativistic calculations at the RPA level of approach (4c/RPA). As in previous works, we divide the LRESC relativistic correcting terms into two groups: core-dependent and ligand-dependent contributions. It is shown here that core-dependent corrections are well-reproduced for the selected DFT functionals, but some differences arise in the ligand-dependent ones. We focus on the performance of different functionals, including the same electron correlation part but containing different amounts of HF exchange. The best results are obtained for the BHandHLYP functional (50% of HF exchange) and the worst for BLYP (0%). When the percentage of HF exchange increases, ligand-dependent contributions are better described, and the final LRESC/DFT results are closer to those obtained with LRESC/HF and 4c/RPA methods. The spin-orbit correction to the shielding constant is one of the main ligand-dependent contributions (there are two more) with total value depending on the amount of HF exchange included in the functional. When the amount of HF exchange decreases, the spin-orbit contribution becomes larger, overestimating the shielding constant even when nonrelativisitc DFT values are much smaller than the nonrelativistic HF ones, as it happens for the heaviest molecular system studied here (SnI₄).

NMR absolute shielding scales and nuclear magnetic dipole moments of transition metal nuclei

This work reports new, accurate nuclear magnetic dipole moments for transition metal nuclei where the long-standing systematic error due to obsolete diamagnetic correction has been eliminated by ab initio calculations of NMR shielding constants.

Is hydrogen electronegativity higher than Pauling’s value? New clues from the13C and29Si NMR chemical shifts of [CHF3] and [SiHF3] molecules

We previously demonstrated that the δ NMR chemical shift of central NMR active atoms (A), in simple halido [AXn] (A=C, Si, Ge, Sn, Pb, Pt; Xn = combination of n halides, n = 4 or 6) derivatives, could be directly related to X radii overall sum, Σ(rL). Further correlation have also been observed for tetrahedral [AX4] (A=C, Si; X4 = combination of four halides) compounds where the X Pauling electronegativities sum,Σ(χLPau),$\Sigma (\chi _L^{{\rm{Pau}}}),$exceeds a specific value (≈12.4). In this work, we focused on these latter systems considering the H vs. X substitution. The analysis of the literature reported δ(13C) and δ(29Si) NMR chemical shift for the mono hydrogenated derivatives and in particular for [CHF3] and [SiHF3], characterized by the lowest Σ(rL) and the highestΣ(χLPau),$\Sigma (\chi _L^{{\rm{Pau}}}),$suggests a revised value for the H electronegativity ranking with respect to Pauling’s.

Spin-orbit effects on the 125Te magnetic-shielding tensor: A cluster-based ZORA/DFT investigation

Cluster-based calculations of ¹²⁵Te magnetic-shielding tensors demonstrate that inclusion of spin-orbit effects is necessary to obtain the best agreement of theoretical predictions with experiment. The spin-orbit contribution to shielding depends on the oxidation state and stereochemistry of the ¹²⁵Te site. Comparison of the performance of various density functionals indicates that GGA functionals behave similarly to each other in predicting NMR magnetic shielding. The use of hybrid functionals improves the predictive ability on average for a large set of ¹²⁵Te-containing materials. The amount of Hartree-Fock exchange affects the predicted parameters. Inclusion of larger Hartree-Fock exchange contributions in hybrid functionals results in larger slopes of the correlation between calculated magnetic-shielding and experimental chemical-shift principal components, by 10-15% from the ideal value.

Influence of the nuclear charge distribution and electron correlation effects on magnetic shieldings and spin-rotation tensors of linear molecules

The nuclear charge distribution effects (NChDE) on two response properties, the NMR magnetic shielding (σ) and the nuclear spin-rotation (SR) constants (M), are analyzed. We do it employing point-like and Gaussian-like models for describing the nuclear charge density of three linear molecules: HBr, HI and HAt. According to our results, both properties are sensitive to the NChDE. We show that the NChDE are almost completely relativistic, i.e., they nearly vanish in the non-relativistic limit of both properties. We calculated the NChDE on σ and M, and analyzed the differences between them in terms of a relativistic relation between these two properties. Using that relation we found that the electronic core mechanisms are the main ones for the NChDE on the shielding of nuclei of both, molecules and free atoms. The NChDE are smaller on SR constants than on shieldings. Nevertheless, within the relativistic polarization propagator formalism at the RPA level of approach they are very important for SR constants of nuclei in heavy-atom-containing compounds. Astatine in HAt has the largest influence: M_At = −9.95 kHz for a point-like model and −50.10 kHz for a Gaussian-like model. Correlation effects must be included and we do it using different DFT schemes. The PBE0 functional gives results that are closest to experiments for Br and I, though the LDA gives the closest for hydrogen. The value of the SR constant of At is reduced among 350 kHz and 500 kHz from its RPA value, when different and usual functionals are applied. Given that the NChDE on M and σ are mostly relativistic in their origin, these effects are also dependent on electron correlation. They have also a nonvanishing dependence with the Gaunt electron–electron interactions.

The origin and the size of both, electron correlation and nuclear charge distribution on shieldings and spin-rotations of heavy-atom containing linear molecules are shown. The analysis is performed on the relativistic and non relativistic regimes.

73Ge, 119Sn and 207Pb: general cooperative effects of single atom ligands on the NMR signals observed in tetrahedral [MXnY4-n] (M = Ge, Sn, Pb; 1 ≤ n ≤ 4; X, Y = Cl, Br, I) coordination compounds of heavier XIV group elements

An inverse linear relationship between ⁷³Ge, ¹¹⁹Sn and ²⁰⁷Pb NMR chemical shifts and the overall sum of ionic radii of coordinated halido ligands has been discovered in tetrahedral [MX_nY_4-n] (M = Ge, Sn, Pb; 1 ≤ n ≤ 4; X, Y = Cl, Br, I) coordination compounds. This finding is consistent with a previously reported correlation found in octahedral, pentacoordinate and square planar platinum complexes. The effect of the coordinated halido ligands acting on the metal as shielding conducting rings is therefore confirmed also by ⁷³Ge, ¹¹⁹Sn and ²⁰⁷Pb NMR spectroscopy.

NMR effective radius of hydrogen in XIV group hydrides evaluated by NMR spectroscopy

In the [ABr_nI_m] (A = C, Si, Ge, Sn; n + m = 4) compounds, with the heavier halido ligands bonded to the central IV group elements, the ¹³C, ²⁹Si, ⁷³Ge and ¹¹⁹Sn NMR chemical shifts were found to be linearly related to the bonded halides ionic radii overall sum, ∑(r_h). The ²⁰⁷Pb NMR chemical shift of the unstable [PbH₄] hydride could be calculated.

On the calculation of second-order magnetic properties using subsystem approaches in a relativistic framework

The implementation of second-order magnetic properties in a frozen density embedding scheme in a four component relativistic framework is outlined and applied to model H₂X–H₂O systems (X = Se, Te, Po).

Role of Spin-Dependent Terms in the Relationship among Nuclear Spin-Rotation and NMR Magnetic Shielding Tensors

The broadly accepted procedure to obtain the experimental absolute scale of NMR magnetic shieldings, σ, is well-known for nonheavy atom-containing molecules. It was uncovered more than 40 years ago by the works of Ramsey and Flygare. They found a quite accurate relationship among σ and the nuclear spin-rotation constants. Its relativistic extension was very recently proposed, although it has an intrinsic weakness because a new SO-S two-component term needs to be considered. We show how to overcome this problem. We found that (ν_Y^S - ν_Y^atom,S) generalizes the SO-S term, where ν_Y^S = ⟨⟨[((r - r_Y) × α)/(|r - r_Y|³)]; S⁽⁴⁾⟩⟩, r - r_Y is the electron position with respect to the position of nucleus Y, and S⁽⁴⁾ is the four-component total electron spin. When including this new term, one finds that the best of our relativistic Flygare-like models fits quite well with the results of the most accurate method available at the moment. We also show that the difference among the parallel component of σ(Xe) in XeF₂ and σ(Xe) of the free atom is almost completely described by that new term.

Spin-orbit effects on the (119)Sn magnetic-shielding tensor in solids: a ZORA/DFT investigation

Periodic-boundary and cluster calculations of the magnetic-shielding tensors of (119)Sn sites in various co-ordination and stereochemical environments are reported. The results indicate a significant difference between the predicted NMR chemical shifts for tin(ii) sites that exhibit stereochemically-active lone pairs and tin(iv) sites that do not have stereochemically-active lone pairs. The predicted magnetic shieldings determined either with the cluster model treated with the ZORA/Scalar Hamiltonian or with the GIPAW formalism are dependent on the oxidation state and the co-ordination geometry of the tin atom. The inclusion of relativistic effects at the spin-orbit level removes systematic differences in computed magnetic-shielding parameters between tin sites of differing stereochemistries, and brings computed NMR shielding parameters into significant agreement with experimentally-determined chemical-shift principal values. Slight improvement in agreement with experiment is noted in calculations using hybrid exchange-correlation functionals.

Relativistic Calculations of Nuclear Magnetic Resonance Parameters

Relativistic effects are important for the accurate evaluation of the observables of nuclear magnetic resonance (NMR) spectroscopy, the nuclear magnetic shielding and the indirect spin–spin coupling tensors. Some of the most notable relativistic effects, in particular for light elements in the vicinity of heavy nuclei, are due to spin–orbit effects, an effect difficult to evaluate when starting from a non-relativistic wavefunction. Two- and four-component relativistic methods include spin–orbit effects variationally, and the recent improvements in the computational efficiency of these methods open new opportunities for accurate calculations of NMR parameters also for molecules with heavy elements. We here present an overview of the different approximations that have been introduced for calculating relativistic effects with two- and four-component methods and how these methods can be used to calculate the NMR parameters. We will also give some examples of systems that have been studied computationally with two- and four-component relativistic methods and discuss the importance of relativistic effects on the shielding and indirect spin–spin coupling constants.

Theoretical analysis of NMR shieldings in XSe and XTe (X = Si, Ge, Sn and Pb): the spin-rotation constant saga

How the electronic contribution to the spin-rotation constant is close to the paramagnetic contribution of the NMR absolute shielding constant?

Toward an absolute NMR shielding scale using the spin-rotation tensor within a relativistic framework

How can one extend Flygare's rule to the relativistic framework? Three models are proposed here. The best of them shows that σ is related with the spin-rotation tensor, the atomic shielding and a new term coined as σ^SO-S.

Big picture of relativistic molecular quantum mechanics

Any quantum mechanical calculation on electronic structure ought to choose first an appropriate Hamiltonian H and then an Ansatz for parameterizing the wave function Ψ, from which the desired energy/property E(λ) can finally be calculated. Therefore, the very first question is: what is the most accurate many-electron Hamiltonian H? It is shown that such a Hamiltonian i.e. effective quantum electrodynamics (eQED) Hamiltonian, can be obtained naturally by incorporating properly the charge conjugation symmetry when normal ordering the second quantized fermion operators. Taking this eQED Hamiltonian as the basis, various approximate relativistic many-electron Hamiltonians can be obtained based entirely on physical arguments. All these Hamiltonians together form a complete and continuous ‘Hamiltonian ladder’, from which one can pick up the right one according to the target physics and accuracy. As for the many-electron wave function Ψ, the most intriguing questions are as follows. (i) How to do relativistic explicit correlation? (ii) How to handle strong correlation? Both general principles and practical strategies are outlined here to handle these issues. Among the electronic properties E(λ) that sample the electronic wave function nearby the nuclear region, nuclear magnetic resonance (NMR) shielding and nuclear spin-rotation (NSR) coupling constant are especially challenging: they require body-fixed molecular Hamiltonians that treat both the electrons and nuclei as relativistic quantum particles. Nevertheless, they have been formulated rigorously. In particular, a very robust ‘relativistic mapping’ between the two properties has been established, which can translate experimentally measured NSR coupling constants to very accurate absolute NMR shielding scales that otherwise cannot be obtained experimentally. Since the most general and fundamental issues pertinent to all the three components of the quantum mechanical equation HΨ = EΨ (i.e. Hamiltonian H, wave function Ψ, and energy/property E(λ)) have fully been understood, the big picture of relativistic molecular quantum mechanics can now be regarded as established.

Performance of Møller-Plesset second-order perturbation theory and density functional theory in predicting the interaction between stannylenes and aromatic molecules

The performances of Møller-Plesset second-order perturbation theory (MP2) and density functional theory (DFT) have been assessed for the purposes of investigating the interaction between stannylenes and aromatic molecules. The complexes between SnX2 (where X = H, F, Cl, Br, and I) and benzene or pyridine are considered. Structural and energetic properties of such complexes are calculated using six MP2-type and 14 DFT methods. The assessment of the above-mentioned methods is based on the comparison of the structures and interaction energies predicted by these methods with reference computational data. A very detailed analysis of the performances of the MP2-type and DFT methods is carried out for two complexes, namely SnH2-benzene and SnH2-pyridine. Of the MP2-type methods, the reference structure of SnH2-benzene is reproduced best by SOS-MP2, whereas SCS-MP2 is capable of mimicking the reference structure of SnH2-pyridine with the greatest accuracy. The latter method performs best in predicting the interaction energy between SnH2 and benzene or pyridine. Among the DFT methods, ωB97X provides the structures and interaction energies of the SnH2-benzene and SnH2-pyridine complexes with good accuracy. However, this density functional is not as effective in reproducing the reference data for the two complexes as the best performing MP2-type methods. Next, the DFT methods are evaluated using the full test set of SnX2-benzene and SnX2-pyridine complexes. It is found that the range-separated hybrid or dispersion-corrected density functionals should be used for describing the interaction in such complexes with reasonable accuracy.

Structure, solvent, and relativistic effects on the NMR chemical shifts in square-planar transition-metal complexes: assessment of DFT approaches

The role of various factors (geometry, solvent, relativistic treatment, functional) influencing the relativistic NMR chemical shift calculations for square-planar transition-metal complexes.

Spin-rotation and NMR shielding constants in XF molecules (X = B, Al, Ga, In, and Tl)

Accurate spin-rotation and absolute shielding constants in a series of XF molecules (X = (11)B, (27)Al, (69)Ga, (115)In, and (205)Tl) determined using high-level ab initio coupled-cluster and four-component relativistic density-functional theory (DFT) calculations are presented. The accuracy of the results is established by comparing the relativistically and vibrationally corrected calculated values with available experimental data; for spin-rotation and shielding constants for which no experimental data exist, we provide new and reliable values. For both properties, our results can be considered as reference values against which more approximate quantum-chemical methods can be benchmarked.

Relativistic and electron-correlation effects on the nuclear magnetic resonance shieldings of molecules containing tin and lead atoms

The reference values for NMR magnetic shieldings, σ(ref), are of the highest importance when theoretical analysis of chemical shifts are envisaged. The fact that the nonrelativistically valid relationship among spin-rotation constants and magnetic shieldings is not any longer valid for heavy atoms requires that the search for σ(ref) for such atoms needs new strategies to follow. We present here results of σ(ref) that were obtained by applying our own simple procedure which mixes accurate experimental chemical shifts (δ) and theoretical magnetic shieldings (σ). We calculated σ(Sn) and σ(Pb) in a family of heavy-halogen-containing molecules. We found out that σ(ref)[Sn;Sn(CH3)4] in gas phase should be close to 3864.11 ± 20.05 ppm (0.5%). For Pb atom, σ(ref)[Pb;Pb(CH3)4] should be close to 14475.1 ± 500.7 ppm. Such theoretical values correspond to calculations with the relativistic polarization propagator method, RelPPA, at the RPA level of approach. They are closer to experimental values as compared to those obtained applying few different functionals such as PBE0, B3LYP, BLYP, BP86, KT2, and KT3 of the density functional theory, DFT. We studied tin and lead shieldings of the XY(4-n)Z(n) (X = Sn, Pb; Y, Z = H, F, Cl, Br, I) and PbH(4-n)I(n) (n = 0, 1, 2, 3, 4) family of compounds with four-component functionals as implemented in the DIRAC code. For these systems results of calculations with RelPPA-RPA are more reliable than DFT ones. We argue about why those DFT functionals must be modified in order to obtain more accurate results of NMR magnetic shieldings within the relativistic regime: first, there is a dependence among both electron-correlation and relativistic effects that should be introduced in some way in the functionals; and second, the DIRAC code uses standard nonrelativistic functionals and the functionals B3LYP and PBE0 were parametrized only with data taken from light elements. It can explain why they are not able to properly introduce relativistic effects on nuclear magnetic shieldings. We finally show that in the analysis of magnetic shieldings for the family of compounds mentioned above, one must consider the newest and so-called heavy-atom effect on vicinal heavy atoms, HAVHA. Such effects are among the most important relativistic effects in these kind of compounds.

Relativistic calculations of magnetic resonance parameters: background and some recent developments

This article outlines some basic concepts of relativistic quantum chemistry and recent developments of relativistic methods for the calculation of the molecular properties that define the basic parameters of magnetic resonance spectroscopic techniques, i.e. nuclear magnetic resonance shielding, indirect nuclear spin-spin coupling and electric field gradients (nuclear quadrupole coupling), as well as with electron paramagnetic resonance g-factors and electron-nucleus hyperfine coupling. Density functional theory (DFT) has been very successful in molecular property calculations, despite a number of problems related to approximations in the functionals. In particular, for heavy-element systems, the large electron count and the need for a relativistic treatment often render the application of correlated wave function ab initio methods impracticable. Selected applications of DFT in relativistic calculation of magnetic resonance parameters are reviewed.

New Experimental NMR Shielding Scales Mapped Relativistically from NSR: Theory and Application

The recently proposed relativistic mapping between nuclear magnetic resonance (NMR) shielding and nuclear spin-rotation (NSR) coupling tensors [J. Chem. Phys. 2013, 138, 134104] is employed to establish new experimental (more precisely, experimentally derived) absolute shielding constants for H and X in HX (X = F, Cl, Br, and I). The results are much more accurate than the old "experimental" values that were based on the well-known nonrelativistic mapping. The relativistic mapping is very robust in the sense that it is rather insensitive to the quality of one-particle basis sets and the treatment of electron correlation. Relativistic effects in the NSR coupling constants are also elucidated in depth.

33S hyperfine interactions in H2S and SO2 and revision of the sulfur nuclear magnetic shielding scale

Using the Lamb-dip technique, the hyperfine structure in the rotational spectra of H2(33)S and (33)SO2 has been resolved and the corresponding parameters--that is, the sulfur quadrupole-coupling and spin-rotation tensors--were determined. The experimental parameters are in good agreement with results from high-level coupled-cluster calculations, provided that up to quadruple excitations are considered in the cluster operator, sufficiently large basis sets are used, and vibrational corrections are accounted for. The (33)S spin-rotation tensor for H2S has been used to establish a new sulfur nuclear magnetic shielding scale, combining the paramagnetic part of the shielding as obtained from the spin-rotation tensor with a calculated value for the diamagnetic part as well as computed vibrational and temperature corrections. The value of 716(5) ppm obtained in this way for the sulfur shielding of H2S is in good agreement with results from high-accuracy quantum-chemical calculations but leads to a shielding scale that is about 28 ppm lower than the one suggested previously in the literature, based on the (33)S spin-rotation constant of OCS.