Abstract:Relativistic magnetic hyperfine interaction Hamiltonian based on the Douglas-Kroll-Hess (DKH) theory up to the second order is implemented within the ab initio multireference methods including spin-orbit coupling in the Molcas/OpenMolcas package. This implementation is applied to calculate relativistic hyperfine coupling (HFC) parameters for atomic systems and diatomic radicals with valence s or d orbitals by systematically varying active space size in the restricted active space self-consistent field (RASSCF) formalism with restricted active space state interaction (RASSI) for spin-orbit coupling. The DKH relativistic treatment of the hyperfine interaction reduces the Fermi contact contribution to the HFC due to the presence of kinetic factors that regularize the singularity of the Dirac delta function in the nonrelativitic Fermi contact operator. This effect is more prominent for heavier nuclei. As the active space size increases, the relativistic correction of the Fermi contact contribution converges well to the experimental data for light and moderately heavy nuclei. The relativistic correction, however, does not significantly affect the spin-dipole contribution to the hyperfine interaction. In addition to the atomic and molecular systems, the implementation is applied to calculate the relativistic HFC parameters for large trivalent and divalent Tb-based single-molecule magnets (SMMs) such as Tb(III)Pc$_2$ and Tb(II)(Cp$^\text{iPr5}$)$_2$ without ligand truncation using well-converged basis sets. In particular, for the divalent SMM which has an unpaired valence 6s/5d hybrid orbital, the relativistic treatment of HFC is crucial for a proper description of the Fermi contact contribution. Even with the relativistic hyperfine Hamiltonian, the divalent SMM is shown to exhibit strong tunability of HFC via an external electric field (i.e., strong hyperfine Stark effect).

On the calculation of second-order magnetic properties using subsystem approaches in the relativistic framework

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

Relativistic Douglas-Kroll-Hess Calculations of Hyperfine Interactions within First Principles Multireference Methods

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

Relativistic magnetic interactions from non-orthogonal basis sets

Fully Relativistic Theories and Methods for NMR Parameters

Calculation of DFT-GIAO NMR shifts with the inclusion of spin-orbit coupling

Ab initio Calculations in Atoms, Molecules, and Solids, Treating Spin-Orbit Coupling and Electron Interaction on Equal Footing

Effective Hamiltonian Crystal Field for Magnetic Interactions in Polynuclear Transition Metal Complexes. Sequential Derivation and Exemplary Numerical Estimates

First-principles calculation of the parameters used by atomistic magnetic simulations

Accurate NMR Shieldings with σ-Functionals

Systematic studies on the computation of nuclear magnetic resonance shielding constants and chemical shifts: the density functional models.

All-electron density functional calculations for electron and nuclear spin interactions in molecules and solids

Investigating solvent effects on the magnetic properties of molybdate ions (MoO$_{4}^{2-}$) with relativistic embedding

Calculation of molecular g-tensors by sampling spin orientations of generalised Hartree-Fock states

Spatial contributions to nuclear magnetic shieldings

Delving into Theoretical and Computational Considerations for Accurate Calculation of Chemical Shifts in Paramagnetic Transition Metal Systems using Quantum Chemical Methods

Non-perturbative Many-Body Treatment of Molecular Magnets

Spin-free formulation of the multireference driven similarity renormalization group: A benchmark study of first-row diatomic molecules and spin-crossover energetics

Calculation of Magnetic Shielding Constants with meta-GGA Functionals Employing the Multipole-Accelerated Resolution of the Identity: Implementation and Assessment of Accuracy and Efficiency

The Synergy between Nuclear Magnetic Resonance and Density Functional Theory Calculations