Abstract
Higher-order topological insulators have been introduced in the precursory Benalcazar-Bernevig-Hughes quadrupole model, but no electronic compound has been proposed to be a quadrupole topological insulator (QTI) yet. In this work, we predict that Ta2M3Te5 (M = Pd, Ni) monolayers can be 2D QTIs with second-order topology due to the double-band inversion. A time-reversal-invariant system with two mirror reflections (Mx and My) can be classified by Stiefel-Whitney numbers (w1, w2) due to the combined symmetry TC2z. Using the Wilson loop method, we compute w1 = 0 and w2 = 1 for Ta2Ni3Te5, indicating a QTI with qxy = e/2. Thus, gapped edge states and localized corner states are obtained. By analyzing atomic band representations, we demonstrate that its unconventional nature with an essential band representation at an empty site, i.e., Ag@4e, is due to the remarkable double-band inversion on Y–Γ. Then, we construct an eight-band quadrupole model with Mx and My successfully for electronic materials. These transition-metal compounds of A2M1,3X5 (A = Ta, Nb; M = Pd, Ni; X = Se, Te) family provide a good platform for realizing the QTI and exploring the interplay between topology and interactions.
Similar content being viewed by others
Introduction
In higher-order topological insulators, the ingap states can be found in (d−n)-dimensional edges (n > 1), such as the corner states of two-dimensional (2D) systems or the hinge states of three-dimensional systems1,2,3,4,5,6,7,8,9,10,11,12. Different from topological insulators with (d−1)-dimensional edge states, the Chern numbers or \({{\mathbb{Z}}}_{2}\) numbers in higher-order topological insulators are zero. The higher-order topology can be captured by topological quantum chemistry13,14,15,16,17,18, nested Wilson-loop method1,8 and second Stiefel–Whitney (SW) class19,20,21,22,23,24,25. Using topological quantum chemistry13, the higher-order topological insulator can be diagnosed by the decomposition of atomic band representations (aBRs) as an unconventional insulator (or obstructed atomic insulator) with mismatching of electronic charge centers and atomic positions16,17,18,26. In contrast to dipoles (Berry phase) for topological insulators, the higher-order topological insulators can be understood by multipole moments1. In a 2D system, the second-order topology corresponds to the quadrupole moment, which can be diagnosed by the nested Wilson-loop method1,8,23. When the system contains space-time inversion symmetries, such as PT and C2zT, where the P and T represent inversion and time-reversal symmetries, the second-order topology can be described by the second SW number (w2)23,27. The second SW number w2 is a well-defined 2D topological invariant of an insulator only when the first SW number w1 = 0. Usually, a 2D quadrupole topological insulator (QTI) with w1 = 0, w2 = 1 has gapped edge states and degenerate localized corner states, which are pinned at zero energy (being topological) in the presence of chiral symmetry. When the degenerate corner states are in the energy gap of bulk and edge states, the fractional corner charge can be maintained due to filling anomaly28.
So far, various 2D systems are proposed to be SW insulators with second-order topology, such as monolayer graphdiyne29,30, liganded Xenes25,31, β-Sb monolayer18 and Bi/EuO32. However, no compound has been proposed to be a QTI with Mx and My symmetries. After considering many-body interactions in transition-metal compounds, superconductivity, exciton condensation and Luttinger liquid could emerge in a transition-metal QTI. In recent years, van der Waals layered materials of A2M1,3X5 (A = Ta, Nb; M = Pd, Ni; X = Se, Te) family have attracted attentions because of their special properties, such as quantum spin Hall effect in Ta2Pd3Te5 monolayer33,34, excitons in Ta2NiSe535,36,37, and superconductivity in Nb2Pd3Te5 and doped Ta2Pd3Te538. In particular, the monolayers of A2M1,3X5 family can be exfoliated easily, serving as a good platform for studying topology and interactions in lower dimensions.
In this work, we predict that based on first-principles calculations, Ta2Ni3Te5 monolayer is a 2D QTI. Using the Wilson-loop method, we show that its SW numbers are w1 = 0 and w2 = 1, corresponding to the second-order topology. We also solve the aBR decomposition for Ta2Ni3Te5 monolayer, and find that it is unconventional with an essential band representation (BR) at an empty Wykoff position (WKP), Ag@4c, which origins from the remarkable double-band inversion on Y–Γ line. To verify the QTI phase, we compute the energy spectrum of Ta2Ni3Te5 monolayer with open boundary conditions in both x and y directions and obtain four degenerate corner states. Then, we construct an eight-band quadrupole model with Mx and My successfully. The double-band-inversion picture widely happens in the band structures of A2M1,3X5 family. The Ta2M3Te5 monolayers are 2D QTI candidates for experimental realization in electronic systems.
Results
Band structures
The band structure of Ta2Ni3Te5 monolayer suggests that it is an insulator with a band gap of 65 meV. We have checked that spin–orbit coupling (SOC) has little effect on the band structure (Supplementary Fig. 2b). We also checked the band structures using GW method and SCAN method, and we find that the band gap remains using these methods (the corresponding band structures are shown in Supplementary Note 1). As shown in the orbital-resolved band structures of Fig. 3a–c, although low-energy bands near the Fermi level (EF) are mainly contributed by Ta-\({d}_{{z}^{2}}\) orbitals (two conduction bands) and NiA-dxz orbitals (two valence bands), the inverted bands of {Y2; GM1−, GM3−} come from Te-px orbitals. The irreducible representations (irreps)39 at Y and Γ are denoted for the inverted bands in Fig. 1b. We notice that the double-band inversion between {Y2; GM1−, GM3−} bands and {Y1; GM1+, GM3+} bands is remarkable, about 1 eV.
Atomic band representations
To analyze the band topology, the decomposition of aBR is performed. In a unit cell of Ta2Ni3Te5 monolayer in Fig. 1a, four Ta atoms, four NiB atoms and eight Te atoms are located at different 4e WKPs. The rest two Te atoms and two NiA atoms are located at 2a and 2b WKPs, respectively. The aBRs are obtained from the crystal structure by pos2aBR16,17,18, and irreps of occupied states are calculated by IRVSP39 at high-symmetry k-points. Then, the aBR decomposition is solved online—http://tm.iphy.ac.cn/UnconvMat.html. The results are listed in Supplementary Table 2 of Supplementary Note 2. Instead of being a sum of aBRs, we find that the aBR decomposition of the occupied bands has to include an essential BR at an empty WKP, i.e., Ag@4c. As illustrated in Fig. 1a, the charge centers of the essential BR are located at the middle of NiB-NiB bonds (i.e., the 4c WKP), indicating that the Ta2Ni3Te5 monolayer is a 2D unconventional insulator with second-order topology.
Double-band inversion
In an ideal atomic limit, Te-p orbitals and Ni-d orbitals are occupied, while Ta-d orbitals are fully unoccupied. Thus, all the occupied bands are supposed to be the aBRs of Te-p and Ni-d orbitals, as shown in the left panel of Fig. 2a. However, in the monolayers of A2M1,3X5 family (see their band structures in Supplementary Note 1), a double-band inversion happens between the occupied aBR A″@4e (Te-px and Ni-d) and unoccupied aBR \(A^{\prime} @4e\) (Ta-\({d}_{{z}^{2}}\)), as shown in the right two panels of Fig. 2a. When the double-band inversion happens between {Y2; GM1−, GM3−} and {Y1; GM2−, GM4−} on Y–Γ line, it results in a semimetal for Ta2NiSe5 monolayer (215-type; Fig. 2b). When it happens between {Y2; GM1−, GM3−} and {Y1; GM1+, GM3+} in Fig. 2c, the system becomes a 2D QTI for Ta2Ni3Te5 monolayer (235-type), resulting in the essential BR of Ag@4c.
Second Stiefel–Whitney class w 2 = 1
To identify the second-order topology of the monolayers, we compute the second SW number by the Wilson-loop method. The first SW class (w1) is
where \({{{{\mathcal{A}}}}}_{mn}({{{\boldsymbol{k}}}})=\left\langle {u}_{m}({{{\boldsymbol{k}}}})\right|{{{\rm{i}}}}{{{{\boldsymbol{\nabla }}}}}_{{{{\boldsymbol{k}}}}}\left|{u}_{n}({{{\boldsymbol{k}}}})\right\rangle\)22. The second SW class (w2) can be computed by the nested Wilson-loop method, or simply by m module 2, where m is the number of crossings of Wilson bands at θ = π. It should be noted that w2 is well-defined only when w1 = 0. With w1 = 0, w2 can be unchanged when choosing the unit cell shifting a half lattice constant. The 1D Wilson-loops are computed along ky. The computed phases of the eigenvalues of Wilson-loop matrices Wy(kx) (Wilson bands) are shown in Fig. 1c as a function of kx. The results show that the first SW class is w1 = 0. In addition, there is one crossing of Wilson bands at θ = π [Fig. 1d], indicating the second SW class w2 = 1. The quadruple moment qxy = e/2 calculated by the nested Wilson-loop method in Supplementary Note 3. Therefore, the Ta2Ni3Te5 monolayer is a QTI with a nontrivial second SW number.
Edge spectrum and corner states
From the orbital-resolved band structures (Fig. 3), the maximally localized Wannier functions of Ta-\({d}_{{z}^{2}}\), NiA-dxz and Te-px orbitals are extracted, to construct a 2D tight-binding (TB) model of Ta2Ni3Te5 monolayer. As shown in Fig. 3d, the obtained TB model fits the density functional theory (DFT) band structure well. First, we compute the (01)-edge spectrum with open boundary condition along y. Instead of gapless edge states for a 2D \({{\mathbb{Z}}}_{2}\)-nontrivial insulator, gapped edge states are obtained for the 2D QTI [Fig. 3e]. Then, we explore corner states as the hallmark of the 2D QTI. We compute the energy spectrum for a nanodisk. For concreteness, we take a rectangular-shaped nanodisk with 50 × 10 unit cells, preserving both Mx and My symmetries in the 0D geometry. The obtained discrete spectrum for this nanodisk is plotted in the inset of Fig. 3f. Remarkably, one observes four degenerate states near EF. The spatial distribution of these four-fold modes can be visualized from their charge distribution, as shown in Fig. 3f. Clearly, they are well localized at the four corners, corresponding to isolated corner states.
Minimum model for the 2D QTI
As shown in Fig. 2, the minimum model for the 2D QTI should be consisted of two BRs of \(A^{\prime} @4e\) and A″@4e. Based on the situation of Ta2Ni3Te5 monolayer in Fig. 2c, the minimum effective model is derived as below:
The terms of HTa(k), HNi(k) and Hint(k) are 4 × 4 matrices, which read
The γ1,2,3(k) matrices are given explicitly in Supplementary Note 4.
We find that ts1 < 0 and tp1, tp2 > 0 for the A2M1,3X5 family (ts3 and tp3 are small). When εs + 2ts1 − 2∣ts2∣ < εp + 2tp1 + 2tp2, the double-band inversion happens in the monolayers of this family. By fitting the DFT bands, we obtain ts2 > 0 for the 215-type, while ts2 < 0 for the 235-type. When εp = − εs and tpi = − tsi(i = 1, 2, 3), the model is chiral symmetric (i.e., δ = 0 in Table 1). Since the second SW insulator or QTI is topological in the presence of chiral symmetry, we would focus on the model (almost respecting chiral symmetry) in the following discussion.
Analytic solution of (01)-edge states
As the remnants of the QTI phase, the localized edge states can be solved analytically for the minimum model. For the (01)-edge, one can treat the model HTB(k) as two parts, H0(k) and \(H^{\prime} ({{{\boldsymbol{k}}}})\),
Note that there is a pair of Dirac points \((\pm {k}_{x}^{D},0)\), with \({k}_{x}^{D}=\arccos \left[\frac{1}{2}{(\frac{{t}_{s3}}{{t}_{s2}})}^{2}-1\right]\). Since kx is still a good quantum number on the (01)-edge, expanding ky to the second order, the zero-mode equation H0(kx, − i∂y)Ψ(kx, y) = 0 can be solved for y ∈ [0, +∞). Taking the trial solution of Ψ(kx, y) = ψ(kx)eλy, we obtain the secular equation and the solution of λ = ±λ±, where
With the boundary conditions Ψ(kx, 0) = Ψ(kx, +∞) = 0, only − λ± are permitted.
In the kx regime of \(\left[-{{{{\boldsymbol{k}}}}}_{x}^{D},{{{{\boldsymbol{k}}}}}_{x}^{D}\right]\), the edge zero-mode states are Ψ(kx, y) = [C1(kx)ϕ1(kx) + C2(kx)ϕ2(kx)]\(({e}^{-{\lambda }_{+}y}-{e}^{-{\lambda }_{-}y})\) with
The edge zero states are Fermi arcs that linking the pair of projected Dirac points \((\!\pm {k}_{x}^{D},0)\). Once \(H^{\prime} ({{{\boldsymbol{k}}}})\) included, the effective (01)-edge Hamiltonian is,
where \(\left|{{\Phi }}\right\rangle \equiv \left|{\phi }_{1}({k}_{x}),{\phi }_{2}({k}_{x})\right\rangle\). Two edge spectra are obtained in Fig. 4a.
Effective Su–Schrieffer–Heeger model on (10)-edge and corner states
Similarly, we derive the (10)-edge modes as \([{F}_{1}({k}_{y}){\varphi }_{1}({k}_{y})+{F}_{2}({k}_{y}){\varphi }_{2}({k}_{y})]\left({{\rm{e}}}^{-{{{\Lambda }}}_{+}x}-{{\rm{e}}}^{-{{{\Lambda }}}_{-}x}\right)\) with
Here, \({{{\Lambda }}}_{\pm }=\frac{2{t}_{sp}\pm \sqrt{4{t}_{sp}^{2}-2(2{t}_{s1}+{t}_{s2})(2{t}_{s1}+2{t}_{s2}+{\varepsilon }_{s})}}{2{t}_{s1}+{t}_{s2}}\), \({{{\Delta }}}_{s}={\left(\frac{{t}_{s3}}{2{t}_{s2}}\right)}^{2}\), and \({{{\Delta }}}_{p}={\left(\frac{{t}_{p3}}{2{t}_{p2}}\right)}^{2}\). Then we obtain the effective Hamiltonian on (10) edge below,
When δ = 0, the minimum QTI model is chiral symmetric and it is gapless on the (10) edge (preserving My symmetry). When the chiral symmetry is slightly broken (δ ≠ 0), the \({H}_{10}^{{\rm{eff}}}\) becomes an Su–Schrieffer–Heeger model (δ < 0 nontrivial; δ > 0 trivial), as presented in Fig. 4b–d. As a result, we obtain a solution state on the end of the edge mode, i.e., the corner. As long as the energy of the corner state is located in the gap of bulk and edge states, the corner state is well localized at the corners, as shown in Fig. 4c, d.
Discussion
In Ta2NiSe5 monolayer, the double-band inversion has also happened between Ta-\({d}_{{z}^{2}}\) and Te-px states, about 0.4 eV, resulting in a semimetal with a pair of nodal lines in the 215-type. The highest valence bands on Y–Γ are from the inverted Ta-\({d}_{{z}^{2}}\) states. However, in Ta2Ni3Te5 monolayer, the double-band inversion strength becomes remarkable, ~1 eV, which is ascribed to the filled B-type voids and more extended Te-p states. On the other hand, the highest valence bands become NiA-dxz states (slightly hybridized with Te-px states). It is insulating with a small gap of 65 meV. When it comes to A2Pd3Te5 monolayer, the remarkable inversion strength is similar to that of Ta2Ni3Te5. But the PdA-dxz states go upwards further due to more expansion of the d orbitals and the energy gap becomes almost zero in Ta2Pd3Te5. In short, the double-band inversion happens in all these monolayers, while the band gap of the 235-type changes from positive (Ta2Ni3Te5), to nearly zero (Ta2Pd3Te5 with a tiny band overlap), to negative (Nb2Pd3Te5), as shown in Supplementary Fig. 1. Although the band structure of Ta2Pd3Te5 bulk is metallic without SOC33,34,38, the monolayer could become a QSH insulator upon including SOC in ref. 33. Since their bulk materials are van der Waals layered compounds, the bulk topology and properties strongly rely on the band structures of the monolayers in the A2M1,3X5 family.
As we find in ref. 33, the band topology of Ta2Pd3Te5 monolayer is lattice sensitive. By applying >1% uniaxial compressive strain along b, it becomes a \({{\mathbb{Z}}}_{2}\)-trivial insulator, being a QTI. On the other hand, due to the quasi-1D crystal structure, the screening effect of carriers is relatively weak and the electron-hole Coulomb interaction may be substantial for exciton condensation. The 1D in-gap edge states as remnants of the QTI are responsible for the observed Luttinger-liquid behavior.
In conclusion, we predict that Ta2M3Te5 monolayers can be QTIs by solving aBR decomposition and computing SW numbers. Through aBR analysis, we conclude that the second-order topology comes from an essential BR at the empty site (Ag@4c), and it origins from the remarkable double-band inversion. The double-band inversion also happens in the band structure of Ta2NiSe5 monolayer. The second SW number of Ta2Ni3Te5 monolayer is w2 = 1, corresponding to a QTI. Therefore, we obtain edge states and corner states of the monolayer. The eight-band quadrupole model with Mx and My has been constructed successfully for electronic materials. With the large double-band inversion and small band energy gap/overlap, these transition-metal materials of A2M1,3X5 family provide a good platform to study the interplay between the topology and interactions (Fig. 4e).
Methods
Calculation method
Our first-principles calculations were performed within the framework of the DFT using the projector augmented wave method40,41, as implemented in Vienna ab-initio simulation package (VASP)42,43. The Perdew–Burke–Ernzerhof (PBE) generalized gradient approximation exchange-correlations functional44 was used. SOC was neglected in the calculations except in Supplementary Fig. 2b. We also used SCAN45 and GW46 method when checking the band gap. In the self-consistent process, 16 × 4 × 1 k-point sampling grids were used, and the cut-off energy for plane wave expansion was 500 eV. The irreps were obtained by the program IRVSP39. The maximally localized Wannier functions were constructed by using the Wannier90 package47. The edge spectra are calculated using surface Green’s function of semi-infinite system48,49.
Data availability
The data sets that support the findings in this study are available from the corresponding author upon request.
Code availability
The codes that support the findings in this study are available from the corresponding author upon request.
References
Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Quantized electric multipole insulators. Science 357, 61–66 (2017).
Schindler, F. et al. Higher-order topological insulators. Sci. Adv. 4, eaat0346 (2018).
Song, Z., Fang, Z. & Fang, C. (d-2)-dimensional edge states of rotation symmetry protected topological states. Phys. Rev. Lett. 119, 246402 (2017).
Langbehn, J., Peng, Y., Trifunovic, L., von Oppen, F. & Brouwer, P. W. Reflection-symmetric second-order topological insulators and superconductors. Phys. Rev. Lett. 119, 246401 (2017).
Schindler, F. et al. Higher-order topology in bismuth. Nat. Phys. 14, 918–924 (2018).
Ezawa, M. Higher-order topological insulators and semimetals on the breathing kagome and pyrochlore lattices. Phys. Rev. Lett. 120, 026801 (2018).
Khalaf, E. Higher-order topological insulators and superconductors protected by inversion symmetry. Phys. Rev. B 97, 205136 (2018).
Wang, Z., Wieder, B. J., Li, J., Yan, B. & Bernevig, B. A. Higher-order topology, monopole nodal lines, and the origin of large Fermi arcs in transition metal dichalcogenides XTe2 (X = Mo,W). Phys. Rev. Lett. 123, 186401 (2019).
Benalcazar, W. A., Li, T. & Hughes, T. L. Quantization of fractional corner charge in Cn-symmetric higher-order topological crystalline insulators. Phys. Rev. B 99, 245151 (2019).
Yue, C. et al. Symmetry-enforced chiral hinge states and surface quantum anomalous Hall effect in the magnetic axion insulator Bi2−xSmxSe3. Nat. Phys. 15, 577–581 (2019).
Xu, Y., Song, Z., Wang, Z., Weng, H. & Dai, X. Higher-order topology of the axion insulator EuIn2As2. Phys. Rev. Lett. 122, 256402 (2019).
Wieder, B. J. et al. Strong and fragile topological Dirac semimetals with higher-order Fermi arcs. Nat. Commun. 11, 627 (2020).
Bradlyn, B. et al. Topological quantum chemistry. Nature 547, 298–305 (2017).
Elcoro, L. et al. Magnetic topological quantum chemistry. Nat. Commun. 12, 5965 (2021).
Gao, J., Guo, Z., Weng, H. & Wang, Z. Magnetic band representations, Fu-Kane-like symmetry indicators, and magnetic topological materials. Phys. Rev. B 106, 035150 (2022).
Nie, S. et al. Application of topological quantum chemistry in electrides. Phys. Rev. B 103, 205133 (2021).
Nie, S., Bernevig, B. A. & Wang, Z. Sixfold excitations in electrides. Phys. Rev. Res. 3, L012028 (2021).
Gao, J. et al. Unconventional materials: the mismatch between electronic charge centers and atomic positions. Sci. Bull. 67, 598–608 (2022).
Fang, C., Chen, Y., Kee, H.-Y. & Fu, L. Topological nodal line semimetals with and without spin-orbital coupling. Phys. Rev. B 92, 081201 (2015).
Zhao, Y. X., Schnyder, A. P. & Wang, Z. D. Unified theory of PT and CP invariant topological metals and nodal superconductors. Phys. Rev. Lett. 116, 156402 (2016).
Zhao, Y. X. & Lu, Y. PT-symmetric real Dirac fermions and semimetals. Phys. Rev. Lett. 118, 056401 (2017).
Ahn, J., Park, S., Kim, D., Kim, Y. & Yang, B.-J. Stiefel-Whitney classes and topological phases in band theory. Chin. Phys. B 28, 117101 (2019).
Ahn, J., Park, S. & Yang, B.-J. Failure of Nielsen-Ninomiya theorem and fragile topology in two-dimensional systems with space-time inversion symmetry: Application to twisted bilayer graphene at magic angle. Phys. Rev. X 9, 021013 (2019).
Wu, Q., Soluyanov, A. A. & Bzdusek, T. Non-abelian band topology in noninteracting metals. Science 365, 1273–1277 (2019).
Pan, M., Li, D., Fan, J. & Huang, H. Two-dimensional Stiefel-Whitney insulators in liganded Xenes. NPJ Comput. Mater. 8, 1–6 (2022).
Xu, Y. et al. Three-dimensional real space invariants, obstructed atomic insulators and a new principle for active catalytic sites. Preprint at http://arxiv.org/abs/2111.02433 (2021).
Wieder, B. J. & Bernevig, B. A. The axion insulator as a pump of fragile topology. Preprint at http://arxiv.org/abs/1810.02373 (2018).
Benalcazar, W. A., Li, T. & Hughes, T. L. Quantization of fractional corner charge in Cn-symmetric higher-order topological crystalline insulators. Phys. Rev. B 99, 245151 (2019).
Sheng, X.-L. et al. Two-dimensional second-order topological insulator in graphdiyne. Phys. Rev. Lett. 123, 256402 (2019).
Lee, E., Kim, R., Ahn, J. & Yang, B.-J. Two-dimensional higher-order topology in monolayer graphdiyne. NPJ Quantum Mater. 5, 1–7 (2020).
Qian, S., Liu, C.-C. & Yao, Y. Second-order topological insulator state in hexagonal lattices and its abundant material candidates. Phys. Rev. B 104, 245427 (2021).
Chen, C. et al. Universal approach to magnetic second-order topological insulator. Phys. Rev. Lett. 125, 056402 (2020).
Guo, Z. et al. Quantum spin Hall effect in Ta2M3Te5 (M = Pd, Ni). Phys. Rev. B 103, 115145 (2021).
Wang, X. et al. Observation of topological edge states in the quantum spin Hall insulator Ta2Pd3Te5. Phys. Rev. B 104, L241408 (2021).
Wakisaka, Y. et al. Excitonic insulator state in Ta2NiSe5 probed by photoemission spectroscopy. Phys. Rev. Lett. 103, 026402 (2009).
Lu, Y. F. et al. Zero-gap semiconductor to excitonic insulator transition in Ta2NiSe5. Nat. Commun. 8, 14408 (2017).
Mazza, G. et al. Nature of symmetry breaking at the excitonic insulator transition: Ta2NiSe5. Phys. Rev. Lett. 124, 197601 (2020).
Higashihara, N. et al. Superconductivity in Nb2Pd3Te5 and chemically-doped Ta2Pd3Te5. J. Phys. Soc. Jpn. 90, 063705 (2021).
Gao, J. C., Wu, Q. S., Persson, C. & Wang, Z. J. Irvsp: to obtain irreducible representations of electronic states in the VASP. Comp. Phys. Commun. 261, 107760 (2021).
Blochl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953 (1994).
Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 59, 1758 (1999).
Kresse, G. & Furthmuller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comp. Mater. Sci. 6, 15–50 (1996).
Kresse, G. & Furthmuller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 54, 11169 (1996).
Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).
Sun, J., Ruzsinszky, A. & Perdew, J. Strongly constrained and appropriately normed semilocal density functional. Phys. Rev. Lett. 115, 036402 (2015).
Shishkin, M. & Kresse, G. Implementation and performance of the frequency-dependent GW method within the PAW framework. Phys. Rev. B 74, 035101 (2006).
Pizzi, G. et al. Wannier90 as a community code: new features and applications. J. Phys. Condens. Matter 32, 165902 (2020).
Sancho, M. P. L., Sancho, J. M. L. & Rubio, J. Quick iterative scheme for the calculation of transfer-matrices - application to Mo(100). J. Phys. F: Met. Phys. 14, 1205–1215 (1984).
Sancho, M. P. L., Sancho, J. M. L. & Rubio, J. Highly convergent schemes for the calculation of bulk and surface Green-functions. J. Phys. F: Met. Phys. 15, 851–858 (1985).
Acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant Nos. 11974395 and 12188101), the Strategic Priority Research Program of Chinese Academy of Sciences (Grant No. XDB33000000), the China Postdoctoral Science Foundation funded project (Grant No. 2021M703461), and the Center for Materials Genome.
Author information
Authors and Affiliations
Contributions
Z.W. conceived and supervised the work. Z.G. performed the first-principles calculation. J.D. and Y.X. contributed to the analytical modeling. All authors contributed to the discussion and wrote 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.
Supplementary information
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
Guo, Z., Deng, J., Xie, Y. et al. Quadrupole topological insulators in Ta2M3Te5 (M = Ni, Pd) monolayers. npj Quantum Mater. 7, 87 (2022). https://doi.org/10.1038/s41535-022-00498-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41535-022-00498-8
This article is cited by
-
A robust and tunable Luttinger liquid in correlated edge of transition-metal second-order topological insulator Ta2Pd3Te5
Nature Communications (2023)
-
Orbital degree of freedom induced multiple sets of second-order topological states in two-dimensional breathing Kagome crystals
npj Quantum Materials (2023)