Giant magnetoresistance, three-dimensional Fermi surface and origin of resistivity plateau in YSb semimetal
Very strong magnetoresistance and a resistivity plateau impeding low temperature divergence due to insulating bulk are hallmarks of topological insulators and are also present in topological semimetals where the plateau is induced by magnetic field, when time-reversal symmetry (protecting surface states in topological insulators) is broken. Similar features were observed in a simple rock-salt-structure LaSb, leading to a suggestion of the possible non-trivial topology of 2D states in this compound. We show that its sister compound YSb is also characterized by giant magnetoresistance exceeding one thousand percent and low-temperature plateau of resistivity. We thus performed in-depth analysis of YSb Fermi surface by band calculations, magnetoresistance, and Shubnikov–de Haas effect measurements, which reveals only three-dimensional Fermi sheets. Kohler scaling applied to magnetoresistance data accounts very well for its low-temperature upturn behavior. The field-angle-dependent magnetoresistance demonstrates a 3D-scaling yielding effective mass anisotropy perfectly agreeing with electronic structure and quantum oscillations analysis, thus providing further support for 3D-Fermi surface scenario of magnetotransport, without necessity of invoking topologically non-trivial 2D states. We discuss data implying that analogous field-induced properties of LaSb can also be well understood in the framework of 3D multiband model.
1]Orest Pavlosiuk 1]Przemysław Swatek 1,*]Piotr Wiśniewski 1]Institute of Low Temperatures and Structure Research, Polish Academy of Sciences, PNr 1410, 50-950 Wrocław, Poland *]Corresponding author: firstname.lastname@example.org
Yttrium monoantimonide has mainly been used as a non-magnetic reference or as a ’solvent’ in monoantimonides of -electron-elements solid solutions with anomalous physical properties such as dense Kondo behavior and complex magnetic ground-states. Within that context it has been characterized as a metal by low-temperature specific heat measurements. Later Hayashi et al. have shown that the first-order phase transition from the NaCl-type to a CsCl-type crystal structure occurs in YSb at 26 GPa . That discovery induced numerous calculations of electronic structure of the compound, among them those by Tütüncü, Bagci and Srivastava, who directly compared electronic structure of YSb with that of LaSb . Results of those calculations were very similar for both compounds, revealing low densities of states at Fermi level and characteristic anti-crossings leading to band inversion at X-points of the Brillouin zone. LaSb has a simple NaCl-type structure without broken inversion symmetry, perfect linear band crossing or perfect electron–hole symmetry, yet it exhibits the exotic magnetotransport properties of complex-structure semimetals like TaAs, NbP (Weyl semimetals)[5, 6], CdAs (Dirac semimetal) and WTe (resonant compensated semimetal)[8, 9]. Recently Tafti et al. discovered in LaSb field-induced resistivity plateau at low temperatures up to K, ultrahigh mobility of carriers in the plateau region, quantum oscillations, and magnetoresistance () of nearly one million percent at 9 T . Their calculations, including spin-orbit coupling (SOC) effect, suggested that LaSb is a topological insulator with a 10 meV gap open near the X-point of the Brillouin zone. They also observed specific angular dependence of frequencies of quantum oscillations and ascribed them to two-dimensional Fermi surface (FS) possibly formed of topologically nontrivial states, and thus proposed LaSb as a model system for understanding the consequences of breaking time-reversal symmetry in topological semimetals.
Tafti et al. also invoked the opening of insulating gap (i.e. metal-insulator transition) as a source of the field-induced resistivity plateau in LaSb but it should be noted that in the case of WTe the existence of a magnetic-field-driven metal-insulator transition has been excluded by means of Kohler scaling analysis of magnetoresistance.
Motivated by these ambiguities in the interpretation of LaSb properties we decided to carry out a comprehensive characterization of magnetotransport properties of a sister compound YSb. We found that YSb displays physical properties in many aspects very similar to those of LaSb. Our results are in accord with those of other groups that appeared during preparation of our article[14, 15]. The interpretation proposed by Ghimire et al. and Yu et al. follows that presented by Tafti et al. for LaSb, implying the role of field-induced metal-insulator transition in YSb.
However, our analysis of magnetoresistance and Shubnikov–de Haas (SdH) effect provides strong support for 3D-Fermi surface scenario of magnetotransport, without necessity of invoking topologically non-trivial 2D states or metal-insulator transition in YSb.
Magnetoresistance and the origin of its plateau
Electrical resistivity () was measured on two samples (denoted as #1 and #2) and its dependence on temperature in zero field is plotted in Supplementary Figure 1(a). Shape of curves is typical for a metal. When measured in different applied fields exhibits universal plateau at temperatures 2–15 K, as shown for sample #1 in Supplementary Figure 1(b). Temperature range of this plateau is very similar to that reported for LaSb.
Magnetoresistance, , is plotted versus magnetic field, , in Figs. 1(a) and 1(b), for samples #1 and #2, respectively. Following the approach of Wang et al. we performed Kohler scaling analysis of magnetoresistance to test for the existence of a magnetic-field-driven metal-insulator transition in YSb. Figures 1(c) and 1(d) show very good Kohler scaling of our data, , yielding at 2.5 K exponents 1.64 and 1.74, respectively, very close for both samples, despite significant difference of their values. Efficiency of this scaling indicates that resistivity plateau is due purely to the magnetoresistance, but not to a field-induced metal-insulator transition.
In order to elucidate the dimensionality of FS we measured magnetoresistance of the sample #2 in fields applied at different angles to its surface, Fig. 2(a). Here denotes the field perpendicular to the sample surface and the current direction, whereas means that the field is parallel to the current. As shown in the inset to Fig. 2(a), at strongest field of 9 T follows . This is typical behavior for materials without magnetic anisotropy, but the change of expressed as anisotropic magnetoresistance, , has an outstanding -80% value.
Moreover, when field strength is scaled by a factor dependent on mass anisotropy and -angle, all data of Fig. 2(a) collapse on single curve, as shown in Fig. 2(b). Inset of Fig. 2(b) shows that values of plotted against field angle can be perfectly fitted with function, shown with red line. Such scaling has initially been proposed for anisotropic superconductors, and recently used to interpret behavior of WTe based on its 3D electronic nature. Parameter represents effective mass anisotropy of carriers mostly contributing to the magnetoresistance. We ascribe this behavior to a strongly anisotropic sheet of 3D-FS revealed by SdH effect data, as shown below.
Hall resistivity of sample #1h (cut from the same single crystal as #1) measured at several temperatures between 2.2 K and 300 K is shown in Fig. 3(a). Nonlinear indicates that at least two types of charge carriers are responsible for the Hall effect observed in YSb. The curves for temperatures from 2.2–15 K range are almost identical, which points to nearly constant carrier concentrations and mobilities, and coincides with the plateau of observed in the same range of . Changes of sign of observed at K, indicate conducting bands of both electrons and holes, at higher temperatures is positive in the whole range of magnetic field (T). Thus, the Hall contributions of holes and electrons nearly compensate, but both depend on temperature in different manner. Clear quantum oscillations are observed in in temperature range 2.2–15 K (cf. Figs. 3(a) and (b)).
Since , the off-diagonal component of conductivity tensor should be used for multiple-band analysis of Hall data. In this case simple Drude model can be used:
summing up conductivities of individual bands, with and denoting respectively concentration and mobility of carriers from the -th band. As Fig. 3(b) shows, for data collected at 300 K accounting for two bands yielded a good fit. On the other hand, fitting with two bands was insufficient for 2.2 K data, but addition of a contribution of another band with small concentration of more mobile holes brought a very satisfactory fit.
For LaSb Tafti et al. estimated uncompensated carrier concentration, , and the Hall mobility , using the relations (with being the zero-temperature limit of ) and . They obtained and .
For our YSb sample the fit with with multiple-band model (Eq. 1), shown in Fig. 3(b), yielded at 2.2 K the concentrations: , , and mobilities , . Third band necessary for that fit consists of more mobile holes [, ]. These results are consistent with band calculations presented below.
The fit to data collected at 300 K yielded: the concentrations: , , and mobilities , , so YSb has similar concentration of carriers but with significantly lower mobility than LaSb. This seems to be the main reason for its significantly smaller magnetoresistance, as it was well demonstrated for WTe.
Overall, Hall effect results are in perfect agreement with characteristics of Fermi surface presented in the next section.
Fermi surface analysis: electronic structure calculations and Shubnikov-de Haas effect
We performed the electronic structure calculations for YSb using a full potential all-electron local orbital code (FPLO) within GGA approximation. Figure 4 shows obtained energy band structure for YSb, with a few bands crossing Fermi level. Near X-point there is an anti-crossing present with a gap of , similar to those reported for lanthanum monopnictides, which led to the proposal of 2D topologically non-trivial states at the origin of extraordinary behavior of their magnetoresistance[19, 10]. Fermi surface (FS) resulting from our calculations is presented in Fig. 5: two electron sheets centered at X-points and three nested hole sheets centered at -point. Our calculations are in good agreement with those of Ref. .
Shubnikov–de Haas oscillations of resistivity are discernible for YSb at temperatures up to 15 K and in fields above 6 T, as seen in Figs. 1(a), 1(b) and 2(a). Since , we may safely assume that the conductivity and analyze directly oscillations of . Points in Fig. 6 represent resistivity of sample #2 measured at 2.5 K, after subtraction of smooth background, plotted versus inverse field (for T). Complex shape of dependence indicates that observed SdH oscillation has several components. Indeed fast Fourier transform (FFT) analysis reveals clearly six well separated frequencies (Fig. 7(a)).
We chose to fit with the multi-frequency Lifshitz-Kosevich function[20, 21, 22] (Eq. 2) because, as the maxima of total do not correspond to the maxima of particular components with different frequencies, it is inadequate to determine phases in a multicomponent SdH oscillation using the so-called Landau-level fan diagram (plot of the values of corresponding to the -th maximum in versus ).
where for each -th component: is the frequency, represents the Dingle factor and comprises the temperature reduction factor and the spin factor . Detailed form of Eq. 2, taking into account harmonic components, and description of , and is presented in Supplementary Material.
The fit including six components was of very good quality, as shown by blue line in Fig. 6. Obtained parameters are collected in Table 1. All frequencies converged almost exactly to values obtained from FFT analysis. We ascribe oscillations with frequencies of 720 and 1072 T to second and third harmonic of the strongest one with T. The phases of fundamental oscillations resulting from the fit are close to Onsager phase factor of expected for free electrons. Thus all components of SdH oscillation can be assigned to 3D-FS sheets predicted by band calculations (shown above in Fig. 5) and no Berry phase of was observed, which could reveal topologically non-trivial charge carriers.
|359||708||753||1066, 1082||1138, 1174||1411, 1452|
We performed FFT analysis for all data sets presented in Fig. 2(a), which allowed us to observe angular behavior of frequencies corresponding to all extreme cross-sections of Fermi sheets and compare them to those derived from our band structure calculations, as shown in Fig. 7(a). Three of observed frequencies were clearly changing upon rotation of the magnetic field: the principal, labeled as (T at ), and its harmonics, and . These frequencies are plotted versus in Fig. 7(b).
It became apparent from the shape of FS obtained from band calculations (cf. Fig. 5) that angular behavior of follows a cross-section area of a prolate ellipsoid (which well approximates the shape of electron sheet centered at X-point shown in Fig. 5) by the ((1- plane passing by X-point. When magnetic field is applied along direction the plane is perpendicular to -X line. When field is tilted from towards by the angle the cross-section increases as . This holds for a two-axial ellipsoid described by the equation: , with . Band structure shown in Fig. 4 yields (with estimated as average size of -sheet of FS along X-U and X-W lines, and as its size along -X line).
After rotation by the frequency meets the one denoted as , initially (i.e. at ) corresponding to the largest cross-section of the same ellipsoidal FS sheet. Two other observed frequencies, and do not change notably with , as expected for almost isotropic hole Fermi sheets centered at -point.
We plotted (for =3.6) in Fig. 7(b) as solid lines. The dependence, expected for two-dimensional FS, behaves similarly and is shown for comparison with dashed lines. The dependence was used by Tafti et al. as a hint of possible topologically-nontrivial states in LaSb. However, Fermi surface of that compound has already been well characterized by band calculations and angle-dependent de Haas–van Alphen measurements[11, 12], revealing FS very similar to the one we found in YSb, namely consisting of one ellipsoidal electron sheet centered at X-point and two isotropic hole pockets centered at -point. Hasegawa and Settai et al. assigned angular dependence of the principal de Haas–van Alphen frequency (identical to the SdH frequency in Ref. ) to the cross-section of the ellipsoidal sheet in accord with our interpretation of behavior for YSb. This underscores the similarities between these two compounds and implies that there is no need to invoke topologically non-trivial states to explain exotic magnetotransport properties neither in YSb nor in LaSb (contrary to Ref. ).
Discussion and conclusions
YSb is another material displaying giant magnetoresistance (1100 % in 9 T), three orders of magnitude smaller than that of sister compound LaSb, thus it cannot be termed ’extreme magnetoresistance’ (XMR). This is due mainly to its lower carrier mobility and weaker electron-hole compensation revealed by our Hall effect measurement.
Kohler scaling analogous to that shown in Figs. 1(c) and (d) has recently been used to explain the remarkable up-turn behavior of in WTe without the field-induced metal-insulator transition or significant contribution of an electronic structure change. The same authors have shown that perfect carrier compensation leads to exponent 2 in this scaling. Kohler scaling for YSb yielded for our samples the exponents 1.64 and 1.74, which seems related to weaker carrier compensation than nearly perfect one in WTe, where 1.92 . Thus, analogously to WTe, Kohler scaling indicates that the field-induced metal-insulator transition is unnecessary to explain up-turn and low-temperature plateau of resistivity in YSb. The origin of the up-turn is a combination of magnetoresistance with the low-temperature resistivity plateau present already at zero field. Given the similarity of YSb and LaSb the same may also be true for the latter compound.
Comparing results of SdH measurements with those of electronic structure calculations we obtained comprehensive description of the Fermi surface of YSb. Presence of both electron and hole sheets of similar volumes provides partial charge compensation responsible for its strong magnetoresistance. Band structures of YSb and LaSb are very similar. All Fermi sheets in YSb but the smallest one centered at -point have their counterparts in LaSb.[11, 12, 10] Our analysis of angular behavior of SdH frequencies in YSb indicates it is related to the three-dimensional FS, in line with Hasegawa and Settai et al. findings for LaSb, but not connected to possible non-trivial topology of electronic structure analogous to that suggested by Tafti et al.
Angular behavior of can also be perfectly explained by anisotropy of 3D-FS. When field strength is scaled by the angle-dependent factor , all data of Fig. 2(a) collapse on single curve. The effective mass anisotropy factor , obtained from the fit of with the expression , is in excellent agreement with we estimated for -sheet of FS. This is not surprising, since the mass anisotropy directly reflects the shape of FS, but it shows that angular behavior of in YSb is mainly governed by anisotropic form of -sheet of FS. That sheet corresponds to the electron band, all other FS-sheets contain holes and are nearly isotropic. The effective mass and mobility of -sheet electrons change significantly with field angle, which strongly modifies the magnetoresistance.
It has been proposed that the magnetic field induces the reconstruction of the FS in a Dirac semimetal by breaking the time reversal invariance[23, 24, 25, 26]. Assisted by the high mobility of carriers such reconstruction has been suggested to induce very large observed in CdAs and NbSb [27, 28]. We also observe features in the electronic structure of YSb, buried under the Fermi level, which may possibly allow the magnetic field to transform this compound into Dirac semimetal. A small gap between inverted bands near the X-point (cf. inset to Fig. 4) might result in topologically non-trivial states. The effect of FS reconstruction could be similar to temperature-induced Lifshitz transition in WTe , whereas its mechanism might be related, for example, to that of Lifshitz transition driven by magnetic field in CeIrIn . Very recently Dirac states have been observed by angle-resolved-photoemission spectroscopy in NbSb, a compound with bulk electronic structure very similar to that of YSb, however topologically protected states were not detected in YSb by this method.
Although a small contribution of topologically non-trivial 2D states cannot be completely excluded our analysis of magnetoresistance and Shubnikov–de Haas effect provides strong support for 3D-Fermi surface scenario of magnetotransport properties in YSb. Analogous field-induced properties of LaSb can most probably be also described in the framework of 3D multiband model.
Measurements were performed using a Physical Property Measurement System (Quantum Design) on two samples cut from one single crystal and labeled as #1 and #1h, and a sample cut from another single crystal and labeled #2. All samples had shapes of rectangular cuboid with all edges along crystallographic directions. Their sizes were: , and , for samples #1, #1h and #2, respectively. The electric current was always flowing along crystallographic direction. Single crystals were grown from Sb flux and their NaCl-type crystal structure was confirmed by powder X-ray diffraction carried out using an Xâpert Pro (PANanalytical) diffractometer with Cu-Ka radiation. No other phases were detected and lattice parameter of 6.163 Å was determined, reasonably close to literature value 6.155 Å. Electronic structure calculations were carried out using FPLO-9.00-34 code within generalized gradient approximation (GGA) method. The full-relativistic Dirac equation was solved self-consistently, treating exactly all relativistic effects, including the spin-orbit interaction without any approximations. The Perdew–Burke–Ernzerhof exchange-correlation potential was applied and the energies were converged on a dense mesh with 24 points. The convergence was set to both the density (10 in code specific units) and the total energy (10 Hartree). For the Fermi surface a 64 mesh was used to ensure accurate determination of the Fermi level.
We are very grateful to Dr. Zhili Xiao for discussion. This work was supported by the National Science Centre of Poland, grant no. 2015/18/A/ST3/00057. The band structure calculations were carried out at the Wrocław Centre for Networking and Supercomputing, grant no. 359.
O.P. conducted all the experiments, O.P. and P.W. conceived the experiments and analyzed their results, P.S. carried out the electronic structure calculations and prepared Figures 5 and 6. All authors reviewed the manuscript.
Supp. Mat. page 1 of 3
Giant magnetoresistance, three-dimensional Fermi surface
and origin of resistivity plateau in YSb semimetal
Orest Pavlosiuk, Przemysław Swatek and Piotr Wiśniewski
Institute of Low Temperatures and Structure Research, Polish Academy of Sciences, Wrocław, Poland
Resistivity and its universal low-temperature plateau
In Fig. S1(a) we show metallic-like temperature dependence of electrical resistivity of samples #1 and #2. The room- to residual-resistivity ratio, , is 13 and 22, respectively. Behavior of of sample #1 in different applied magnetic fields is shown in Fig. 1(b). Plateau at K is observed independent of strength of applied field.
Supp. Mat. page 2 of 3
Determination of effective masses
Observation of clear SdH oscillations at different temperatures allowed us to determine effective masses of charge carriers responsible for their strongest components. Fitting of FTT amplitudes for components denoted as and with function with set at 9 T and constant T/K), as shown in Fig. S2, yielded effective masses: and .
In their recent arXiv preprint Yu et al. reported and .
Multi-component Lifshitz-Kosevich function
where for -th SdH component:
is frequency, – the Dingle reduction factor, – the temperature reduction factor and – the spin factor.
is the phase, denotes harmonic number , cyclotron mass (in units of free electron mass ), effective g-factor and the Dingle temperature. T/K) is constant.
In order to fit complex data shown in Fig. 6 without preassuming the harmonic numbers of components we used a simplified version of the Equation S1:
Supp. Mat. page 3 of 3
Here represents , whereas
comprises and (with and ).
Initial fit revealed that and , thus we could identify second and fourth component as second and third harmonic, respectively, of the strongest oscillation with T.
Accordingly, we put constraints on parameters: and .
For K (temperature of our measurement), and (since and , as it was shown above) one obtains T and T.
When and were fixed at these values, the final fit yielded parameters collected in the Table below:
Supplementary Table Parameters obtained from fit of the multi-component Lifshitz-Kosevich function to data-points shown in Fig. 6 of the main text. Effective masses and Dingle temperatures calculated from and parameters are also shown.
1 2 3 4 5 6 (T) 360 720 740 1072 1160 1430 FS-sheet assignment: 0.74(3) 0.50(4) 0.74(2) 0.60(2) 0.68(3) 0.31(2) 0.41 0.41 0.11 0.41 0.005 0.0037 11.51 23.02 11.91 34.53 30 40 () 0.29 – 0.3 – 0.76 1.00 4.56 9.11 22.93 13.67 54.24 64.00 (K) 1.07 – 5.2 – 4.9 4.35
The ratio is the anisotropy of effective mass of electrons on the -sheet of FS. This value is in excellent agreement with values of 3.4 and 3.6 derived from angular behavior of and SdH oscillations, respectively.
- Vogt, O. & Mattenberger, K. Magnetic measurements on rare earth and actinide monopnictides and monochalcogenides. in Handbook on the Physics and Chemistry of Rare Earths, vol. 17 (eds. Choppin, G. R., Eyring, L. & Lander, G.) chapter 114 (Elsevier, 1993).
- Gambino, R. J., Eastman, D. E., McGuire, T. R., Moruzzi, V. L. & Grobman, W. D. Evidence of the intrinsic metallic character of some rare-earth pnictides. J. Appl. Phys. 42, 1468 (1971).
- Hayashi, J. et al. Structural phase transition of ScSb and YSb with a NaCl-type structure at high pressures. Solid State Commun. 125, 543–546 (2003).
- Tütüncü, H. M., Bagci, S. & Srivastava, G. P. Electronic, elastic and phonon properties of the rock-salt LaSb and YSb. J. Phys. Condens. Matter 19, 156207 (2007).
- Huang, X. et al. Observation of the Chiral-Anomaly-Induced Negative Magnetoresistance in 3D Weyl Semimetal TaAs. Phys. Rev. X 5, 031023 (2015).
- Shekhar, C. et al. Extremely large magnetoresistance and ultrahigh mobility in the topological Weyl semimetal candidate NbP. Nat. Phys. 11, 645–649 (2015).
- Liang, T. et al. Ultrahigh mobility and giant magnetoresistance in the Dirac semimetal CdAs. Nat. Mater. 14, 280–284 (2014).
- Ali, M. N. et al. Large, non-saturating magnetoresistance in WTe. Nature 514, 205 (2014).
- Ali, M. N. et al. Correlation of crystal quality and extreme magnetoresistance of WTe. Europhys. Lett. 110, 67002 (2015).
- Tafti, F. F., Gibson, Q. D., Kushwaha, S. K., Haldolaarachchige, N. & Cava, R. J. Resistivity plateau and extreme magnetoresistance in LaSb. Nat. Phys. 12, 272–277 (2016).
- Hasegawa, A. Fermi surface of LaSb and LaBi. J. Phys. Soc. Japan 54, 677–684 (1985).
- Settai, R. et al. Acoustic de Haas-Van Alphen effect in LaSb and CeSb. Physica B 186-188, 176–178 (1993).
- Wang, Y. L. et al. Origin of the turn-on temperature behavior in WTe. Phys. Rev. B 92, 180402 (2015).
- Ghimire, N. J., Botana, A. S., Phelan, D., Zheng, H. & Mitchell, J. F. Magnetotransport of single crystalline YSb. J. Phys. Condens. Matter 28, 235601 (2016).
- Yu, Q.-H., Wang, Y.-Y., Xu, S. & Xia, T.-L. Magnetoresistance and Shubnikov-de Hass oscillation in YSb. arXiv preprint no. 1604.05912 http://arxiv.org/abs/1604.05912 (2016).
- Blatter, G., Geshkenbein, V. B. & Larkin, A. I. From isotropic to anisotropic superconductors: A scaling approach. Phys. Rev. Lett. 68, 875–878 (1992).
- Thoutam, L. R. et al. Temperature-Dependent Three-Dimensional Anisotropy of the Magnetoresistance in WTe. Phys. Rev. Lett. 115, 046602 (2015).
- Hurd, C. M. The Hall Effect in Metals and Alloys (Springer US, Boston, MA, 1972).
- Zeng, M. et al. Topological semimetals and topological insulators in rare earth monopnictides. arXiv preprint no. 1504.03492 http://arxiv.org/abs/1504.03492 (2015).
- Lifshitz, I. M. & Kosevich, L. M. On the Theory of the Shubnikov-De Haas Effect. Sov. Phys. JETP 6, 67–77 (1958).
- Shoenberg, D. Magnetic Oscillations in Metals (Cambridge University Press, Cambridge, 1984).
- Seiler, D. G. The Shubnikov-de Haas Effect: A Powerful Tool for Characterizing Semiconductors. In Landwehr, G. (ed.) High Magn. Fields Semicond. Phys. II, 578–587 (Springer, 1989).
- Young, S. M. et al. Dirac Semimetal in Three Dimensions. Phys. Rev. Lett. 108, 140405 (2012).
- Wang, Z. et al. Dirac semimetal and topological phase transitions in Bi (=Na, K, Rb). Phys. Rev. B 85, 195320 (2012).
- Burkov, A. A., Hook, M. D. & Balents, L. Topological nodal semimetals. Phys. Rev. B 84, 235126 (2011).
- Liu, Z. K. et al. Discovery of a three-dimensional topological Dirac semimetal, NaBi. Science 343, 864–867 (2014).
- Wang, Z., Weng, H., Wu, Q., Dai, X. & Fang, Z. Three-dimensional Dirac semimetal and quantum transport in CdAs. Phys. Rev. B 88, 125427 (2013).
- Wang, K., Graf, D., Li, L., Wang, L. & Petrovic C. Anisotropic giant magnetoresistance in NbSb. Sci. Rep. 4, 7328 (2014).
- Wu, Y. et al. Temperature-Induced Lifshitz Transition in WTe. Phys. Rev. Lett. 115, 166602 (2015).
- Aoki, D. et al. Field-induced Lifshitz transition without metamagnetism in CeIrIn. Phys. Rev. Lett. 116, 037202 (2016).
- Neupane, M. et al. Observation of Dirac-like semi-metallic phase in NdSb. J. Phys. Condens. Matter 28, 23LT02 (2016).
- He, J. et al. Distinct Electronic Structure for the Extreme Magnetoresistance in YSb. arXiv preprint no. 1607.07855 http://arxiv.org/abs/1607.07855 (2016).
- Koepernik, K. & Eschrig, H. Full-potential nonorthogonal local-orbital minimum-basis band-structure scheme. Phys. Rev. B 59, 1743–1757 (1999).
- Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865–3868 (1996).