Parity nonconservation contribution to the nuclear magnetic resonance shielding constants of chiral molecules: a four-component relativistic study

A systematic four-component relativistic study of the parity nonconservation (PNC) contribution to the (isotropic) NMR shielding constants of chiral molecules is presented for the P enantiomers of the series H(2)X(2) (X=(17)O,(33)S,(77)Se,(125)Te,(209)Po). The PNC contributions are obtained within a linear response approach at the Hartree-Fock level. A careful design of the basis sets is necessary. The four-component relativistic results based on the Dirac-Coulomb Hamiltonian are compared with the nonrelativistic Levy-Leblond results and those obtained by the spin-free modified Dirac Hamiltonian. The calculations confirm the nonrelativistic scaling law Z(2.4) of the PNC contribution with respect to nuclear charge Z. However, the calculations also show that the overall scaling is significantly modified by relativistic effects. The scalar relativistic effect scales as Z(4.7) for the selected set of molecules, whereas the spin-orbit effect, of opposite sign, scales better than Z(6) and completely dominates the PNC contribution for the heaviest elements. This opens up the intriguing possibility of the experimental observation of PNC effects on NMR parameters of molecules containing heavy atoms. The presented formalism is expected to be valuable in assisting the search for suitable candidate molecules.     

Use of locally dense basis sets for nuclear magnetic resonance shielding calculations

The locally dense basis set approach to the calculation of nuclear magnetic resonance shieldings is one in which a sufficiently large or dense set of basis functions is used for an atom or molecular fragment containing the resonant nucleus or nuclei of interest and fewer or attenuated sets of basis functions employed elsewhere. Provided the dense set is of sufficient size, this approach is capable of determining chemical shieldings nearly as well as a calculation with a balanced basis set of quality equal to the locally dense set, but with considerable savings of CPU time. Detailed comparisons are provided of locally dense and balanced calculations in the gauge including atomic orbital (GIAO) method for the individual principal values, the isotropic shieldings, and the tensor orientations for hydrogen, carbon, nitrogen, oxygen, fluorine, and phosphorus nuclei. It is seen that chemical functional groups can often define the appropriate molecular fragment to be taken locally dense. While the present test cases are for the most part small molecules, the value of the method is that it will allow calculations on systems that would otherwise presently be computationally expensive or inaccessible. © John Wiley & Sons, Inc.     

Leading-order relativistic effects on nuclear magnetic resonance shielding tensors

We present perturbational ab initio calculations of the nuclear-spin-dependent relativistic corrections to the nuclear magnetic resonance shielding tensors that constitute, together with the other relativistic terms reported by us earlier, the full leading-order perturbational set of results for the one-electron relativistic contributions to this observable, based on the (Breit-)Pauli Hamiltonian. These contributions are considered for the H(2)X (X = O,S,Se,Te,Po) and HX (X = F,Cl,Br,I,At) molecules, as well as the noble gas (Ne, Ar, Kr, Xe, Rn) atoms. The corrections are evaluated using the relativistic and magnetic operators as perturbations on an equal footing, calculated using analytical linear and quadratic response theory applied on top of a nonrelativistic reference state provided by self-consistent field calculations. The (1)H and heavy-atom nuclear magnetic shielding tensors are compared with four component, nearly basis-set-limit Dirac-Hartree-Fock calculations that include positronic excitations, as well as available literature data. Besides the easy interpretability of the different contributions in terms of familiar nonrelativistic concepts, the accuracy of the present perturbational scheme is striking for the isotropic part of the shielding tensor, for systems including elements up to Xe.     

Four-component relativistic theory for nuclear magnetic shielding constants: the orbital decomposition approach

The authors present a scheme to simplify four-component relativistic calculations of nuclear magnetic shielding constants. The central idea is to decompose each first order orbital into two terms, one is magnetically balanced and directly leads to the diamagnetic term, and the other is, to leading order of relativity, kinetically balanced and can therefore simply be represented in the basis of unperturbed positive energy states. As a matrix formulation, the present approach is far simpler than other operator theories. Combined with the Dirac-Kohn-Sham ansatz, the nuclear magnetic shielding constants for the Kr, Xe, and Rn atoms as well as the HBr and HI molecules are calculated, and the results compare favorably with those of other schemes.     

Four-component relativistic theory for nuclear magnetic shielding constants: critical assessments of different approaches

Both formal and numerical analyses have been carried out on various exact and approximate variants of the four-component relativistic theory for nuclear magnetic shielding constants. These include the standard linear response theory (LRT), the full or external field-dependent unitary transformations of the Dirac operator, as well as the orbital decomposition approach. In contrast with LRT, the latter schemes take explicitly into account both the kinetic and magnetic balances between the large and small components of the Dirac spinors, and are therefore much less demanding on the basis sets. In addition, the diamagnetic contributions, which are otherwise "missing" in LRT, appear naturally in the latter schemes. Nevertheless, the definitions of paramagnetic and diamagnetic terms are not the same in the different schemes, but the difference is only of O(c(-2)) and thus vanishes in the nonrelativistic limit. It is shown that, as an operator theory, the full field-dependent unitary transformation approach cannot be applied to singular magnetic fields such as that due to the magnetic point dipole moment of a nucleus. However, the inherent singularities can be avoided by the corresponding matrix formulation (with a partial closed summation). All the schemes are combined with the Dirac-Kohn-Sham ansatz for ground state calculations, and by using virtually complete basis sets a new and more accurate set of absolute nuclear magnetic resonance shielding scales for the rare gases He-Rn have been established.     

Density-functional calculations of relativistic spin-orbit effects on nuclear magnetic shielding in paramagnetic molecules

Terms arising from the relativistic spin-orbit effect on both hyperfine and Zeeman interactions are introduced to density-functional theory calculation of nuclear magnetic shielding in paramagnetic molecules. The theory is a generalization of the former nonrelativistic formulation for doublet systems and is consistent to O(alpha4), the fourth power of the fine structure constant, for the spin-orbit terms. The new temperature-dependent terms arise from the deviation of the electronic g tensor from the free-electron g value as well as spin-orbit corrections to hyperfine coupling tensor A, the latter introduced in the present work. In particular, the new contributions include a redefined isotropic pseudocontact contribution that consists of effects due to both the g tensor and spin-orbit corrections to hyperfine coupling. The implementation of the spin-orbit terms makes use of all-electron atomic mean-field operators and/or spin-orbit pseudopotentials. Sample results are given for group-9 metallocenes and a nitroxide radical. The new O(alpha4) corrections are found significant for the metallocene systems while they obtain small values for the nitroxide radical. For the isotropic shifts, none of the three beyond-leading-order hyperfine contributions are negligible.     

Effect of perturbation of the atomic basis on nuclear magnetic shielding constants for small molecules

We have calculated the components of the paramagnetic part of the magnetic shielding tensor for nuclei in molecules of LiH, HF, and H₂O within the uncoupled variant of Hartree-Fock-Roothaan perturbation theory, taking into account the dependence of the original basis set of Slater-type AO's (STO's) on the perturbation parameter. We have shown that it is necessary to take into account such a dependence when calculating the components of the magnetic shielding tensor in minimal basis sets of STO's. We have carried out a comparative analysis of the data obtained with results of other approaches.Translated from Teoreticheskaya i Éksperimental'naya Khimiya, Vol. 24, No. 5, pp. 527–532, September–October, 1988     

Electric field effects on the shielding constants of noble gases: a four-component relativistic Hartree-Fock study

Second derivatives of nuclear shielding constants with respect to an electric field, i.e., shielding polarizabilities, have been calculated for the noble gas atoms from helium to xenon. The calculations have been carried out using the four-component relativistic Hartree-Fock method. In order to assess the importance of the individual relativistic corrections, the shielding polarizabilities have also been calculated at the nonrelativistic Hartree-Fock level, with spin-orbit and scalar (Darwin and mass-velocity) effects having been established by perturbative methods. Electron correlation effects have been estimated using the second-order polarization propagator approach. The relativistic effects on the tensor components of the shielding polarizabilities are found to be larger and changing less regularly with the atomic number than for the shielding constant itself. However, there is a partial cancellation of the contributions to the parallel and perpendicular components of the shielding polarizability and as a consequence the mean shielding polarizability is far less affected than the individual components.     

A comparison of two-component and four-component approaches for calculations of spin-spin coupling constants and NMR shielding constants of transition metal cyanides

Relativistic density functional theory (DFT) calculations of nuclear spin-spin coupling constants and shielding constants have been performed for selected transition metal (11th and 12th group of periodic table) and thallium cyanides. The calculations have been carried out using zeroth-order regular approximation (ZORA) Hamiltonian and four-component Dirac-Kohn-Sham (DKS) theory with different nonrelativistic exchange-correlation functionals. Two recent approaches for representing the magnetic balance (MB) between the large and small components of four-component spinors, namely, mDKS-RMB and sMB, have been employed for shielding tensor calculations and their results have been compared. Relativistic effects have also been analysed in terms of scalar and spin-orbit contributions at the two-component level of theory, including discussion of heavy-atom-on-light-atom effects for (1)J(CN), σ(C), and σ(N). The results for molecules containing metals from 4th row of periodic table show that relativistic effects for them are small (especially for spin-spin coupling constants). The biggest effects are observed for the 6th row where nonrelativistic theory reproduces only about 50%-70% of the two-component ZORA results for (1)J(MeC) and about 75% for heavy metal shielding constants. It is important to employ a full Dirac picture for calculations of heavy metal shielding constants, since ZORA reproduces only 75%-90% of the DKS results. Smaller discrepancies between ZORA-DFT and DKS are observed for nuclear spin-spin coupling constants. No significant differences are observed between the results obtained using mDKS-RMB and sMB approaches for magnetic balance in four-component calculations of the shielding constants.     

Self-consistent perturbation-INDO calculations of nuclear magnetic resonance carbon–carbon coupling constants for methylcycloalkanes

One bond ¹³C,¹³C coupling constants have been calculated for some methylcycloalkanes, as well as for 2-methylbutane, using the self-consistent perturbation theory as formulated by Blizzard and Santry at the INDO (intermediate neglect of differential overlap) level of approximation. Together with previously published experimental data, the results can be interpreted satisfactorily on the basis of the Walsh model for cyclopropane and the known s character relationship ¹J(¹³C¹³C) = 575·s(i)· s(j).     

Nuclear magnetic resonance shielding constants and chemical shifts in linear 199Hg compounds: a comparison of three relativistic computational methods

We investigate the importance of relativistic effects on NMR shielding constants and chemical shifts of linear HgL(2) (L = Cl, Br, I, CH(3)) compounds using three different relativistic methods: the fully relativistic four-component approach and the two-component approximations, linear response elimination of small component (LR-ESC) and zeroth-order regular approximation (ZORA). LR-ESC reproduces successfully the four-component results for the C shielding constant in Hg(CH(3))(2) within 6 ppm, but fails to reproduce the Hg shielding constants and chemical shifts. The latter is mainly due to an underestimation of the change in spin-orbit contribution. Even though ZORA underestimates the absolute Hg NMR shielding constants by ～2100 ppm, the differences between Hg chemical shift values obtained using ZORA and the four-component approach without spin-density contribution to the exchange-correlation (XC) kernel are less than 60 ppm for all compounds using three different functionals, BP86, B3LYP, and PBE0. However, larger deviations (up to 366 ppm) occur for Hg chemical shifts in HgBr(2) and HgI(2) when ZORA results are compared with four-component calculations with non-collinear spin-density contribution to the XC kernel. For the ZORA calculations it is necessary to use large basis sets (QZ4P) and the TZ2P basis set may give errors of ～500 ppm for the Hg chemical shifts, despite deceivingly good agreement with experimental data. A Gaussian nucleus model for the Coulomb potential reduces the Hg shielding constants by ～100-500 ppm and the Hg chemical shifts by 1-143 ppm compared to the point nucleus model depending on the atomic number Z of the coordinating atom and the level of theory. The effect on the shielding constants of the lighter nuclei (C, Cl, Br, I) is, however, negligible.     

Quantitative prediction of gas-phase 19F nuclear magnetic shielding constants

Benchmark calculations of (19)F nuclear magnetic shielding constants are presented for a set of 28 molecules. Near-quantitative accuracy (ca. 2 ppm deviation from experiment) is achieved if (1) electron correlation is adequately treated by employing the coupled-cluster singles and doubles (CCSD) model augmented by a perturbative correction for triple excitations [CCSD(T)], (2) large (uncontracted) basis sets are used, (3) gauge-including atomic orbitals are used to ensure gauge-origin independence, (4) calculations are performed at accurate equilibrium geometries [obtained from CCSD(T)/cc-pVTZ calculations correlating all electrons], and (5) vibrational averaging and temperature corrections via second-order vibrational perturbation theory (VPT2) are included. For the CCSD(T)/13s9p4d3f calculations corrected for vibrational effects, mean and standard deviation from experiment are -1.9 and 1.6 ppm, respectively. Less elaborate theoretical treatments result in larger errors. Consideration of relative shifts can reduce the mean deviation (through an appropriately chosen reference compound), but does not change the standard deviation. Density-functional theory calculations of absolute and relative (19)F nuclear magnetic shielding constants are found to be, at best, as accurate as the corresponding Hartree-Fock self-consistent-field calculations and are not improved by consideration of vibrational effects. Molecular systems containing fluorine-oxygen, fluorine-nitrogen, and fluorine-fluorine bonds are found to be more challenging than the other investigated molecules for the considered theoretical methods.     

A fully relativistic method for calculation of nuclear magnetic shielding tensors with a restricted magnetically balanced basis in the framework of the matrix Dirac-Kohn-Sham equation

A new relativistic four-component density functional approach for calculations of NMR shielding tensors has been developed and implemented. It is founded on the matrix formulation of the Dirac-Kohn-Sham (DKS) method. Initially, unperturbed equations are solved with the use of a restricted kinetically balanced basis set for the small component. The second-order coupled perturbed DKS method is then based on the use of restricted magnetically balanced basis sets for the small component. Benchmark relativistic calculations have been carried out for the (1)H and heavy-atom nuclear shielding tensors of the HX series (X=F,Cl,Br,I), where spin-orbit effects are known to be very pronounced. The restricted magnetically balanced basis set allows us to avoid additional approximations and/or strong basis set dependence which arises in some related approaches. The method provides an attractive alternative to existing approximate two-component methods with transformed Hamiltonians for relativistic calculations of chemical shifts and spin-spin coupling constants of heavy-atom systems. In particular, no picture-change effects arise in property calculations.     

Improved description of nuclear magnetic resonance chemical shielding constants using the M06-L meta-generalized-gradient-approximation density functional

The performance of the M06-L density functional has been tested for four databases of NMR isotropic chemical shielding constants. Comparison with the B3LYP, BLYP, HCTH, KT1, KT2, LSDA, OPBE, OLYP, PBE, TPSS, and VSXC functionals shows that M06-L has improved performance for calculating NMR chemical shielding constants, especially for highly correlated systems. We also found that VSXC and M06-L have encouraging accuracy for calculating (13)C chemical shielding constants, and both functionals perform very well for the chemical shielding constants in the o-benzyne molecule.     

Calculation of nuclear magnetic shielding constants in molecules by the extended Hückel method with expansion by Gaussian functions

Theoretical and Experimental Chemistry -     

A simple method for correcting baseline artifacts in nuclear magnetic resonance spectroscopy

A simple baseline correction program is described. It depends on an adjustable threshold for initiation of a data-sorting process between the signals and baseline of a computer-stored n.m.r. spectrum. An application is discussed.     

Calculation of nuclear magnetic shieldings using an analytically differentiated relativistic shielding formula

Two expressions for nuclear-magnetic-shielding tensor components based on analytically differentiating the electronic energy of a system are presented. The first is based on a second-order Douglas-Kroll-Hess approach, in which the off-diagonal block terms of the transformed Dirac Hamiltonian are diminished to second order with respect to both the electrostatic nuclear attraction potential V and the magnetic vector potential A. The second expression is based on the method of Barysz-Sadlej-Snijders, in which the off-diagonal block terms in the transformed Dirac Hamiltonian are completely eliminated with respect to purely V terms, while they are diminished to second order with respect to terms including A. The two approaches are applied to the calculation of nuclear magnetic shieldings of HX (X=F, Cl, Br, I), H2X (X=O, S, Se, Te), and noble gas X (X =He,Ne,Ar,Kr,Xe) systems with common gauge origins. The results show that relativistic corrections of higher than second order are negligibly small, except for the paramagnetic parts of I, Te, and Xe shieldings. The present calculations yield very large positive values for the anisotropy of proton shielding, deltasigma(H) = sigmaparallel(H)-sigmaperpendicular(H), of HI compared to previous reports. Unfortunately, no experimental values for the anisotropy of proton shielding in HI are available for verification.     

Iglo calculations of the antisymmetric components of nuclear magnetic shielding tensors

Calculations of the antisymmetric components of ¹³C shielding tensors have been performed in several compounds using the IGLO (individual gauge for localized orbitals) method. The results indicate chemical situations in which the detection of the antisymmetric components should be easier than in the cases analyzed by Iwai and Saika.     

Quasirelativistic study of 125Te nuclear magnetic shielding constants and chemical shifts

Calculations for ¹²⁵Te magnetic shielding constants and chemical shifts were carried out using a quasirelativistic Hamiltonian including the spin‐free relativistic, one‐ and two‐electron spin–orbit, and relativistic magnetic interaction terms. For the tellurium‐containing series Te(CH₃)₂, TeH₂, TeF₆, Te(CH₃)₄, and Te(CH₃)₂Cl₂, the relativistic effects amounted to as much as 1300 ppm and were very important for qualitatively reproducing the absolute value of the ¹²⁵Te shielding constants obtained experimentally. On the other hand, for the ¹²⁵Te chemical shifts the relativistic effects were less important, because they cancelled each other between the sample and reference compounds. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1502–1508, 2001     

Highly accurate relativistic universal Gaussian basis set for Dirac–Fock–Breit calculations

Ion mobilities of H₂O⁺ drifting in helium are calculated and compared with experiment. These calculations employ global potential energy surfaces of the H₂O⁺–He complex, which in the present case were calculated ab initio at the unrestricted MP2 level of theory using a basis set of aug‐cc‐pVTZ quality, and treating the ion as a rigid body. Details are presented of the general characteristics of both the ground and first‐excited electronic states of the complex. Although only the ground‐state surface was used for the mobility calculations, the ab initio determination of the ground state necessitated the inclusion of the first‐excited state owing to the presence of a crossing between the two. This crossing is also described. Mobilities calculated from the global surfaces are in good agreement with experiment. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005