Role of spin currents on electron-electron interaction in the quantum spin Hall phase

The spin current density functional theory (SCDFT) is the generalization of the standard DFT to treat a fermionic system embedded in the effective external field produced by the spin-orbit coupling interaction. Even in the absence of a spin polarization, the SCDFT requires the electron-electron potential to depend on the spin currents ${\mathbf{J}}^{x}, {\mathbf{J}}^{y}$, and ${\mathbf{J}}^{z}$, which only recently was made possible for practical relativistic quantum-mechanical simulations, through the inclusion of a fraction of nonlocal Fock exchange into local-density or generalized-gradient density functional approximations [J. K. Desmarais, J.-P. Flament, and A. Erba, Phys. Rev. B 102, 235118 (2020)]. Here, we apply the SCDFT to the quantum spin Hall phase and show how it improves (even qualitatively) the description of its electronic features relative to the DFT. We study the Bi (001) 2D bilayer and its band insulator to topological insulator phase transition (via $s+{p}_{z}\ensuremath{\leftrightarrow}{p}_{x}+i{p}_{y}$ band inversion) as a function of mechanical strain. We show that the explicit account of spin currents in the electron-electron potential of the SCDFT is key to the appearance of a Dirac cone at the $\mathrm{\ensuremath{\Gamma}}$ point in the valence band structure at the onset of the topological phase transition. Finally, the valence band structure of this system is rationalized using a simple first-order $\mathbf{k}\ifmmode\cdot\else\textperiodcentered\fi{}\mathbf{p}$ quasidegenerate perturbation theory model.