Search

A previously proposed non-canonical coupled-perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment for spin-orbit coupling is here generalized to infinite periodic systems. The scalar-relativistic periodic KS-DFT/HF solution, obtained with a relativistic effective core potential, is taken as the zeroth-order approximation. Explicit expressions are given for the total energy through 3rd-order, which satisfy the 2N + 1 rule (i.e. requiring only the 1st-order perturbed wave function for determining the energy through 3rd-order). Expressions for additional 2nd-order corrections to the perturbed wave function (as well as related one-electron properties) are worked out at the uncoupled-perturbed level of theory. The approach is implemented in the \textsc{Crystal} program and validated with calculations of the total energy, electronic band structure, and density variables of spin-current DFT on the tungsten dichalcogenide hexagonal bilayer series (i.e. WSe$_2$, WTe$_2$, WPo$_2$, WLv$_2$), including 6p and 7p elements as a stress test. The computed properties through second- or third-order match well with those from reference two-component self-consistent field (2c-SCF) calculations. For total energies, $E^{(3)}$ was found to consistently improve the agreement against the 2c-SCF reference values. For electronic band structures, visible differences w.r.t. 2c-SCF remained through second-order in only the single-most difficult case of WLv$_2$. As for density variables of spin-current DFT, the perturbed electron density, being vanishing in first-order, is the most challenging for the perturbation theory approach. The visible differences in the electron densities are, however, largest close to the core region of atoms and smaller in the valence region. Perturbed spin-current densities, on the other hand, are well reproduced in all tested cases.

An analysis of orbital magnetization in band insulators is provided. It is shown that a previously proposed electronic orbital angular-momentum operator generalizes the ``modern theory of orbital magnetization'' to include non-local Hamiltonians. Expressions for magnetic transition dipole moments needed for the calculation of optical rotation (OR) and other properties are developed. A variety of issues that arise in this context are critically analyzed. These issues include periodicity of the operators, previously proposed band dispersion terms as well as, if and where needed, evaluation of reciprocal space derivatives of orbital coefficients. Our treatment is used to determine the optical rotatory power of band insulators employing a formulation that accounts for electric dipole - electric quadrupole (DQ), as well as electric dipole-magnetic dipole, contributions. An implementation in the public \textsc{Crystal} program is validated against a model finite system and applied to the $\alpha$-quartz mineral through linear-response time-dependent density functional theory with a hybrid functional. The latter calculations confirmed the importance of DQ terms. Agreement against experiment was only possible with i) use of a high quality basis set, ii) inclusion of a fraction of non-local Fock exchange, and iii) account of orbital-relaxation terms in the calculation of response functions.