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

A second-order perturbation theory treatment of spin-orbit corrections to hyperfine coupling tensors has been implemented within a density-functional framework. The method uses the all-electron atomic mean-field approximation and/or spin-orbit pseudopotentials in incorporating one- and two-electron spin-orbit interaction within a first-principles framework. Validation of the approach on a set of main-group radicals and transition metal complexes indicates good agreement between all-electron and pseudopotential results for hyperfine coupling constants of the lighter nuclei in the system, except for cases in which scalar relativistic effects become important. The nonrelativistic Fermi contact part of the isotropic hyperfine coupling constants is not always accurately reproduced by the exchange-correlation functionals employed, particularly for the triplet and pi-type doublet radicals in the present work. For this reason, ab initio coupled-cluster singles and doubles with perturbative triples results for the first-order contributions have been combined in the validation calculations with the density-functional results for the second-order spin-orbit contributions. In the cases where spin-orbit corrections are of significant magnitude relative to the nonrelativistic first-order terms, they improve the agreement with experiment. Antisymmetric contributions to the hyperfine tensor arise from the spin-orbit contributions and are discussed for the IO2 radical, whereas rovibrational effects have been evaluated for RhC, NBr, and NI.     

Finite perturbation MCSCF and CI calculations of nuclear spin-spin coupling constants for some molecules with multiple bonds

Non-empirical finite perturbation calculations at the Hartree-Fock, multiconfigurational self-consistent field and configuration interaction levels of approximation are presented for the Fermi contact contribution in multiply-bonded molecules ethene, formimide, formaldehyde, ethyne and hydrogen cyanide. The finite perturbation-multiconfigurational SCF (FPMC) method (with few configurations) is free from the UHF triplet instabilities normally present in single configuration coupled Hartree-Fock (CHF) calculations of the Fermi contact (and spin-dipolar) contribution for π-bonded systems. The behaviour of coupling constants, calculated using FPMC for multiply-bonded systems parallels the behaviour of the CHF coupling constants in comparable systems with single bonds only. The effects of dynamic electron correlations are important and are obtained using the Cl method. After accounting for the orbital contribution by means of the single configuration CHF method, agreement with experiment is excellent for systems containing only carbon and hydrogen, when a double-zeta quality basis set is used. For systems containing nitrogen and oxygen agreement is still reasonable, but the use of larger basis sets seems to be necessary if good agreement with experiment is to be obtained.     

Calculations of the indirect nuclear spin–spin coupling constants of PbH4

We report ab initio calculations of the indirect nuclear spin–spin coupling constants of PbH₄ using a basis set which was specially optimized for correlated calculations of spin–spin coupling constants. All nonrelativistic contributions and the most important part of the spin–orbit correction were evaluated at the level of the random phase approximation. Electron correlation corrections to the coupling constants were calculated with the multiconfigurational linear-response method using extended complete and restricted active space wavefunctions as well as with the second-order polarization propagator approximation and the second-order polarization propagator approximation with coupled-cluster singles and doubles amplitudes. The effects of nuclear motion were investigated by calculating the coupling constants as a function of the totally symmetric stretching coordinate. We find that the Fermi contact term dominates the Pb‐H coupling, whereas for the H‐H coupling it is not more important than the orbital paramagnetic and diamagnetic contributions. Correlation affects mainly the Fermi contact term. Its contribution to the one-bond coupling constant is reduced by correlation, independent of the method used; however, the different correlated methods give ambiguous results for the Fermi contact contribution to the H‐H couplings. The dependence of both coupling constants on the Pb‐H bond length is dominated by the change in the Fermi contact term. The geometry dependence is, however, overestimated in the random phase approximation. Received: 16 November 1998 / Accepted: 30 March 1999 / Published online: 14 July 1999     

The calculation of indirect nuclear spin-spin coupling constants in large molecules

We present calculations of indirect nuclear spin-spin coupling constants in large molecular systems, performed using density functional theory. Such calculations, which have become possible because of the use of linear-scaling techniques in the evaluation of the Coulomb and exchange-correlation contributions to the electronic energy, allow us to study indirect spin-spin couplings in molecules of biological interest, without having to construct artificial model systems. In addition to presenting a statistical analysis of the large number of short-range coupling constants in large molecular systems, we analyse the asymptotic dependence of the indirect nuclear spin-spin coupling constants on the internuclear separation. In particular, we demonstrate that, in a sufficiently large one-electron basis set, the indirect spin-spin coupling constants become proportional to the inverse cube of the internuclear separation, even though the diamagnetic and paramagnetic spin-orbit contributions to the spin-spin coupling constants separately decay as the inverse square of this separation. By contrast, the triplet Fermi contact and spin-dipole contributions to the indirect spin-spin coupling constants decay exponentially and as the inverse cube of the internuclear separation, respectively. Thus, whereas short-range indirect spin-spin coupling constants are usually dominated by the Fermi contact contribution, long-range coupling constants are always dominated by the negative diamagnetic spin-orbit contribution and by the positive paramagnetic spin-orbit contribution, with small spin-dipole and negligible Fermi contact contributions.     

Computational studies of electron paramagnetic resonance parameters for paramagnetic molybdenum complexes. 1. Method validation on small and medium-sized systems

A variety of density functional methods have been evaluated in the computation of electronic g-tensors and molybdenum hyperfine couplings for systems ranging from the Mo atom through MoIIIN, [MoVOCl4]-, and [MoVOF5]2- to two larger MoV complexes MoXLCl2 (X=O, S; L=tris(3,5-dimethylpyrazolyl)hydroborate anion). In particular, the influence of the molybdenum basis set and of various exchange-correlation functionals with variable admixtures of Hartree-Fock exchange on the computed EPR parameters have been evaluated in detail. Careful basis-set studies have provided a moderate-sized 12s6p5d all-electron basis on molybdenum that gives hyperfine tensors in excellent agreement with much larger basis sets and that will be useful for calculations on larger systems. The best agreement with experimental data for both hyperfine and g-tensors is obtained with hybrid functionals containing approximately 30-40% Hartree-Fock exchange. Only for MoSLCl2 does increasing spin contamination with increasing exact-exchange admixture restrict the achievable computational accuracy. In all cases, spin-orbit corrections to the hyperfine tensors are sizable and have to be included in accurate calculations. Scalar relativistic effects enhance the isotropic Mo hyperfine coupling by approximately 15-20%. Two-component g-tensor calculations with variational inclusion of spin-orbit coupling show that the Deltag parallel components in [MoVOCl4]- and [MoVOF5]2- depend on higher-order spin-orbit contributions and are thus described insufficiently by the usual second-order perturbation approaches. Computed orientations of g- and hyperfine tensors relative to each other and to the molecular framework for the MoXLCl2 complexes provide good agreement between theory and single-crystal electron paramagnetic resonance experiments. In these cases, the hyperfine tensor orientations are influenced only slightly by spin-orbit effects.     

Unrestricted Hartree-Fock instabilities in semiempirical CNDO/S and INDO/S calculations of spin-spin coupling constants

It is shown that non-convergent calculations of the Fermi contact term of spin-spin coupling constants within the self-consistent and finite perturbation schemes used to solve the coupled Hartree-Fock equations, are originated in non-singlet Hartree-Fock instabilities of the closed-shell restricted Hartree-Fock wavefunction. In CNDO/S and INDO/S wavefunctions, where the electronic system response has been successfully reproduced, all investigated molecules containing MOs were found to be unstable. Results of spin-spin coupling constants are given and compared with experimental as well as FP and SOS INDO values.Part of a Ph.D. thesis (G.E.S.) to be presented to the University of Buenos Aires.Comisión de Investigaciones Científicas (CIC, Pcia. de Bs. As.) fellow.     

Density functional theory study of indirect nuclear spin-spin coupling constants with spin-orbit corrections

This work outlines the calculation of indirect nuclear spin-spin coupling constants with spin-orbit corrections using density functional response theory. The nonrelativistic indirect nuclear spin-spin couplings are evaluated using the linear response method, whereas the relativistic spin-orbit corrections are computed using quadratic response theory. The formalism is applied to the homologous systems H2X (X=O,S,Se,Te) and XH4 (X=C,Si,Ge,Sn,Pb) to calculate the indirect nuclear spin-spin coupling constants between the protons. The results confirm that spin-orbit corrections are important for compounds of the H2X series, for which the electronic structure allows for an efficient coupling between the nuclei mediated by the spin-orbit interaction, whereas in the case of the XH4 series the opposite situation is encountered and the spin-orbit corrections are negligible for all compounds of this series. In addition we analyze the performance of the density functional theory in the calculations of nonrelativistic indirect nuclear spin-spin coupling constants.     

Two-component relativistic density-functional calculations of the dimers of the halogens from bromine through element 117 using effective core potential and all-electron methods

A two-component quasirelativistic Hamiltonian based on spin-dependent effective core potentials is used to calculate ionization energies and electron affinities of the heavy halogen atom bromine through the superheavy element 117 (eka-astatine) as well as spectroscopic constants of the homonuclear dimers of these atoms. We describe a two-component Hartree-Fock and density-functional program that treats spin-orbit coupling self-consistently within the orbital optimization procedure. A comparison with results from high-order Douglas-Kroll calculations--for the superheavy systems also with zeroth-order regular approximation and four-component Dirac results--demonstrates the validity of the pseudopotential approximation. The density-functional (but not the Hartree-Fock) results show very satisfactory agreement with theoretical coupled cluster as well as experimental data where available, such that the theoretical results can serve as an estimate for the hitherto unknown properties of astatine, element 117, and their dimers.     

Improving the calculation of Electron Paramagnetic Resonance hyperfine coupling tensors for d-block metals

Calculation of hyperfine coupling constants (HFCs) of Electron Paramagnetic Resonance from first principles can be a beneficial complement to experimental data in cases where the molecular structure is unknown. We have recently investigated basis set convergence of HFCs in d-block complexes and obtained a set of basis functions for the elements Sc-Zn, which were saturated with respect to both the Fermi contact and spin-dipolar components of the hyperfine coupling tensor [Hedeg?rd et al., J. Chem. Theory Comput., 2011, 7, 4077-4087]. Furthermore, a contraction scheme was proposed leading to very accurate, yet efficient basis sets for the elements Sc-Zn. Here this scheme is tested against a larger test set of molecules and a wider range of DFT functionals. We further investigate the regular aug-cc-pVTZ and core-valence correlation aug-cc-pCVTZ basis sets as well as another core-property basis set, CP(PPP). While aug-cc-pVTZ-J provides hyperfine coupling constants that are almost identical to the converged series (aug-cc-pVTZ-Juc), we observe that not only the regular but also the core-valence correlation basis sets provide results far from the converged results. The usage of specialized core-basis sets leads to a large and highly significant improvement of the calculated hyperfine couplings in comparison with experimental data.     

Comparative density-functional study of the electron paramagnetic resonance parameters of amavadin

Electronic g tensors and hyperfine coupling tensors have been calculated for amavadin, an unusual eight-coordinate vanadium(IV) complex isolated from Amanita muscaria mushrooms. Different density-functional methods have been compared, ranging from local via gradient-corrected to hybrid functionals with a variable Hartree-Fock exchange admixture. For both electron paramagnetic resonance (EPR) properties, hybrid functionals with an appreciable exact-exchange admixture provide the closest agreement with experimental data. Second-order spin-orbit corrections provide non-negligible contributions to the 51V hyperfine tensor. The orientation of g and A tensors relative to each other also depends on spin-orbit corrections to the A tensor. A rationalization for the close resemblance of the EPR parameters of amavadin to those of the structurally rather different vanadyl complexes is provided, based on the nature of the relevant frontier orbitals.     

The permanent dipole moment of gas-phase para-amino benzoic acid revisited

The permanent dipole moment of para-amino benzoic acid has been calculated at various theoretical levels, including Hartree-Fock, second-order M?ller-Plesset perturbation (MP2), coupled cluster singles and doubles (CCSD), and triples corrections CCSD(T), and hybrid density functional theory at B3LYP level. It is found that the B3LYP method fails to provide correct results for the geometry and the permanent dipole moment. These results are significantly improved by MP2 calculations. Our best estimated dipole moments obtained at CCSD and CCSD(T) levels are in good agreement with experiment.     

Unexpected geometrical effects on paramagnetic spin-orbit and spin-dipolar 2J(FF) couplings

The second-rank tensor character of the paramagnetic spin-orbit and spin-dipolar contributions to nuclear spin-spin coupling constants is usually ignored when NMR measurements are carried out in the isotropic phase. However, in this study it is shown that isotropic (2)J(FF) couplings strongly depend on the relative orientation of the C-F bonds containing the coupling nuclei and the eigenvectors of such tensors. Predictions about such effect are obtained using a qualitative approach based on the polarization propagator formalism at the RPA, and results are corroborated performing high-level ab initio spin-spin coupling calculations at the SOPPA(CCSD)/EPR-III//MP2/EPR-III level in a model system. It is highlighted that no calculations at the RPA level were carried out in this work. The quite promising results reported in this paper suggest that similar properties are expected to hold for the second-rank nuclear magnetic shielding tensor.     

Elucidation of the electronic structure of molecules with the help of NMR spin-spin coupling constants: the FH molecule

It is demonstrated how the one-bond NMR spin-spin coupling constant (SSCC) (1)J(FH) can be used as a source of information on the electronic structure of the FH molecule. For this purpose, the best possible agreement between measured and calculated SSCC is achieved by large basis set coupled perturbed density functional theory calculations. Then, the calculated value is dissected into its four Ramsey terms: Fermi contact, the paramagnetic spin-orbit term, the diamagnetic spin-orbit term, and the spin dipole term, which in turn are decomposed into orbital contributions and then described by their spin densities and orbital current densities. In this way, the SSCC gives detailed information about the electronegativity of F, the bond polarity, the bond polarizability, the volume and the polarizability of sigma and pi lone pair orbitals, the s- or p-character of the bond orbital, the nature of the LUMO, and the density distribution around F.     

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.     

The effect of spin-orbit coupling on the surface dynamical properties and electron-phonon interaction of Tl(0001)

We present an ab initio study of the effect of spin-orbit coupling on the dynamical properties of the Tl(0001) surface as well as on the electron-phonon interaction at the surface. The calculations based on density-functional theory were carried out using a linear response approach and a mixed-basis pseudopotential method. It is shown that the spin-orbit effects on the phonon spectrum and the electron-phonon interaction at the Fermi level of the surface are weak but conspire to a reduction in the electron-phonon coupling strength by 16%.     

DFT study of the NMR properties of xenon in covalent compounds and van der waals complexes-implications for the use of 129Xe as a molecular probe

The NMR properties (chemical shift and spin-spin coupling constants) of (129)Xe in covalent compounds and weakly bound complexes have been investigated by DFT methods including relativistic effects. For covalent species, a good agreement between experimental and calculated results is achieved without scalar relativistic effects, but their inclusion (with a triple-zeta, double-polarization basis set) leads to some improvement in the quality of the correlation. The spin-orbit coupling term has a significant effect on the shielding constant, but makes a small contribution to the chemical shift. Coupling constants contain substantial contributions from the Fermi contact and paramagnetic spin-orbit terms; unlike light nuclei the spin-dipole term is also large, whereas the diamagnetic spin-orbit term is negligible. For van der Waals dimers, the dependence of the xenon chemical shift and anisotropy is calculated as a function of the distance. Small (<1 Hz) but non-negligible through-space coupling constants between (129)Xe and (13)C or (1)H are predicted. Much larger couplings, of the order of few Hz, are calculated between xenon and (17)O in a model silicate residue.     

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.     

Vicinal fluorine-fluorine coupling constants: Fourier analysis

Stereochemical dependences of vicinal fluorine-fluorine nuclear magnetic resonance coupling constants (3JFF) have been studied with the multiconfigurational self-consistent field in the restricted active space approach, with the second-order polarization propagator approximation (SOPPA), and with density functional theory. The SOPPA results show the best overall agreement with experimental couplings. The relationship with the dihedral angle between the coupled fluorines has been studied by Fourier analysis, the result is very different from that of proton-proton couplings. The Fourier coefficients do not resemble those of a typical Karplus equation. The four nonrelativistic contributions to the coupling constants of 1,2-difluoroethane configurations have been studied separately showing that up to six Fourier coefficients are required to reproduce the calculated values satisfactorily. Comparison with Fourier coefficients for matching hydrogen fluoride dimer configurations suggests that the higher order Fourier coefficients (Cn> or =3) originate mainly from through-space Fermi contact interaction. The through-space interaction is the main reason 3JFF do not follow the Karplus equation.     

Experimental, SOPPA(CCSD), and DFT analysis of substitutent effects on NMR 1JCF coupling constants in fluorobenzene derivatives

Interesting insight into the electronic molecular structure changes associated with substituent effects on the Fermi contact (FC) and paramagnetic spin-orbit (PSO) terms of (1)J(CF) NMR coupling constants (SSCCs) in o-X-, m-X-, and p-X-fluorobenzenes (X = NH(2); NO(2)) is presented. The formulation of this approach is based on the influence of different conjugative and hyperconjugative interactions on a second-order property, which can be qualitatively predicted if it is known how they affect the main virtual excitations entering into that second-order property. A set of consistent approximations are introduced in order to analyze the behavior of occupied and virtual orbitals, which define some experimental trends for (1)J(CF) spin-spin coupling constants. In addition, DFT hybrid functionals were used, and a similar degree of confidence to compute the (1)J(CF) with those observed for the SOPPA(CCSD) method was obtained. The (1)J(CF) SSCCs for ezetimibe, a commercially fluorinated drug used to reduce cholesterol levels, were measured and DFT-calculated, and the qualitative approach quoted above was applied. As a byproduct, a possible method to determine experimentally a significant PSO contribution to (1)J(CF) SSCCs is discussed.     

Basis set and electron correlation effects on the polarizability and second hyperpolarizability of model open-shell pi-conjugated systems

The basis set and electron correlation effects on the static polarizability (alpha) and second hyperpolarizability (gamma) are investigated ab initio for two model open-shell pi-conjugated systems, the C(5)H(7) radical and the C(6)H(8) radical cation in their doublet state. Basis set investigations evidence that the linear and nonlinear responses of the radical cation necessitate the use of a less extended basis set than its neutral analog. Indeed, double-zeta-type basis sets supplemented by a set of d polarization functions but no diffuse functions already provide accurate (hyper)polarizabilities for C(6)H(8) whereas diffuse functions are compulsory for C(5)H(7), in particular, p diffuse functions. In addition to the 6-31G(*)+pd basis set, basis sets resulting from removing not necessary diffuse functions from the augmented correlation consistent polarized valence double zeta basis set have been shown to provide (hyper)polarizability values of similar quality as more extended basis sets such as augmented correlation consistent polarized valence triple zeta and doubly augmented correlation consistent polarized valence double zeta. Using the selected atomic basis sets, the (hyper)polarizabilities of these two model compounds are calculated at different levels of approximation in order to assess the impact of including electron correlation. As a function of the method of calculation antiparallel and parallel variations have been demonstrated for alpha and gamma of the two model compounds, respectively. For the polarizability, the unrestricted Hartree-Fock and unrestricted second-order M?ller-Plesset methods bracket the reference value obtained at the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples level whereas the projected unrestricted second-order M?ller-Plesset results are in much closer agreement with the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples values than the projected unrestricted Hartree-Fock results. Moreover, the differences between the restricted open-shell Hartree-Fock and restricted open-shell second-order M?ller-Plesset methods are small. In what concerns the second hyperpolarizability, the unrestricted Hartree-Fock and unrestricted second-order M?ller-Plesset values remain of similar quality while using spin-projected schemes fails for the charged system but performs nicely for the neutral one. The restricted open-shell schemes, and especially the restricted open-shell second-order M?ller-Plesset method, provide for both compounds gamma values close to the results obtained at the unrestricted coupled cluster level including singles and doubles with a perturbative inclusion of the triples. Thus, to obtain well-converged alpha and gamma values at low-order electron correlation levels, the removal of spin contamination is a necessary but not a sufficient condition. Density-functional theory calculations of alpha and gamma have also been carried out using several exchange-correlation functionals. Those employing hybrid exchange-correlation functionals have been shown to reproduce fairly well the reference coupled cluster polarizability and second hyperpolarizability values. In addition, inclusion of Hartree-Fock exchange is of major importance for determining accurate polarizability whereas for the second hyperpolarizability the gradient corrections are large.