## Topological Metal of NaBi with Ultralow Lattice Thermal Conductivity and Electron-phonon Superconductivity

Xing-Qiu Chen, Ronghan Li, Yan Sun, Xiyue Cheng, Dianzhong Li, and Yiyi Li Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Sciences, Shenyang 110016, China (Dated: August 5, 2014)

∗

By means of first-principles and ab initio tight-binding calculations, we found that the compound of NaBi is a three-dimensional non-trivial topological metal. Its topological feature can be confirmed by the presence of band inversion, the derived e ff ective Z2 invariant and the non-trivial surface states with the presence of Dirac cones. Interestingly, our calculations further demonstrated that NaBi exhibits the uniquely combined properties between the electron-phonon coupling superconductivity in nice agreement with recent experimental measurements and the obviously anisotropic but extremely low thermal conductivity. The spin-orbit coupling e ff ects greatly a ff ect those properties. NaBi may provide a rich platform to study the relationship among metal, topology, superconductivity and thermal conductivity.

PACS numbers: 71.20.-b, 71.15.-m, 71.38.-k, 73.20.-r

Because the topological concept was successfully introduced into insulators, various insulators can be classified into topological trivial and non-trivial states[1-4], in which topological insulators are highlighting an exciting scientific frontier of the topological electronic states. In analog of insulators, semimetals can also be classified from topological points as trivial semimetals and topological non-trivial semimetals (TSMs). Among TSMs, there are two classes of peculiar materials, topological Dirac semimetals (TDSs)[5-8] and topological Weyl semimetals (TWSs)[9-18], in which Fermi surfaces are consisted of isolated Fermi points in lattice momentum space. In general, the TDSs are predicted to exist at the critical phase transition point from a normal insulator and a topological one through the spin-orbit coupling e ff ect or by tuning the chemical composition [19, 20]. However, such bulk Dirac points are occasionally degeneracies and not stable. Interestingly, very recently the systems of the P 63 / mmc -Na3Bi [6, 7, 21, 22] and β -BiO2 [23] and Cd3As2 [8, 24-27] have been predicted theoretically and then Na3Bi and Cd3As2 have been experimentally confirmed to be robust TDSs protected by crystal symmetry. TWSs have been theoretically suggested to appear in skutterudite-structure pnictides[14], pyrochlore iridates[16], doped compound Hg1 - -x y Cd Mn Te[15] x y and some constructed heterostructures[17], but to date no experimental verification has been achieved.

Certainly, there is no doubt that the topological concept can be also introduced into metals. Hence, metals would be also classified into two typical types of trivial metals (Ms) and nontrivial topological metals (TMs). In fact, many studies have been focused on the realization and the properties of TMs [2833]. In general, the TMs can be achieved just by the e ff ects of imperfections ( i.e. , chemical doping, strain engineering, heterostructure, etc) on topological insulators. However, to date for this search of native TMs (without any doping and strain applications) no material in reality has been reported successfully. It needs to be emphasized that TMs would indeed extensively exist. Nevertheless, because the topological non-trivial states of specified TMs' surface crossing the Fermi level can be easily mixed by trivial metallic bands, the real realization of TMs indeed poses a challenge.

FIG. 1: Structure, Fermi surface and densities of states (DOSs) of NaBi. (a) The unit cell, (b) the Fermi surface with the SOC e ff ect, (c) and (d) the derived total and projected DOSs without and with the SOC e ff ect, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0491v1-figure-2c104554-0f35-4dcd-ab94-6aacd225d857.jpg)

Within this context, through first-principles calculations with the framework of Density Functional Theory (DFT) by employing the VASP code [39, 40], here we reported a native 3D TM, NaBi, which exhibits the combined interesting properties of the electron-phonon induced superconductivity and the obviously anisotropic but extremely low bulk thermal conductivity. Its topological feature has been analyzed according to the band inversion occurrence between Nas and Bip orbits at the Γ point, the Z2 number based on the derived parities, and the two selected surface non-trivial helical states. Without (with) the spin-orbit coupling (SOC) effect the superconducting transition temperature of Tc is derived to be 1.82-2.59 (2.92-3.75) Kelvin from the electronphonon coupling strength λ = 0.71 (0.84) and the average velocity &lt; ω &gt; ln = 40.8 (38.7) cm -1 , agreeing well with the

experimental findings [34, 35]. In addition, by considering phonon vibrational eigenvalues in the whole of Brillioun zone (BZ) and the phonon relaxation time derived from third-order force constants, we have further revealed that NaBi exhibits an extremely low lattice thermal conductivity but an obviously anisotropic feature of κ a -axis ω = 3.98 Wm -1 K -1 along the a -axis and κ c -axis ω = 1.53 Wm -1 K -1 along the c -axis at room temperature, respectively.

As early as 1932, the compound of NaBi was synthesized to crystallize in a body-centered tetragonal CuAu-type structure (the space group of P 4 / mmm , No.123, see Fig. 1(a)) with Na at the 1 d (1 2, / 1 2, / 1 2) site / and Bi at the 1 a (0, 0, 0) site [36]. The optimized DFT lattice constants [39] of NaBi at the ground state, a = 3.4116 Å and c = 4.9530 Å , are in nice agreement with the experimental lattice constants [36] ( a = 3.46 Å and c = 4.80 Å ). As illustrated in Fig. 1(c and d), NaBi is a typical metal. The SOC inclusion results in several apparent features. Without the SOC inclusion, the Fermi level lies in the declining shoulder of the densities of states (DOS), indicating a relatively high state of N(E F ) = 0.85 states eV -1 f.u. -1 . In contrast, the SOC inclusion significantly reduces the N(E ) to 0.52 states eV F -1 f.u. -1 , due to the fact that the Fermi level now stays at the valley of the pseudogap. In addition, from Fig. 1(d) in the occupied states of the DOS profile the SOC e ff ect even induces the appearance of two obvious peaks dominated by Bip -like states at about -3 eV to -1 eV below the Fermi level, respectively. The presence of those features indicates the significance of the SOC e ff ect for NaBi.

The SOC e ff ect is even more obvious from the electronic band structures in Fig. 2(a and b). Firstly, without the SOC inclusion the three bands (as marked by No.1-3 in Fig. 2(a)) around the Fermi level heavily overlap each other along some high-symmetry lines. The large SOC e ff ect results in their separations, as evidenced in Fig. 2(b). It interprets well as to why two main peaks (corresponding to No.1 and No.2 bands) occur in the occupied states of the DOS profile (Fig. 1(d)). In addition, due to the SOC separation between No.2 and No.3 the Fermi level now locates at the valley of the pesudogap. Secondly, from Fig. 2(a), because of the tetragonal symmetry with c &gt; a and Bi atoms separated by the body-centered Na atom, at Γ the Bi pz orbital is lower in energy than both the degenerated Bi px y , orbital and the Nas orbital. In particular, the band inversion and the anti-crossing feature between Nas and Bipx y , orbitals occur around Γ , showing a nontrivial gap of about 2.5 eV even without the SOC e ff ect. It uncovers that this feature is indeed induced by both the crystal symmetry and the crystal field e ff ect. Furthermore, under the SOC e ff ect and the D 1 4 h symmetry the doubly degenerated Bipx y , orbitals are further split into | Px y -, , ± 3 2 &gt; (as marked by P 1 in Fig. 2(b)) and | Px y -, , ± 1 2 &gt; (as marked by P 2 in Fig. 2b) states. This leads to a reduced non-trivial gap of about 1.1 eV between | Pz -, ± 1 2 &gt; and Na s ( | s + 1 2 , ± 1 2 &gt; ) states at Γ . Despite of the existence of the non-trivial gap, it is intrinsically di ff erent from topological insulators and topological semimetal since NaBi is a typical metal. Therefore, it would be extremely interesting to see

FIG. 2: Bulk and surface electronic band structures of NaBi. (a) and (b) the DFT electronic band structures along the high-symmetry points without and with SOC, respectively. The signs ' + ' and '-' denote the parities of bands at the time-reversal invariant momenta (TRIMs). (c) and (d) corresponds to the (001) and (100) surface electronic band structures derived from the tight-binding model (supplementary materials). The red dots show the helical spin-resolved metallic states on the surfaces. Inset of the panel (d) displays the surface Dirac cone with the enlarged scale around the Fermi level at the Γ point.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0491v1-figure-976f9f67-ac18-4559-a765-cc9a85ca0bd3.jpg)

whether or not NaBi is a non-trivial TM.

To answer this problem, the most important aspect is to elucidate whether or not the continuous energy gap exists between No.2 and No. 3 bands in the whole BZ in Fig. 2(b). On the one hand, we have performed the band structure calculations using a very dense k -mesh set (in total 187836 k -point number in the whole BZ, see Ref.[39]) and the results demonstrated these two bands never touch each other at any k -point and, on the other hand, the calculations even uncovered that the smallest energy gap between No.2 and No.3 bands is about 0.08 eV at the four equivalent ( ± 0.3461 × 2 π a , ± 0.3494 × 2 π a , 0.5 × 2 π c ) points in the k -space kz = 0.5 × 2 π c plane (here, a and c are the lattice constants)[39]. We have also constructed the tight-binding (TB) model Hamilton according to the DFT band structure with the SOC inclusion to further calculate the Berry phase of each energy band in the kz = 0.5 × 2 π c plane. The result uncovers that the Berry phase of the No.2 band is zero, thereby evidencing that the No.2 band never touches the No.1 and No.3 bands. These results fully evidence the existence of a continuous gap between No.2 and No.3 bands in the whole BZ. Given the fact that the non-trivial gap exists between these two bands at Γ , we can further derive the topological invariant, according to the Berry curvature and connection [37]. Interestingly, for the center-symmetric structure (with the inversion symmetry) that NaBi crystallizes in, the e ff ective Z2 invariant can be obtained in terms of the

method proposed by Fu and Kane [38]. Because the bands below No.1 band are fully filled and far away in energy, the topological order just depends on the No.1 and No.2 bands starting from No.1 band around the Fermi level. As shown in Fig. 2(a and b), the product of the parities at the eight timereversal invariant momentums (TRIMs) is -1, corresponding to Z 2 of (1; 0 0 0). It indicates that NaBi is a strong 3D TM with the presence of the topological non-trivial states.

We have further examined the intrinsic surface properties of NaBi. In principles, in similarity to topological insulators, TMs would have an odd number of Dirac cones to appear at any surface orientation because the topological order exists. However, for TMs the behaviors can be highly complex, mainly because the surface Dirac cones perhaps submerge into the bulk metallic states. Therefore, in some orientations it would have no chance to see the presence of surface Dirac cones for TMs. To prove these expectations, we shall now compute the band dispersions for the (001) and (100) surfaces using the ab initio TB model. The ab initio TB model is constructed by downfolding the bulk energy bands, obtained by first-principles calculations using maximally localized Wannier functions (MLWFs). The MLWFs are derived from atomic p -like and s -like states. The surface slab models (with the terminations of Bi atoms) for the (001) and (100) surfaces have been constructed with the thickness of 199 and 399 atomic layers, respectively. The results of the TB calculations are summarized in see Fig. 2(c and d). For the (001) surface, the surface electronic bands (as marked by the solid red circles) connecting the bulk electronic states derived from the No.2 and No.3 bands in Fig. 2(b) cross the Fermi level only once (odd number) for both M -Γ and Γ -X . In addition, for this surface no Dirac cone appears because the surface electronic bands at Γ mix totally with the bulk electronic bands stemmed from the No. 3 band. However, the di ff erent behavior has been observed for the (100) surface [see Fig. 2b]. At Γ the clear Dirac cone appears with surface non-trivial states (as marked by solid red circles) which only once cut the Fermi level in the R -Γ direction. In the Γ -Z direction, there is no crossing at the Fermi level because the surface non-trivial states in this direction submerges into the bulk band states derived from the No.2 band. From the viewpoint of the topology, the cutting number of the Fermi level can be adjusted in di ff erent odd number just by shifting Fermi energy (such as chemical electronic and hole doping treatments). All these facts further evidence that NaBi is a 3D non-trivial TM.

We have utilized the linear response theory and fine k and q meshes [39, 41] to calculate the phonon dispersion, phonon density of states (PHDOS), Eliashberg function ( α 2 F( ω )), and the strength of the electron-phonon ( e -ph) coupling ( λ ω ( )) with and without the SOC inclusion. The phonon spectrum and the phonon densities of states in Fig. 3(a and b) can be divided into two main regions with mostly Bi (but also slightly mixed with Na) modes (0 - 60 cm -1 for SOC and 0 - 85 cm -1 for non-SOC) and highly pure Na modes (80 - 155 cm -1 ) for SOC and 85 - 170 cm -1 for non-SOC). As can be inferred for the phonon DOSs in Fig. 3(b), the SOC inclusion results in the

FIG. 3: Phonon dispersion and electron-phonon coupling strength of NaBi with the SOC inclusion. (a) Phonon dispersion curves along the high symmetry lines of the BZ with the SOC inclusion, (b) total and projected phonon density of states (PHDOS) in NaBi with without / the SOC inclusion, (c) Eliashberg function and the strength of the electron-phonon coupling with without the SOC inclusion. /

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0491v1-figure-bc6b339e-4b9e-4e85-9895-8a5f90b0ca40.jpg)

average softening of over 15% for the transverse modes, and about 10% for the longitudinal ones. The Eliashberg function integrates to a large e -ph coupling strength λ = 0.72 (0.84) without (with) the SOC inclusion but gives the highly low logarithmic average &lt; ω &gt; ln = 40.9 (38.7 for SOC) cm -1 . Although λ is very close to the value of λ ∼ 0.8 for MgB2 which mainly comes from high-frequency boron modes [44], from Fig. 3(c) it is very clear that nearly over 95% of λ in NaBi is generated by the dominated Bi modes in the low-frequency acoustic branches. Strikingly, &lt; ω &gt; ln in NaBi is found to be only one tenth of the MgB2 value of ∼ 450 cm -1 [44]. Using the Allen-Dynes formula [45] and typical µ of 0.14-0.10 we further estimate the Tc in NaBi to be 1.82-2.59 K (2.923.75 K for SOC) (see Table I). Although the estimated data without the SOC inclusion yields a perfect agreement with the experimental data [34], the SOC inclusion indeed exhibits a significant e ff ect on these superconducting parameters.

In particular, it needs to be emphasized that the compound of NaBi have two types of Fermi surfaces: one is a 2D hole Fermi surface (Fig. 1(b)) processing such a shape of quite tetragonal prism centered at the zone center Γ and the other one is a 3D electron Fermi surface (Fig. 1(b)) centered at the zone corner A . These two 2D and 3D Fermi surfaces are obviously originated from the No.2 and No.3 bands (see Fig. 2(b)), respectively. In addition, we also noted that, from Fig. 1(d) since the large gradient of the DOS around the Fermi level which just locates at the valley, the superconducting properties of NaBi may be highly a ff ected by chemical impurities, vacancies, and external strains.

Furthermore, utilizing the frequency ( ω ) of phonon and the phonon group velocity ( v ) in the given transport direction and the derived relaxation time ( τ ) at wave vector q and polarization j within the framework of linear Boltzman's equation, we

TABLE I: Superconducting parameters ( λ -electron-phonon coupling strength, &lt; ω &gt; ln -logarithmic average in cm -1 , Tc -superconducting transition temperature in K , and Θ D -Debye temperature in K ) of NaBi without and with the SOC inclusion. Note that Debye temperature has been derived according to the elastic constants of NaBi with without the SOC inclusion. /

|          |    λ | < ω > ln   | T c       |   Θ D |
|----------|------|------------|-----------|-------|
| non-SOC  | 0.72 | 40.9       | 1.82-2.59 | 147.4 |
| SOC      | 0.84 | 38.7       | 2.92-3.75 | 151.2 |
| Expt[34] | 0.62 |            | 2.15      | 140   |

FIG. 4: Upper panel: the derived lattice thermal conductivities along the a - and c -axis with and without the SOC e ff ect, respectively. Lower panel: the phonon group velocities along the a and c -axis per phonon mode as a function of the distances between the given k point and the centered Γ point (in total, 125 000 k -point number in the BZ) with and without the SOC e ff ect.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0491v1-figure-65908725-c3fb-49c1-95ef-1249a1f058af.jpg)

further derived the bulk lattice thermal conductivity as a function of temperatures as follows,

$$\kappa _ { \omega } = \frac { 1 } { k _ { B } T ^ { 2 } V } \sum _ { q, j } n ( q, j ) [ n ( q, j ) + 1 ] \hslash ^ { 2 } \omega ( q, j ) ^ { 2 } v ( q, j ) _ { z } ^ { 2 } \tau ( q, j ) _ { z } & \text{array} \text{ \quad \text{our ca} \\ \text{contril} \\ \text{throws $b_{-}$ $V$ is and $zero$ bottom count $other$ the central $\iota$ to}$}$$

where kB , V , n and z are Boltzman constant, the crystal volume, and Bose-Einstein distribution as well as the direction of the thermal transportation. Specifically, the phonon's group velocity and specific heat per mode have been derived according to the second-order interatomic force constants obtained by the Phonopy code [42]. The relaxation time is in general determined by third-order force constants, which are derivatives of the total energy with respect to the atomic displacements in any three atoms i , j , and k in directions a , b , and c within a large supercell [46]. We have employed a real-space supercell approach to anharmonic third-order force constant calculations using the script of thirdorder.py [43], which analyzes the symmetries of the crystal and significantly reduce the enormous number of DFT runs that would be required to characterize all relevant third-order derivatives of the energy. This method has been successfully applied to calculate the lattice thermal conductivity for a number of materials (such as, Si, diamond, InAs, and lonsdaleite, etc ) [43, 47, 48].

Currently, our derived temperature-dependent lattice thermal conductivities of NaBi have been compiled in Fig. 4. To our surprising, it can been seen that NaBi exhibits apparent anisotropic but extremely low lattice thermal conductivities. It has been also noted that the SOC e ff ect plays an important role in a ff ecting the κω . In comparison with the non-SOC case in Fig. 4, the SOC e ff ect heavily reduced the thermal conductivities of both a -and c -axes directions. The mechanism is mainly attributed to the that the SOC inclusion indeed results in the softer phonon modes, as compared with those without the SOC e ff ect (Fig. 3b). At room temperature, along the a -axis direction the κω is found to be 4.40 Wm -1 K 1 -(3.98 Wm -1 K 1 ), whereas along the -c -axis direction the κω is extremely low, only about 1.98 Wm -1 K 1 -(1.53 Wm -1 K 1 ) -with (without) the SOC e ff ect. In particular, this low thermal conductivity is indeed comparable to those of widely known materials with ultralow thermal conductivities [49-52], such as PbS, PbSe, PbTe, PtLaSb, and SnSe.

Interestingly, the low lattice thermal conductivity of NaBi exhibits an obviously anisotropic ratio of κ a -axis ω / κ c -axis ω ≈ 2.2 (2.6) without (with) the SOC e ff ect. Its anisotropy can be interpreted well, according to the group velocities of the phonon as illustrated in Fig. 4(b-m) in which the group velocities, vx along the a -axis (Fig. 4(a-f)) and vz along the c -axis (Fig. 4(b-m)), have been visualized as a function of the k -space distances between any phonon mode in the whole BZ and the zone centered Γ point. No matter whether the SOC e ff ect is included, both vx and vz show the quite similar character, as evidenced in Fig.4(b-m). The acoustic modes play a main role in determining the lattice thermal conductivities. In particular, along the a -axis direction the phonon group velocities are overall larger than those along the c -axis direction (Fig. 4(bd) for a-axis and Fig. 4(h-j) for c-axis), thereby resulting in a higher κω along the a -axis. In general, the optical modes nearly makes no contributions to the κω . However, based on our calculations for NaBi, the optical modes make a certain contribution to the κω along the a -axis direction. For instance, at room temperature along the a -axis direction the κω that the optical modes contributed to is about 0.822 (0.742) Wm -1 K 1 -without (with) the SOC e ff ect, being about 18% of the whole κω . This is mainly because its optical modes exhibit very large group velocities along the a -axis direction (Fig. 4(e-g)) and, in the meanwhile, the low frequencies (Fig. 3b) which are comparable to those of the acoustic modes. In contrast, from Fig. 4(k-m) the optical modes almost contribute nothing to the κω along the c -axis direction. This fact further enlarges the anisotropic ratio of the κω .

In summary, through first-principles calculations we have

found that NaBi is an intrinsic 3D TM with the combined properties of electron-phonon induced superconducting and obviously anisotropic but extremely low lattice thermal conductivity. The SOC e ff ect has been demonstrated to have the significant impacts on those properties. Compared to topological insulators and topological semimetals, our results for NaBi suggest that the topological metal can be realized in a simple body-centred tetragonal structure without any doping or strain treatments.

Acknowledgements We thank Z. Fang for useful discussions. This work was supported by the 'Hundred Talents Project' of the Chinese Academy of Sciences and from the Major Research Plan (Grand Number: 91226204) of the NSFC of China (Grand Numbers:51174188 and 51074151) and Beijing Supercomputing Center of CAS (including its Shenyang branch) as well as Vienna Supercomputing Center (VSC cluster).

- ∗ Corresponding author:xingqiu.chen@imr.ac.cn
- [1] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95 , 226801 (2005).
- [2] B. A. Bernevig, Taylor L. Hughes, and Shou-Cheng Zhang, Science 314 ,1757 (2006).
- [3] X.-L. Qi, and S.-C. Zhang, Rew. Mod. Phys. 83 , 1057 (2011).
- [4] M. Z. Hasan, and C. L. Kane, Rew. Mod. Phys. 82 , 3945 (2010).
- [5] S. M. Young, S. Zaheer, J. C.Y. Teo, C. L. Kane, E. J. Mele, and A. M. Rappe, Phys. Rev. Lett. 108 , 140405 (2012).
- [6] Z. Wang, Y. Sun, X.-Q. Chen, C. Franchini, G. Xu, H. Weng, X. Dai, Z. Fang, Phys. Rev. B 85 , 195320 (2012).
- [7] X. Y. Cheng, R. H. Li, Y. Sun, X.-Q. Chen, D. Z. Li and Y. Y. Li, Phys. Rev. B, 89 , 245201 (2014).
- [8] Z. Wang, H. Weng, Q. Wu, X. Dai, and Z. Fang, Phys. Rev. B, 88 , 205101 (2013).
- [9] H. Weyl, Z. Phys. 56 , 330 (1929).
- [10] H. B. Nielsen and M. Ninomiya, Phys. Lett. B 130 , 389 (1983).
- [11] Y. Xu, R.-L. Chu, and C.W. Zhang Phys. Rev. Lett. 112 , 136402 (2014).
- [12] R. R. Biswas, and S. Ryu, Phys. Rev. B 89 , 014205 (2014).
- [13] T. Ojanen, Phys. Rev. B 87 , 245112 (2013).
- [14] V. Pardo, J. C. Smith and W. E. Pickett, Phys. Rev. B 85 , 214531 (2012).
- [15] D. Bulmash, C.-X. Liu, and X.-L. Qi, Phys. Rev. B 89 , 081106(R) (2014).
- [16] X.G. Wan, Ari M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B 83 , 205101 (2011).
- [17] T. Das, Phys. Rev. B 88 , 035444 (2013).
- [18] J.-H. Zhou, H. Jiang, Q. Niu, J.-R. Shi, Chin. Phys. Lett. 30 , 027101 (2013).
- [19] S. Murakami, New J. Phys. 9 , 356 (2007).
- [20] S. M. Young, S. Chowdhury, E. J. Walter, E. J. Mele, C. L. Kane, and A. M. Rappe, Phys. Rev. B 84 , 085106 (2011).
- [21] Z. K. Liu, et al. Science 343 , 864 (2014).
- [22] S.-Y. Xu, et al. Preprint at http: // arXiv: 1312.7624 (2013).
- [23] S. M. Young, S. Zaheer, J. C. Y. Teo, C. L. Kane, E. J. Mele, and A. M. Rappe, Phys. Rev. Lett. 108 , 140405 (2012).
- [24] M. Neupane, et al. Nature Communications, 5, 3786 (2014)
- [25] S. Borisenko, Q. Gibson, D. Evtushinsky, V. Zabolotnyy, B.

Buechner, R. J. Cava, Phys. Rev. Lett., 113 , 027603 (2014).

- [26] Z. K. Liu, et al. Nature Mater. 13 , 677 (2014)
- [27] S. J. Seon, et al. , Nature Mater. 13 , (2014), DOI: 10.1038 NMAT4023 /
- [28] J.W. Wells, J. H. Dil, F. Meier, J. Lobo-Checa, V. N. Petrov, J. Osterwalder, M. M. Ugeda, I. Fernandez-Torrente, J. I. Pascual, E. D. L. Rienks, M. F. Jensen, and P. Hofmann, Phys. Rev. Lett. 102 , 096802 (2009).
- [29] D. Xiao, Y.G. Yao, W.X. Feng, J. Wen, W.G. Zhu, X.-Q. Chen, G. Malcolm Stocks, and Z.Y. Zhang, Phys. Rev. Lett. 105 , 096404 (2010).
- [30] H. Lin, L. AndrewWray, Y. Xia, S.Y. Xu, S. Jia, R. J. Cava, A. Bansil and M. Z. Hasan, Nature Mater. 9 , 546 (2010).
- [31] D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, and M. Z. Hasan, Nature (London), 452 ,970 (4522008).
- [32] M. Brahlek, N. Bansal, N. Koirala, S.-Y. Xu, M. Neupane, C. Liu, M. Z. Hasan, and O. Seongshik, Phys. Rev. Lett. 109 , 186403 (2012).
- [33] Z. Y. Zhu, Y.C. Cheng, and U. Schwingenschlogl, Phys. Rev. Lett. 110 , 077202 (2013).
- [34] S. K. Kushwaha, J. W. Krizan, J. Xiong, T.Klimczuk, Q. D. Gibson, T. Liang, N. P. Ong, and R. J. Cava, J. Phys.: Condens. Matter 26 212201 (2014).
- [35] J. M. Reynolds, and C. T. Lane, Phys. Rev. B, 79 , 405 (1950).
- [36] E. Zintl and W. Dullenkopf, Z. Phys. Chem. Abt. 16 , 183 (1932).
- [37] L. Fu and C. L. Kane, Phys. Rev. B, 74 , 195312 (2006).
- [38] L. Fu and C. L. Kane, Phys. Rev. B, 76 , 045302 (2007).
- [39] See Supplemental Material http: // link.aps.org supplemental xxx. / / The supplemental materials contains (1) the computational methods of electronic structure (VASP [40]), vibrational phonon dispersions (Phonopy [42] and QE [41]), electronphonon coupling coe ffi cients (QE [41]), and lattice thermal conductivities (ShengBTE [43]) and (2) the detailed analysis of the band gap between No.2 and No.3 bands in the whole BZ.
- [40] G. Kresse and J. Furthmu¨ uller, Phys. Rev. B 54 , 11169 (1996).
- [41] P. Giannozzi, S. Baroni, N. Bonini, et al. , J. Phys.: Condens. Matter 21 395502 (2009).
- [42] L. Chaput, A. Togo, I. Tanaka, and G. Hug, Phys. Rev. B, 84 , 094302 (2011).
- [43] W. Li, J. Carrete, N. A. Katcho, N. Mingo, Comp. Phys. Comm., 185 , 1747-1758 (2014).
- [44] A. Y. Liu, I. I. Mazin, and J. Kortus, Phys. Rev. Lett. 87 , 087005 (2001).
- [45] P. B. Allen and R. C. Dynes, Phys. Rev. B 12 , 905 (1975).
- [46] A. Ward, D. A. Broido, D. A. Stwwart and G. Deinzer, Phys. Rev. B, 80 , 125203 (2009).
- [47] L. Lindsay, D. A. Broido, and T. L. Reinecke, Phys. Rev. Lett., 111 , 025901 (2013).
- [48] L. Chaput, Phys. Rev. Lett., 110 , 265506 (2013).
- [49] J. M. Skelton, S. C. Parker, A. Togo, I. Tanaka, A. Walsh, Phys. Rev. B, 89 , 205223 (2014).
- [50] K. Biswas, J.Q. He, I. D. Blum, C.-I. Wu, T. P. Hogan, D. N. Seidman, V. P. Dravid, M. G. Kanatzidis, Nature, 489 , 414-418 (2012).
- [51] J. Carrete, W. Li, N. Mingo, S. D. Wang and S. Curtarolo, Phys. Rev X, 4 , 011019 (2014).
- [52] L.-D. Zhao, S.-H. Lo, Y.S. Zhang, H. Sun, G. J. Tan, C. Uher, C. Wolverton, V. P. Dravid and M. G. Kanatzidis, Nature, 508 , 373-377 (2014).
Xing-Qiu Chen Ronghan Li Yan Sun Xiyue Cheng Dianzhong Li Yiyi Li
Topological Metal of NaBi with Ultralow Lattice Thermal Conductivity and Electron-phonon Superconductivity
By means of first-principles and \emph{ab initio} tight-binding calculations,
we found that the compound of NaBi is a three-dimensional non-trivial
topological metal. Its topological feature can be confirmed by the presence of
band inversion, the derived effective Z$_2$ invariant and the non-trivial
surface states with the presence of Dirac cones. Interestingly, our
calculations further demonstrated that NaBi exhibits the uniquely combined
properties between the electron-phonon coupling superconductivity in nice
agreement with recent experimental measurements and the obviously anisotropic
but extremely low thermal conductivity. The spin-orbit coupling effects greatly
affect those properties. NaBi may provide a rich platform to study the
relationship among metal, topology, superconductivity and thermal conductivity.