Perturbation Theory Treatment of Spin-Orbit Coupling, Part I: Double Perturbation Theory Based on a Single-Reference Initial Approximation

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

Perturbation Theory Treatment of Spin-Orbit Coupling. III: Coupled Perturbed Method for Solids

A previously proposed noncanonical 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 third-order, which satisfy the 2N + 1 rule (i.e., requiring only the first-order perturbed wave function for determining the energy through third-order). Expressions for additional second-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 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₂, WTe₂, WPo₂, WLv₂), 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⁽³⁾ 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₂. 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.