Topological electronic structure of YbMg2Bi2 and CaMg2Bi2

Zintl compounds have been extensively studied for their outstanding thermoelectric properties, but their electronic structure remains largely unexplored. Here, we present a detailed investigation of the electronic structure of the isostructural thermopower materials YbMg2Bi2 and CaMg2Bi2 using angle-resolved photoemission spectroscopy (ARPES) and density functional theory (DFT). The ARPES results show a significantly smaller Fermi surface and Fermi velocity in CaMg2Bi2 than in YbMg2Bi2. Our ARPES results also reveal that in the case of YbMg2Bi2, Yb-4f states reside well below the Fermi level and likely have a negligible impact on transport properties. To properly model the position of 4f-states, as well as the overall electronic structure, a Hubbard U at the Yb sites and spin-orbit coupling (SOC) have to be included in the DFT calculations. The theoretical results reveal that both materials belong to a Z2 topological class and host topological surface states around EF. Due to the intrinsic hole doping, the topological states reside above the Fermi level, inaccessible by ARPES. Our results also suggest that in addition to SOC, vacancies and the resulting hole doping play an important role in the transport properties of these materials.

Among Zintl families, AB 2 X 2 (A = Ca, Yb, Eu, Sr; B = Zn, Mn, Cd, Mg; X = Sb, Bi) type compounds gain more attention as many of them show high thermoelectric performance [19][20][21][22][23][24][25][26][27] . Recent works have demonstrated the coexistence of intrinsic magnetism, topological Dirac electronic states, and moderate thermopower efficiency in EuMg 2 Bi 2 5,12,13,21,28 , which could provide an interesting playground to investigate the interplay between magnetism, topology, and thermoelectricity. Interestingly, it has been found that when the A-sites of these compounds contain rare-earth elements or are partially substituted by rare-earth elements, they show enhanced carrier mobility and carrier concentration compared to the alkaline-earth compounds 26,29,30 . For example, going from CaMg 2 Bi 2 to YbMg 2 Bi 2 , both the mobility and hole concentration are increased by a factor of~2 30 . This is quite unexpected, as for chemically-doped semiconductors, the mobility generally decreases with increasing carrier concentration 31 .
The stronger hybridization between Bi and Yb/Eu was considered to be partially responsible for higher mobilities in Yb/Eu-based compounds 21 . On the other hand, if the 4f states were near the valence-band maximum, one could expect lower mobilities due to heavy bands and stronger scattering (impurity, electron-phonon, electron-magnon), exactly opposite to the experimental results 21,30 . Flage et al. 32 have reported that in the prototype system YbZn 2 Sb 2 , the Yb-4f states make a nonnegligible contribution to the valence band edge. The exact effects of heavy 4f electrons on carrier concentrations and transport properties in rare-earth-based compounds are still unclear due to the lack of experimental studies.
To achieve a better TE performance, high carrier mobility, heavy effective mass, and low lattice thermal conductivity are highly desirable 1,21 . Various strategies have been proposed to achieve this goal, such as band convergence and resonant states for heavy effective mass [33][34][35] , band alignment to achieve high carrier mobility 21,36 and introducing microstructural defects (point defects and nanostructures) and alloying to enhance phonon scattering 21,34,37 . However, many of the theoretical calculations 25,38-40 do not take into account the effects of spin-orbit coupling (SOC) in the prediction of the transport properties. The effect of SOC could change the band degeneracy (which is related to carrier concentration), band hybridization, and band gap. Thus, the consideration of the effects of SOC is necessary for accurate predictions of new relevant materials and their associated properties. In addition to the more realistic modeling, the experimental measurement of the electronic structure is an irreplaceable component of studies of new materials.
Here, we report comprehensive experimental and theoretical studies of the electronic structures of YbMg 2 Bi 2 and CaMg 2 Bi 2 using ARPES and first-principles calculations. Our results show the importance of inclusion of the SOC in describing these materials. We also show that the Yb-4f states are far away from E F , implying that they do not play a role in transport. Our studies show that both materials are narrow band-gap topological insulators with topological surface states. Figure 1 shows the schematics of the crystal structure, x-ray diffraction (XRD) and the low-energy electron diffraction (LEED) patterns of the AMg 2 Bi 2 (A = Yb, Ca) samples. Figure 1a, b show the trigonal unit cell and the top views of the (001) surface (a-b plane) of the crystal structure, respectively 41 . Using the Zintl concept, the trigonal structure of AMg 2 Bi 2 compounds can be viewed as polyanionic [Mg 2 Bi 2 ] 2− layers stacked along the c-axis and separated by the trigonal layers of A 2+ . The roomtemperature powder XRD patterns taken on crushed YbMg 2 Bi 2 and CaMg 2 Bi 2 single crystals are shown in Fig. 1c, d, respectively. The experimental data fit well with the CaAl 2 Si 2 -type crystal structure having a trigonal lattice with space group P3m1 (No. 164) 42 . The refined lattice parameters are a = b = 4.7258(4) Å and c = 7.6453(14) Å for YbMg 2 Bi 2 and a = b = 4.7236(3) Å and c = 7.6512(10) Å for CaMg 2 Bi 2 , in good agreement with the earlierreported values 30 . LEED patterns taken from in-situ UHV-cleaved YbMg 2 Bi 2 and CaMg 2 Bi 2 crystals are shown in Fig. 1e, f, respectively. Hexagonal patterns are obtained for both samples, confirming that the cleaved surface is the (001) plane. The schematics of the hexagonal bulk Brillouin zone and its surface projection, i.e. surface Brillouin zone (SBZ), are shown in Fig. 1g. Figure 2 shows the bulk electronic structure along various symmetry paths for AMg 2 Bi 2 . Figure 2a, b illustrate the band dispersions of YbMg 2 Bi 2 , without and with spin-orbit coupling (SOC), respectively. A Hubbard U parameter of 4 eV is included to describe the localized 4f states of Yb and to match the observed splitting. Multiple Dirac-like band crossings are observed (e.g., at K and between K and M symmetry points) in the absence of SOC; however, most of the Dirac points become gapped in the presence of SOC due to lifting of band degeneracy. The Yb-4f derived flat band that was initially at about −1.3 eV (Fig. 2a) splits into two levels, 4f 5/2 (−2.85 eV) and 4f 7/2 levels (−1.25 eV) in the presence of SOC. We also see the strong hybridization between the 4f states and highly-dispersive bands when they cross each other. To better understand the role of 4f states, we have also treated the 4f electrons as core electrons and the resulting band structure is shown in Fig. 2c. Note that the low-energy electronic structure is very similar to the one in which the 4f electrons are considered as valence electrons (Fig. 2b), except for a slight downward shift of the hole-bands around the Γ point. To understand whether this apparent energy shift is related to the intrinsic properties of the 4f electrons, we performed the calculation considering higher U values (see Supplementary note 1). By increasing the U value from 4 eV to 8 eV, we essentially move the 4f states much deeper in energy, but still observe a similar low-energy electronic structure without any energy shift of the hole-bands. This suggests that the energy shift is related to inclusion of the external parameter U in the DFT calculation. Without SOC a spectral gap of about 0.1 eV can be observed between the valence and conduction bands, and the system behaves like a semiconductor while it is nearly gapless in the presence of SOC. Interestingly, a closer look around the Γ point (inset of Fig. 2c) reveals that a small energy gap of about~15 meV is still present and the shapes of the valence band maximum (VBM) and conduction band minimum (CBM) suggest that there might be a band inversion. The parity of the bands at Γ also reverses when SOC is turned on, suggesting a band inversion scenario. Similar electronic structures have also been demonstrated in previous theoretical studies; however, they were conducted without SOC 21 .

Theoretical bulk electronic structures
Similarly, Fig. 2d, e illustrate the electronic structure of CaMg 2 Bi 2 without and with SOC, respectively. A spectral gap of about 0.5 eV can be seen between the VBM and CBM in the absence of SOC while the gap is extremely reduced when SOC is turned on. This is consistent with the previous reports 21, 27 Fig. 1 Crystal structure of AMg 2 Bi 2 (A = Yb, Ca). a Unit cell of trigonal AMg 2 Bi 2 . b Projection of the crystal structure onto the a-b plane with a unit cell indicated by a rhombus. The crystal structures were generated using the VESTA software 41 . c, d XRD patterns of powdered YbMg 2 Bi 2 and CaMg 2 Bi 2 crystals, along with the refinement of the data using FULLPROF software 42 . e, f Hexagonal LEED patterns of YbMg 2 Bi 2 and CaMg 2 Bi 2 single crystals obtained after cleaving the samples in-situ under UHV, respectively. g Schematics of hexagonal bulk Brillouin zone and its surface projection (SBZ). To correctly verify the nontrivial topology of these materials hinted by the shape of the valence/conduction bands near the zone center, we have calculated the Z 2 topological numbers using the Wilson loop (Wannier charge center) method 43 for the six time-reversal invariant momentum planes. The obtained results show Z 2 topological numbers v 0 ;(v 1 v 2 v 3 ) = 1;(000) for both systems, which indicates that both materials are strong topological insulators. Here, we want to point out that although the lowenergy electronic structure is very similar to the type-II nodal-line semimetal Mg 3 Bi 2 11 , the present materials are not nodal-line semimetals due to the presence of a band gap~0.10-0.5 eV without SOC. In contrast, nodal-line semimetals should show conduction and valence-band crossing in the absence of SOC. Figure 2f, g show the partial densities of states (PDOS) of YbMg 2 Bi 2 and CaMg 2 Bi 2 , respectively, in the presence of SOC. Both systems show very similar PDOS for Bi and Mg atoms, while they differ significantly for Ca and Yb atoms. Yb contributes more to the valence and conduction bands than Ca. The contribution of Yb-f states is strongly localized around −1 eV, and it diminishes around the band edges. This implies a negligible role of 4f electrons in the transport properties in this material. In order to obtain a better understanding of the orbital character of the bands and their hybridization, orbital-resolved band dispersions are shown with and without SOC for YbMg 2 Bi 2 ( Fig. 3a-f) and CaMg 2 Bi 2 ( Fig. 3g-l) systems, respectively. The contributions from Ca/Yb-s, Bi-p, and Mg-s are shown. It can be seen that the outer hole-like valence band has dominant Bi-p character while the inner one has mixed Bi-p, Mg-s, and Ca-s/Yb-s character. It appears that all orbitals are strongly hybridized with each other, both in the presence and absence of SOC. We note that some contributions of other orbitals such as Bi-s, Mg-p, Ca-p, Yb-d, and Yb-p are also present around E F (not shown).
Comparison between the experimental and theoretical band dispersions: YbMg 2 Bi 2 Figure 4 shows the comparison between the experimental and the calculated electronic structure of YbMg 2 Bi 2 . Panel 4 a represents the measured ARPES spectrum along the Γ À M line, which, at the photon energy used in the experiments, is very close to the bulk Γ-M line. The calculated bulk electronic structure along the Γ-M line is superimposed for comparison. To match the experimental dispersions, the theoretical spectrum is shifted up by~0.3 eV, implying that the measured sample is heavily hole-doped. According to the previous reports, intrinsic hole doping in these samples is predominantly due to A-site vacancies 21,30 . In addition to numerous dispersing bands, two flat bands at energies −1.15 eV and −2.5 eV are visible in the ARPES spectrum. These flat bands are better resolved when probed with the He-II photon energy (40.8 eV) as shown in Supplementary note 2. They arise from spin-orbit-split Yb-4f states, matching nicely with the calculations when SOC and U = 4 eV are taken into account. It is also visible that the dispersive bands and the Yb-4f bands hybridize, in accord with our theoretical results. We note that the position of the Yb-4f states varies among different Ybcompounds 26,44 , implying that U is not unique. In previous studies, a reasonable value of U has been found to be 5.3 eV for YbZn 2 Sb 2 26 and 8 eV for YbB 12 44 . High-resolution ARPES spectra in the vicinity of E F , near the zone center in the first and second SBZ, are shown in panels b and c, respectively. Both spectra show linearly-dispersive hole-like bands crossing the Fermi level without any measurable renormalization, indicative of weak electron correlation in YbMg 2 Bi 2 . A visible change in the Fermi wave vector (k F ) of both the inner and outer bands between the first and second SBZ originates from the difference in probed k Z for the two spectra and points to their bulk electronic origin. This change of k Z can be better realized from the momentum-dispersion curves (MDC) around E F as shown in Supplementary note 3. The   bulk origin of these states is in good agreement with the recent photon-energy-dependent ARPES studies on similar materials 10,11 . These linearly-dispersing bands form nearly-circular Fermi surfaces, Fig. 4d. The constant-energy contours at −0.11 eV and −0.39 eV are also shown in Fig. 4e, f, respectively. The inner two contours are very closely spaced and can only be resolved away from the Fermi level. The outer contour exhibits hexagonal warping that strengthens with increasing binding energy. The 3-fold symmetry of the intensity of the outer contour also suggests its significant bulk-band character. Similar contours are obtained from the bulk electronic structure near the Γ point of the 3D BZ, as shown in Fig. 4g. Since ARPES is a surface-sensitive probe and the bulk bandstructure calculations imply the existence of topological surface states, we have also calculated the surface electronic structure for various possible terminations of the (001) surface. The Yb-4f states are not included in these calculations for simplicity. Figure 4h-j show the calculated surface electronic spectra from a (001) semi-infinite slab terminated by Yb, Mg, and Bi, respectively, where the enlarged views near the Fermi energy around Γ are shown in the left insets of each figure. The spectral brightness indicates the integrated charge density over the top six atomic layers, roughly equivalent to the ARPES probing depth.
It is clear that a topological surface state (TSS), spanning the bulk band gap, is present on every termination. The Dirac point is buried inside the valence band, but should be exposed on Yb-and Bi-terminations. A similar TSS was also reported in other topological materials where the Dirac point lies inside the valence band, such as p-type Bi 2 Te 3 and Mg 3 Bi 2 11,18,45 . The topological character of TSS is further supported by its helical spin-texture as shown in Fig. 4n. Overall, the calculated spectra for different terminations are very similar, resembling the bulk band features. A few noticeable differences are: (1) in the case of Mg termination, the outer state seems sharper, while for Yb and Bi terminations the inner state is more pronounced and (2) the Yb and Bi terminations show an intense surface state around Γ at −0.85 eV, while the Mg termination does not. A direct side-by-side comparison of the ARPES spectra with the calculated surface state spectra (Fig. 4k, l) indicate somewhat better agreement with the Mg termination. However, the theoretical spectra for the Mg termination show only two states forming the Fermi surface (Fig. 4m), whereas three states are observed experimentally. On the other hand, the Yb terminations display an additional surface-resonance band (SR) (inset of Fig. 4h).
The coexistence of different terminations would probably give the best agreement with the observed spectra, as their contributions would be averaged over the macroscopic size of the excitation spot. To get further insight about the surface structure, we performed atomic force microscopy (AFM) measurements and the results are shown in Supplementary note 6. On the opticallyflat surface regions, we observe different step heights, suggesting different surface terminations. The coexistence of different surface terminations has been observed for various materials such as Pt 3 Te 4 , YBa 2 Cu 4 O 8 and Bi 4 Se 3 using micro/nano ARPES and scanning-tunneling microscopy (STM) [46][47][48] . Similar microscopic studies would be needed to disentangle the contributions from different terminations in the present materials.
Potassium-deposited YbMg 2 Bi 2 Due to high hole doping of pristine samples, the VBM, TSS and conduction band are unoccupied and cannot be probed in conventional ARPES. Therefore, we have performed in-situ electron-doping of cleaved surfaces by potassium deposition. The obtained results are shown in Fig. 5. Figure 5a shows the spectrum of the pristine YbMg 2 Bi 2 surface in the second SBZ. The second SBZ is chosen because bands near E F show a higher photoemission intensity. Figure 5b-d show the spectra after 3, 7, and 12 min of potassium deposition, respectively. It can be seen that the states shift towards higher binding energy and become blurred due to the increase in scattering rates. In the last sequence, an intensity (V-shaped-like feature) very close to E F can be seen in Fig. 5d, suggesting that the CBM just becomes occupied. We note that in Fig. 5d, the ARPES intensity was divided by the Fermi-Dirac function to better resolve the spectral feature around E F . The raw ARPES data shown in Supplementary Fig. 4, are consistent with that conclusion. To better understand the change of electronic structure due to the potassium deposition, we have plotted the second derivatives of the spectra in Fig. 5e-h and the corresponding energy distribution curves (EDCs) in Fig. 5i. The energy shift of the 4f 7/2 peak with potassium deposition time is shown in Fig.  5j. A very similar energy shift is also observed for the Bi 5d core-levels (Supplementary note 5), suggesting a nearly rigidband-shift scenario. Thus, on top of the second derivative plots (Fig. 5e-h), the calculated bulk bands (red curves) are superimposed and shifted to match the experimental spectra by using the energies from Fig. 5i. Even though there is a relatively good agreement between theory and experiment, we could not resolve the TSS and conduction band edge individually as the electronic states become too diffuse. We note that the deposition of potassium on the surface alters the surface potential, causing band-bending in the relatively thick range from the surface (~10-50 nm) 49 . Previous ARPES studies have shown that energy shifts of valence and conduction states are nearly equal in topological insulators 49,50 . In case of YbMg 2 Bi 2 , the observed band structure shift is nearly-rigid, similar to that observed in other topological insulators 49,50 .
Comparison between experimental and theoretical band dispersions: CaMg 2 Bi 2 Figure 6 shows a comparison between the experimental and the calculated electronic structure of CaMg 2 Bi 2 . Figure 6a represents the ARPES spectrum close to the K À Γ À K line. Theoretical electronic structures along the Γ-K and A-H lines are superimposed on Fig. 6a. It is clear that most of the spectral features are well reproduced by the calculation, except for the states indicated by an arrow. For better visualization of the states in the vicinity of E F , a high-resolution ARPES spectrum is shown in Fig. 6b and its second derivative in Fig. 6c. It is evident that two sets of hole-like bands cross the Fermi level. The FS and a constant energy contour (at −0.1 eV) formed by these bands are shown in Fig. 6d, e, respectively. The shape of these contours is similar to that for YbMg 2 Bi 2 and well reproduced by the theoretical calculation, as shown in Fig. 6f. Figure 6g-i show the calculated surface electronic spectra from a (001) semi-infinite slab terminated by Ca, Mg, and Bi, respectively. As in the case of YbMg 2 Bi 2 , a TSS (indicated by arrows) is present on every termination. In addition to the TSS, the Bi-terminated surface shows trivial surface states near the VBM that are absent on the Ca-and Mg-terminated surfaces. These surface states originate from the unsaturated surface dangling bonds 51,52 . The absence of such states in the measured spectra (Fig. 6a-c) indicates that the surface is not Bi-terminated. Direct comparisons with the calculated spectra for Ca and Mg terminations are shown in Figs. 6j, k, respectively. Spectral features around Γ are in somewhat better agreement with the Ca-termination, while near K the agreement is better with the Mg termination. Therefore, the coexistence of Mg and Ca terminations would probably give the best agreement with the observed spectra. The AFM measurements (Supplementary note 7) showing different step heights support that scenario. We note that according to the Zintl formalism both YbMg 2 Bi 2 and CaMg 2 Bi 2 should be nominally charge balanced and predicted to be narrow-gap semiconductors. However, samples are usually highly p-doped, suggesting that holes are generated due to intrinsic defects 30 . According to the previous Hall studies 21,30 , the estimated values of carrier concentrations and mobilities are lower in CaMg 2 Bi 2 than in the rare-earth-based compounds YbMg 2 Bi 2 and EuMg 2 Bi 2 . From ARPES measurements it is also possible to estimate the carrier concentrations by measuring the volume of the FS. In Fig. 7a, the MDC-derived bands crossing the Fermi level in YbMg 2 Bi 2 and CaMg 2 Bi 2 are plotted together. The Fermi wave vectors (k F ) of both inner and outer bands are larger in YbMg 2 Bi 2 than in CaMg 2 Bi 2 , implying higher carrier concentrations. By assuming an isotropic spherical shape of the two FSs, formed by the two hole-like bands, we estimate the hole concentrations to be~1.4 × 10 20 cm −3 and~5 × 10 19 cm −3 for YbMg 2 Bi 2 and CaMg 2 Bi 2 , respectively. The values estimated from the Hall coefficient are 4.6 × 10 19 cm −3 and 1.7 × 10 19 cm −3 , respectively 30 . Also, CaMg 2 Bi 2 has a lower Fermi velocity than YbMg 2 Bi 2 . This behavior can be realized by putting the calculated valence bands of both materials on the same graph with their observed fillings (Fig. 7b).
Based on the shape and sharpness of the ARPES spectra, it seems that the effective masses and scattering rates are similar in these systems, which could lead to similar mobilities. However, the exact quantification of mobility is not possible from ARPES as the transport scattering rates are, in general, different from quasiparticle ones and the relationship between the effective mass and mobility for non-parabolic bands might be quite complex 53 . Moreover, the origin of higher mobilities in YbMg 2 Bi 2 is still unclear and remains an open question.
We note that in predicting/interpreting the thermoelectric properties of these classes of materials, some studies disregarded SOC 25,[38][39][40] , leading to overestimated gaps, while the experiments and inclusion of SOC in calculations led to better agreement.
In summary, we have probed the electronic structure of the Zintl-phase isostructural analogs YbMg 2 Bi 2 and CaMg 2 Bi 2 using ARPES studies complemented with first-principles calculations. It is shown that the different carrier concentrations in these materials appear primarily due to a relative shift in E F with respect to the valence band edge. Our results confirmed that the localized 4f electrons do not play a considerable role in the low-energy electronic structure of YbMg 2 Bi 2 . We have also shown the importance of the SOC effect in these materials. Calculations suggest that these materials are narrow-band-gap topological insulators that support topological surface states with helical spin-texture. Due to the intrinsic hole doping, the actual E F is located well inside the valence bands in these materials, which further restricts the effect of nontrivial band topology to be reflected in the transport properties. By tuning the E F to the bulk band gap through chemical substitution or strain could lead to the observation of anomalous surface transport properties such as Shubnikov-de Haas oscillations, opening the possibility for their use in spintronics applications.

Single crystal growth
Single crystals of CaMg 2 Bi 2 and YbMg 2 Bi 2 were grown with the self-flux solution-growth technique. The high-purity starting materials Ca (99.999%), Mg (99.98%), Bi (99.9999%) from Alfa Aesar and Yb from Ames Laboratory were used in nominal compositions CaMg 4 Bi 6 and YbMg 4 Bi 6 . The elements were placed into alumina crucibles and sealed under ≈ 1/4 atm of high-purity inert Ar gas inside silica tubes. The assembly was then heated to 900°C at a rate of 50°C/h where the temperature was maintained for 12 h. The tubes were then cooled to 850°C in 1 h followed by cooling to 750°C in 10 h and finally to 650°C in 24 h. Bulk three-dimensional-shaped single crystals with trigonal facets were obtained by removing the excess flux using a centrifuge at the last temperature.

Experimental details
Bulk structural characterizations were carried out by room-temperature powder x-ray diffraction (XRD) measurements on crushed single crystals using a Rigaku Geigerflex x-ray diffractometer with Cu-K α radiation. The structural refinement was performed through Rietveld refinement using FULLPROF software package 42 . For the ARPES measurements, the samples were mounted on a sample holder and cleaved in-situ just before the measurements. The surface structure of the cleaved surface was checked using low-energy electron diffraction (LEED) and an atomic-force microscope (AFM) (Digital Instruments, NanoScope III). The ARPES experiments were carried out at OASIS-laboratory at Brookhaven National Laboratory using a Scienta SES-R4000 electron spectrometer with monochromatized He I α (21.22 eV) and He II α (40.8 eV) radiation (VUV-5k) 54 . The total instrumental energy resolution was~10 meV for He I α and 20 meV for He II α . The angular resolution was better than 0.15°a nd 0.4°along and perpendicular to the slit of the analyzer, respectively. Most of the data were taken at~20 K.

Theoretical details
First-principles calculations have been carried out using Quantum ESPRESSO [55][56][57] . For the exchange and correlation energy/potential we used the PBEsol functional 58 . The projected-augmented-wave 59 method has been used to represent the core electrons. All the calculations have Fig. 7 Comparison of valence bands in YbMg 2 Bi 2 and CaMg 2 Bi 2 . a MDC-derived bands of YbMg 2 Bi 2 and CaMg 2 Bi 2 near the Fermi level. The black (red) circles correspond to YbMg 2 Bi 2 (CaMg 2 Bi 2 ). Band slopes at E F are highlighted by light-blue lines. b Shapes of the valence bands from our calculations. The CaMg 2 Bi 2 spectra is shifted downward by 0.15 eV to reproduce the experimental features. Experimental E F position is marked by dotted-line.
been carried out using experimental lattice parameter, as noted in the 'structural details of samples' subsection. The cut-off energy for the plane waves and charge density were set to 60 Ry and 520 Ry, respectively. A k-mesh of 12 × 12 × 10 has been used for the Brillouinzone integration. The electronic structure calculations have been carried out both with and without SOC. Pw2wannier interface and WANNIER90 were used for the construction of the first-principle tight binding Hamiltonian 60 . Surface-states spectra have been calculated using the WANNIERTOOLS package 61 .

DATA AVAILABILITY
The data that support the findings of this study are available from the corresponding author upon reasonable request.