Abstract
We present a novel methodology to identify topologically non-trivial materials based on band inversion induced by spin-orbit coupling (SOC) effect. Specifically, we compare the density functional theory (DFT) based wavefunctions with and without spin-orbit coupling and compute the âspin-orbit-spillageâ as a measure of band-inversion. Due to its ease of calculation, without any need for symmetry analysis or dense k-point interpolation, the spillage is an excellent tool for identifying topologically non-trivial materials. Out of 30000 materials available in the JARVIS-DFT database, we applied this methodology to more than 4835 non-magnetic materials consisting of heavy atoms and low bandgaps. We found 1868 candidate materials with high-spillage (using 0.5 as a threshold). We validated our methodology by carrying out conventional Wannier-interpolation calculations for 289 candidate materials. We demonstrate that in addition to Z2 topological insulators, this screening method successfully identified many semimetals and topological crystalline insulators. Importantly, our approach is applicable to the investigation of disordered or distorted as well as magnetic materials, because it is not based on symmetry considerations. We discuss some individual example materials, as well as trends throughout our dataset, which is available at the websites: https://www.ctcms.nist.gov/~knc6/JVASP.html and https://jarvis.nist.gov/.
Similar content being viewed by others
Introduction
Topological materials are a new class of quantum materials, which are characterized by surface states induced by the topology of the bulk band structure1. In the past decade, there has been a huge interest in the field of topological materials, which has rapidly expanded to include Chern insulators2, topological insulators (TI)3,4, crystalline topological insulators (CTI)5, Dirac semimetals (DS)6, Weyl semimetals (WS)7, quadratic band-crossing materials8 with many other related subclasses and variations. While some topological materials classes exist without spin-orbit coupling, many topological behaviors are related to spin-orbit coupling (SOC) induced band inversion9. Despite their huge potential for technological applications such as robust quantum computation10, only a few examples of most classes of TMs are known, and a workflow/metrology for discovering, characterizing and finally developing a repository of such materials is still in the developing phase. In this work, we use the spin-orbit spillage criteria11 (discussed later) and Wannier interpolation12,13,14 to quickly screen topological materials. The materials of interest in the present work include topological insulators3,15, Dirac semimetals16, Weyl-semimetals17, quadratic band-crossing materials8, crystalline topological insulators5,18 and nodal materials19,20.
Historically, individual TMs have been initially identified21,22,23 mainly by intuition, trial and error, heuristic screening24, and analogy to other materials25, often with first principles calculations leading the way. Very recently, there have been several efforts to perform a systematic analysis of band structure symmetries to help identify TMs in systems with enough crystalline symmetries such as by Slager et al., Bradlyn et al., Tang et al. and Zhang et al.26,27,28,29,30,31,32,33,34. In this work, we take a complementary approach to the symmetry-based analysis by directly searching for materials with large changes in the occupied wave functions due to including SOC in a first principles calculation. This approach can identify topological band features and phase transitions in materials with low symmetry, or âaccidentalâ band crossings at k-points with low symmetry. We can apply our search technique to materials with disordered or distorted structures, which are difficult to analyze with symmetry-based techniques, and our search method can easily be extended to systems with broken time-reversal symmetry, like Chern insulators or magnetic Weyl semimetals. As mentioned above, most classes of TMs require SOC, which leads to band-inversion or splitting of otherwise degenerate bands11,35. We calculate the k-dependent spin-orbit spillage, and then we screen for materials with large spillage at least at one k-point, as an indicator of band inversion. To classify the previously identified materials, we then use Wannier-interpolation to efficiently calculate topological indices. Using this workflow, we discover hundreds of potential TMs. We added SOC and non-SOC band structures as well as their spillage results to the JARVIS-DFT repository. The JARVIS-DFT website contains about 30000 bulk and 800 two-dimensional materials with their density functional theory (DFT)-computed structural, energetics36, elastic37, optoelectronic38 and solar-cell efficiency39 properties. Finally, we show examples of several specific TMs identified in this search, and we provide a full list of TMs we have identified in the Supplementary Materials.
Changes in bands structures due to SOC are generally â¤1âeV and are larger in materials with heavy atoms. Hence, we limit our initial set of materials to those having heavy atoms (atomic weightââ¥â65) with bandgapââ¤â0.6âeV and zero magnetic moment value, which results in 4835 materials out of 30000 materials. We calculate the spillage of these materials by comparing wavefunctions from DFT calculations with and without SOC for each material. As introduced by Liu et al.11, the spin-orbit spillage is a method for comparing wave functions at a given k-point with and without SOC, with the goal of identifying band inversion. We calculate the spin-orbit spillage, \(\eta ({\bf{k}})\), given by the following equation:
where \(P({\bf{k}})={\sum }_{n=1}^{{n}_{occ}({\bf{k}})}\,|{\psi }_{n{\bf{k}}}\rangle \langle {\psi }_{n{\bf{k}}}|\) is the projector onto the occupied wavefunctions without SOC, and \(\tilde{P}\) is the same projector with SOC. Here, âTrâ denotes trace over the occupied bands. We can write the spillage equivalently as:
where \({M}_{mn}({\bf{k}})=\langle {\psi }_{m{\bf{k}}}|{\tilde{\psi }}_{n{\bf{k}}}\rangle \) is the overlap between occupied Bloch functions with and without SOC at the same wave vector k. If the SOC does not change the character of the occupied wavefunctions, the spillage will be near zero, while a band inversion will result in a large spillage.
We have extended the original definition of the spillage, which was only defined for insulators, for use in metals by allowing the number of occupied electrons, \({n}_{occ}({\bf{k}})\), to vary at different k-points by fixing it to the zero-temperature occupation of the non-SOC calculation. As argued in ref.11, in the case where the non-SOC calculation gives a trivial insulator and the SOC calculation gives a Z2 topological insulator, the spin-orbit spillage will be at least 2.0 at the k-point(s) where band inversion occurs. At other k-points, or in cases where both materials/calculations are trivial, the spillage is generally much lower. We find that a similar relation holds when SOC drives a crystalline topological phase transition as well. Due to its ease of calculation, without any need for symmetry analysis or dense k-point interpolation, the spillage is an excellent tool for identifying topologically non-trivial insulators and understanding where in k-space the band inversion occurs; although, it is not sufficient to classify the type of TI in general.
In this work, we extend the use of the spillage to cases where one or both of the non-SOC or SOC calculations are metallic, which significantly broadens the class of TMs that we can search for. In the case where either the non-SOC or SOC calculations are metallic, there is no specific spillage value that guarantees a band inversion, but we find empirically that spillages above 0.5 are often due to topologically non-trivial features in metals. For example, in the well-known Dirac material Na3Bi16, both the SOC and non-SOC calculations have a 3D Dirac cone in the bulk, but the addition of SOC causes the Dirac point to shift off the high symmetry point and onto a symmetry line. Due to this shift, the symmetry of the occupied bands in this region of k-space will be different in SOC and non-SOC calculations, which results in a large spillage. Similarly, Weyl points and nodes require SOC, so Weyl materials will always have large spillage near these features. In addition, there are many cases where the addition of SOC can lift a degeneracy, turning a metal into a non-trivial insulator. This gap opening will also generate a large spillage, as again the occupied wavefunctions will change due to SOC. We note that we do not expect the spillage to be useful in identifying TMs where the SOC does not play an important role.
Results and Discussion
To further elucidate how spillage helps in the screening of TMs, we show an example of Ba2HgPb in Fig. 1. Looking at the band structures shown in Fig. 1a (without SOC) and 1b (with SOC), it is unclear if there as a topological phase transition due to SOC. However, the band structures projected onto the Pb [p-orbital] states (Fig. 1c,d, without and with SOC, respectively), suggest that there is a SOC-induced band-inversion at the L-point. While this band inversion is a strong sign that Ba2HgPb is topologically non-trivial, automating this type of visual analysis based on orbital projections is difficult, as many orbitals can contribute to bands near the Fermi level in complex and low symmetry materials. However, the non-trivial band inversion becomes immediately clear in Fig. 1e, which shows the spillage at each k-point, including a peak above 2.0 at the L point. This band inversion at the L point is consistent with Wannier charge center (Figs S1 and S2) calculations that show that Ba2HgPb is a crystalline topological insulator. In the JARVIS-DFT database (https://www.ctcms.nist.gov/~knc6/JVASP.html and https://jarvis.nist.gov) we have provided plots similar to Fig. 1, including projected band structures, the spillage, and information on band gap, for all materials investigated in this work. In the Supplementary Materials, we provide examples of these plots for several trivial and non-trivial materials classes (Figs S3 and S7) to elaborate on the spillage criteria. We expect the most promising TMs to have a small number of high peaks in the spillage. We also encounter some metals with high spillage values across the Brillouin zone, which indicates a strong change in the band structure at the Fermi level due to SOC. We find that these materials are less likely to have simple topological features, and instead display many complicated band crossings near the Fermi level. Therefore, after using the maximum spillage as a screening criterion, it is necessary to perform additional analysis to determine the source of the high spillage and classify any topological behavior.
As discussed earlier, we calculate the spillage values for a large set of low bandgap and high atomic mass materials, with a spillage value of 0.5 as a threshold to screen for potential topologically non-trivial materials. In Fig. 2a, we show the obtained distribution of SOC spillage, with most materials having low spillage even if they had high atomic mass. In Fig. 2bâf) results are shown only for high-spillage materials (spillageââ¥â0.5). Overall the spillage ranges from zero to four.
In Fig. 2b we analyze the SOC bandgap distribution of all the materials with spillage greater than or equal to 0.5. We find that there are many metallic/semi-metallic materials (about 40% of the high-spillage materials), and wide distribution of band gaps amongst the insulating materials. In Fig. 2c we show the SOC band gap as a function of spillage, and we find clusters of data at spillage values 1 and 2. Our data indicate no correlation between SOC band gap and spillage. It is to be noted that both NSOC and SOC band gaps have the well-known problem of being underestimated as in the case of conventional DFT. However, as this is a systematic effect, it should not affect the comparison between NSOC and SOC values. In Fig. 2d we plot the SOC vs NSOC bandgaps of both high-spillage (spillageââ¥â0.5, blue stars) and low-spillage (red dots) materials. For low spillage materials, there is a fairly strong linear relationship between the NSOC and SOC gaps, indicating that in most cases spin-orbit coupling has little effect on band gaps of topologically trivial materials. However, there are also materials where a large change in band gap occurs, even among low spillage materials. For high spillage materials, there is no relationship between SOC and NSOC band gaps, as expected for materials which have a different band ordering near the Fermi level. This suggests that NSOC band structures should be used with care when studying small gap semiconductors, even if the topological behavior is not specifically desired. In Fig. 2e we show the crystal system distribution of all the high spillage materials. We observe that many high spillage materials are cubic, which is in part due to a large number of both half-Heusler alloys and antiperovskites in our DFT database, which are known to display topologically non-trivial behavior40,41,42. There are also many hexagonal materials, including both materials related to the topological insulator Bi2Se3, as well as hexagonal ABC materials43,44,45,46. Similar behavior is observed in the space-group distribution in Fig. 2f. Results in Fig. 2e,f imply that high symmetry in a crystal favors the rise of topologically non-trivial behavior. Moreover, in Fig. 2g we see that most of these materials are ternary; however, there are elemental, binary, and multicomponent non-trivial materials as well. Furthermore, we classify the non-trivial materials in terms of the compositional prototype in Fig. 2h. We also identify the dimensionality of the materials based on the lattice-constant36 and data-mining approaches47. We found 89.61% materials to be 3D-bulk, 10.02% to be 2D-bulk, 0.32% to be 1D-bulk and 0.05% to be 0D materials with high spillage values.
Next, we analyze the distribution of various high spillage materials as a function of composition. For each element, we calculate the probability that a material in our dataset containing that element has a spillage above 0.5 and is, therefore, a TM candidate. We present the results in the form of a periodic table in Fig. 3. Consistent with known TMs, we observe that materials containing the heavy elements Bi and Pb are by far the likeliest ones to have high spillage. Other high spillage elements are Ta and the 5d elements Re-Pt. To contribute to SOC-induced band inversion, an element must both have significant SOC and contribute to bands located near the Fermi level, which favors heavy elements with moderate electronegativity.
Using the spillage as a screening tool, we greatly reduce the number of candidate TMs to consider, but further analysis is necessary to understand the class of each potential TM. For a few hundred promising candidates, we use Wannier-based tight-binding models to perform efficient Wannier charge center calculations and node finding calculations. We present the results of these calculations in the Supplementary Information. We characterize 40 strong topological insulators (TI), 7 weak TIs, 24 crystalline topological insulators (CTI), 30 Dirac semimetals, 3 Weyl semimetals, 47 quadratic band-crossing materials. We identify many previously known topological materials of various classes. In addition, we find a variety of promising TMs that, as far as we know, have not been previously reported in the literature. As examples of various types, we present the band structures of a variety of TMs in Fig. 4. As shown in Fig. 4a, we find that PbS in space group P63/mmc (JVASP-35680) is a strong Z2 topological insulator due to band inversion at Î. This topology is reflected in the (001) surface band structure (Fig. 4e), which shows a Dirac cone feature at Î. In Fig. 4b, we show the band structure of LiBiS2 (P4/mmm, JVASP-52332), which is a weak topological insulator due to band inversion at the R point. Figure 4f shows the surface band structure, including Dirac cone features at M and Z points. In Fig. 4c, we show the band structure of KHgAs (P63/mmc, JVASP-15801), which is a crystalline topological insulator48. The surface band structure (Fig. 4g) shows non-trivial surface bands. In Fig. 4d,h, we show the band structure of InSb (JVASP-36123) and GaSb (JVASP-35711) in the P63mc space group. InSb has a simple bulk Dirac cone feature along the Î â A line. In contrast, GaSb has a slightly different band ordering that leads to two Dirac cones along Î â A and complicated Weyl crossings along Î â K and Î â M. We find that materials with several crossings near the Fermi level, like GaSb, are relatively common, but simple semimetals like InSb are relatively rare.
Conclusions
In summary, we have used the spin-orbit spillage criteria to quickly screen a large material database for SOC-related TMs, which we then classified using Wannier charge sheets and band crossing searches. This approach for high-throughput TM discovery has the advantages of being applicable to materials with any symmetry, computationally efficient, and applicable to both strong and weak TIâs, as well as many crystalline TIs and semimetals. We found at least 1868 materials to be candidate TMs out 4835 investigated materials, which we present in a database, with the specific topology of hundreds already classified. These new TMs include materials from many crystal structure and chemistries can significantly expand the range of known TMs. Additionally, we show a unique way of extracting information from DFT wavefunctions. We believe that the database can be useful to accelerate the discovery and characterization for TMs and other properties tied to SOC.
Methods
We perform our spillage-based screening using Vienna Ab-initio Simulation Package (VASP)49,50 and the PBE functional51 (with and without SOC). All input structures in this work are relaxed using the OptB88vdW functional52, which is known to give accurate results for both vdW and not-vdW-bonded materials, as discussed in ref.37. However, OptB88vdWâ+âSOC is not publicly available in VASP-package, so we computed spillage using PBE53. This approach is reliable because recently Cao et al.54 have shown that band-structures obtained using PBE on vdW-functional-relaxed structure display all the correct features. We also added a comparison of non-spin-orbit coupling bandstructures on same OptB88vdW relaxed geometry using PBE and OptB88vdW functionals for Bi2Te3 and Ge as examples (Fig. S10).
We calculate the spillage using 1000/atom k-points, as well as by analyzing the spillage along a high-symmetry Brillouin-zone path, with a 600âeV plane-wave cut-off. High values of the spillage are generally restricted to regions of k-space where band inversion occurs, so as long as the k-point grid is dense enough to locate these regions, the spillage can be calculated with a sparse k-point grid (see Table S1). For topological insulators, band inversion generally happens at high-symmetry points, which are included in any typical k-point grid. For semimetals, where topological features can occur at arbitrary k-points, it is necessary to search on a grid; however, we allow tolerance for what we consider high spillage (spillageââ¥â0.5) in part to avoid the need to hit the highest spillage point exactly in order to identify a potentially interesting material.
Analysis of the spillage can in some cases help identify the type of topological phase. For example, strong topological insulators with Bi2Se3-like structures have band inversion at the Ð-point only, and the spillage reflects this pattern, with spillage above 2 at Ð-point, but dropping to near zero when moving along away from Ð. In contrast, for the weak topological insulator BiSe (JVASP-5410), which has a strongly two-dimensional band structure, the bands are inverted along the entire Ð-A line, which is also reflected in the spillage. However, the spillage cannot be used to identify the topological phase in all cases. For example, Tl8Sn2Te6 (JVASP-9377) does have a SOC-induced band inversion at Ð, however, more detailed analysis (as discussed below) reveals that the resulting band structure is topologically trivial. False positive examples like this one are very rare. These examples are presented in the Supplementary Materials (Figs S3, S8 and S9).
Therefore, in order to investigate 289 candidate materials further, we used Quantum Espresso (QE)55 with norm-conserving pseudopotentials56,57 and Wannier9058,59 to generate first principles tight-binding models. Using WannierTools12, we then performed Wannier change center (Wilson loop)60,61 calculations of topological invariants, a numerical search for band crossings, and surface state calculations.
Data Availability
The electronic structure data is available at the JARVIS-DFT website: https://www.ctcms.nist.gov/~knc6/JVASP.html and http://jarvis.nist.gov. We also provide a JSON file with the spillage, bandgap (with and without spin-orbit coupling) and topological material type information at the figshare link: https://doi.org/10.6084/m9.figshare.7594571.v1.
References
Yan, B. & Felser, C. Topological materials: Weyl semimetals. Ann. Rev. Cond. Mat. Phys. 8, 337â354 (2017).
Thonhauser, T. & Vanderbilt, D. Insulator/Chern-insulator transition in the Haldane model. Physical Review B 74, 235111 (2006).
Hasan, M. Z. & Kane, C. L. Colloquium: topological insulators. Reviews of Modern Physics 82, 3045 (2010).
Ortmann, F., Roche, S. & Valenzuela, S. O. Topological insulators: Fundamentals and perspectives. (Wiley-VCH Verlag, Weinheim, Germany, 2015).
Fu, L. Topological crystalline insulators. Physical Review Letters 106, 106802 (2011).
Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197 (2005).
Huang, S.-M. et al. A Weyl Fermion semimetal with surface Fermi arcs in the transition metal monopnictide TaAs class. Nature communications 6, 7373 (2015).
Sun, K., Yao, H., Fradkin, E. & Kivelson, S. A. Topological insulators and nematic phases from spontaneous symmetry breaking in 2D Fermi systems with a quadratic band crossing. Physical Review Letters 103, 046811 (2009).
Moore, J. E. The birth of topological insulators. Nature 464, 194 (2010).
Nayak, C., Simon, S. H., Stern, A., Freedman, M. & Sarma, S. D. Non-Abelian anyons and topological quantum computation. Reviews of Modern Physics 80, 1083 (2008).
Liu, J. & Vanderbilt, D. Spin-orbit spillage as a measure of band inversion in insulators. Physical Review B 90, 125133 (2014).
Wu, Q., Zhang, S., Song, H.-F., Troyer, M. & Soluyanov, A. A. WannierTools: An open-source software package for novel topological materials. Computer Physics Communications 224, 405â416 (2018).
Kane, C. L. & Mele, E. J. Z2 topological order and the quantum spin Hall effect. Physical review letters 95, 146802 (2005).
Gresch, D. et al. Z2 Pack: Numerical implementation of hybrid Wannier centers for identifying topological materials. Physical Review B 95, 075146 (2017).
Liu, Y. et al. Tuning Dirac states by strain in the topological insulator Bi2Se3. Nature Physics 10, 294 (2014).
Liu, Z. et al. Discovery of a three-dimensional topological Dirac semimetal, Na3Bi. Science 343, 864â867 (2014).
Burkov, A. & Balents, L. Weyl semimetal in a topological insulator multilayer. Physical review letters 107, 127205 (2011).
Ma, J. et al. Experimental evidence of hourglass fermion in the candidate nonsymmorphic topological insulator KHgSb. Science Advances 3, e1602415 (2017).
Bian, G. et al. Topological nodal-line fermions in spin-orbit metal PbTaSe2. Nature communications 7, 10556 (2016).
Burkov, A., Hook, M. & Balents, L. Topological nodal semimetals. Physical Review B 84, 235126 (2011).
Fu, L. & Kane, C. L. Topological insulators with inversion symmetry. Physical Review B 76, 045302 (2007).
Fang, C., Gilbert, M. J. & Bernevig, B. A. Bulk topological invariants in noninteracting point group symmetric insulators. Physical Review B 86, 115112 (2012).
Hughes, T. L., Prodan, E. & Bernevig, B. A. Inversion-symmetric topological insulators. Physical Review B 83, 245132 (2011).
Yang, K., Setyawan, W., Wang, S., Nardelli, M. B. & Curtarolo, S. A search model for topological insulators with high-throughput robustness descriptors. Nature materials 11, 614 (2012).
Lin, H. et al. Adiabatic transformation as a search tool for new topological insulators: Distorted ternary Li2 AgSb-class semiconductors and related compounds. Physical Review B 87, 121202 (2013).
Slager, R.-J., Mesaros, A., JuriÄiÄ, V. & Zaanen, J. The space group classification of topological band-insulators. Nature Physics 9, 98 (2013).
Bradlyn, B. et al. Topological quantum chemistry. Nature 547, 298 (2017).
Tang, F., Po, H. C., Vishwanath, A. & Wan, X. Towards ideal topological materials: Comprehensive database searches using symmetry indicators. arXiv preprint arXiv:1807.09744 (2018).
Tang, F., Po, H. C., Vishwanath, A. & Wan, X. Topological Materials Discovery By Large-order symmetry indicators. arXiv preprint arXiv:1806.04128 (2018).
Tang, F., Po, H. C., Vishwanath, A. & Wan, X. Efficient Topological Materials Discovery Using Symmetry Indicators. arXiv preprint arXiv:1805.07314 (2018).
Zhou, X. et al. Topological crystalline insulator states in the Ca $ _2 $ As family. arXiv preprint arXiv:1805.05215 (2018).
Vergniory, M., Elcoro, L., Felser, C., Bernevig, B. & Wang, Z. The (High Quality) Topological Materials In The World. arXiv preprint arXiv:1807.10271 (2018).
Po, H. C., Vishwanath, A. & Watanabe, H. Symmetry-based indicators of band topology in the 230 space groups. Nature Communications 8, 50 (2017).
Zhang, T. et al. Catalogue of Topological Electronic Materials. arXiv preprint arXiv:1807.08756 (2018).
Zhu, Z., Cheng, Y. & Schwingenschlögl, U. Band inversion mechanism in topological insulators: A guideline for materials design. Physical Review B 85, 235401 (2012).
Choudhary, K., Kalish, I., Beams, R. & Tavazza, F. High-throughput Identification and Characterization of Two-dimensional Materials using Density functional theory. Scientific Reports 7, 5179 (2017).
Choudhary, K., Cheon, G., Reed, E. & Tavazza, F. Elastic properties of bulk and low-dimensional materials using van der Waals density functional. Physical Review B 98, 014107 (2018).
Choudhary, K. et al. Computational screening of high-performance optoelectronic materials using OptB88vdW and TB-mBJ formalisms. Scientific data 5, 180082 (2018).
Choudhary, K. et al. Accelerated Discovery of Efficient Solar-cell Materials using Quantum and Machine-learning Methods. arXiv preprint arXiv:1903.06651 (2019).
Chadov, S. et al. Tunable multifunctional topological insulators in ternary Heusler compounds. Nature materials 9, 541 (2010).
Yu, R., Weng, H., Fang, Z., Dai, X. & Hu, X. Topological node-line semimetal and Dirac semimetal state in antiperovskite Cu3PdN. Physical review letters 115, 036807 (2015).
Hsieh, T. H., Liu, J. & Fu, L. Topological crystalline insulators and Dirac octets in antiperovskites. Physical Review B 90, 081112 (2014).
Zhang, J. et al. Topological band crossings in hexagonal materials. arXiv preprint arXiv:1805.05120 (2018).
Zhang, H. et al. Topological insulators in Bi2Se3, Bi2Te3 and Sb2Te3 with a single Dirac cone on the surface. Nature physics 5, 438 (2009).
Monserrat, B., Bennett, J. W., Rabe, K. M. & Vanderbilt, D. Antiferroelectric Topological Insulators in Orthorhombic A MgBi Compounds (A = Li, Na, K). Physical review letters 119, 036802 (2017).
Di Sante, D. et al. Intertwined Rashba, Dirac, and Weyl fermions in hexagonal hyperferroelectrics. Physical review letters 117, 076401 (2016).
Cheon, G. et al. Data mining for new two-and one-dimensional weakly bonded solids and lattice-commensurate heterostructures. Nano letters 17, 1915â1923 (2017).
Wang, Z., Alexandradinata, A., Cava, R. J. & Bernevig, B. A. Hourglass fermions. Nature 532, 189 (2016).
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Computational Materials Science 6, 15â50 (1996).
Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Physical Review B 54, 11169â11186 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Physical review letters 77, 3865 (1996).
Klimeš, J., Bowler, D. R. & Michaelides, A. Chemical accuracy for the van der Waals density functional. Journal of Physics: Condensed Matter 22, 022201 (2009).
Li, W. F., Fang, C. & van Huis, MarjinA. Strong spin-orbit splitting and magnetism of point defect states in monolayer WS2. Physical Review B 94, 195425 (2016).
Cao, G. et al. Rhombohedral Sb2Se3 as an intrinsic topological insulator due to strong van der Waals interlayer coupling. Physical Review B 97, 075147 (2018).
Giannozzi, P. et al. Advanced capabilities for materials modelling with Quantum ESPRESSO. Journal of Physics: Condensed Matter 29, 465901 (2017).
Hamann, D. Optimized norm-conserving Vanderbilt pseudopotentials. Physical Review B 88, 085117 (2013).
Scherpelz, P., Govoni, M., Hamada, I. & Galli, G. Implementation and validation of fully relativistic GW calculations: spinâorbit coupling in molecules, nanocrystals, and solids. Journal of chemical theory and computation 12, 3523â3544 (2016).
Mostofi, A. A. et al. wannier90: A tool for obtaining maximally-localised Wannier functions. Computer physics communications 178, 685â699 (2008).
Marzari, N., Mostofi, A. A., Yates, J. R., Souza, I. & Vanderbilt, D. Maximally localized Wannier functions: Theory and applications. Reviews of Modern Physics 84, 1419 (2012).
Yu, R., Qi, X. L., Bernevig, A., Fang, Z. & Dai, X. Equivalent expression of Z2 topological invariant for band insulators using the non-abelian Berry connection. Physical Review B 84, 075119 (2011).
Soluyanov, A. A. & Vanderbilt, D. Computing topological invariants without inversion symmetry. Physical Review B 83, 235401 (2011).
Author information
Authors and Affiliations
Contributions
K.C. and K.G. jointly developed the method. K.C. and K.G. carried out the high-throughput DFT calculations using VASP and QE packages respectively. F.T. helped in writing the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisherâs note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Choudhary, K., Garrity, K.F. & Tavazza, F. High-throughput Discovery of Topologically Non-trivial Materials using Spin-orbit Spillage. Sci Rep 9, 8534 (2019). https://doi.org/10.1038/s41598-019-45028-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-019-45028-y