The four-component DFT method for the calculation of the EPR g-tensor using a restricted magnetically balanced basis and London atomic orbitals

ABSTRACT Four-component relativistic treatments of the electron paramagnetic resonance g-tensor have so far been based on a common gauge origin and a restricted kinetically balanced basis. The results of such calculations are prone to exhibit a dependence on the choice of the gauge origin for the vector potential associated with uniform magnetic field and a related dependence on the basis set quality. In this work, this gauge problem is addressed by a distributed-origin scheme based on the London atomic orbitals, also called gauge-including atomic orbitals (GIAOs), which have proven to be a practical approach for calculations of other magnetic properties. Furthermore, in the four-component relativistic domain, it has previously been shown that a restricted magnetically balanced (RMB) basis for the small component of the four-component wavefunctions is necessary for achieving robust convergence with regard to the basis set size. We present the implementation of a four-component density functional theory (DFT) method for calculating the g-tensor, incorporating both the GIAOs and RMB basis and based on the Dirac–Coulomb Hamiltonian. The approach utilizes the state-of-the-art noncollinear Kramers-unrestricted DFT methodology to achieve rotationally invariant results and inclusion of spin-polarization effects in the calculation. We also show that the gauge dependence of the results obtained is connected to the nonvanishing integral of the current density in a finite basis, explain why the results of cluster calculations exhibit surprisingly low gauge dependence, and demonstrate that the gauge problem disappears for systems with certain point-group symmetries

Linking the Character of the Metal-Ligand Bond to the Ligand NMR Shielding in Transition-Metal Complexes: NMR Contributions from Spin-Orbit Coupling

Relativistic effects significantly affect various spectroscopic properties of compounds containing heavy elements. Particularly in Nuclear Magnetic Resonance (NMR) spectroscopy, the heavy atoms strongly influence the NMR shielding constants of neighboring light atoms. In this account we analyze paramagnetic contributions to NMR shielding constants and their modulation by relativistic spin-orbit effects in a series of transition-metal complexes of Pt(II), Au(I), Au(III), and Hg(II). We show how the paramagnetic NMR shielding and spin-orbit effects relate to the character of the metal-ligand (M-L) bond. A correlation between the (back)-donation character of the M-L bond in d10 Au(I) complexes and the propagation of the spin-orbit (SO) effects from M to L through the M-L bond influencing the ligand NMR shielding via the Fermi-contact mechanism is found and rationalized by using third-order perturbation theory. The SO effects on the ligand NMR shielding are demonstrated to be driven by both the electronic structure of M and the nature of the trans ligand, sharing the σ-bonding metal orbital with the NMR spectator atom L. The deshielding paramagnetic contribution is linked to the σ-type M-L bonding orbitals, which are notably affected by the trans ligand. The SO deshielding role of σ-type orbitals is enhanced in d10 Hg(II) complexes with the Hg 6p atomic orbital involved in the M-L bonding. In contrast, in d8 Pt(II) complexes, occupied π-type orbitals play a dominant role in the SO-altered magnetic couplings due to the accessibility of vacant antibonding σ-type MOs in formally open 5d-shell (d8). This results in a significant SO shielding at the light atom. The energy- and composition-modulation of σ- vs π-type orbitals by spin-orbit coupling is rationalized and supported by visualizing the SO-induced changes in the electron density around the metal and light atoms (spin-orbit electron deformation density, SO-EDD). © 2017 American Chemical Society.Czech Science Foundation [16-05961S, 15-09381S]; Ministry of Education, Youth and Sports of the Czech Republic [LQ1601, LO1504]; multilateral cooperation project [8X17009]; SASPRO Program [1563/03/02]; European Union; Slovak Academy of Sciences; Grant Agency of the Ministry of Education of the Slovak Republic; Slovak Academy of Sciences VEGA [2/0116/17]; Research Council of Norway [179568]; Norwegian supercomputing program NOTUR [NN4654K

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, DdSO, which roughly corresponds to an increase in DdSO 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.Czech Science Foundation [15-09381S, 14-03564S]; European Regional Development Fund [CZ.1.05/1.1.00/02.0068]; Research Council of Norway through a Centre of Excellence [179568, 214095, 177558]; Czech-Norway mobility grant from Norway Funds [NF-CZ07-MOP-3-245-2015]; program Center CERIT Scientific Cloud, part of the Operational Program Research and Development for Innovations [CZ.1.05/3.2.00/08.0144

The existence and unambiguity of the principal axis system of the EPR tensors

Although the role of the electron paramagnetic resonance (EPR) g-tensor and hyperfine coupling tensor in the EPR effective spin Hamiltonian is discussed extensively in many textbooks, certain aspects of the theory are missing. In this text we will cover those gaps and thus provide a comprehensive theory about the existence of principal axes of the EPR tensors. However, an important observation is that both g- and a-tensors have two sets of principal axes -- one in the real and one in the fictitious spin space -- and, in fact, are not tensors. Moreover, we present arguments based on the group theory why only eigenvalues of the G-tensor, \mb{G} = \mb{g}\mb{g}^{\!\mathsf{T}}, and the sign of the determinant of the g-tensor are observable quantities (an analogical situation also holds for the hyperfine coupling tensor). We keep the number of assumptions to a minimum and thus the theory is applicable in the framework of the Dirac--Coulomb--Breit Hamiltonian and for any spatial symmetry of the system

Four-component relativistic time-dependent density-functional theory using a stable noncollinear DFT ansatz applicable to both closed- and open-shell systems

We present a formulation of relativistic linear response time-dependent density functional theory for the calculation of electronic excitation energies in the framework of the four-component Dirac-Coulomb Hamiltonian. This approach is based on the noncollinear ansatz originally developed by Scalmani and Frisch [J. Chem. Theory Comput. 8, 2193 (2012)] and improves upon the past treatment of the limit cases in which the spin density approaches zero. As a result of these improvements, the presented approach is capable of treating both closed- and open-shell reference states. Robust convergence of the Davidson-Olsen eigenproblem algorithm for open-shell reference states was achieved through the use of a solver which considers both left and right eigenvectors. The applicability of the present methodology on both closed- and open-shell reference states is demonstrated on calculations of low-lying excitation energies for Group 3 atomic systems (Sc3+–Ac3+) with nondegenerate ground states, as well as for Group 11 atomic systems (Cu–Rg) and octahedral actinide complexes (PaCl2−6, UCl−6, and NpF6) with effective doublet ground states

Relativistic four-component linear damped response TDDFT for electronic absorption and circular dichroism calculations

We present a detailed theory, implementation, and a benchmark study of a linear damped response time-dependent density functional theory (TDDFT) based on the relativistic four-component (4c) Dirac–Kohn–Sham formalism using the restricted kinetic balance condition for the small-component basis and a noncollinear exchange–correlation kernel. The damped response equations are solved by means of a multifrequency iterative subspace solver utilizing decomposition of the equations according to Hermitian and time-reversal symmetry. This partitioning leads to robust convergence, and the detailed algorithm of the solver for relativistic multicomponent wavefunctions is also presented. The solutions are then used to calculate the linear electric- and magnetic-dipole responses of molecular systems to an electric perturbation, leading to frequency-dependent dipole polarizabilities, electronic absorption, circular dichroism (ECD), and optical rotatory dispersion (ORD) spectra. The methodology has been implemented in the relativistic spectroscopy DFT program ReSpect, and its performance was assessed on a model series of dimethylchalcogeniranes, C4H8X (X = O, S, Se, Te, Po, Lv), and on larger transition metal complexes that had been studied experimentally, [M(phen)3]3+ (M = Fe, Ru, Os). These are the first 4c damped linear response TDDFT calculations of ECD and ORD presented in the literature

How does relativity affect magnetically induced currents?

Magnetically induced probability currents in molecules are studied in relativistic theory. Spin–orbit coupling (SOC) enhances the curvature and gives rise to a previously unobserved current cusp in AuH or small bulge-like distortions in HgH2 at the proton positions. The origin of this curvature is magnetically induced spin-density arising from SOC in the relativistic description