Interactions between individual atoms underpin the structure and behavior of matter. These interactions govern atomic positions and dynamics, as well as the organization of electronsparticularly in the frontier region. Because electrons lie at the core of chemical phenomena, numerous theoretical frameworks have been developed to rationalize the molecular structure and properties. Electronic motion within molecules and the resulting induced currents provide powerful probes of the molecular or supramolecular structure, building on and going beyond molecular orbital and valence bond theories. In particular, current density offers a spatially-resolved description of the electronic response to external perturbations, enabling direct analysis of electron delocalization and magnetic response in molecular systems. In this work, the effect of relativistic spin-orbit (SO) coupling on the strength and topology of the magnetically induced current density (MICD) is analyzed in depth for a series of model heavy-atom hydrides at the four-component Dirac-Kohn-Sham level. For the most simple molecules, TlH, HAt, and AuH, we demonstrate a connection between the SO effects on the molecular geometry, strength and topology of MICDs, and ligand 1H NMR shielding. For model HMX molecules, where M = AuI, HgII; X = F, Cl, Ph, CH3, H, SiH3, BH2, the hydride deshielding due to the slight elongation of the M-H bond upon increasing the trans-ligand influence (TLI) of X is shown to be marginal when compared to that originating from the electronic SO effect. In particular, the inclusion of SO effects gives rise to highly localized paratropic MICD vortices on the hydride position of those complexes bearing strong TLI ligands. Our results disprove the previously proposed governing role of the current around the metal atom (similar to the classical Buckingham-Stephens model for transition metal hydrides) associated with TLI-induced variations in the metal-hydrogen bond length in determining the characteristic ligand 1H NMR shifts.

Quaternion algebra provides a natural representation of time-reversal-symmetric matrix structures in relativistic electronic-structure theory, whereas complementary time-reversal-antisymmetric structures extend this representation to complex quaternions, or biquaternions. Here, we explicitly exploit biquaternion algebra for fundamental objects, including operators, kinetically and magnetically balanced basis functions, and their expectation values, within a unified framework encompassing real, complex, and real quaternion subalgebras as special cases. To realize this framework computationally, we have developed HMATLIB, a biquaternion matrix library implemented within the ReSpect package for pure CPU and hybrid CPU/GPU execution. A key development is a bilinear algorithm for matrix-valued biquaternion multiplication that reduces the number of real matrix–matrix multiplications from 64 to 24 compared with the conventional component-wise approach, while the biquaternion representation effectively doubles the maximum accessible matrix dimension under the same 32-bit indexing constraint. Numerical benchmarks on modern CPU and GPU architectures demonstrate that the biquaternion formulation consistently outperforms its isomorphic complex-algebra counterpart. For the largest matrices reported, hybrid CPU/GPU execution accelerates bilinear matrix multiplication by approximately 5-8 times over pure CPU execution and matrix diagonalization by approximately 12 and 61 times relative to CPU oneAPI MKL and NVHPC OpenBLAS, respectively. These results establish biquaternion algebra as a general and computationally advantageous framework for modern relativistic electronic-structure calculations.

Rovibrational energies, wave functions, and Raman transition moments are reported for the lowest-energy states of the H_{3}^{+} molecular ion including the magnetic couplings of the proton spins and molecular rotation in the presence of a weak external magnetic field. The rovibrational-hyperfine-Zeeman Hamiltonian matrix is constructed and diagonalized using the rovibrational eigenstates and the proton spin functions. The developed methodology can be used to compute hyperfine-Zeeman effects also for higher-energy rovibrational excitations of H_{3}^{+} and other closed-shell polyatomic molecules. These developments will guide future experiments extending quantum logic spectroscopy to polyatomic systems.

Since its inception, the ReSpect program has been evolving to provide powerful tools for simulating spectroscopic processes and exploring emerging research areas, all while incorporating relativistic effects, particularly spin-orbit interactions, in a fully variational manner. Recent developments have focused on exact two-component (X2C) Hamiltonian models that go beyond the standard one-electron X2C approach by incorporating two-electron picture-change corrections. This paper presents the theoretical foundations of two distinct atomic mean-field X2C models, amfX2C and extended eamfX2C, which offer computationally efficient and accurate alternatives to fully relativistic four-component methods. These models enable simulations of complex phenomena, such as time-resolved pump-probe spectroscopies and cavity-modified molecular properties, which would otherwise be computationally prohibitive. ReSpect continues to evolve, providing state-of-the-art quantum chemical methods and post-processing tools, all available free of charge through our website, www.respectprogram.org, to support researchers exploring relativistic effects across various scientific disciplines.

A summary of current state-of-the-art relativistic theory for the calculation of the EPR and (paramagnetic) NMR parameters is presented. The theory is valid within the framework of the Dirac–Coulomb–Breit Hamiltonian, as only the time-reversal symmetry of the relativistic operators is assumed. It is therefore applicable to every approximate relativistic theory that is derived from the Dirac–Coulomb–Breit Hamiltonian. This chapter looks at both the exact-state theory, i.e., theory applicable within ab initio wavefunction methodologies, and DFT approaches. Most of the theory in this chapter has been described elsewhere. However, some of the material represents theoretical advances that are presented here for the first time: (1) it is argued that the Helmholtz free energy and its derivatives are the correct starting point for any approach to the calculation of molecular properties of any order at finite temperature; (2) the proposed definition of the orbital contribution to the pNMR tensor is now consistent with both the exact-state pNMR theory and the pNMR theory based on the parameters of the NMR and EPR effective Hamiltonians; (3) a comprehensive description of the theory of using noncollinear DFT for the calculation of the g-tensor and hyperfine coupling tensor is presented, including its connection to the exact-state EPR theory.

Magnetic resonance techniques represent a fundamental class of spectroscopic methods used in physics, chemistry, biology, and medicine. Electron paramagnetic resonance (EPR) is an extremely powerful technique for characterizing systems with an open-shell electronic nature, whereas nuclear magnetic resonance (NMR) has traditionally been used to investigate diamagnetic (closed-shell) systems. However, these two techniques are tightly connected by the electron-nucleus hyperfine interaction operating in paramagnetic (open-shell) systems. Hyperfine interaction of the nuclear spin with unpaired electron(s) induces large temperature-dependent shifts of nuclear resonance frequencies that are designated as hyperfine NMR shifts (delta(HF)). Three fundamental physical mechanisms shape the total hyperfine interaction: Fermi-contact, paramagnetic spin-orbit, and spin-dipolar. The corresponding hyperfine NMR contributions can be interpreted in terms of through-bond and through-space effects. In this Account, we provide an elemental theory behind the hyperfine interaction and NMR shifts and describe recent progress in understanding the structural and electronic principles underlying individual hyperfine terms. The Fermi-contact (FC) mechanism reflects the propagation of electron-spin density throughout the molecule and is proportional to the spin density at the nuclear position. As the imbalance in spin density can be thought of as originating at the paramagnetic metal center and being propagated to the observed nucleus via chemical bonds, FC is an excellent indicator of the bond character. The paramagnetic spin-orbit (PSO) mechanism originates in the orbital current density generated by the spin-orbit coupling interaction at the metal center. The PSO mechanism of the ligand NMR shift then reflects the transmission of the spin polarization through bonds, similar to the FC mechanism, but it also makes a substantial through-space contribution in long-range situations. In contrast, the spin-dipolar (SD) mechanism is relatively unimportant at short-range with significant spin polarization on the spectator atom. The PSO and SD mechanisms combine at long-range to form the so-called pseudocontact shift, traditionally used as a structural and dynamics probe in paramagnetic NMR (pNMR). Note that the PSO and SD terms both contribute to the isotropic NMR shift only at the relativistic spin-orbit level of theory. We demonstrate the advantages of calculating and analyzing the NMR shifts at relativistic two- and four-component levels of theory and present analytical tools and approaches based on perturbation theory. We show that paramagnetic NMR effects can be interpreted by spin-delocalization and spin-polarization mechanisms related to chemical bond concepts of electron conjugation in pi-space and hyperconjugation in sigma-space in the framework of the molecular orbital (MO) theory. Further, we discuss the effects of environment (supramolecular interactions, solvent, and crystal packing) and demonstrate applications of hyperfine shifts in determining the structure of paramagnetic Ru(III) compounds and their supramolecular host-guest complexes with macrocycles.In conclusion, we provide a short overview of possible pNMR applications in the analysis of spectra and electronic structure and perspectives in this field for a general chemical audience.

Understanding how the electronic g-tensor is linked to the electronic structure is desirable for the correct interpretation of electron paramagnetic resonance spectra. For heavy-element compounds with large spin–orbit (SO) effects, this is still not completely clear. We report our investigation of quadratic SO contributions to the g-shift in heavy transition metal complexes. We implemented third-order perturbation theory in order to analyze the contributions arising from frontier molecular spin orbitals (MSOs). We show that the dominant quadratic SO term—spin-Zeeman (SO2/SZ)—generally makes a negative contribution to the g-shift, irrespective of the particular electronic configuration or molecular symmetry. We further analyze how the SO2/SZ contribution adds to or subtracts from the linear orbital-Zeeman (SO/OZ) contribution to the individual principal components of the g-tensor. Our study suggests that the SO2/SZ mechanism decreases the anisotropy of the g-tensor in early transition metal complexes and increases it in late transition metal complexes. Finally, we apply MSO analysis to the investigation of g-tensor trends in a set of closely related Ir and Rh pincer complexes and evaluate the influence of different chemical factors (the nuclear charge of the central atom and the terminal ligand) on the magnitudes of the g-shifts. We expect our conclusions to aid the understanding of spectra in magnetic resonance investigations of heavy transition metal compounds.

X-ray absorption spectroscopy (XAS) has gained popularity in recent years as it probes matter with high spatial and elemental sensitivities. However, the theoretical modeling of XAS is a challenging task since XAS spectra feature a fine structure due to scalar (SC) and spin-orbit (SO) relativistic effects, in particular near L and M absorption edges. While full four-component (4c) calculations of XAS are nowadays feasible, there is still interest in developing approximate relativistic methods that enable XAS calculations at the two-component (2c) level while maintaining the accuracy of the parent 4c approach. In this article we present theoretical and numerical insights into two simple yet accurate 2c approaches based on an (extended) atomic mean-field exact two-component Hamiltonian framework, (e)amfX2C, for the calculation of XAS using linear eigenvalue and damped response time-dependent density functional theory (TDDFT). In contrast to the commonly used one-electron X2C (1eX2C) Hamiltonian, both amfX2C and eamfX2C account for the SC and SO two-electron and exchange-correlation picture-change (PC) effects that arise from the X2C transformation. As we demonstrate on L- and M-edge XAS spectra of transition metal and actinide compounds, the absence of PC corrections in the 1eX2C approximation results in a substantial overestimation of SO splittings, whereas (e)amfX2C Hamiltonians reproduce all essential spectral features such as shape, position, and SO splitting of the 4c references in excellent agreement, while offering significant computational savings. Therefore, the (e)amfX2C PC correction models presented here constitute reliable relativistic 2c quantum-chemical approaches for modeling XAS.

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.

Relativistic effects are known to alter the chemical bonds and spectroscopic properties of heavy-element compounds. In this work, we introduce the concept of spin-orbit (SO) electronegativity of a heavy atom, as reflected by an SO-induced change in the interatomic distance between the heavy atom (HA) and a neighboring light atom (LA). We provide a transparent interpretation of these SO effects by using the concept of spin-orbit electron deformation density (SO-EDD). Spin-orbit coupling at the HA induces rearrangement of the electron density for the scalar-relativistically optimized geometry that, in turn, exerts a new force on the LA. The resulting expansion or contraction of the HA-LA bond depends on the nature and electron configuration of the HA. In addition, we quantify the change in atomic electronegativity induced by SO coupling for a series of hydrides, thereby complementing the SO-EDD picture. The trends in the SO-induced electronegativity and the HA-LA bond length across the periodic table of elements are demonstrated and interpreted, and also linked, intuitively, with the SO-induced NMR shielding at the LA.

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.

We revise formal and numerical aspects of collinear and non-collinear density functional theories in the context of a two-component self-consistent treatment of spin-orbit coupling. Theoretical and numerical analyses of the non-collinear approaches confirm their ability to yield the proper collinear limit and provide rotational invariance of the total energy for functionals in the local-density or generalized-gradient approximations (GGAs). Calculations on simple molecules corroborate the formal considerations and highlight the importance of an effective screening algorithm to provide the sufficient level of numerical stability required for a rotationally invariant implementation of non-collinear GGA functionals. The illustrative calculations provide a first numerical comparison of both previously proposed non-collinear formulations for GGA functionals. The proposed screening procedure allows us to effectively deal with points of small magnetization, which would otherwise be problematic for the evaluation of the exchange-correlation energy and/or potential for non-collinear GGA functionals. Both previously suggested formulations for the non-collinear GGA are confirmed to be adequate for total energy calculations, provided that the screening is achieved on a sufficiently fine grid. All methods are implemented in the Crystal program.

Nuclear magnetic resonance (NMR) spectroscopy of paramagnetic molecules provides detailed information about their molecular and electron-spin structure. The paramagnetic NMR spectrum is a very rich source of information about the hyperfine interaction between the atomic nuclei and the unpaired electron density. The Fermi-contact contribution to ligand hyperfine NMR shifts is particularly informative about the nature of the metal-ligand bonding and the structural arrangements of the ligands coordinated to the metal center. In this account, we provide a detailed experimental and theoretical NMR study of compounds of Cr(III) and Cu(II) coordinated with substituted acetylacetonate (acac) ligands in the solid state. For the first time, we report the experimental observation of extremely paramagnetically deshielded 13C NMR resonances for these compounds in the range of 900-1200 ppm. We demonstrate an excellent agreement between the experimental NMR shifts and those calculated using relativistic density-functional theory. Crystal packing is shown to significantly influence the NMR shifts in the solid state, as demonstrated by theoretical calculations of various supramolecular clusters. The resonances are assigned to individual atoms in octahedral Cr(acac)3 and square-planar Cu(acac)2 compounds and interpreted by different electron configurations and magnetizations at the central metal atoms resulting in different spin delocalizations and polarizations of the ligand atoms. Further, effects of substituents on the 13C NMR resonance of the ipso carbon atom reaching almost 700 ppm for Cr(acac)3 compounds are interpreted based on the analysis of Fermi-contact hyperfine contributions.

The simulation of X-ray absorption spectra requires both scalar and spin-orbit (SO) relativistic effects to be taken into account, particularly near L- and M-edges where the SO splitting of core p and d orbitals dominates. Four-component Dirac-Coulomb Hamiltonian-based linear damped response time-dependent density functional theory (4c-DR-TDDFT) calculates spectra directly for a selected frequency region while including the relativistic effects variationally, making the method well suited for X-ray applications. In this work, we show that accurate X-ray absorption spectra near L2,3- and M4,5-edges of closed-shell transition metal and actinide compounds with different central atoms, ligands, and oxidation states can be obtained by means of 4c-DR-TDDFT. While the main absorption lines do not change noticeably with the basis set and geometry, the exchange-correlation functional has a strong influence with hybrid functionals performing the best. The energy shift compared to the experiment is shown to depend linearly on the amount of Hartee-Fock exchange with the optimal value being 60% for spectral regions above 1000 eV, providing relative errors below 0.2% and 2% for edge energies and SO splittings, respectively. Finally, the methodology calibrated in this work is used to reproduce the experimental L2,3-edge X-ray absorption spectra of [RuCl2(DMSO)2(Im)2] and [WCl4(PMePh2)2], and resolve the broad bands into separated lines, allowing an interpretation based on ligand field theory and double point groups. These results support 4c-DR-TDDFT as a reliable method for calculating and analyzing X-ray absorption spectra of chemically interesting systems, advance the accuracy of state-of-the art relativistic DFT approaches, and provide a reference for benchmarking more approximate techniques.

Chemical shifts present crucial information about an NMR spectrum. They show the influence of the chemical environment on the nuclei being probed. Relativistic effects caused by the presence of an atom of a heavy element in a compound can appreciably, even drastically, alter the NMR shifts of the nearby nuclei. A fundamental understanding of such relativistic effects on NMR shifts is important in many branches of chemical and physical science. This review provides a comprehensive overview of the tools, concepts, and periodic trends pertaining to the shielding effects by a neighboring heavy atom in diamagnetic systems, with particular emphasis on the "spin-orbit heavy-atom effect on the light-atom" NMR shift (SO-HALA effect). The analyses and tools described in this review provide guidelines to help NMR spectroscopists and computational chemists estimate the ranges of the NMR shifts for an unknown compound, identify intermediates in catalytic and other processes, analyze conformational aspects and intermolecular interactions, and predict trends in series of compounds throughout the Periodic Table. The present review provides a current snapshot of this important subfield of NMR spectroscopy and a basis and framework for including future findings in the field.

With the increasing interest in compounds containing heavier elements, the experimental and theoretical community requires computationally efficient approaches capable of simultaneous non-perturbative treatment of relativistic, spin-polarization, and electron correlation effects. The ReSpect program has been designed with this goal in mind and developed to perform relativistic density functional theory (DFT) calculations on molecules and solids at the quasirelativistic two-component (X2C Hamiltonian) and fully relativistic four-component (Dirac-Coulomb Hamiltonian) level of theory, including the effects of spin polarization in open-shell systems at the Kramers-unrestricted self-consistent field level. Through efficient algorithms exploiting time-reversal symmetry, biquaternion algebra, and the locality of atom-centered Gaussian-type orbitals, a significant reduction of the methodological complexity and computational cost has been achieved. This article summarizes the essential theoretical and technical advances made in the program, supplemented by example calculations. ReSpect allows molecules with >100 atoms to be efficiently handled at the four-component level of theory on standard central processing unit-based commodity clusters, at computational costs that rarely exceed a factor of 10 when compared to the non-relativistic realm. In addition to the prediction of band structures in solids, ReSpect offers a growing list of molecular spectroscopic parameters that range from electron paramagnetic resonance parameters (g-tensor, A-tensor, and zero-field splitting), via (p)NMR chemical shifts and nuclear spin-spin couplings, to various linear response properties using either conventional or damped-response time-dependent DFT (TDDFT): excitation energies, frequency-dependent polarizabilities, and natural chiroptical properties (electronic circular dichroism and optical rotatory dispersion). In addition, relativistic real-time TDDFT electron dynamics is another unique feature of the program. Documentation, including user manuals and tutorials, is available at the program's website http://www.respectprogram.org.

An unambiguous assignment of coupling pathways plays an important role in the description and rationalization of NMR indirect spin-spin coupling constants (SSCCs). Unfortunately, the SSCC analysis and visualization tools currently available to quantum chemists are restricted to nonrelativistic theory. Here, we present the theoretical foundation for novel relativistic SSCC visualization techniques based on analysis of the SSCC densities and the first-order current densities induced by the nuclear magnetic dipole moments. Details of the implementation of these techniques in the ReSpect program package are discussed. Numerical assessments are performed on through-space SSCCs, and we choose as our examples the heavy-atom Se-Se, Se-Te, and Te-Te coupling constants in three similar molecules for which experimental data are available. SSCCs were calculated at the nonrelativistic, scalar relativistic, and four-component relativistic density functional levels of theory. Furthermore, with the aid of different visualization methods, we discuss the interpretation of the relativistic effects, which are sizable for Se-Se, very significant for Se-Te, and cannot be neglected for Te-Te couplings. A substantial improvement of the theoretical SSCC values is obtained by also considering the molecular properties of a second conformation.

We report a combined experimental-theoretical study on the 31P NMR chemical shift for a number of trans-platinum(ii) complexes. Validity and reliability of the 31P NMR chemical shift calculations are examined by comparing with the experimental data. A successful computational protocol for the accurate prediction of the 31P NMR chemical shifts was established for trans-[PtCl2(dma)PPh3] (dma = dimethylamine) complexes. The reliability of the computed values is shown to be critically dependent on the level of relativistic effects (two-component vs. four component), choice of density functionals, dynamical averaging, and solvation effects. Snapshots obtained from ab initio molecular dynamics simulations were used to identify those solvent molecules which show the largest interactions with the platinum complex, through inspection by using the non-covalent interaction program. We observe satisfactory accuracy from the full four-component matrix Dirac-Kohn-Sham method (mDKS) based on the Dirac-Coulomb Hamiltonian, in conjunction with the KT2 density functional, and dynamical averaging with explicit solvent molecules.

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.