Electronic phase transitions of bismuth under strain from relativistic self-consistent GW calculations

Electronic phase transitions of bismuth under strain from relativistic self-consistent Gw calculations

Irene Aguilera, Christoph Friedrich, Stefan Blügel Peter Grünberg Institute and Institute for Advanced Simulation, Forschungszentrum Jülich and JARA, 52425 Jülich, Germany
July 20, 2019

We present quasiparticle self-consistent (QS) calculations of semimetallic bulk Bi. We go beyond the conventional QS method by including the spin-orbit coupling throughout the self-consistency cycle. This approach improves the description of the electron and the hole pockets considerably with respect to standard density functional theory (DFT), leading to excellent agreement with experiment. We employ this relativistic QS approach to conduct a study of the semimetal-to-semiconductor and the trivial-to-topological transitions that Bi experiences under strain. DFT predicts that an unphysically large strain is needed for such transitions. We show, by means of the relativistic QS description of the electronic structure, that an in-plane tensile strain of only 0.3% and a compressive strain of 0.4% are sufficient to cause the semimetal-to-semiconductor and the trivial-to-topological phase transitions, respectively. Thus, the required strain moves into a regime that is likely to be realizable in experiment, which opens up the possibility to explore bulklike topological behavior of pure Bi.


Bismuth exhibits a series of peculiarities that has made it the subject of experimental and theoretical interest for decades. The very low density of carriers with high carrier mobility and a very long mean free path make bismuth very interesting for electronic transport studies. With the recent introduction of the classification of solids according to topology classes of the electronic structure,Moore (2010) Bi has again moved into the limelight of current discussions. With a lattice structure that bears similarities to graphene and a very large spin-orbit interaction, it has all the potential to be a topological semimetal or semiconductor. From the theoretical point of view, the large spin-orbit coupling (SOC), the semimetallic character, and the very small local direct band gap at the L point of the Brillouin zone make ab initio calculations of the electronic structure of Bi very challenging.

As for its topological classification, some important aspects are not yet fully understood. It is widely accepted that bulk Bi is topologically trivial in contrast to Sb. This is in agreement with angle-resolved photoemission spectroscopy (ARPES) measurements of BiAst and Höchst (2003); Koroteev et al. (2004); Hofmann (2006) and BiSb alloys,Hsieh et al. (2008); Roushan et al. (2009); Hsieh et al. (2009) as well as transport studies, in which a semimetal-to-semiconductor (SMSC)Tichovolsky and Mavroides (1969) and a topological transitionRed’ko and Rodionov (1985) are observed as the Sb concentration is varied from the topologically trivial Bi to the topological semimetal Sb, resulting in three-dimensional (3D) topological insulators for some of the alloys. However, controversies concerning the topological properties of bismuth persist. For example, it has recently been claimedOhtsubo et al. (2013) on the basis of ARPES measurements that bulk bismuth is topologically nontrivial, a claim that is in conflict with most previous findings including ab initio calculations based on density functional theory (DFT). The authors of Ref. Ohtsubo et al., 2013 attributed these discrepancies to the systematic underestimation of the band gaps obtained in the local-density (LDA) or generalized gradient (GGA) approximations of DFT. Given the very small experimental value (11-15 meVMaltz and Dresselhaus (1970); Vecchi and Dresselhaus (1974); Isaacson and Williams (1969); Brown et al. (1963); Smith et al. (1964)) for the direct band gap at the L point (), whose “sign” (order of the states) controls the topological nature of Bi, it is very conceivable that DFT might incorrectly predict the sign of the band gap, and not just its magnitude. Our parity analysis within LDA is in agreement with results in the literature (see, e.g., Refs. Golin, 1968; Timrov et al., 2012; Teo et al., 2008) and confirms the trivial character of Bi. But surprisingly, LDA presents an unusual overestimation of the band gap at L (86 meV; see Table 1), which is an effect usually characteristic of inverted gaps, as with those in topological insulators.Yazyev et al. (2012); Aguilera et al. (2013a) The gap at L is indeed inverted, as can be seen by varying the SOC strength from zero to 100%: the two states at L exchange order for a certain SOC strength. But the gap at T is inverted as well. There is, hence, an even number of inverted gaps at the time-reversal invariant momenta, giving rise to a topologically trivial value of the invariant.Kane and Mele (2005); Fu et al. (2007) The GGA approximation hardly improves the LDA results (Table 1). The poor description of the band gaps casts doubt on the LDA and GGA studies of the trivial-to-topological (TT)Hirahara et al. (2012); Liu et al. (2011a) or the SMSCGutierrez et al. (2006); Hirahara et al. (2012) transitions of Bi, as well as the TT transition in BiSb alloys, as these studies critically depend on the sign and the value of the band gap at L.

The cause of the band-gap problem of DFT is well known. Strictly, the single-particle energies cannot be interpreted as quasiparticle energies of the interacting electron system. The approximationHedin (1965) for the electronic self-energy, applied as a one-shot correction to DFT or, more rarely, in the quasiparticle self-consistent (QS) approach,Faleev et al. (2004) remedies the aforementioned problem and is often used to correct the band gaps obtained from LDA or GGA. In most cases it increases the band gap leading to a much better agreement with experiments.van Schilfgaarde et al. (2006) However, we show here (cf. Table 1) that in the case of bismuth, the quasiparticle correction is negative, i.e., it correctly reduces the band gap to 32 meV in one-shot (which corresponds to the first iteration of QS) and to 13 meV within the self-consistent QS scheme. The latter value lies in the range of experimental values. For these calculations, we have devised and implemented a method that combines the SOC with the many-body renormalization in a consistent way (see below). The sign of the band gap is not changed by the quasiparticle correction, i.e., bismuth remains trivial, in accordance with the vast majority of experimental and theoretical studies. However, the phase boundary at which bismuth becomes topologically nontrivial is considerably closer in QS than in LDA. We will show that only a relatively modest in-plane tensile strain of 0.3% and compressive strain of 0.4% are sufficient to cause the SMSC and TT phase transition, respectively, which opens up the possibility to explore the 3D topological properties of pure Bi.

direct gap at L () 186 74 32 13 11 – 1511111.0 [Maltz and Dresselhaus, 1970], 13.6 [Vecchi and Dresselhaus, 1974], 15.0 [Isaacson and Williams, 1969; Brown et al., 1963], 15.3 [Smith et al., 1964]
indirect gap T-L () 105 72 66 33 32 – 39222-32.0 [Edel’man, 1975], -36.0 [Isaacson and Williams, 1969], -38.0 [Dinger and Lawson, 1973], -38.2 [Isaacson and Williams, 1969],-38.5 [Smith et al., 1964]

Table 1: Values of the direct band gap at L () and the indirect band gap () [see Fig. 2(b)] for bulk bismuth (in meV) calculated within LDA, GGA, , and QS, compared with experimental results.

Bulk bismuth crystallizes in the A7 rhombohedral structure with Rm space group and two atoms per unit cell. The cell and the atomic positions are determined by the internal parameter , the rhombohedral lattice parameter , and the rhombohedral angle [see Fig. 1(a)]. The layered structure consists of Bi bilayers separated by a van der Waals gap. For the undistorted reference lattice, we employ experimental parameters from Ref. Schiferl and Barrett, 1969: =4.7458 Å, =57.230, =0.23389. Equivalently, the structure can be described as hexagonal with in-plane and out-of-plane lattice constants =4.5460 Å and =11.862 Å (=2.6093). These lattice parameters were measured at 298 K. The gray vertical areas in Fig. 3 show the range of experimental values of for temperatures between 4 and 298 K.Schiferl and Barrett (1969) All of our calculations are carried out with the DFT code fleurfle () and the code spex,Friedrich et al. (2010) realizations in the all-electron full-potential linearized augmented-plane-wave (FLAPW) formalism. The convergence parameters of the DFT, , and QS calculations are given in the Appendix. To include relativistic effects in the DFT calculations, we employ the scalar-relativisticKoelling and Harmon (1977) approximation inside the muffin-tin spheres, and the SOC is included in a self-consistent manner.Li et al. (1990)

Figure 1: (Color online) (a) Rhombohedral (green) unit cell and hexagonal crystal structure of bulk Bi. The blue and purple atoms distinguish the atoms in the two layers of each bilayer. (b) Bulk Brillouin zone of Bi expressed in the rhombohedral unit cell.

Most of the one-shot calculations that include SOC and all of the QS ones published so far, employ the SOC as an a posteriori correction after the quasiparticle correction to DFT has been carried out without SOC. This is, of course, an approximation. For example, many-body renormalization effects of spin-orbit split bands are then not taken into account (see Ref. Aguilera et al., 2013b for a detailed discussion). In our relativistic calculations, the SOC is already incorporated in the noninteracting system that serves as the starting point for the quasiparticle calculation.Aguilera et al. (2013b); Sakuma et al. (2011) The one-shot self-energy thus acquires terms that couple the two spin channels, enabling a many-body renormalization of the SOC itself. In this work, we propose to extend this principle to the self-consistent QS calculations, so that the self-energy contains spin off-diagonal blocks that it inherits from the SOC throughout the whole self-consistent process. (In analogy to the notation of Ref. Aguilera et al., 2013b, one could denote this approach by QS, but for simplicity we will simply write and QS in the following.) Our results show that this treatment is crucial for Bi. The a posteriori correction mentioned above, on the other hand, leads to unphysical results for Bi (see Appendix), which then appears as a topological insulator with large band gaps (83 and 259 meV in one-shot and QS, respectively) instead of a trivial semimetal.

In addition to the value of the band gap at the L point, the properties of Bi are also determined by the states at the symmetry point T. In particular, the semimetallic character of bismuth is caused by the overlap between the highest valence band at T and the lowest conduction band at L [see Fig. 2(b)]. This creates three very small electron pockets (at the three L points in the Brillouin zone) and one hole pocket at T. Several attempts to describe these pockets theoretically have been published in recent years. For a review and comparison of them see, e.g., Table II in Ref. Timrov et al., 2012. The approaches used include LDA and GGA,Timrov et al. (2012); Gonze et al. (1988); Shick et al. (1999) tight-binding models,Liu and Allen (1995) and empirical pseudopotential calculations,Golin (1968) but no results have been published so far. The results from the models and DFT calculations have not been satisfactory: the absolute value of the overlap, or in other words the indirect band gap , came out too large in general compared to experiment.

Figure 2: (Color online) LDA, , and QS band structures of bulk Bi. The inset in (b) shows the electron and hole pockets together with the and gaps corresponding to the LDA bands.

Our values of are shown in Table 1. Both the LDA and GGA overestimate indeed this indirect band gap, in line with previous publications. While the one-shot approach does correct in the right direction, it uses LDA as a starting point and, thus, inherits some of this overestimation. The self-consistent procedure of QS, on the other hand, makes the result independent of the starting point and moves inside the range of experimental values.

Figure 2 shows the band structure of bulk bismuth obtained with LDA, , and QS along the path connecting the high-symmetry k points represented in Fig. 1(b). Whereas the overall shape of the band structure remains basically unaffected by the many-body corrections, there are important changes in the details. In particular, while in most parts the valence and conduction bands are shifted apart by the many-body corrections as commonly seen in calculations, the changes of the dispersion in the vicinity of the Fermi energy are such that the absolute values of and get smaller and reach values of and  meV, respectively. As has already been discussed before,Sakuma et al. (2011); Yazyev et al. (2012); Aguilera et al. (2013a, b) this is the result of a delicate interplay between the SOC and the many-body renormalization described by the approximation. We conclude that the present relativistic QS is the method of choice to study the band structure of Bi and its topological character as well as the TT and SMSC transitions that Bi experiences under strain.

A TT transition was indeed recently observed experimentally for strained Bi(111) ultrathin films,Hirahara et al. (2012) and a similar transition was found for strained Bi nanowires elongated about 2%.Condrea et al. (2013) In the case of thin films, a strain can be introduced naturally by a lattice mismatch between film and substrate. To simulate the strain in our calculations, we apply a volume-conserving distortion of the structure by varying the ratio, and keeping the internal parameter constant. A compression in the direction is accompanied by a corresponding increase of , and vice versa. This is equivalent to varying the rhombohedral angle , see Fig. 1(a). As the ratio increases, the nearest-neighbor distance reduces and the van der Waals gap expands. Both effects lead to a strengthening of the in-plane interactions, while the interbilayer coupling weakens. The altered crystal field reduces the band gap at L until it reaches a zero value at a critical ratio and becomes negative. This undoes the band inversion and causes a sign change in the calculation of the topological invariants,Kane and Mele (2005); Fu et al. (2007) which, for systems with inversion symmetry, derive from the parities of the occupied states at the time-reversal invariant momenta. The parity of a wave function remains unchanged upon applying the quasiparticle correction. The topological invariants can, thus, only change if there is an interchange of valence and conduction states with respect to LDA. A topological-to-trivial transition in the topological insulator BiSe as a result of strain has also been discussed in the literature. It has been arguedYoung et al. (2011); Liu et al. (2011b) that tensile out-of-plane strain reduces the spin-orbit strength leading to a trivial phase. We do not find such an effect for bulk Bi: no substantial increase or decrease in the spin-orbit strength is observed upon application of strain.

Figure 3: (Color online) (a) LDA, GGA, , and QS direct band gap of bismuth at L for different strains (varying the ratio and keeping the volume constant). (b) Three electronic phases are found for Bi as a function of the ratio within all the approaches. (c) LDA, GGA, , and QS indirect band gap of bismuth for different values of . On the upper axis of (a) and (c), we show the percentage of in-plane strain. The range of experimental values of ,Schiferl and Barrett (1969) the experimental gap at L,Maltz and Dresselhaus (1970); Vecchi and Dresselhaus (1974); Isaacson and Williams (1969); Brown et al. (1963); Smith et al. (1964) as well as the experimental indirect gapEdel’man (1975); Dinger and Lawson (1973); Isaacson and Williams (1969); Smith et al. (1964) are represented as gray shaded areas in (a) and (c). Small differences to the QS values in Table 1 are due to the use of a reduced k-point set; see Appendix.

Figure 3(a) shows the direct band gap as a function of the ratio with the range of equilibrium experimental values of bulk Bi represented as gray shaded areas. Positive values of the gap denote a trivial phase, whereas negative values correspond to a topologically nontrivial case. In addition to the semimetallic state (solid symbols) in Fig. 3(a), we also observe a transition to a semiconducting state (open symbols) for tensile in-plane strain. This SMSC transition can also be seen in Fig. 3(c) showing the variation of the indirect band gap, which changes sign and becomes positive at the point of transition.

At the SMSC transition, the and QS curves exhibit abnormal behavior; there is a jump (dashed lines) where we expect the curve to be continuous. An analogous steplike behavior can be seen in the QS curve at the TT transition. As shown in the Appendix, these anomalies are an artifact of the k-point convergence: The steps get systematically smaller when one improves the k-point sampling, and convergence is achieved much more readily in the trivial semimetallic phase than in the other two phases, so that we can take the points between the transitions in Fig. 3 as converged, and the fitted curves (solid lines) show the expected behavior in the sampling limit. We have employed a quadratic function for the direct gap and a linear function for the indirect gap as a function of , which for QS give and in meV.

(Å) (Å) (meV)
(a) TT:
LDA 4.461(1.9%) 12.315 2.761(+5.8%) 15.8
GGA 4.475(1.6%) 12.238 2.735(+4.8%) 11.7
4.513(0.7%) 12.036 2.667(+2.2%) 4.1
QS 4.527(0.4%) 11.960 2.642(+1.2%) 1.9
(b) SMSC:
LDA 4.579(+0.7%) 11.690 2.553(2.2%) 0.5
GGA 4.580(+0.7%) 11.685 2.551(2.2%) 0.6
4.569(+0.5%) 11.744 2.570(1.5%) 0.2
QS 4.562(+0.3%) 11.776 2.581(1.1%) 0.1

Table 2: Values of the lattice parameters and and the ratio for which (a) the TT and (b) the SMSC transitions occur for the different theoretical approaches. The values in brackets indicate the compressive () or tensile (+) strain in percentage. The difference in the total energies () between the structure for which the transition occurs and the unperturbed reference lattice has been calculated within LDA. Note that we used =4.5460 Å, =11.862 Å, and =2.6093Schiferl and Barrett (1969) for the undistorted lattice.

Figure 3(b) shows the three electronic phases that we can distinguish in bulk Bi under strain: a trivial semiconducting phase, a trivial semimetallic phase, and a topological semimetallic phase. The overlap between the electron and hole pockets, or the indirect gap, is caused by the distortion of the cell with respect to a cubic structure. It seems, therefore, that bismuth cannot be in a topological insulating phase because the stress necessary to make it topological distorts the cell rhombohedrally, distancing it even further from the cubic symmetry so that Bi becomes more and more metallic. Figure 3 demonstrates that, whereas the two transitions are qualitatively described correctly by the four approaches, the LDA and GGA predict a large critical strain that would be very difficult to achieve experimentally. In particular, the TT transition is predicted to occur at a critical in-plane compressive strain of 1.9% and 1.6%, respectively (see Table 2). On the contrary, the strain needed to cause the transition is significantly smaller in (0.7%) and, in particular, in QS (0.4%). In other words, taking into account structural relaxation, significantly thicker samples could be grown under the relatively modest critical strain that is predicted by QS than under that predicted by DFT.

In conclusion, we propose a QS scheme in which the SOC is included throughout the whole self-consistent process. Our results show that such relativistic QS calculations are the right ab initio approach to determine the electronic structure and the topological classification of bismuth. They predict its two most significant and delicate electronic properties in quantitative agreement with experiments: The direct gap at L, which determines whether Bi is topological or trivial, and the indirect gap between L and T, which describes the overlap between the electron and hole pockets and thus determines the semimetallic character of Bi. In response to a recent controversy,Ohtsubo et al. (2013) we conclude the following from our calculations: bulk Bi is a trivial semimetal.

We showed that Bi can undergo electronic phase transitions under rather small lattice strain: a trivial SMSC transition under 0.3% in-plane tensile strain and a TT semimetal transition under 0.4% compressive strain. In contrast to LDA and GGA, the critical strains predicted by QS are attainable by experiments. The presented results motivate additional experimental efforts to prepare Bi as a topological insulator, by opening band gaps in slightly strained Bi. For example, the small strain needed to cause the TT transition concluded by this work, together with the SMSC transition driven by quantum-size effects in relatively thick films,Xiao et al. (2012) reveals the potential to observe a topological insulating behavior in bulk-like films of pure Bi. In fact, first experimental evidences of 3D-like topological thin films of Bi are currently under discussion.Zhu et al. (2014)

The present results demonstrate that LDA and GGA lack the quantitative accuracy to predict the critical strain for which the electronic phase transitions occur. We speculate that standard DFT also yields wrong predictions for the critical concentration at which Bi experiences a TT transition upon alloying with other elements (such as in BiSb alloysTichovolsky and Mavroides (1969); Red’ko and Rodionov (1985); Teo et al. (2008); Hsieh et al. (2008, 2009)) or for the critical thickness at which a SMSC transition caused by quantum-size effects takes place. Bi is one example, but we believe it is much more general. Other examples could include the very large strains necessary to induce topological phase transitions in BiSe Young et al. (2011); Liu et al. (2011b) and in TlBiS and TlSbS.Zhang et al. (2015) It would thus be desirable to reinvestigate the critical points for these transitions within the relativistic QS method proposed in this work.

We thank Xiaofeng Jin for the very interesting discussions that motivated this work, Gregor Mussler for his comments on the experimental growth of Bi, and Gustav Bihlmayer for valuable discussions and a critical reading of the manuscript. This work was supported by the Alexander von Humboldt Foundation through a postdoctoral fellowship, and by the Helmholtz Association through the Virtual Institute for Topological Insulators (VITI).


The calculations are carried out with the DFT code fleurfle () and the code spex,Friedrich et al. (2010) which are based on the all-electron FLAPW formalism. For the valence electrons, space is partitioned into spherical regions around the atoms (muffin-tin spheres) and interstitial region between the spheres. We use an angular momentum cutoff in the muffin-tin spheres and a plane-wave cutoff of 4.5 bohr in the interstitial region. For the DFT calculations, we employ either the Perdew-ZungerPerdew and Zunger (1981) parametrization of the LDA exchange-correlation functional or the Perdew-Burke-ErnzerhofPerdew et al. (1997) parametrization of the GGA. The one-shot calculations are always performed using the LDA mean-field system as a starting point.

Due to the very small band gap of bulk Bi (11 to 15 meV for the direct oneMaltz and Dresselhaus (1970); Vecchi and Dresselhaus (1974); Isaacson and Williams (1969); Smith et al. (1964) and 32 to 39 for the indirect oneEdel’man (1975); Dinger and Lawson (1973); Isaacson and Williams (1969); Smith et al. (1964)), all computational parameters entering the and QS calculations are chosen after a systematic study in order to obtain thoroughly converged results. An angular momentum cutoff of and a linear momentum cutoff of 2.9 bohr are employed to construct the mixed product basis,Kotani and van Schilfgaarde (2002); Friedrich et al. (2010) used to represent the dielectric matrix and the screened interaction.

The one-shot and self-consistent QS results for the undistorted structure are obtained with a 666 k-point mesh sampling the Brillouin zone (Table 1 and Fig. 2). Fortunately, we found that the QS results converge faster with respect to the number of k points than the results of one-shot so that we could afford to reduce the k-point mesh to 444 for the QS calculations of bismuth under strain. For more details about the k-point convergence, see below.

To compute the Green function and the polarization function, 500 bands are used. This corresponds to approximately 130 eV above the Fermi energy. The complete fifth shell of Bi is treated as valence states by the use of local orbitals. In this way, the contribution of the 5s5p5d states to the Green function and the electronic screening is fully taken into account. Despite their low energetic position, inclusion of the 5s and 5p states effects a change of 5 meV in the band gap at L. To describe high-lying states accurately and to avoid linearization errors,Friedrich et al. (2006, 2011); Michalicek et al. (2013) we complement the basis for the valence electrons for each atom by two local orbitals per angular momentum up to with energy parameters far up in the unoccupied states. For the interpolation of the band structures in Fig. 2, maximally localized Wannier functions obtained by the wannier90 libraryMostofi et al. (2008) and a 666 k-point mesh were employed.

We employ the Dirac equation for the core electrons (up to and including the fourth shell) so that relativistic effects are fully accounted for. The core states are thus represented by four-component spinors. This representation is retained for the evaluation of the core-valence contribution to the exchange self-energy. (Using averaged core states that are represented in terms of non-relativistic spherical harmonics instead of the spinors introduces an error of 4 and 7 meV in the direct band calculated with one-shot and QS, respectively.)

As a semimetal, bismuth exhibits metallic screening, which is described technically by the so-called Drude term in the screened interaction . This term stems from virtual intraband transitions across the Fermi surface. It can be formulated in a functional form being proportional to the square of the plasma frequency, which in turn is evaluated by an integration over the Fermi surface. The Drude term can be treated analyticallyFriedrich et al. (2010) and, as long as the Fermi surface is sufficiently big, it normally does not pose any numerical problem. However, bismuth has a very small Fermi surface due to the tiny electron and hole pockets, which eventually leads to a very sharp “Drude peak” in the self-energy, impeding a straight-forward numerical solution of the non-linear quasiparticle equation. One could also say that, while the Drude term is actually treated correctly, it gets too much weight because of the finite k-point mesh. Therefore, we neglect the Drude term in our calculations and, instead, simply scale the head element of in the limit to enforce metallic screening. This hardly changes , and we have found that it leads to a very favorable k-point convergence; see below. (In the limit of dense k-point sets, the treatment of the Drude term, which affects only a single k point, i.e., the point, becomes immaterial.)

Figure 4: (Color online) Convergence of the QS direct band gap at the L point: without SOC (circles), with SOC added a posteriori (triangles), and with the relativistic approach used in the main text (squares). Iteration 0, 1, and 10, correspond to the LDA, one-shot , and QS values, respectively.

In Fig. 4, we show the direct band gap at the L point for each iteration in the QS self-consistency procedure, in particular, we distinguish between calculations (i) without SOC, (ii) with spin-orbit coupling added a posteriori to a result obtained without SOC (QS+SOC), which is a correction scheme used in most calculations so far, and (iii) the relativistic approach used in the present work in which spin-orbit coupling is included in the calculation of both and throughout the whole self-consistent procedure. (In an earlier publication,Aguilera et al. (2013b) we have used the notation for the corresponding one-shot calculation. By analogy, one could use the notation QS for the present self-consistent approach.) The a posteriori SOC correction leads to unphysical results for Bi, which appears as a topological insulator instead of a trivial semimetal. We attribute the lack of published results of bulk Bi to the crucial and delicate treatment of SOC, but also to the difficulties encountered with the semimetallicity of the system exhibiting very tiny electron and hole pockets (leading to a very sharp Drude peak in the self-energy, as discussed above) and to the smallness of the gaps requiring the calculations to be converged to within a few meV.

We have also performed relativistic QS calculations, i.e., we keep the at the LDA level while is updated self-consistently. The converged QS result for the band gap at L differs only by 0.3 meV from the QS one. This indicates that the screening obtained with LDA is a good approximation and that self-consistency induces changes mainly in . Finally, to quantify the effects of dynamics, we have also performed self-consistent calculations within the static Coulomb-hole screened-exchange (COHSEX)Hedin (1965) approximation to the self-energy (including SOC). In this approximation, the changes are more pronounced. With respect to QS, the direct gap at L is reduced by 12 meV, nearly making the system topological, and the indirect band gap shrinks to -4 meV, nearly making the system semiconducting.

To explain the discontinuities in Figs. 3(a) and 3(c) and to prove their origin as a k-point convergence issue, we have performed and QS calculations with increasing number of k points (Fig. 5). To allow for these calculations, we had to reduce some of the parameters in the calculation such as the plane-wave cutoff in the interstitial region for the LDA calculations (now 4.0 bohr), and both the angular momentum cutoff () and the linear momentum cutoff (2.6 bohr) to construct the mixed product basis for the calculations. We also reduced the number of bands for the calculations to 100 and only 5d local orbitals were used in the calculations. This obviously reduces the accuracy of the band gaps obtained—and therefore they differ from those of the main text— but it is not our intention in this appendix to analyze the quantitative results, only the qualitative behavior due to the k-point convergence.

Figure 5: (Color online) LDA, , and QS values of the direct and indirect band gaps of bulk Bi as a function of the ratio for different k-point sets, but for a set of reduced convergence parameters, see text.

From Fig. 5, we can conclude the following about the k-point convergence of one-shot and self-consistent results: (1) The calculations of the trivial semimetallic phases can be considered converged with a 666 k-point set for the one-shot calculations and 444 for the QS ones. (2) The semiconducting and topological phases require much larger k-point meshes, but (3) the corresponding results converge systematically towards an extrapolated quadratic (linear) fit to the converged trivial semimetallic phases for the direct (indirect) band gap. We note that in the case of one-shot calculations, the discontinuity in the curves due to the k-point convergence of the topological phases does not appear at the TT transition of but at that of LDA, because of the use of the LDA mean-field system as the starting point for the one-shot quasiparticle correction.


  • Moore (2010) J. E. Moore, Nature (London) 464, 194 (2010).
  • Ast and Höchst (2003) C. R. Ast and H. Höchst, Phys. Rev. B 67, 113102 (2003).
  • Koroteev et al. (2004) Y. M. Koroteev, G. Bihlmayer, J. E. Gayone, E. V. Chulkov, S. Blügel, P. M. Echenique,  and P. Hofmann, Phys. Rev. Lett. 93, 046403 (2004).
  • Hofmann (2006) P. Hofmann, Prog. Surf. Sci. 81, 191 (2006).
  • Hsieh et al. (2008) D. Hsieh, D. Qian, L.Wray, Y. Xia, Y. S. Hor, R. J. Cava,  and M. Z. Hasan, Nature (London) 452, 970 (2008).
  • Roushan et al. (2009) P. Roushan, J. Seo, C. V. Parker, Y. S. Hor, D. Hsieh, D. Qian, A. Richardella, M. Z. Hasan, R. J. Cava,  and A. Yazdani, Nature (London) 460, 1106 (2009).
  • Hsieh et al. (2009) D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, J. H. Dil, J. Osterwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. S. Hor, R. J. Cava,  and M. Z. Hasan, Science 323, 919 (2009).
  • Tichovolsky and Mavroides (1969) E. Tichovolsky and J. Mavroides, Solid State Commun. 7, 927 (1969).
  • Red’ko and Rodionov (1985) N. A. Red’ko and N. A. Rodionov, Pis’ma Zh. Eksp. Teor. Fiz. 42, 246 (1985).
  • Ohtsubo et al. (2013) Y. Ohtsubo, L. Perfetti, M. O. Goerbig, P. L. Fèvre, F. Bertran,  and A. Taleb-Ibrahimi, New J. Phys. 15, 033041 (2013).
  • Maltz and Dresselhaus (1970) M. Maltz and M. S. Dresselhaus, Phys. Rev. B 2, 2877 (1970).
  • Vecchi and Dresselhaus (1974) M. P. Vecchi and M. S. Dresselhaus, Phys. Rev. B 10, 771 (1974).
  • Isaacson and Williams (1969) R. T. Isaacson and G. A. Williams, Phys. Rev. 185, 682 (1969).
  • Brown et al. (1963) R. N. Brown, J. G. Mavroides,  and B. Lax, Phys. Rev. 129, 2055 (1963).
  • Smith et al. (1964) G. E. Smith, G. A. Baraff,  and J. M. Rowell, Phys. Rev. 135, A1118 (1964).
  • Golin (1968) S. Golin, Phys. Rev. 166, 643 (1968).
  • Timrov et al. (2012) I. Timrov, T. Kampfrath, J. Faure, N. Vast, C. R. Ast, C. Frischkorn, M. Wolf, P. Gava,  and L. Perfetti, Phys. Rev. B 85, 155139 (2012).
  • Teo et al. (2008) J. C. Y. Teo, L. Fu,  and C. L. Kane, Phys. Rev. B 78, 045426 (2008).
  • Yazyev et al. (2012) O. V. Yazyev, E. Kioupakis, J. E. Moore,  and S. G. Louie, Phys. Rev. B 85, 161101(R) (2012).
  • Aguilera et al. (2013a) I. Aguilera, C. Friedrich, G. Bihlmayer,  and S. Blügel, Phys. Rev. B 88, 045206 (2013a).
  • Kane and Mele (2005) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005).
  • Fu et al. (2007) L. Fu, C. L. Kane,  and E. J. Mele, Phys. Rev. Lett. 98, 106803 (2007).
  • Hirahara et al. (2012) T. Hirahara, N. Fukui, T. Shirasawa, M. Yamada, M. Aitani, H. Miyazaki, M. Matsunami, S. Kimura, T. Takahashi, S. Hasegawa,  and K. Kobayashi, Phys. Rev. Lett. 109, 227401 (2012).
  • Liu et al. (2011a) Z. Liu, C.-X. Liu, Y.-S. Wu, W.-H. Duan, F. Liu,  and J. Wu, Phys. Rev. Lett. 107, 136805 (2011a).
  • Gutierrez et al. (2006) G. Gutierrez, E. Menendez-Proupin,  and A. K. Singh, J. Appl. Phys. 99, 103504 (2006).
  • Hedin (1965) L. Hedin, Phys. Rev. 139, A796 (1965).
  • Faleev et al. (2004) S. V. Faleev, M. van Schilfgaarde,  and T. Kotani, Phys. Rev. Lett. 93, 126406 (2004).
  • van Schilfgaarde et al. (2006) M. van Schilfgaarde, T. Kotani,  and S. Faleev, Phys. Rev. Lett. 96, 226402 (2006).
  • Edel’man (1975) V. S. Edel’man, J. Exp. Theor. Phys. 41, 125 (1975).
  • Dinger and Lawson (1973) R. J. Dinger and A. W. Lawson, Phys. Rev. B 7, 5215 (1973).
  • Schiferl and Barrett (1969) D. Schiferl and C. S. Barrett, J. Appl. Crystallogr. 2, 30 (1969).
  • (32) www.flapw.de.
  • Friedrich et al. (2010) C. Friedrich, S. Blügel,  and A. Schindlmayr, Phys. Rev. B 81, 125102 (2010).
  • Koelling and Harmon (1977) D. D. Koelling and B. N. Harmon, J. Phys. C 10, 3107 (1977).
  • Li et al. (1990) C. Li, A. J. Freeman, H. J. F. Jansen,  and C. L. Fu, Phys. Rev. B 42, 5433 (1990).
  • Aguilera et al. (2013b) I. Aguilera, C. Friedrich,  and S. Blügel, Phys. Rev. B 88, 165136 (2013b).
  • Sakuma et al. (2011) R. Sakuma, C. Friedrich, T. Miyake, S. Blügel,  and F. Aryasetiawan, Phys. Rev. B 84, 085144 (2011).
  • Gonze et al. (1988) X. Gonze, J.-P. Michenaud,  and J.-P. Vigneron, Phys. Scr. 37, 785 (1988).
  • Shick et al. (1999) A. B. Shick, J. B. Ketterson, D. L. Novikov,  and A. J. Freeman, Phys. Rev. B 60, 15484 (1999).
  • Liu and Allen (1995) Y. Liu and R. E. Allen, Phys. Rev. B 52, 1566 (1995).
  • Condrea et al. (2013) E. Condrea, A. Gilewski,  and A. Nicorici, J. Phys.: Condens. Matter 25, 205303 (2013).
  • Young et al. (2011) S. M. Young, S. Chowdhury, E. J. Walter, E. J. Mele, C. L. Kane,  and A. M. Rappe, Phys. Rev. B 84, 085106 (2011).
  • Liu et al. (2011b) W. Liu, X. Peng, C. Tang, L. Sun, K. Zhang,  and J. Zhong, Phys. Rev. B 84, 245105 (2011b).
  • Xiao et al. (2012) S. Xiao, D. Wei,  and X. Jin, Phys. Rev. Lett. 109, 166805 (2012).
  • Zhu et al. (2014) K. Zhu, L. Wu, X. Gong, S. Xiao, S. Li, X. Jin, M. Yao, D. Qian, M. Wu, J. Feng, Q. Niu, F. de Juan,  and D.-H. Lee, arXiv:1403.0066 (2014).
  • Zhang et al. (2015) Q. Zhang, Y. Cheng,  and U. Schwingenschlögl, Scientific Reports 5, 8379 (2015).
  • Perdew and Zunger (1981) J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981).
  • Perdew et al. (1997) J. Perdew, K. Burke,  and M. Ernzerhof, Phys. Rev. Lett. 78, 1396(E) (1997).
  • Kotani and van Schilfgaarde (2002) T. Kotani and M. van Schilfgaarde, Solid State Commun. 121, 461 (2002).
  • Friedrich et al. (2006) C. Friedrich, A. Schindlmayr, S. Blügel,  and T. Kotani, Phys. Rev. B 74, 045104 (2006).
  • Friedrich et al. (2011) C. Friedrich, M. C. Müller,  and S. Blügel, Phys. Rev. B 83, 081101(R) (2011).
  • Michalicek et al. (2013) G. Michalicek, M. Betzinger, C. Friedrich,  and S. Blügel, Comput. Phys. Commun. 184, 2670 (2013).
  • Mostofi et al. (2008) A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Vanderbilt,  and N. Marzari, Comput. Phys. Commun. 178, 685 (2008).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description