# Hybrid Dirac Semimetal in CaAgBi Materials Family

###### Abstract

Based on their formation mechanisms, Dirac points in three-dimensional systems can be classified as accidental or essential. The former can be further distinguished into type-I and type-II, depending on whether the Dirac cone spectrum is completely tipped over along certain direction. Here, we predict the coexistence of all three kinds of Dirac points in the low-energy band structure of CaAgBi-family materials with a stuffed Wurtzite structure. Two pairs of accidental Dirac points reside on the rotational axis, with one pair being type-I and the other pair type-II; while another essential Dirac point is pinned at the high symmetry point on the Brillouin zone boundary. Due to broken inversion symmetry, the band degeneracy around accidental Dirac points is completely lifted except along the rotational axis, which may enable the splitting of chiral carriers at a ballistic p-n junction with a double negative refraction effect. We clarify their symmetry protections, and find both the Dirac-cone and Fermi arc topological surface states.

Topological semimetals have been attracting tremendous interest in current research Chiu et al. (2016); Dai (2016); Zhao and Wang (2013), partly because they offer a convenient platform to explore the intriguing physics of high-energy elementary particles. For example, Weyl semimetals possess linearly-dispersing two-fold-degenerate Weyl points close to Fermi energy Wan et al. (2011); Murakami (2007); Burkov and Balents (2011); Weng et al. (2015); Huang et al. (2015); Lv et al. (2015); Xu et al. (2015a). Each Weyl point has a definite chirality of , around which the quasiparticle excitations mimic the relativistic Weyl fermions Volovik (2003); Nielsen and Ninomiya (1983). Two Weyl points with opposite chirality would be unstable towards gap-opening when they meet at the same -point, unless there exists additional symmetry protection, such as the case in so-called Dirac semimetals Young et al. (2012), in which stable Dirac points with four-fold-degeneracy make it possible to simulate massless Dirac fermions.

Dirac points can be classified as accidental or essential Yang and Nagaosa (2014) (see Fig. 1). Accidental Dirac points require band inversion (which is in a sense accidental) and are stablized by certain symmorphic crystalline symmetries such as rotation, so such Dirac points typically reside on high-symmetry lines. Examples include the experimentally confirmed Dirac semimetals NaBi Wang et al. (2012); Liu et al. (2014a) and CdAs Wang et al. (2013); Liu et al. (2014b); Neupane et al. (2014); Borisenko et al. (2014). Essential Dirac points do not need band inversion, and their presence is solely determined by certain nonsymmorphic symmtries at high-symmetry points on the boundary of the Brillouzin zone (BZ). Examples include the first few proposals such as -cristobalite BiO Young et al. (2012) and several Bi-containing distorted spinels Steinberg et al. (2014). Accidental Dirac points can be removed by reverting the band ordering without changing the symmetry, whereas essential Dirac points cannot. In addition, since accidental Dirac points are located at -points with lower symmetry (compared with essential ones), the dispersions around them are less constrained. It is possible to have the Dirac cone completely tipped over along certain direction Soluyanov et al. (2015); Xu et al. (2015b), realizing a so-called type-II Dirac point Chang et al. (), which has been identified in VAl Chang et al. (), PdTe Fei et al. (); Xiao et al. (); Noh et al. (), PtSe family Zhang et al. (); Huang et al. (2016); Yan et al. () and others Guo et al. (2017); Le et al. (); Wang et al. (). Fascinating yet distinct physics have been proposed for type-I and type-II points, and it has been theoretically argued in the context of Weyl points that when both types coexist in a single hybrid material, even more interesting effects could appear Li et al. (2016); Kong et al. (2017).

In this work, we show that all the above-mentioned three kinds of Dirac points are simultaneously present in the low-energy band structure of CaAgBi-family materials. Two pairs of accidental Dirac points are realized on the -axis and are protected by the point group symmetry, with one pair being type-I and the other pair type-II. Meanwhile, another single essential Dirac point occurs at the A point on the boundary of BZ, as dictated by the nonsymmorphic space group symmetry. The inversion symmetry, preserved in most Dirac semimetals proposed so far, is broken in the CaAgBi-family materials. Consequently, the dispersions around the Dirac points exhibit distinct features. While each band is two-fold degenerate around the essential Dirac point, the band degeneracy around the accidental Dirac points is lifted due to the broken inversion symmetry (except along the rotational axis), making it possible to split the chiral carriers at a ballistic p-n junction with a double negative refraction effect. We clarify the symmetry protection mechanisms, and find interesting surface-orientation-dependent topological surface states including both the Dirac-cone surface bands and the surface Fermi arcs. This class of hybrid Dirac materials offers a promising platform to explore the intriguing physics of Dirac fermions.

The CaAgBi-family materials take a stuffed Wurtzite-type structure, with space group No. 186 () Sun et al. (2007), as shown in Fig. 2. The structure may be viewed as a zinc-blende structure compressed along the [111]-direction. Since members of this family show similar low-energy band features, we shall focus on CaAgBi as a representative in the following discussion. It is found that the low-energy states are mainly from the orbitals of Ag and Bi atoms, which lie in the (110) mirror plane. Notably, the structure does not preserve an inversion center. Besides the (110) mirror plane (here we take [110] as -direction), two other important crystalline symmetries are the three-fold rotation and the two-fold screw rotation .

We performed first-principles calculations based on the density functional theory (DFT). The calculation details are presented in the Supplemental Material SI (). Figure 3(a) shows the plot of calculated band structure of CaAgBi with spin-orbit coupling (SOC) included. One observes that the material is a semimetal, and the bands around Fermi level are dominated by the Ag- and Bi- orbitals. Normally, Ag- orbitals have a higher energy than Bi- orbitals. However, around -point, one observes a band inversion feature with Ag- lower than Bi-. Importantly, along the -axis (-A), the three low-energy bands cross pairwisely and linearly at three discrete points, which we label as , , and , as shown in Fig. 3(b). and are on the rotational axis, which are accidental crossing points originated from the band inversion at . In contrast, is pinned at A-point on the BZ boundary, which is an essential crossing point dictated by the space group symmetry. One notes that each of the three crossing bands are two-fold degenerate (with SOC) on -A, which is caused by the non-commutativity between (or ) and symmetries along this path. Thus each crossing point () is of four-fold degeneracy, corresponding to a Dirac point.

Let’s first consider the two accidental Dirac points and . On -A path, the point group symmetry is of , and the three doubly-degenerate crossing bands are found to form three distinct two-dimensional irreducible representations , and of the double group (see Fig. 3(b)). Therefore, these bands cannot hybridize along this line, and their crossing points are protected and of linear type. The symmetry protection is lost for -points deviating from the rotational axis, hence each crossing should be an isolated point rather than a line. Another important observation is that while the slopes of the two crossing bands at have opposite signs, the two bands at have the same sign. Thus and are type-I and type-II Dirac points, respectively.

The dispersions around and in the - plane also show interesting features. As shown in Fig. 3(c,d), the bands completely split in directions different from the -axis, which can be understood by noticing that the two-fold band degeneracy on -A cannot be maintained at generic -points in the - plane because such points are not invariant under the rotational symmetries and . This is in contrast to most proposed Dirac semimetals, where the two-fold band degeneracy is maintained at any -point due to the presence of inversion symmetry (in combination with time reversal symmetry).

To characterize the low-energy Dirac fermions, we construct the effective models around each Dirac point, subjected to the symmetry constraints Bradley and Cracknell (1972). For point , it is at the intersection between band and band on -A. Using the four-states with such symmetries at as basis, the Hamiltonian around up to linear order in takes the form

(1) |

where the wave-vector and the energy are measured from , , the model parameters , , are real, and is complex. The first term in Eq. (1) represents a tilt of spectrum. Figure 3(c) shows the fitting of this model to the DFT band structure. Particularly, we find that , such that the point is of type-I with dispersion along -direction. Similarly, we can find the effective model around point , which is at the intersection between and bands. We find that the obtained Hamiltonian has the same form as , and only the values of the parameters are different. Most importantly, now . The dispersion along -direction is dominated by the tilt term, so that is a type-II Dirac point (see Fig. 3(d)).

Next, we come to the essential Dirac point at A. Its essential character can be argued as follows. The bands along -A can be chosen as eigenstates. Since (the minus sign is from spin rotation, and is measured in unit of ), the eigenvalues are given by

(2) |

Meanwhile, we have , so the two states and form a degenerate pair with opposite eigenvalues, where denotes an eigenstate of with eigenvalue . This explicitly demonstrates the two-fold band degeneracy along -A, as we claimed before. Note that A is invariant under time reversal operation (), such that at A, each state has a Kramers degenerate partner . Since at point A (), shares the same eigenvalue as . Consequently, the four states at A are linearly independent and must be degenerate with the same energy. Thus the nonsymmorphic space group symmetry necessitates the four-fold degeneracy at A-point. The effective model up to linear order in expanded around is obtained as

(3) |

Here and energy are measured from , and are the Pauli matrices, and model parameters and correspond to the Fermi velocities. Figure 3(e) shows the fitting of DFT band structure using model (3). Compared with model (1) for the two accidental Dirac points, model (3) has less free parameters due to the additional time reversal symmetry at A. Consequently, the dispersion around is different from the other two in that: (i) there is no energy tilt term (hence it must be type-I); and (ii) the bands in plane are still two-fold degenerate without splitting.

We further clarify the distinction between the accidental Dirac points (, ) and the essential Dirac point (). The realization of accidental Dirac points requires band inversion. In the current system, this can be characterized by two indices and . Here , with , where ( or , and , A) denotes the energy of the band at point . Evidently, corresponds to the nontrivial case with band inversion, leading to the Dirac point. Similarly, is defined with replaced by , and characterizes the band inversion between and bands for the formation of . The nontriviality of () is accidental, depending on the material. For example, in CaAgAs, which is isostructural to CaAgBi, we find that (see Fig. 4), indicating that the band inversion between and is removed. Therefore, point is eliminated. Meanwhile, remains , so is kept, but it may ultimately be eliminated by a symmetry-preserving strain. In contrast, is invariably present in all these cases, since it is only determined by symmetry.

Band topology and intrinsic anisotropy lead to interesting surface spectra of CaAgBi-family materials, with distinct topological surface states on different crystal surfaces. Here the Dirac points are all residing on -axis. They are projected onto the same image point on (001) surface. Since the band is inverted at -point in the bulk, we expect to have Dirac-cone type surface states similar to the topological insulator case. In Fig. 5(a), we indeed find such surface states buried in the bulk valence bands. Figure 5(b) shows the constant energy slice above the surface Dirac point, at which the surface states exhibit a left-handed spin helicity, similar to the pattern for strong topological insulators Sheng et al. (2014, ). Meanwhile, the Dirac points get projected to different image points on the side surfaces. For (010) surface, we find Fermi arcs connecting a pair of projected Dirac points, as shown in Fig. 5(c,d) for CaAgAs (for CaAgBi, the arcs are buried in the bulk bands SI ()). These Fermi arcs are in fact protected by a bulk invariant. Consider the -invariant planes and . The two planes are fully gapped hence each has a well-defined 2D invariant. Using the Wilson loop method Yu et al. (2011); Soluyanov and Vanderbilt (2011); Wu et al. (), we find that for plane, whereas for plane. Thus a Kramers pair of surface states must exist on the path for the side surfaces, ensuring the presence of a pair of surface Fermi arcs.

As mentioned, there exist a large family of ternary compounds with stuffed Wurtzite structure. Many members of this family share the same band features as discussed above (see Supplemental Material SI ()), exhibiting a hybrid Dirac semimetal phase with multiple Dirac points. Since the different kind of Dirac points are at different energies, it is possible to control the type of Dirac fermions in a single material by tuning the doping level, probing each kind separately. Alternatively, one may choose the proper material from the family which has the desired Dirac point most close to the Fermi level.

A salient feature of these materials is the broken inversion symmetry, which makes the Dirac bands completely split except along the -axis. This could lead to a new double negative refraction effect at a p-n junction of such material. Consider the case when the Fermi level is most close to the point. Consider ballistic transport across a p-n junction along -direction, and we focus on the carrier dynamics in the - plane (the transverse momenta and are conserved). As illustrated in Fig. 6(a), with a constant , there are two Fermi circles around the Dirac point, due to the band splitting. Carriers on each Fermi circle is chiral and undergo refractions at the junction interface, similar to that in graphene Cheianov et al. (2007). For p-n junction, this leads to negative refraction effect (Fig. 6(b)), and for incident state with two transmitted states, there are two negative refraction paths (see Fig. 6(c)). In addition, by controlling the bias voltage and the incident angle of an incident electron beam, one can have transmission of carriers from only one of the Fermi circle (see Fig. 6(d)), whereas carriers from the other circle are totally reflected. This could be useful in electron optics to spatially separate the chiral carriers.

In conclusion, we have demonstrated a novel hybrid Dirac semimetal phase in CaAgBi-family materials. The low-energy band structure features three kinds of Dirac points: two pairs of accidental Dirac points with type-I and type-II dispersions, and a single essential Dirac point dictated by nonsymmorphic symmetry. We clarify their symmetry protections and characterize their individual low-energy Dirac fermions. Distinct from centrosymmetric Dirac semimetals, the Dirac bands here completely split due to the lack of inversion center. We show that this feature leads to interesting transport features at a p-n junction. We also find that different crystalline surfaces possess different topological surface states. The bulk Dirac points and the surface states can be experimentally probed by the angle-resolved photoemission spectroscopy (ARPES). Our findings make it possible to study multiple types of Dirac fermions and their interplays within a single material system.

## References

- Chiu et al. (2016) C.-K. Chiu, J. C. Y. Teo, A. P. Schnyder, and S. Ryu, Rev. Mod. Phys. 88, 035005 (2016).
- Dai (2016) X. Dai, Nature Physics 12, 727 (2016).
- Zhao and Wang (2013) Y. X. Zhao and Z. D. Wang, Phys. Rev. Lett. 110, 240404 (2013).
- Wan et al. (2011) X. Wan, A. M. Turner, A. Vishwanath, and S. Y. Savrasov, Phys. Rev. B 83, 205101 (2011).
- Murakami (2007) S. Murakami, New Journal of Physics 9, 356 (2007).
- Burkov and Balents (2011) A. A. Burkov and L. Balents, Phys. Rev. Lett. 107, 127205 (2011).
- Weng et al. (2015) H. Weng, C. Fang, Z. Fang, B. A. Bernevig, and X. Dai, Phys. Rev. X 5, 011029 (2015).
- Huang et al. (2015) S.-M. Huang, S.-Y. Xu, I. Belopolski, C.-C. Lee, G. Chang, B. Wang, N. Alidoust, G. Bian, M. Neupane, C. Zhang, S. Jia, A. Bansil, H. Lin, and M. Z. Hasan, Nature Communications 6, 7373 (2015).
- Lv et al. (2015) B. Q. Lv, H. M. Weng, B. B. Fu, X. P. Wang, H. Miao, J. Ma, P. Richard, X. C. Huang, L. X. Zhao, G. F. Chen, Z. Fang, X. Dai, T. Qian, and H. Ding, Phys. Rev. X 5, 031013 (2015).
- Xu et al. (2015a) S.-Y. Xu, I. Belopolski, N. Alidoust, M. Neupane, G. Bian, C. Zhang, R. Sankar, G. Chang, Z. Yuan, C.-C. Lee, S.-M. Huang, H. Zheng, J. Ma, D. S. Sanchez, B. Wang, A. Bansil, F. Chou, P. P. Shibayev, H. Lin, S. Jia, and M. Z. Hasan, Science 349, 613 (2015a).
- Volovik (2003) G. E. Volovik, The Universe in a Helium Droplet (Clarendon Press, Oxford, 2003).
- Nielsen and Ninomiya (1983) H. Nielsen and M. Ninomiya, Physics Letters B 130, 389 (1983).
- Young et al. (2012) 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).
- Yang and Nagaosa (2014) B.-J. Yang and N. Nagaosa, Nature Communications 5, 4898 (2014).
- Wang et al. (2012) Z. Wang, Y. Sun, X.-Q. Chen, C. Franchini, G. Xu, H. Weng, X. Dai, and Z. Fang, Phys. Rev. B 85, 195320 (2012).
- Liu et al. (2014a) Z. K. Liu, B. Zhou, Y. Zhang, Z. J. Wang, H. M. Weng, D. Prabhakaran, S.-K. Mo, Z. X. Shen, Z. Fang, X. Dai, Z. Hussain, and Y. L. Chen, Science 343, 864 (2014a).
- Wang et al. (2013) Z. Wang, H. Weng, Q. Wu, X. Dai, and Z. Fang, Phys. Rev. B 88, 125427 (2013).
- Liu et al. (2014b) 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, Nature Materials 13, 677 (2014b).
- Neupane et al. (2014) M. Neupane, S.-Y. Xu, R. Sankar, N. Alidoust, G. Bian, C. Liu, I. Belopolski, T.-R. Chang, H.-T. Jeng, H. Lin, A. Bansil, F. Chou, and M. Z. Hasan, Nature Communications 5, 386 (2014).
- Borisenko et al. (2014) S. Borisenko, Q. Gibson, D. Evtushinsky, V. Zabolotnyy, B. Büchner, and R. J. Cava, Phys. Rev. Lett. 113, 027603 (2014).
- Steinberg et al. (2014) J. A. Steinberg, S. M. Young, S. Zaheer, C. L. Kane, E. J. Mele, and A. M. Rappe, Phys. Rev. Lett. 112, 036403 (2014).
- Soluyanov et al. (2015) A. A. Soluyanov, D. Gresch, Z. Wang, Q. Wu, M. Troyer, X. Dai, and B. A. Bernevig, Nature 527, 495 (2015).
- Xu et al. (2015b) Y. Xu, F. Zhang, and C. Zhang, Phys. Rev. Lett. 115, 265304 (2015b).
- (24) T.-R. Chang, S.-Y. Xu, D. S. Sanchez, S.-M. Huang, G. Chang, C.-H. Hsu, G. Bian, I. Belopolski, Z.-M. Yu, X. Xu, C. Xiang, S. A. Yang, T. Neupert, H.-T. Jeng, H. Lin, and M. Zahid Hasan, arXiv:1606.07555 .
- (25) F. Fei, X. Bo, R. Wang, B. Wu, J. Jiang, D. Fu, M. Gao, H. Zheng, Y. Chen, F. Song, X. Wan, B. Wang, X. Wang, and G. Wang, arXiv:1611.08112 .
- (26) R. C. Xiao, P. L. Gong, Q. S. Wu, W. J. Lu, M. J. Wei, J. Y. Li, H. Y. Lv, X. Luo, P. Tong, X. B. Zhu, and Y. P. Sun, arXiv:1705.05708 .
- (27) H.-J. Noh, J. Jeong, E.-J. Cho, K. Kim, B. I. Min, and B.-G. Park, arXiv:1612.06946 .
- (28) K. Zhang, M. Yan, H. Zhang, H. Huang, M. Arita, Z. Sun, W. Duan, Y. Wu, and S. Zhou, arXiv:1703.04242 .
- Huang et al. (2016) H. Huang, S. Zhou, and W. Duan, Phys. Rev. B 94, 121117 (2016).
- (30) M. Yan, H. Huang, K. Zhang, E. Wang, W. Yao, K. Deng, G. Wan, H. Zhang, M. Arita, H. Yang, Z. Sun, H. Yao, Y. Wu, S. Fan, W. Duan, and S. Zhou, arXiv:1607.03643 .
- Guo et al. (2017) P.-J. Guo, H.-C. Yang, K. Liu, and Z.-Y. Lu, Phys. Rev. B 95, 155112 (2017).
- (32) C. Le, S. Qin, X. Wu, X. Dai, P. Fu, and J. Hu, arXiv:1606.05042 .
- (33) H.-X. Wang, Y. Chen, Z. H. Hang, H.-Y. Kee, and J.-H. Jiang, arXiv:1703.09899 .
- Li et al. (2016) F.-Y. Li, X. Luo, X. Dai, Y. Yu, F. Zhang, and G. Chen, Phys. Rev. B 94, 121105 (2016).
- Kong et al. (2017) X. Kong, J. He, Y. Liang, and S.-P. Kou, Phys. Rev. A 95, 033629 (2017).
- Sun et al. (2007) Z.-M. Sun, J.-Y. Xie, D.-C. Pan, and J.-G. Mao, Journal of Alloys and Compounds 430, 71 (2007).
- (37) See Supplemental Material.
- Bradley and Cracknell (1972) C. J. Bradley and A. P. Cracknell, The Mathematical Theory of Symmetry in Solids (Clarendon Press, Oxford, 1972).
- Sheng et al. (2014) X.-L. Sheng, Z. Wang, R. Yu, H. Weng, Z. Fang, and X. Dai, Phys. Rev. B 90, 245308 (2014).
- (40) X.-L. Sheng, Z.-M. Yu, R. Yu, H. Weng, and S. A. Yang, arXiv:1703.09040 .
- Yu et al. (2011) R. Yu, X. L. Qi, A. Bernevig, Z. Fang, and X. Dai, Phys. Rev. B 84, 075119 (2011).
- Soluyanov and Vanderbilt (2011) A. A. Soluyanov and D. Vanderbilt, Phys. Rev. B 83, 035108 (2011).
- (43) Q. Wu, S. Zhang, H.-F. Song, M. Troyer, and A. A. Soluyanov, arXiv:1703.07789 .
- Cheianov et al. (2007) V. V. Cheianov, V. Fal’ko, and B. L. Altshuler, Science 315, 1252 (2007).