Stanislav Komorovsky

Relativistic four-component DFT calculations of 1H NMR chemical shifts in transition-metal hydride complexes: unusual high-field shifts beyond the Buckingham-Stephens model.

State-of-the-art relativistic four-component DFT-GIAO-based calculations of (1)H NMR chemical shifts of a series of 3d, 4d, and 5d transition-metal hydrides have revealed significant spin-orbit-induced heavy atom effects on the hydride shifts, in particular for several 4d and 5d complexes. The spin-orbit (SO) effects provide substantial, in some cases even the dominant, contributions to the well-known characteristic high-field hydride shifts of complexes with a partially filled d-shell, and thereby augment the Buckingham-Stephens model of off-center paramagnetic ring currents. In contrast, complexes with a 4d(10) and 5d(10) configuration exhibit large deshielding SO effects on their hydride (1)H NMR shifts. The differences between the two classes of complexes are attributed to the dominance of π-type d-orbitals for the true transition-metal systems compared to σ-type orbitals for the d(10) systems.

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

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.

Four-Component Relativistic Density Functional Theory Calculations of NMR Shielding Tensors for Paramagnetic Systems

A four-component relativistic method for the calculation of NMR shielding constants of paramagnetic doublet systems has been developed and implemented in the ReSpect program package. The method uses a Kramer unrestricted noncollinear formulation of density functional theory (DFT), providing the best DFT framework for property calculations of open-shell species. The evaluation of paramagnetic nuclear magnetic resonance (pNMR) tensors reduces to the calculation of electronic g tensors, hyperfine coupling tensors, and NMR shielding tensors. For all properties, modern four-component formulations were adopted. The use of both restricted kinetically and magnetically balanced basis sets along with gauge-including atomic orbitals ensures rapid basis-set convergence. These approaches are exact in the framework of the Dirac-Coulomb Hamiltonian, thus providing useful reference data for more approximate methods. Benchmark calculations on Ru(III) complexes demonstrate good performance of the method in reproducing experimental data and also its applicability to chemically relevant medium-sized systems. Decomposition of the temperature-dependent part of the pNMR tensor into the traditional contact and pseudocontact terms is proposed.

Excitation Energies from Real-Time Propagation of the Four-Component Dirac–Kohn–Sham Equation

We report the first implementation of real-time time-dependent density functional theory (RT-TDDFT) at the relativistic four-component level of theory. In contrast to the perturbative linear-response TDDFT approach (LR-TDDFT), the RT-TDDFT approach performs an explicit time propagation of the Dirac-Kohn-Sham density matrix, offering the possibility to simulate molecular spectroscopies involving strong electromagnetic fields while, at the same time, treating relativistic scalar and spin-orbit corrections variationally. The implementation is based on the matrix representation of the Dirac-Coulomb Hamiltonian in the basis of restricted kinetically balanced Gaussian-type functions, exploiting the noncollinear Kramers unrestricted formalism implemented in the program ReSpect. We also present an analytic form for the delta-type impulse commonly used in RT-TDDFT calculations, as well as a dipole-weighted transition matrix analysis, facilitating the interpretation of spectral transitions in terms of ground-state molecular orbitals. The possibilities offered by the methodology are illustrated by investigating vertical excitation energies and oscillator strengths for ground-state to excited-state transitions in the Group 12 atoms and in heavy-element hydrides. The accuracy of the method is assessed by comparing the excitation energies obtained with earlier relativistic linear response TDDFT results and available experimental data.

Four-Component Relativistic Density-Functional Theory Calculations of Nuclear Spin–Rotation Constants: Relativistic Effects in p-Block Hydrides

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.

Spin-rotation and NMR shielding constants in HCl

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.

Communication: The absolute shielding scales of oxygen and sulfur revisited

We present an updated semi-experimental absolute shielding scale for the (17)O and (33)S nuclei. These new shielding scales are based on accurate rotational microwave data for the spin-rotation constants of H2(17)O [Puzzarini et al., J. Chem. Phys. 131, 234304 (2009)], C(17)O [Cazzoli et al., Phys. Chem. Chem. Phys. 4, 3575 (2002)], and H2(33)S [Helgaker et al., J. Chem. Phys. 139, 244308 (2013)] corrected both for vibrational and temperature effects estimated at the CCSD(T) level of theory as well as for the relativistic corrections to the relation between the spin-rotation constant and the absolute shielding constant. Our best estimate for the oxygen shielding constants of H2(17)O is 328.4(3) ppm and for C(17)O -59.05(59) ppm. The relativistic correction for the sulfur shielding of H2(33)S amounts to 3.3%, and the new sulfur shielding constant for this molecule is 742.9(4.6) ppm.

Four-Component Relativistic DFT Calculations of (13)C Chemical Shifts of Halogenated Natural Substances.

We have calculated the (13)C NMR chemical shifts of a large ensemble of halogenated organic molecules (81 molecules for a total of 250 experimental (13)C NMR data at four different levels of theory), ranging from small rigid organic compounds, used to benchmark the performance of various levels of theory, to natural substances of marine origin with conformational degrees of freedom. Carbon atoms bonded to heavy halogen atoms, particularly bromine and iodine, are known to be rather challenging when it comes to the prediction of their chemical shifts by quantum methods, due to relativistic effects. In this paper, we have applied the state-of-the-art four-component relativistic density functional theory for the prediction of such NMR properties and compared the performance with two-component and nonrelativistic methods. Our results highlight the necessity to include relativistic corrections within a four-component description for the most accurate prediction of the NMR properties of halogenated organic substances.

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 (structure, solvent, and relativistic treatment) was evaluated for square-planar 4d and 5d transition-metal complexes. The DFT method for calculating the structures was calibrated using a cluster approach and compared to X-ray geometries, with the PBE0 functional (def2-TZVPP basis set) providing the best results, followed closely by the hybrid TPSSH and the MN12SX functionals. Calculations of the NMR chemical shifts using the two-component (2c, Zeroth-Order Regular Approximation as implemented in the ADF package) and four-component (4c, Dirac-Coulomb as implemented in the ReSpect code) relativistic approaches were performed to analyze and demonstrate the importance of solvent corrections (2c) as well as a proper treatment of relativistic effects (4c). The importance of increased exact-exchange admixture in the functional (here PBE0) for reproducing the experimental data using the current implementation of the 2c approach is partly rationalized as a compensation for the missing exchange-correlation response kernel. The kernel contribution was identified to be about 15-20% of the spin-orbit-induced NMR chemical shift, ΔδSO, which roughly corresponds to an increase in ΔδSO introduced by the artificially increased exact-exchange admixture in the functional. Finally, the role of individual effects (geometry, solvent, relativity) in the NMR chemical shift is discussed in selected complexes. Although a fully relativistic DFT approach is still awaiting the implementation of GIAOs for hybrid functionals and an implicit solvent model, it nevertheless provides reliable NMR chemical shift data at an affordable computational cost. It is expected to outperform the 2c approach, in particular for the calculation of NMR parameters in heavy-element compounds.

Absolute NMR shielding scales and nuclear spin–rotation constants in 175LuX and 197AuX (X = 19F, 35Cl, 79Br and 127I)

We present nuclear spin-rotation constants, absolute nuclear magnetic resonance (NMR) shielding constants, and shielding spans of all the nuclei in (175)LuX and (197)AuX (X = (19)F, (35)Cl, (79)Br, (127)I), calculated using coupled-cluster singles-and-doubles with a perturbative triples (CCSD(T)) correction theory, four-component relativistic density functional theory (relativistic DFT), and non-relativistic DFT. The total nuclear spin-rotation constants determined by adding the relativistic corrections obtained from DFT calculations to the CCSD(T) values are in general in agreement with available experimental data, indicating that the computational approach followed in this study allows us to predict reliable results for the unknown spin-rotation constants in these molecules. The total NMR absolute shielding constants are determined for all the nuclei following the same approach as that applied for the nuclear spin-rotation constants. In most of the molecules, relativistic effects significantly change the computed shielding constants, demonstrating that straightforward application of the non-relativistic formula relating the electronic contribution to the nuclear spin-rotation constants and the paramagnetic contribution to the shielding constants does not yield correct results. We also analyze the origin of the unusually large absolute shielding constant and its relativistic correction of gold in AuF compared to the other gold monohalides.