We quantify the accuracy of different non-self-consistent and self-consistent spin-orbit coupling (SOC) treatments in Kohn-Sham and hybrid density-functional theory by providing a band structure benchmark set for the valence and low-lying conduction energy bands of 103 inorganic compounds, covering chemical elements up to Po. Reference energy band structures for the PBE density functional are obtained using the full-potential (linearized) augmented plane wave code Wien2k, employing its self-consistent treatment of SOC including Dirac-like p1/2 orbitals in the basis set. We use this benchmark set to benchmark a computationally simpler, non-self-consistent all-electron treatment of SOC based on scalar-relativistic orbitals and numeric atom-centered orbital basis functions. For elements up to Z≈50, both treatments agree virtually exactly. For the heaviest elements considered (Tl, Pb, Bi, Po), the band structure changes due to SOC are captured with a relative deviation of 11% or less. For different density functionals (PBE vs. the hybrid HSE06), we show that the effect of spin-orbit coupling is usually similar but can be dissimilar if the qualitative features of the predicted underlying scalar-relativistic band structures do not agree. All band structures considered in this work are available online via the NOMAD Repository to aid in future benchmark studies and methods development.
@article{arxiv.1705.01804,
title = {103-Compound Band Structure Benchmark of Post-SCF Spin-Orbit Coupling Treatments in Density-Functional Theory},
author = {William P. Huhn and Volker Blum},
journal= {arXiv preprint arXiv:1705.01804},
year = {2017}
}
Comments
20 pages, 9 figures, 5 tables, supplemental material included in ancillary files