Theoretical investigation of the evolution of the topological phase of Bi$_{2}$Se$_{3}$ under mechanical strain

The topological insulating phase results from inversion of the band gap due to spin-orbit coupling at an odd number of time-reversal symmetric points. In Bi$_2$Se$_3$, this inversion occurs at the $Gamma$ point. For bulk Bi$_2$Se$_3$, we have analyzed the effect of arbitrary strain on the $Gamma$ point band gap using Density Functional Theory. By computing the band structure both with and without spin-orbit interactions, we consider the effects of strain on the gap via Coulombic interaction and spin-orbit interaction separately. While compressive strain acts to decrease the Coulombic gap, it also increases the strength of the spin-orbit interaction, increasing the inverted gap. Comparison with Bi$_2$Te$_3$ supports the conclusion that effects on both Coulombic and spin-orbit interactions are critical to understanding the behavior of topological insulators under strain, and we propose that the topological insulating phase can be effectively manipulated by inducing strain through chemical substitution.

Optimized norm-conserving Hartree-Fock pseudopotentials for plane-wave calculations

We report Hartree-Fock (HF) based pseudopotentials suitable for plane-wave calculations. Unlike typical effective core potentials, the present pseudopotentials are finite at the origin and exhibit rapid convergence in a plane-wave basis; the optimized pseudopotential method [A. M. Rappe et. al, Phys. Rev. B 41 1227--30 (1990)] improves plane-wave convergence. Norm-conserving HF pseudopotentials are found to develop long-range non-Coulombic behavior which does not decay faster than 1/r, and is non-local. This behavior, which stems from the nonlocality of the exchange potential, is remedied using a recently developed self-consistent procedure [J. R. Trail and R. J. Needs, J. Chem. Phys. 122, 014112 (2005)]. The resulting pseudopotentials slightly violate the norm conservation of the core charge. We calculated several atomic properties using these pseudopotentials, and the results are in good agreement with all-electron HF values. The dissociation energies, equilibrium bond lengths, and frequency of vibrations of several dimers obtained with these HF pseudopotentials and plane waves are also in good agreement with all-electron results.