# Full text of "Combined Fast Reversible Liquid-like Elastic Deformation with Topological Phase Transition in Na$_3$Bi"

## See other formats

arXiv:1505.03936vl [cond-mat.mtrl-sci] 15 May 2015 Combined Fast Reversible Liquid-like Elastic Deformation with Topological Phase Transition in Na 3 Bi Xiyue Cheng, 1 Ronghan Li, 1 Dianzhong Li, 1 Yiyi Li, 1 and Xing-Qiu Chen' ll 1 Shenyang National Laboratory for Materials Science, Institute of Metal Research, Chinese Academy of Science, 110016 Shenyang, Liaoning, China (Dated: May 18, 2015) By means of first-principles calculations, we identified the structural phase transition of Na 3 Bi from hexag¬ onal ground state to cubic cF 16 phase above 0.8 GPa, in agreement with the experimental findings. Upon the releasing of pressure, cF 16 phase of Na 3 Bi is mechanically stable at ambient condition. The calculations re¬ vealed that the cF 16 phase is topological semimetal, in similarity to well-known HgTe and it even exhibits an unusually low C' modulus (only about 1.9 GPa) and a huge anisotropy, A" of as high as 11, the third highest value among all known cubic crystals in their elastic behaviors. These facts render cF16-type Na 3 Bi very soft with a liquid-like elastic deformation in the (110)< 110> slip system. Importantly, as accompanied with this deformation, Na 3 Bi shows a topological phase transition from a topological semimetal state at its strain-free cubic phase to a topological insulating state at its distorted phase. Because the C' elastic deformation almost costs no energy in a reversible and liquid-like soft manner, cF16-type Na 3 Bi would potentially provide a fast on/off switching way between topological insulator and topological semimetal, which would be beneficial to the quantum electronic devices for practical applications. PACS numbers: 71.70.Ej, 73.20.At, 62.20.Dc The topological material have been extensively studied and its number and types have increased dramatically over the past decades. As the first member of the topolog¬ ical family, the topological insulators (TIs) are materials with a bulk band gap but have protected metallic states on their edge/surface, originating from the combination of spin- orbit interactions and time-reversal symmetry (l]-0, includ¬ ing Bi 2 Se 3 0, 0 and the topological crystalline insulators (TCIs) such as Pbi_ A Sn A Te(Se), SnSe, and SnS|)6|-0. In an- glogy to TIs, the non-trivial topology of the band structure in semimetal gives rise to other new states of matters, topologi¬ cal Dirac semimetals (TDSs) and topological Weyl semimet¬ als (TWSs). For both TDSs and TWSs, the bulk conduction and valence bands touch only at discrete (Dirac) points and disperse linearly along all momentum directions, while their low-energy bulk excitations are described by the Dirac and Weyl equations, respectively. Compared to Dirac semimetal. inversion or time reversal symmetry must be broken in Weyl semimetal. For example, Na 3 Bi|i4H and Cd 3 As 2 0G0 were theoretically and experimentally demonstrated to be 3D TDSs protected by crystal symmetry and the noncentrosym- metric materials including TaAs, TaP, NbAs, and NbP jig] was recently predicted as natural TWSs. Certainly, the topo¬ logical concept can be also introduced into metals leading to topological metals (TMs). Generally, most of TMs can be manually achieved from TIs or semimetals, yet so far several compounds, such as HgTe IB 01 and NaBi 10 1. were found to be a native 3D TM. Although the flourishing topological material has inspired great interest, one has to admit the fact that some existed challenges are restricting its development, no matter from the aspects of scientific researches or techno¬ logical application. Besides two often facing challenges (small band gap and small quantized conductance of the edge state) 0000, another main challenge is the lacking of effective ways 0 to control the phase transition between topological materials and regular insulators or metals. Traditionally, the control of topological phase transitions was realized mainly by pressure, lattice strain, or chemical substitution. For instance, hP%- type Na 3 Bi undergoes a topological phase transition from a Dirac semimental to a trivial insulator upon doping Sb or P through varying the concentrations of Na 3 Bii_ A (Sb,P) A 0. However, via such a way it is extremely difficult to make the transition reversible once the material is fabricated. But, is there any way to rapidly and effectively switch on or off these phase transitions, rather than traditional means? This defi¬ nitely poses a challenging issue. As advanced functional materials, the elastic or mechani¬ cal properties of topological materials were often neglected and concealed under their wonderful and shining electronic properties. It is true that a vast majority of topological ma¬ terials only have ordinary elastic or mechanical properties which may just indicate the structural mechanical stability 1 2211 . Thus, it would be highly interesting to find a topologi¬ cal material with extinct elastic or mechanical properties. Re¬ cently, Na 3 Bi was proposed and demonstrated to be a 3D TDS I9l 4l2ll , and it is found that Na 3 Bi would undergo a structural phase transition of from its hexagonal ground state (lrP24) to a cubic cFl6 phase at 0.8 G Pa 12311 , in well agreement with the experimental findings 12442611. To one’s surprise, this cF 16- Na 3 Bi exhibits an unusually low C' modulus (only about 1.9 GPa) and a huge anisotropy. A" of as high as 11, the third high¬ est value among all known cubic crystals in its elastic prop¬ erties. These facts render cFl 6-type Na 3 Bi very soft with a liquid-like elastic deformation in the (110)< 110> slip system. Meanwhile, as the C' elastic deformation induced a topologi¬ cal phase transition from TS to TI and it almost cost no energy in a reversible and liquid-like elastic manner, this cubic Na 3 Bi 2 (a) hP 8 (Jual O Na2 O Na 1 • Bi L, c- kh a 60 120 160 200 Piassura (GPa) FIG. 1: (Color online) The DFT-derived pressure-dependent phase transitions of NajBi. (a) hP8 hexagonal phase, which was experimentally claimed to be stable at the ground state under ambient pressure, (b) hP24 hexagonal phase, which was theoretically found to be stable at the ground state under ambient pressure, (c) cF\6 cubic phase, which is more stable in energy above 0.8 GPa than both /?P8 and liP24 phases, and (d) oC16 orthorhombic phase, which is more stable above 118 GPa. (f,g,h,i) the derived phonon dispersions of hP8 (0 GPa), hP24 (0 GPa), cF 16 (0 GPa) and oC 16 (120 GPa), respectively, (e and j) the pressure-dependent enthalpies of hP 8, hP24, cF 16 and oC 16 phases to show the pressures of phase transitions. The red circles represent the Bi atoms while the Na atoms are shown in different color according to their Wyckoff sites. Here, yellow, blue and green circles represents Nal, Na2 and Na3 atoms, respectively. phase would potentially provide a fast on/ofF switching way between TI and TS. The structural optimization, electronic properties and elas¬ tic behaviors of the Na 3 Bi phases were calculated within the framework of density functional theory (DFT) 0 Bj using the Vienna ab initio Simulation Package (VASPl fe W with the projector augmented wave (PAW) method 113111 and gen¬ eralized gradient approximation (GGA) within the Perdew- Burke-Ernzerhof (PBE) exchange-correlation functional 13211 . The cutoff energy for the expansion of the wavefunction into plane waves was set at 350 eV and the tetrahedron method with Blochl corrections was ultilized. All the Brillouin zone integrations were performed on Monkhorst-Pack k-meshes and was sampled with a resolution of 2 7r x 0.07 A -1 , which showed excellent convergence of the energy differences and stress tensors. For hexagonal structures, the Brillouin zone integrations were performed on the T-centered symmetry. To check the dynamical stability, we further derived the phonon dispersion curves using the finite-displacement approach as implemented in the Phonopy code [33fl- The phonon frequen¬ cies are constructed from forces, resulting from displacements of certain atoms in a supercell containing typically 80-100 atoms for each Na 3 Bi phases, respectively. The electronic lo¬ calized function (ELF) 1341 - 136 1 was done using the grid-based algorithm with 100 x 100 x 100 grids. In addition, all crys¬ tal structures and ELF diagrams were generated using VESTA & The latest discoveries that native Na 3 Bi is a three dimen¬ sional (3D) Dirac semimetal represent a significant advance in topological Dirac materials lM-ll ill, because this material enables the study of a new type of quantum state. At ambi¬ ent condition, Na 3 Bi possesses bulk Dirac fermions in 3D at the Fermi level, which disperse linearly along all three mo¬ mentum directions [liiL lTTll . in contrast to the two-dimensional Dirac fermions present on the surfaces of 3D topological in¬ sulators and in graphene |381. Recently, the first-principles calculations revealed that the previously experimentally characterized /zP8 phase (Fig. m is unstable at the ground state due to the presence of a negative phonon branch around the high-symmetry X point, as illus¬ trated in Fig. |T}f). Instead, the calculation suggested that the real ground-state phase would be the lrP24 phase (Fig. [Q b)), which is indeed a distorted superlattice version of the /zP8 phase 0. Importantly, this stable hP24 structure exhibits a semblable 3D TDS feature as what the /zP8 phase presents & Furthermore, the previous experiments demonstrated that Na 3 Bi undergoes a phase transition above the pressures of 0.7-1.0 GPal24-26ll. In nice agreement with the experimental measurements, our first-principles calculations successfully reproduced this phase transition from the ground-state liP24 (or /zP8) phase to the cubic cFl6 phase (Fig. [He)) above 0.8 GPa [231], as shown in Fig. 0e). The calculations demon¬ strated that, even upon the releasing pressure, the cF 16 phase is still mechanically and dynamically stable at zero pressure, as elucidated by its phonon dispersions (Fig. |T]h)). This fact indicates that the cF 16-type phase is quenchable at ambient condition once it is synthesized above 0.8 GPa. Specifically, this cubic Na 3 Bi crystalizes in a BiF 3 -type structure (space group of Fm3m) with a- 7.550 A at ambient zero pressure. The Bi atoms lie at 4 b site (0.5, 0.5, 0.5) and the Na atoms occupy two inequivalent sites, 4 a (0, 0, 0) and 8c (3/4, 3/4, 3/4). The calculation also revealed that, by further increasing the pressure up to above 118 GPa, the cubic cF 16 phase can 3 be transformed into an orthorhombic 0 CI 6 phase (Fig. Ed) and Fig. Ej)X which is a wide-gap insulator 11231. Surprisingly, this cubic cF 16 phase exhibits several unusual elastic behaviors in its mechanical properties. At 0 GPa, the DFT-derived single crystal elastic constants are C n = 22.3 GPa, C \2 = 18.4 GPa and C 44 =21.9 GPa. Firstly, Cn is very close to C 44 , which gives nearly coinciding longitudinal and transverse acoustic phonons along the direction of < 001 >. This fact is reflected well by that two phonon branches along this I -X direction (Fig. Eh)) are nearly degenerated at small q values close to the F point. Secondly, it has been noted that the difference of C 44 -C 12 is only about 3.5 GPa, reveal¬ ing a very small deviation from the Cauchy relation with C 44 = C 12 , which measures the importance of the angular depen¬ dence in atomic forces as compared to the central force de¬ scription. Thirdly, cFlb-NaiBi has an rather low C' = r " 3 Cl ’ = 1.9 GPa. Note that C' represents the resistance to shear deformation by a shear stress applied across the ( 110 ) plane in the < 110 > direction and it indeed measures the rigidity against the volume-conserving tetragonal deformation. We further derived its universal anisotropic ratio to be as high as A u = 11.07, according to the recentl y pr oposed definition of the universal elastic anisotropy index [39], highlighting a huge anisotropy in its elastic properties. Furthermore, we plot the C' versus A 1 ' for a variety of cubic crystals (330 compounds, elastic data from the supplementary table of ref. 140]) as il¬ lustrated in Fig. E a and b). Among those compounds, we have observed from the plot that the cubic Na 3 Bi almost ex¬ hibits the lowest C' value and the third highest anisotropic ratio A" (which are only lower than two known superelastic compounds of CuZn and CuAuZn->) lf4l[ El . These unique features render cFlb-Na^Bi special: it is very soft, showing a liquid-like behaviour along the C' deformation corresponding to the shear slip ( 110 )< 110 > system within the elastic regime. This fact has also been evidenced by Fig. E C X from which it can be seen that the C' deformation almost costs no en¬ ergy, based on the deformation energy-versus-strain relation. In particular, it has been also found that, the C' and universal anisotropic ratio A u show unusual trends with increasing pres¬ sure, as evidenced in Fig. Ed)- In the first range below 1 GPa, both C' and A 1 ' remain almost constant and, in the pressure range from 1 GPa to 4 GPa C' rapidly increases from 2 GPa to 6 GPa whereas A“ shows an apparent drop from about 11 to 6 . Above about 4 GPa, they again almost become constant. In particular, the liquid-like and soft elastic behavior of Na 3 Bi remains robust in a certain pressure range at least from 0 GPa to 1 GPa. Importantly, the huge anisotropic ratio of rFI 6 -type Na 3 Bi has been further reflected by the stereographic projections at ambient pressure in Fig. E a an d c), which are used in order to show the particular crystal orientation and the cor¬ responding Young’s and shear moduli more clearly, while the 3D representation surfaces (Fig. Eb and d)) visualized the crystallographic-orientation dependent Young’s and shear moduli which, however, presents the strong and weak direc¬ tions in one crystal structure directly [[43;]. Note that the neg¬ 500 400 300 U diamond , . (a) 20 CuAuZn,° □ - 15 _□ Zener anisotropy, A CuZn 0 Universal anisotropy, A u t<> J 10 a Na s Bi / < ifl* 5 EL * CuAuZn 2 0 ^ Na 3 Bi CuZn \ , 0 « « 1 0 , ^1 , 0- o. (b) nun 1 r( 1 1 1 1 1 1 0 5 10 15 20 Universal anisotropy, A u 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 Chung anisotropy, A c FIG. 2: (Color online) The liquid-like elastic behavior of Na 3 Bi. (a) C' versus universal anisotropic ratio A", showing that Na 3 Bi has an extremely low C' and a highly large A“ among a variety of cubic crystals, (b) anisotropic ratio A (A“) versus Chung anisotropy A c for a variety of cubic crystals, (c) the DFT-derived deformation energy as a functional of strain according to the definition of elastic energies of bulk modulus, C 44 , and C' at 0 GPa, (d) the pressure-dependent universal anisotropic ratio A “ and C'. ative sign only denotes the negative direction corresponding to the positive one. In its stereographic projections, each blue point denotes one crystal orientation and the lines with differ¬ ent colour and density represent the counter lines of Young’s or shear moduli. It can be seen that the elastic properties un¬ der different crystal orientations highly vary from each other. For Young’s modulus, the lowest and highest orientations are [001] (£[ 001 ] = 6.42 GPa) and [111] (£ [m] = 53.9 GPa), re¬ spectively. In contrast, its shear modulus exhibits an opposite trend. The [001] orientation has the highest shear modulus of G[ooi] = 24.5 GPa and the [111] one with the lowest value of G[in] =3.17 GPa. There is no doubt that their opposite trends between both [ 001 ] and [ 111 ] directions substantially con¬ tribute to the huge anisotropy of this cubic crystal. Mechan¬ ically, the huge anisotropy can be attributed to the electronic structures and atomic arrangements. As illustrated in the elec¬ tronic localization function (ELF) images (Fig. Ee,f,g)), the (001) plane consists of Bi and Na atoms within the ionic¬ bonding manner (Fig. E e )X whereas the ( 111 ) plane is indeed comprised by the densest metallic Na atoms with pure metal¬ lic bonds (Fig. Eg))- It is clear that the dense ( 111 ) plane is hard to be compressed, but it can be easily sheared under deformation due to the pure metallic bonds. However, for the (001) plane, all bonds are ionic Na-Bi bonds with low atomic density. This interprets as to why the resistance to the com¬ pression of the ( 001 ) plane is low, whereas its shear modulus is much higher. Similar interpretation also holds for the (110) 4 plane (Fig. [3jf)), from which we can see that both Na-Na metallic bonds and Na-Bi ionic bonds co-exist ( c.f and the metallic bonds form a pure channel along the < 110 > direc¬ tion, Fig. [3j), leading to the intermediate Young’s and shear moduli. FIG. 3: (Color online) Huge elastic anisotropy of Na^Bi. (a and c) The stereographic projections of the Young’s and shear mod¬ uli, respectively, (b and d) The representation surfaces for the crystallographic-orientation depended Young’s and shear moduli, re¬ spectively. (e, f, and g) The electron localization function (ELF) isosurface maps for the (001), (110) and (111) planes, respectively. Na 3 Bi naturally processes the inverted band ordering, which is a prerequisite for the occurrence of topological electronic properties. Interestingly, when the SOC is switch on, the in¬ verted band ordering is further enhanced. The Na-3,v state (I 'jt in Fig. He)) is significantly reduced in energy by increasing the energy difference up to 0.74 eV (Fig. 4(c)). Therefore, in its strain-free cubic case, cf lO-NaiBi is a TS at 0 GPa, no matter whether the spin-orbit coupling effect (SOC) is con¬ sidered. Note that this topological semimetal feature of cF 16- NagBi remains unchanged up to 3.65 GPa before it transforms to a regular insulator 112311 . Certainly, it needs to be empha¬ sized that the band inversion at F is caused by the crystal field effect via the protection of lattice symmetry, rather than the chemical doping, pressure or strains. This feature is similar to the famous case of HgTe with a nearly zero direct band gap at F point in the BZ and an analogous inverted band ordering & 1.0 0.5 > 0 U- LU I LU 0 - 0.5 - 1.0 1.0 > <D li- lu ■ LU 0.5 0 - 0.5 - 1.0 rxL r wtxl r w In addition to its unusual elastic properties, the cubic cF 16 phase of Na 3 Bi exhibits another type of attractive electronic properties, in different from the ground-state hP24 phase with unique feature of 3D TDS. From Fig. Ha), it has been seen that the Fermi level exactly crosses the threefold degenerated Bi- P( X ,y,z ) states (r^ at I ’ for cF 16-type Na 3 Bi at 0 GPa. If the spin-orbit coupling (SOC) effect is included, this threefold degenerated F 4 p-states would further split into the doubly de¬ generated |/7 3 ) and nondegenerated |pi) states, as shown in Fig. He). With this electronic feature, it is clear that in the strain-free cubic case Na 3 Bi is a semimetal because its con¬ duction and valence bands touch only at the F point. Another typical feature is that even in the case without SOC, the Na- 3s states (T|) are energetically lower by about 0.4 eV than the Bi- 6 p states ^ 4 ) at the F point (Fig. 4(a)). More impor¬ tantly, T 1 and 1 4 exhibits the opposite parities. However, it is worth noting that no any other band inversions present for the rest high symmetric points in the BZ. These facts imply that FIG. 4: (Color online) Topological phase transition from topological semimetal to topological insulator under the C' tetragonal deforma¬ tion. Band structures for cubic cF16-type NajBi (a and c) and for its distorted structure via the C deformation with a 4% strain (b and d). a and b: no SOC inclusion; c and d: with SOC inclusion. The solid circles denote the j-like states. The doubly degenerated |p 3 ) states crossing the Fermi level in Fig. He) can be certainly broken by the engineering of strains. Importantly, via the liquid-like C' soft deformation of the ( 110 )< 110 > slip system, our calculation demonstrated that the cF lh-NatBi can be indeed transformed into a TI from its intrinsic feature of a TS. From Fig. Hb and d) the band structure at around the Fermi level depends sensitively on the lattice distortion along the soft ( 110 )< 110 > tetragonal elas¬ tic deformation. Without the SOC inclusion any slight lattice distortion along the C' deformation breaks up the triply de¬ generated Bi-p( A ViZ) states ^ 4 ) at F into doubly degenerated Bi-p VJ and nondegenerated Bi-p, states (see Fig. Hb) with 5 a 4% distortion). Furthermore, with the SOC inclusion, the doubly degenerated Bi -p xy states at the Fermi level is further broken to open a band gap, as evidenced in Fig. |4]d). In the meanwhile, such a distortion does not affect its band inversion. From Fig. [4] d), the nondegenerated Na-3,v state (T*) is ener¬ getically lower by about 0.76 eV than the Bi -p state (TL,) at the T point. It needs to be emphasized that for cubic structure the product of the parities of the Bloch wave function for the occupied bands at time-reversal invariant momenta (TRIMs) of three equivalent X (0.5,0,0), three equivalent M (0.5,0.5,0), and single R (0.5,0.5,0.5) has been derived to be -1, whereas at r it is +1. Via the C' deformation, although that these three equivalent momenta of X (or M) can be lifted, the product of the parities for all occupied bands at the TRIMs of 1', R, X and M still remains -1. Therefore, these results evidence that the distorted cubic Na 3 Bi by the strain engineering along the C' deformation is a three-dimensional topological insula¬ tor. Importantly, because of the reversible liquid-like elastic behavior along the C' deformation which almost cost no en¬ ergy, the cF 16-type Na 3 Bi becomes a highly attractive com¬ pound. Through such a way, one would easily realize rapid on/off switching from a TS state in its cubic phase to a TI state in its distorted lattice. In summary, we have demonstrated that a soft liquid-like elastic control of topological phase transition in cubic cF 16- type Na 3 Bi can be achieved through the C' tetragonal defor¬ mation along the (110)< 110> slip system, highlighting a real¬ ization of a fast and reversible on/off switching way between topological semimetal and topological insulator. This would be beneficial to the quantum electronic devices for practical applications. Acknowledgments This work was supported by the ’’Hun¬ dred Talents Project” of the Chinese Academy of Sciences and by the Key Research Program of Chinese Academy of Sciences (Grant No. KGZD-EW-T06) and NSFC of China (Grand Numbers: 51074151). The computational resource is using the local HPC cluster of the Materials Process Modeling Division in the IMR as well as the National Supercomputing Center in Tianjin (TH-1A system). * Corresponding author:xingqiu.chen@imr.ac.cn [1] M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). [2] X.-L. Qi and S.-C. Zhang. Rev. Mod. Phys. 83, 1057 (2011). [3] B. Yan and S. C. Zhang, Rep. Prog. Phys. 75, 096501 (2012). [4] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang and S.-C. Zhang, Nat. Phys., 5, 438 (2009). [5] Y. Xia, D. Qian, D. Hsieh. L. Wray, A. Pal, H. Lin. A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava and M. Z. Hasan, Nat. Phys. 5, 398(2009) [6] T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil and L. Fu, Nat. Commun. 3, 982 (2013). [7] S. Y. Xu, C. Liu, N. Alidoust, M. Neupane, D. Qian, I. Belopol- ski, J. D. Denlinger, Y. J. Wang, H. Lin, L. A. Wray, G. Landolt, B. Slomski, J. H. Dil, A. Marcinkova, E. Morosan, Q. Gibson, R. Sankar, F. C. Chou, R. J. Cava, A. Bansil and M. Z. Hasan, Nat. Commun. 3, 1192 (2013). [ 8 ] Y. Sun, Z.C. Zhong, T. Shirakawa, C. Franchini, D. Z. Li, Y.Y. Li, S. Yunoki, and X.-Q. Chen, Phys. Rev. B 88 . 235122 (2013). [9] Z. Wang, Y. Sun, X.-Q. Chen, C. Franchini, G. Xu, H. Weng, X. Dai, and Z. Fang, Phys. Rev. B 85, 205101 (2012). [10] Z. K. Liu, B. Zhou, Y. Zhang, Z. J. Wang, H. M. Weng, D. Prab- hakaran, S.-K. Mo, Z. X. Shen, Z. Fang, X. Dai, Z. Hussain, and Y. L. Chen, Science 343, 864 (2014). [11] X. Y. Cheng, R. H. Li, X.-Q. Chen, D. Z. Li, and Y. Y. Li, Phys. Rev. B 89, 245201 (2014). [12] S. Y. Xu, C. Liu, S. K. Kushwaha, R. Sankar, J. W. Krizan, I. Belopolski, M. Neupane, G. Bian, N. Alidoust, T.-R. Chang. H.-T. Jeng, C.-Y. Huang, W.-F. Tsai, H. Lin, P. P. Shibayev, E- C. Chou, R. J. Cava and M. Z. Hasan, Science 347, 294 (2015). [13] Z. Wang, H. Weng, Q. Wu, X. Dai, and Z. Fang, Phys. Rev. B 88 , 125427 (2013). [14] M. N. Ali, Q. Gibson, S. Jeon, B. B. Zhou, A. Yazdani and R. J. Cava, Inorg. Chem. 53, 4062 (2014). [15] Z. K. Liu, J. Jiang, B. Zhou, Z. J. Wang, Y. Zhang, H. M. Weng, D. Prabhakaran, S. K. Mo, H. Peng, P. Dudin, T. Kim, M. Hoesch, Z. Fang, X. Dai, Z. X. Shen, D. L. Feng, Z. Hussain and Y. L. Chen, Nat. Mater. 13, 677 (2014). [16] H. Weng, C. Fang, Z. Fang, A. Bernevig and X. Dai, Phys. Rev. X 5, 011029 (2015). [17] B. A. Bernevig, T. L. Hughes and S. C. Zhang, Science 314, 1757 (2006). [18] M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, L. Wiedmann, X. L. Qi and S. C. Zhang, Science, 318, 766 (2007). [19] X. F. Qian, J. W. Liu, L. Fu and J. Li, Science, 346, 1344 (2014). [20] R. H. Li, X. Y. Cheng, Q. Xie, Y. Sun, D. Z. Li, Y. Y. Li and X.-Q. Chen, Sci. Rep., 5, 8446 (2015). [21] A. Narayan, D. Disante, S. Picozzi, and S. Sanvito, Phys. Rev. Lett., 113, 256403 (2014). [22] B. Sa, J. Zhou, R. Ahuja and Z. Sun, Comput. Mater. Sci. 82, 66 (2014). [23] X. Y. Cheng. R. H. Li, D. Z. Li, Y. Y. Li and X.-Q.Chen, Phys. Chem. Chem. Phys. 17, 6933 (2015). [24] M. E. Leonova, S. A. Kulinich, and L. G. Sevast’yanova, Exp. Geosci. 7, 55 (1998). [25] M.E. Leonova, I. K. Bdikin, S.A. Kulinich, O.K. Gulish, L.G. Sevast’yanova, K.P. Burdina, Inorg. Mater. 39, 266 (2003). [26] S. A. Kulinich. M. E. Leonova, and L. G. Sevast’yanova, Zh. Obshch. Khim. 69, 681 (1999). [27] P. Hohenberg, Phys. Rev. B 136, 864 (1964). [28] W. Kohn and L. J. Sham, Phys. Rev. A 140, 1133 (1965). [29] G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). [30] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 (1996). [31] P. E. Blochl, Phys. Rev. B 50, 17953 (1994). [32] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). [33] A. Togo, F. Oba, I. Tanaka, Phys. Rev. B 78, 134106 (2008). [34] A. D. Becke and K. E. Edgecombe, J Chem. Phys. 92, 5397 (1990). [35] B. Silvi and A. Savin, Nature 371, 683 (1994). [36] A. Savin, R. Nesper, S. Wengert, and T. F. Fassler, Angew. Chem. Int. Ed. Engl. 36, 1808 (1977). [37] K. Momma and F. Izumi, J. Appl. Crystallogr. 44, 1272 (2011). [38] A. H. Castro Neto, N. M. R. Peres, K. S. Novoselov and A. K. Geim, Rev. Mod. Phys., 81, 109 (2009). [39] S. I. Ranganathan and M. Ostojia-Starzewski, Phys. Rev. Lett., 101, 055504 (2008). [40] H. Niu, X.-Q. Chen, P. Liu, W. Xing, X. Cheng, D. Li and Y. Li, Sci. Rep. 2, 718 (2012). 6 [41] G. M. McManus, Phys. Rev. 129, 2004 (1963). [43] X. Y. Cheng, X.-Q. Chen, D. Z. Li and Y. Y. Li, Acta Crystal- [42] M. Meyers and K. K. Chawla, in Mechanical Behavior of Ma- logr. C Struct. Chem. 70, 85 (2014). terials, Cambridge University Press, (2009).