Note: This article is part of the JCP Special Topic on Electronic Structure Software.

The effect of a finite size model for both the nuclear charge and magnetic moment distributions on calculated EPR hyperfine structure have been studied using a relativistic four-component method based on density functional theory. This approach employs a restricted kinetically balanced basis (mDKS-RKB) and includes spin-polarization using noncollinear spin-density exchange-correlation functionals in the unrestricted fashion. Benchmark calculations have been carried out for a number of small molecules containing Zn, Cd, Ag, and Hg. The present results are compared with those obtained at the Douglas-Kroll-Hess second order (DKH-2) method. The dependence of the results on the quality of the orbital and auxiliary basis sets has been studied. It was found that some basis sets contain irregularities that deteriorate the results. Especial care has to be taken also on the construction of the auxiliary basis for fitting the total electron and spin-densities.

A scalar relativistic method to calculate hyperfine coupling tensors at the Douglas-Kroll-Hess level has been extended to incorporate a finite-size nucleus model using a Gaussian charge and magnetic moment distribution. Density functional calculations at gradient-corrected and hybrid functional levels have been carried out for the group 11 atoms and for a set of small group 12 molecules, comparing nonrelativistic as well as scalar relativistic results at second-order Douglas-Kroll-Hess level with and without finite-size nucleus. While nonrelativistic calculations underestimate isotropic hyperfine couplings increasingly with increasing nuclear charge, scalar relativistic calculations with point nucleus provide somewhat overestimated values. Inclusion of the finite-size nuclear model in the calculation of the wavefunction, and in the transformed hyperfine operators both decrease the magnitude of the hyperfine couplings. The effects, which are cumulative, improve agreement with experiment.

No abstract

We demonstrate that the apparent disagreement between experimental determinations and four-component relativistic calculations of the absolute shielding constants of heavy nuclei is due to the breakdown of the commonly assumed relation between the electronic contribution to the nuclear spin-rotation constants and the paramagnetic contribution to the NMR shielding constants. We demonstrate that this breakdown has significant consequences for the absolute shielding constant of (119)Sn, leading to errors of about 1000 ppm. As a consequence, we expect that many absolute shielding constants of heavy nuclei will be in need of revision.

We present an implementation of the nuclear spin-rotation (SR) constants based on the relativistic four-component Dirac-Coulomb Hamiltonian. This formalism has been implemented in the framework of the Hartree-Fock and Kohn-Sham theory, allowing assessment of both pure and hybrid exchange-correlation functionals. In the density-functional theory (DFT) implementation of the response equations, a noncollinear generalized gradient approximation (GGA) has been used. The present approach enforces a restricted kinetic balance condition for the small-component basis at the integral level, leading to very efficient calculations of the property. We apply the methodology to study relativistic effects on the spin-rotation constants by performing calculations on XHn (n = 1-4) for all elements X in the p-block of the periodic table and comparing the effects of relativity on the nuclear SR tensors to that observed for the nuclear magnetic shielding tensors. Correlation effects as described by the density-functional theory are shown to be significant for the spin-rotation constants, whereas the differences between the use of GGA and hybrid density functionals are much smaller. Our calculated relativistic spin-rotation constants at the DFT level of theory are only in fair agreement with available experimental data. It is shown that the scaling of the relativistic effects for the spin-rotation constants (varying between Z(3.8) and Z(4.5)) is as strong as for the chemical shieldings but with a much smaller prefactor.

The spin-rotation and nuclear magnetic shielding constants are analysed for both nuclei in the HCl molecule. Nonrelativistic ab initio calculations at the CCSD(T) level of approximation show that it is essential to include relativistic effects to obtain spin-rotation constants consistent with accurate experimental data. Our best estimates for the spin-rotation constants of (1)H(35)Cl are CCl = -53.914 kHz and C(H) = 42.672 kHz (for the lowest rovibrational level). For the chlorine shielding constant, the ab initio value computed including the relativistic corrections, σ(Cl) = 976.202 ppm, provides a new absolute shielding scale; for hydrogen we find σ(H) = 31.403 ppm (both at 300 K). Combining the theoretical results with our new gas-phase NMR experimental data allows us to improve the accuracy of the magnetic dipole moments of both chlorine isotopes. For the hydrogen shielding constant, including relativistic effects yields better agreement between experimental and computed values.

By using complementary experimental techniques and first-principles theoretical calculations, magnetic anisotropy in a series of five hexacoordinated nickel(II) complexes possessing a symmetry close to C , has been investigated. Four complexes have the general formula [Ni(bpy)X ] (bpy=2,2'-bipyridine; X =bpy (1), (NCS ) (2), C O (3), NO (4)). In the fifth complex, [Ni(HIM -py) (NO )] (5; HIM -py=2-(2-pyridyl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazolyl-1-hydroxy), which was reported previously, the two bpy bidentate ligands were replaced by HIM -py. Analysis of the high-field, high-frequency electronic paramagnetic resonance (HF-HFEPR) spectra and magnetization data leads to the determination of the spin Hamiltonian parameters. The D parameter, corresponding to the axial magnetic anisotropy, was negative (Ising type) for the five compounds and ranged from -1 to -10 cm . First-principles SO-CASPT2 calculations have been performed to estimate these parameters and rationalize the experimental values. From calculations, the easy axis of magnetization is in two different directions for complexes 2 and 3, on one hand, and 4 and 5, on the other hand. A new method is proposed to calculate the g tensor for systems with S=1. The spin Hamiltonian parameters (D (axial), E (rhombic), and g ) are rationalized in terms of ordering of the 3 d orbitals. According to this orbital model, it can be shown that 1) the large magnetic anisotropy of 4 and 5 arises from splitting of the e -like orbitals and is due to the difference in the σ-donor strength of NO and bpy or HIM -py, whereas the difference in anisotropy between the two compounds is due to splitting of the t -like orbitals; and 2) the anisotropy of complexes 1-3 arises from the small splitting of the t -like orbitals. The direction of the anisotropy axis can be rationalized by the proposed orbital model.