Abstract
Topological lattice defects, such as dislocations and grain boundaries (GBs), are ubiquitously present in the bulk of quantum materials and externally tunable in metamaterials. In terms of robust modes, localized near the defect cores, they are instrumental in identifying topological crystals, featuring the hallmark band inversion at a finite momentum (translationally active type). Here we show that the GB superlattices in both two-dimensional and three-dimensional translationally active higher-order topological insulators harbor a myriad of dispersive modes that are typically placed at finite energies, but always well-separated from the bulk states. However, when the Burgers vector of the constituting edge dislocations points toward the gapless corners or hinges, both second-order and third-order topological insulators accommodate self-organized emergent topological metals near the zero energy (half-filling) in the GB mini Brillouin zone. We discuss possible material platforms where our proposed scenarios can be realized through the band-structure and defect engineering.
Similar content being viewed by others
Introduction
Topological lattice defects are ubiquitous in crystalline materials and are of central importance for their structural properties. In the last decade, they have emerged as viable platforms for probing topological phases of matter through a subtle interplay of the topology of the electronic wavefunctions and lattice geometry. In particular, dislocations, the defects associated with lattice translations, can probe a wide range of topological crystals featuring the band-inversion at a finite momentum directly in the bulk via the symmetry and topology protected localized defect modes either at zero1,2,3,4,5,6,7,8,9,10,11,12 or finite13 energy. As such, these modes are immune to interface contamination and surface termination, and have been experimentally observed in both topological crystals14,15 and their metamaterial analogues16,17.
These developments boosted the exploration of extended lattice defects on topological platforms, among which grain boundaries (GBs) are the most prominent ones. A GB develops at the interface between two misoriented crystalline grains due to the accumulated elastic stress, and at low angles, it consists of an array of dislocations18,19. See Fig. 1. In turn, by virtue of the hybridization between the localized zero-energy dislocation modes, the GBs can host a wide range of quantum phases in both static20,21 and dynamic22 settings, when a parent topological insulator (TI) is first-order in nature, featuring gapless modes on the edges or surfaces.
In higher-order topological insulators (HOTIs), accommodating localized states on lower-dimensional boundaries, such as hinges and corners23,24,25,26,27,28,29,30,31,32,33, the dislocation modes on the other hand typically move to finite energies, controlled by the relative orientation between its Burgers vector (\(\textbf{b}\)) and the axis of inversion (domain wall) of the discrete symmetry breaking Wilson-Dirac (WD) mass, responsible for the higher-order topology13. Such a nontrivial interplay between the real-space geometry of lattice defects and the momentum-space topology of WD mass propels the current pursuit to unveil its signatures on the emergent electronic bands along the GBs in HOTIs.
Key results
We show that the GB superlattice in both two-dimensional (2D) and three-dimensional (3D) HOTIs harbors a variety of dispersive bands, which are typically placed at finite energies. However, as the Burgers vector of the constituting edge dislocations and the domain wall directions of the WD mass approach each other these dispersive bands come closer, hybridize and finally touch when the \(\textbf{b}\)-vector points toward gapless corners or hinges. The GB then fosters an emergent self-organized metallicity near the zero energy, stemming from an interplay of lattice geometry and momentum space topology. These outcomes are anchored from the Fourier transformation of the defect modes in the GB mini Brillouin zone (BZ), their local density of states (LDOS) in real space lattice with GB defects and the spectral flow of the defect modes with varying WD mass domain wall in 2D (Fig. 2) and 3D (Fig. 3) second-order TIs, and 3D third-order TI (Fig. 4). Our results are thus consequential for the quantum crystals ubiquitously hosting GB defects34, and metamaterials, where such extended defects can be engineered externally35,36,37,38.
Lattice model
The universal Bloch Hamiltonian for HOTIs in two (\(D=2\)) and three (\(D=3\)) dimensions can be split as \(\hat{h}_\text{HOTI} = \hat{h}_{1} + \hat{h}_{\Delta }\). The Hamiltonian for the parent first-order TIs takes the form
where \(d_i(\textbf{k})=t \sin (k_i a)\), t is the hopping amplitude set to be unity, and a is the lattice spacing. The first-order WD mass, preserving all non-spatial and crystal symmetries, and featuring band inversion (and thus TIs) within the parameter regime \(0< \Delta _1/B<4 D\), reads
A tower of HOTIs can now be constructed by adding discrete symmetry breaking WD masses (\(\hat{h}_\Delta\)) to \(\hat{h}_1\), which can be decomposed as \(\hat{h}_\Delta =\hat{h}_{2} + \hat{h}_{3}\). The second-order WD mass in \(D=2\) and \(D=3\) takes the general form13
where \(0 \le \theta \le \pi /2\) (about which more in a moment), \(d^\text{lat}_{x^2-y^2}=\cos (k_x a) - \cos (k_y a)\), \(d^\text{lat}_{xy}=\sin (k_x a)\sin (k_y a)\). The third-order WD mass (only in \(D \ge 3\)) reads as39
Here \(\Gamma _1, \ldots ,\Gamma _{D+2},\Gamma _{D+3}\) are mutually anticommuting Hermitian matrices each of which squares to unity. Results are independent of their explicit representation, which are shown in the Supplementary Information.
The Bloch Hamiltonian \(\hat{h}_1 + \hat{h}_2\) describes a second-order TI in both \(D=2\) and 3, since \(\hat{h}_{2}\) gaps out the topological edge (surface) states accommodated by the first-order phase, and leaves only four corners (four z-directional hinges and two xy surfaces) gapless in D=2 (D=3). As such, for each value of the parameter \(\theta\), \(\hat{h}_2\) features a domain wall, which lies along the principal axes \(k_x=0\) and \(k_y=0\) (the diagonals \(k_y=\pm k_x\)) for \(\theta =\pi /2\) (\(\theta =0\)). Sharp corner or hinge modes then appear only when they lie on the axes of inversion (domain wall) for the WD mass. For additional details on the role of \(\theta\), see Ref.13.
Since \(\hat{h}_1 + \hat{h}_2\) involves four (five) mutually anticommuting \(\Gamma\) matrices in \(D=2\) (\(D=3\)), their dimensionality is four. Consequently, the corner modes are pinned at zero energy due to a unitary (generated by \(\Gamma _5\)) and an antiunitary particle-hole (PH) symmetry in \(D=2\), while in \(D=3\) only an antiunitary PH symmetry pins hinge modes to zero energy40, as the maximal number of mutually anticommuting four-dimensional Hermitian \(\Gamma\) matrices is five. Notice that \(\hat{h}_1 + \hat{h}_2\) also enjoys the antiunitary composite \(C_4 {\mathcal T}\) symmetry, a product of the four-fold rotation about the z-axis (\(C_4\)), generated by \(i \Gamma _1 \Gamma _2\) and under which \((k_x,k_y) \rightarrow (-k_y,k_x)\), and the time-reversal (\({\mathcal T}\)) in both \(D=2\) and \(D=3\). The explicit forms of the antiunitary PH and \({\mathcal T}\) symmetry generators however depend on the \(\Gamma\) matrix representation. The model Hamiltonian \(\hat{h}_1 + \hat{h}_2\) also breaks the unitary parity (inversion) symmetry (\({\mathcal P}\)) under which \(\varvec{k} \rightarrow -\varvec{k}\). But, its explicit form is dimension and \(\Gamma\) matrix representation dependent. Therefore, 2D and 3D second-order TIs additionally preserve composite \(C_4 {\mathcal {P}}\) and \({\mathcal {P}} {\mathcal {T}}\) symmetries.
When we turn on the mass term \(\hat{h}_3\), the Hamiltonian \(\hat{h}_1+\hat{h}_2+\hat{h}_3\) (with finite \(\Delta _2\) and \(\Delta _3\)) describes a third-order topological insulator in \(D=3\). Since it requires six mutually anticommuting \(\Gamma\) matrices, they must therefore be at least eight-dimensional. The additional mass term (\(\hat{h}_3\)) gaps out otherwise gapless hinge and xy (top and bottom) surface states, yielding the eight zero-energy corner-localized modes in the cubic geometry. For the sharp corner localization, eight corners of the cubic lattice must coincide with the directions along which \(\hat{h}_2=0=\hat{h}_3\)41. The corner modes in a third-order TI are pinned at zero energy due to a unitary (generated by \(\Gamma _7\)) as well as an antiunitary PH symmetry. The 3D third-order TI also breaks the individual \({\mathcal T}\), \(C_4\), and \({\mathcal P}\) symmetries, and preserves the composite \(C_4 {\mathcal T}\), \(C_4 {\mathcal P}\), and \({\mathcal P} {\mathcal T}\) symmetries.
Since HOTIs are obtained as descendants of the first-order TIs upon systematically switching on the discrete symmetry breaking WD masses (\(\hat{h}_2\) and \(\hat{h}_3\)), in what follows, we consider the translationally active \(\textrm{M}\) and \(\textrm{R}\) phases of the 2D and 3D first-order TIs, respectively, featuring band inversion at finite momentum (\(\textbf{K}_\text{inv}\)) \(\textrm{M}=(1,1)\pi /a\) and \(\textrm{R}=(1,1,1)\pi /a\) points in the corresponding BZ of the parent square and cubic lattices. For numerical calculations in \(D=2\) and \(D=3\), we throughout set \(\Delta _1/B=6\) and \(\Delta _1/B=11\), respectively [Eq. (2)].
GB: construction
When introduced in a TI, the elementary building block of a GB defect, a single dislocation, sources an effective hopping phase across the defect \(\Phi _\text{dis}=\textbf{K}_\text{inv}\cdot \textbf{b}\,\,(\textrm{mod}\,\, 2\pi )\)1,3. Therefore, electrons encircling the defect can pick up a nontrivial hopping phase \(\Phi _\text{dis}=\pi\) in the \(\textrm{M}\) (\(\textrm{R}\)) phase on the square (cubic) lattice, while \(\Phi _\text{dis}=0\) always in the \(\Gamma\) phase as \(\textbf{K}_\text{inv}=0\) therein. As a result, in a translationally-active first-order phase, a single dislocation hosts localized modes at zero energy. By contrast, in HOTIs the additional higher-order WD masses typically gap out the dislocation modes and place them at finite energies, unless the Burgers vector pierces a gapless lower-dimensional boundary, while throughout retaining their topological and symmetry protection13.
When a GB is immersed in a parent translationally-active first-order TI, the localized zero modes at individual dislocations hybridize, giving rise to an emergent topological metal near the half-filling, or zero energy, along the GB superlattice20. The emergence of such a self-organized topological metal can be corroborated from the Fourier transform of the GB modes with respect to the superlattice periodicity d, yielding a band structure in the superlattice mini BZ. When the dislocation modes hybridize, they also develop a comparable weight at the middle of two successive defect cores. Thus, we denote the distance between them by 2d (Fig. 1). Next, we report the nature of emergent metallic states in the GB defect consisting of edge dislocations in both 2D (Fig. 2) and 3D (Fig. 3) second-order TIs, and 3D third-order TI (Fig. 4), showcasing an intriguing interplay between the defect geometry and the topology of discrete symmetry breaking WD masses (\(\hat{h}_2\) and \(\hat{h}_3\)). For concreteness, throughout the Burgers vector of constituting dislocation in the GB is taken to be \(\textbf{b}=a\textbf{e}_x\). Here we consider only small angle GBs, characterized by \(\varphi =\sin ^{-1}(a/(2d)) \approx 14.47^\circ\) (\(< 15^\circ\)) for \(d=2a\).
GB in 2D second-order TI
The evolution of the GB band structure as the orientation (\(\theta\)) of the WD mass domain wall is varied with respect to \(\textbf{b}\) is shown in Fig. 2a. When \(\theta =0\), the axes of inversion for the second-order WD mass (\(\hat{h}_2\)) is perpendicular to \(\textbf{b}\), and a single dislocation features PH symmetric localized modes at finite energies13. Hybridization among them then yields isolated dispersive bands centered at finite energies along the GB superlattice, separated by a gap G (Fig. 2b). It measures the energy difference between the bottom and top of such dispersive bands living within the conduction and valence bands, respectively. As \(\theta \rightarrow \pi /2\), these bands come close to each other and start to hybridize. For \(\theta \approx \pi /2\), when the domain wall of \(\hat{h}_2\) is parallel to \(\textbf{b}\), this gap eventually closes, yielding an emergent metal around the zero energy. For any \(\theta\), the dispersive or metallic states are always separated from the bulk states, like their parent dislocation modes. For \(\theta =\pi /2\), besides the dispersive metallic states a pair of zero-energy modes also appear in the spectra (Fig. 2a and c). All of them are delocalized along the GB, but confined in its close vicinity, as shown from their LDOS in Fig. 2d. The near-zero energy modes, as well as the rest of the GB modes, do not show any spectral weight in the region between the GB and anti-GB, separated by a large distance of 10a, such that the modes bound to them do not overlap with each other.
GB in 3D second-order TI
In \(D=3\), a GB consists of an array of edge dislocations, which forms a plane spanned by the directions of the dislocation line and the dislocation array, respectively in the z-axis and y-axis for \(\textbf{b}=a \textbf{e}_x\). Each dislocation line then hosts localized modes along the z direction, which are typically at finite energies for \(\theta \ne \pi /2\), unless \(\textbf{b}\) points toward gapless hinges (\(\theta =\pi /2\)), when the dislocation modes become gapless13. These modes, when hybridize, yield a plethora of 2D dispersive bands along the GB defect, manifesting a gap depending on the parameter \(\theta\), as shown in the upper panel of Fig. 3a, with the metallic band structure centered around the zero energy emerging for \(\theta =\pi /2\). Furthermore, as shown in the lower panel of Fig. 3a, the closest to zero-energy modes feature almost the entire spectral weight at momenta \((\pm \pi /(2d),\pi /a)\), which is independent of the parameter \(\theta\). The emergent metallicity around the zero energy is further corroborated from the \(\theta\)-dependence of the gap (G) between bottom and top of the dispersive bands respectively residing within the conduction and valence bands, showing that G goes to zero as \(\theta \rightarrow \pi /2\) (Fig. 3b). The GB modes are always well separated from the bulk states for any \(\theta\) (Fig. 3a), which we also explicitly show for \(\theta =\pi /2\) in Fig. 3c. Finally, these modes are highly localized in the yz-plane constituted by the GB (for any \(\theta\)), as explicitly shown from their LDOS for \(\theta =\pi /2\) in Fig. 3d. Appearance of two gapless Dirac points at \((\pm \pi /(2d),\pi /a)\) when \(\theta =\pi /2\) in a 2D GB mini BZ conforms to the Nielsen-Ninomiya Fermion doubling theorem42.
GB in 3D third-order TI
An edge dislocation in a 3D third-order TI harbors modes localized near its ends on the top and bottom surfaces, for example, when \(\textbf{b}=a \textbf{e}_x\), which are at finite energies, unless the Burgers vector points toward eight gapless corners, when they become gapless13. Consequently, a GB in a 3D third-order TI is expected to host gapped dispersive bands localized near the top and the bottom edges of the defect, which become gapless only when \(\theta =\pi /2\). Indeed Fig. 4a confirms it, showing that the gap between the dispersive GB bands decreases as the WD mass domain wall approaches the principal crystallographic axes (\(\theta \rightarrow \pi /2\)), when it should become gapless. A small residual gap between these dispersive modes at \(\theta =\pi /2\) is purely due to a finite thickness of the system in the z-direction (\(L_z\)), which approaches zero as \(L_z \rightarrow \infty\). See Fig. 4b. The GB modes (gapped dispersive or gapless metallic about the half-filling) are always well-separated from the bulk states (Fig. 4a), as explicitly shown in Fig. 4c and localized near the top and bottom surfaces, as shown in Fig. 4d for \(\theta =\pi /2\). The fact that these modes are maximally localized just below (above) the top (bottom) surfaces possibly stems from their parent corner modes, also localized slightly away from the terminal surfaces.
Discussion and outlook
Here we show that 2D and 3D HOTIs foster a variety of dispersive bands within the conduction and valence bands of the emergent BZ constituted by the GB superlattice, with the gap between them tunable by the relative orientation of the WD mass domain wall and the Burgers vector of individual dislocations. Especially, when the Burgers vector pierces gapless corners or hinges, these dispersive bands touch each other at zero energy, giving birth to self-organized topological metals around zero energy confined to the GB defect. See Figs. 2, 3, 4. The GB modes (dispersive and metallic) are robust against weak on-site disorder, the dominant source of elastic scattering in any real material. Explicit results are shown in the Supplementary Information.
Note that the bulk gap is determined by the first-order [Eq. (2)] and higher-order [Eqs. (3) and (4)] masses, while the energy scale of the GB modes is set by the latter ones and \(\theta\) [Eq. (3)]. As all the \(\Gamma\) matrices appearing in the universal model Hamiltonian for HOTIs mutually anticommute with each other, there exists a finite energy separation between the bulk and GB modes (\(\Delta _\mathrm{bulk-GB}\)). However, as \(\Delta _2\) and/or \(\Delta _3\) become sufficiently large, the higher-order mass overwhelms the first-order mass, and \(\Delta _\mathrm{bulk-GB} \rightarrow 0\), but it never vanishes. The explicit dependence of \(\Delta _\mathrm{bulk-GB}\) on the amplitude of higher-order masses is shown in the Supplementary Information.
While here we focus on small angle GBs, described by an array of dislocations (Fig. 1), at large opening angles (\(\varphi > 15^\circ\)) a GB defect may be described as an array of elementary disclination defects43 or in terms of a lattice of partial dislocations with stacking faults44. As a single disclination45,46 and a partial dislocation47 defect can also harbor localized topological modes, the possibility of emergent metallicity on high-angle GBs stands as a fascinating avenue for future investigation.
The prerequisite for the realization of the proposed emergent metallicity along the GB in HOTIs is a band-inversion at a finite momentum in the BZ (translationally active3,5) such that \(\textbf{K}_\text{inv} \cdot \textbf{b}=\pi\) (modulo \(2\pi\)), as GB defects are rather ubiquitous in quantum crystals. We emphasize that a metallic behavior has been recently observed at GBs in 1T’-MoTe\(_2\)34, while some of the realized HOTIs are of translationally-active type48,49, which should motivate further investigation of the materials prospects for the realization of our proposal. On the other hand, in metamaterials the predicted GB band structures can be engineered by artificially tuning the tunneling processes and manipulating defects therein. Among them designer metamaterials35, topolectric circuits36,37, photonic38 and mechanical50 lattices are the most promising platforms. While the metallic nature of the GB modes can be probed via electronic transport measurements in quantum and designer crystals, classical metamaterials can only reveal their dispersive nature from the energy-conserved momentum relation of the associated classical modes, such as the vibrational ones in mechanical lattices. By contrast, the local density of states of the GB modes can be probed in quantum and designer crystals via scanning tunneling spectroscopy, as well as in classical metamaterials, through the measurements of local electric (in topolectric circuits) or mechanical (in mechanical lattice) impedance or two-point pump-probe spectroscopy (in photonic lattices). Our findings should therefore motivate experimental efforts for the realization of the defect-based emergent band structures in a wide range of both quantum and classical topological materials.
Data availability
The datasets used and/or analysed during the current study available from the corresponding authors on reasonable request. Main codes and the data for generating the figures presented in the main text and Supplementary Information are already available at https://doi.org/10.5281/zenodo.8341312.
References
Ran, Y., Zhang, Y. & Vishwanath, A. One-dimensional topologically protected modes in topological insulators with lattice dislocations. Nat. Phys. 5, 298. https://doi.org/10.1038/nphys1220 (2009).
Teo, J. C. Y. & Kane, C. L. Topological defects and gapless modes in insulators and superconductors. Phys. Rev. B 82, 115120. https://doi.org/10.1103/PhysRevB.82.115120 (2010).
Juričić, V., Mesaros, A., Slager, R.-J. & Zaanen, J. Universal probes of two-dimensional topological insulators: Dislocation and \(\pi\) flux. Phys. Rev. Lett. 108, 106403. https://doi.org/10.1103/PhysRevLett.108.106403 (2012).
Asahi, D. & Nagaosa, N. Topological indices, defects, and Majorana fermions in chiral superconductors. Phys. Rev. B 86, 100504. https://doi.org/10.1103/PhysRevB.86.100504 (2012).
Slager, R.-J., Mesaros, A., Juričić, V. & Zaanen, J. The space group classification of topological band-insulators. Nat. Phys. 9, 98. https://doi.org/10.1038/nphys2513 (2013).
Hughes, T. L., Yao, H. & Qi, X.-L. Majorana zero modes in dislocations of Sr2RuO4. Phys. Rev. B 90, 235123. https://doi.org/10.1103/PhysRevB.90.235123 (2014).
Slager, R.-J., Mesaros, A., Juričić, V. & Zaanen, J. Interplay between electronic topology and crystal symmetry: Dislocation-line modes in topological band insulators. Phys. Rev. B 90, 241403. https://doi.org/10.1103/PhysRevB.90.241403 (2014).
Nag, T. & Roy, B. Anomalous and normal dislocation modes in Floquet topological insulators. Commun. Phys. 4, 157. https://doi.org/10.1038/s42005-021-00659-4 (2021).
S.K. Das & B. Roy, Dynamic melting and condensation of topological dislocation modes, arXiv:2210.15661.
Panigrahi, A., Moessner, R. & Roy, B. Non-Hermitian dislocation modes: Stability and melting across exceptional points. Phys. Rev. B 106, L041302. https://doi.org/10.1103/PhysRevB.106.L041302 (2022).
Panigrahi, A., Juričić, V. & Roy, B. Projected topological branes. Commun. Phys. 5, 230. https://doi.org/10.1038/s42005-022-01006-x (2022).
L.-H. Hu & R.-X. Zhang, Dislocation Majorana Bound States in Iron-based Superconductors, arXiv:2207.10113.
Roy, B. & Juričić, V. Dislocation as a bulk probe of higher-order topological insulators. Phys. Rev. Res. 3, 033107. https://doi.org/10.1103/PhysRevResearch.3.033107 (2021).
Hamasaki, H., Tokumoto, Y. & Edagawa, K. Dislocation conduction in Bi-Sb topological insulators. Appl. Phys. Lett. 110, 092105. https://doi.org/10.1063/1.4977839 (2017).
Nayak, A. K. et al. Resolving the topological classification of bismuth with topological defects. Sci. Adv. 5, eaax6996. https://doi.org/10.1126/sciadv.aax6996 (2019).
Xue, H. et al. Observation of dislocation-induced topological modes in a three-dimensional acoustic topological insulator. Phys. Rev. Lett. 127, 214301. https://doi.org/10.1103/PhysRevLett.127.214301 (2021).
Ye, L. et al. Topological dislocation modes in three-dimensional acoustic topological insulators. Nat. Commun. 13, 508. https://doi.org/10.1038/s41467-022-28182-2 (2022).
Sutton, A. & Balluffi, R. Interfaces in Crystalline Materials (Oxford University Press, UK, 2007).
Han, J., Thomas, S. L. & Srolovitz, D. J. Grain-boundary kinetics: A unified approach. Prog. Mater. Sci. 98, 386. https://doi.org/10.1016/j.pmatsci.2018.05.004 (2018).
Slager, R.-J., Juričić, V., Lahtinen, V. & Zaanen, J. Self-organized pseudo-graphene on grain boundaries in topological band insulators. Phys. Rev. B 93, 245406. https://doi.org/10.1103/PhysRevB.93.245406 (2016).
Amundsen, M. & Juričić, V. Grain-boundary topological superconductor. Commun. Phys. 6, 232. https://doi.org/10.1038/s42005-023-01351-5 (2023).
D. J. Salib & B. Roy, Dynamic metal along grain boundaries in Floquet topological crystals, arXiv:2212.08060.
Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Quantized electric multipole insulators. Science 357, 61. https://doi.org/10.1126/science.aah6442 (2017).
Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Electric multipole moments, topological multipole moment pumping, and chiral hinge states in crystalline insulators. Phys. Rev. B 96, 245115. https://doi.org/10.1103/PhysRevB.96.245115 (2017).
Song, Z., Fang, Z. & Fang, C. \((d-2)\)-dimensional edge states of rotation symmetry protected topological states. Phys. Rev. Lett. 119, 246402. https://doi.org/10.1103/PhysRevLett.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. https://doi.org/10.1103/PhysRevLett.119.246401 (2017).
Schindler, F. et al. Higher-order topological insulators. Sci. Adv. 4, eaat0346. https://doi.org/10.1126/sciadv.aat0346 (2018).
Khalaf, E. Higher-order topological insulators and superconductors protected by inversion symmetry. Phys. Rev. B 97, 205136. https://doi.org/10.1103/PhysRevB.97.205136 (2018).
Hsu, C.-H., Stano, P., Klinovaja, J. & Loss, D. Majorana Kramers pairs in higher-order topological insulators. Phys. Rev. Lett. 121, 196801. https://doi.org/10.1103/PhysRevLett.121.196801 (2018).
Matsugatani, A. & Watanabe, H. Connecting higher-order topological insulators to lower-dimensional topological insulators. Phys. Rev. B 98, 205129. https://doi.org/10.1103/PhysRevB.98.205129 (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 \(X{\rm Te}_{2}\) (\(X={\rm Mo},{\rm W}\)). Phys. Rev. Lett. 123, 186401. https://doi.org/10.1103/PhysRevLett.123.186401 (2019).
Trifunovic, L. & Brouwer, P. W. Higher-order bulk-boundary correspondence for topological crystalline phases. Phys. Rev. X 9, 011012. https://doi.org/10.1103/PhysRevX.9.011012 (2019).
Călugăru, D., Juričić, V. & Roy, B. Higher-order topological phases: A general principle of construction. Phys. Rev. B 99, 041301. https://doi.org/10.1103/PhysRevB.99.041301 (2019).
Kim, H. W. et al. Symmetry dictated grain boundary state in a two-dimensional topological insulator. Nano Lett. 20, 5837. https://doi.org/10.1021/acs.nanolett.0c01756 (2020).
Kempkes, S. N. et al. Robust zero-energy modes in an electronic higher-order topological insulator. Nat. Mater. 18, 1292. https://doi.org/10.1038/s41563-019-0483-4 (2019).
Imhof, S. et al. Topolectrical-circuit realization of topological corner modes. Nat. Phys. 14, 925. https://doi.org/10.1038/s41567-018-0246-1 (2018).
Dong, J., Juričić, V. & Roy, B. Topolectric circuits: Theory and construction. Phys. Rev. Res. 3, 023056. https://doi.org/10.1103/PhysRevResearch.3.023056 (2021).
Li, F.-F. et al. Topological light-trapping on a dislocation. Nat. Commun. 9, 2462. https://doi.org/10.1038/s41467-018-04861-x (2018).
Nag, T., Juričić, V. & Roy, B. Hierarchy of higher-order Floquet topological phases in three dimensions. Phys. Rev. B 103, 115308. https://doi.org/10.1103/PhysRevB.103.115308 (2021).
Roy, B. Antiunitary symmetry protected higher-order topological phases. Phys. Rev. Res. 1, 032048. https://doi.org/10.1103/PhysRevResearch.1.032048 (2019).
Roy, B. & Juričić, V. Mixed-parity octupolar pairing and corner Majorana modes in three dimensions. Phys. Rev. B 104, L180503. https://doi.org/10.1103/PhysRevB.104.L180503 (2021).
Nielsen, H. & Ninomiya, M. Absence of neutrinos on a lattice: (I). Proof by homotopy theory. Nucl. Phys. B 185, 20. https://doi.org/10.1016/0550-3213(81)90361-8 (1981).
Nazarov, A. A., Shenderova, O. A. & Brenner, D. W. Elastic models of symmetrical \(\langle 001 \rangle\) and \(\langle 011 \rangle\) tilt grain boundaries in diamond. Phys. Rev. B 61, 928. https://doi.org/10.1103/PhysRevB.61.928 (2000).
King, A. H. & Chen, F.-R. Interactions between lattice partial dislocations and grain boundaries. Mater. Sci. Eng. 66, 227. https://doi.org/10.1016/0025-5416(84)90183-6 (1984).
Benalcazar, W. A., Li, T. & Hughes, T. L. Quantization of fractional corner charge in \({C}_{n}\)-symmetric higher-order topological crystalline insulators. Phys. Rev. B 99, 245151. https://doi.org/10.1103/PhysRevB.99.245151 (2019).
Geier, M., Fulga, I. C. & Lau, A. Bulk-boundary-defect correspondence at disclinations in rotation-symmetric topological insulators and superconductors. SciPost Phys. 10, 092. https://doi.org/10.21468/SciPostPhys.10.4.092 (2021).
Queiroz, R., Fulga, I. C., Avraham, N., Beidenkopf, H. & Cano, J. Partial lattice defects in higher-order topological insulators. Phys. Rev. Lett. 123, 266802. https://doi.org/10.1103/PhysRevLett.123.266802 (2019).
Schindler, F. et al. Higher-order topology in bismuth. Nat. Phys. 14, 918. https://doi.org/10.1038/s41567-018-0224-7 (2018).
Noguchi, R. et al. Evidence for a higher-order topological insulator in a three-dimensional material built from van der Waals stacking of bismuth-halide chains. Nat. Mater. 20, 473. https://doi.org/10.1038/s41563-020-00871-7 (2021).
Grinberg, I. H., Lin, M., Benalcazar, W. A., Hughes, T. L. & Bahl, G. Trapped state at a dislocation in a weak magnetomechanical topological insulator. Phys. Rev. Appl. 14, 064042. https://doi.org/10.1103/PhysRevApplied.14.064042 (2020).
Acknowledgements
D.J.S. and B.R. were supported by NSF CAREER Grant No. DMR- 2238679 of B.R. V.J. acknowledges support of the Swedish Research Council (VR 2019-04735) and Fondecyt (Chile), Grant No. 1230933. Nordita is partially supported by Nordforsk.
Funding
Open access funding provided by Stockholm University.
Author information
Authors and Affiliations
Contributions
D.J.S. performed all the numerical calculations. V.J. and B.R. conceived and structured the project, and wrote the manuscript. B.R. supervised the project.
Corresponding authors
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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Salib, D.J., Juričić, V. & Roy, B. Emergent metallicity at the grain boundaries of higher-order topological insulators. Sci Rep 13, 15308 (2023). https://doi.org/10.1038/s41598-023-42279-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-42279-8
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.