Abstract
We discover three-dimensional intertwined Weyl phases, by developing a theory to create topological phases. The theory is based on intertwining existing topological gapped and gapless phases protected by the same crystalline symmetry. The intertwined Weyl phases feature both unconventional Weyl semimetallic (monopole charge>1) and higher-order topological phases, and more importantly, their exotic intertwining. While the two phases are independently stabilized by the same symmetry, their intertwining results in the specific distribution of them in the bulk. The construction mechanism allows us to combine different kinds of unconventional Weyl semimetallic and higher-order topological phases to generate distinct phases. Remarkably, on 2D surfaces, the intertwining causes the Fermi-arc topology to change in a periodic pattern against surface orientation. This feature provides a characteristic and feasible signature to probe the intertwining Weyl phases. Moreover, we provide guidelines for searching candidate materials, and elaborate on emulating the intertwined double-Weyl phase in cold-atom experiments.
Similar content being viewed by others
Introduction
Owing to fertile ground of crystalline symmetries, topological gapless phases, characterized by nontrivial band degeneracies, have been undergoing rapid development in condensed matter physics1,2,3. A representative example is the diversification of prototypical Weyl phases. Beyond conventional Weyl fermions with monopole chargeâ±â14,5,6,7,8,9, unconventional Weyl fermions, which possess a higher monopole charge due to crystalline symmetry, were discovered10,11,12,13,14,15,16,17,18,19. For instance, threefold Weyl fermions with chargeâ±â2 can be stabilized by a nonsymmorphic symmetry10,11,12,13, and double(triple)-Weyl fermions with a quadratic(cubic) twofold degeneracy can be protected by a rotation symmetry14,15,16,17,18,19.
On the other hand, crystalline symmetries have largely enriched topological gapped phases20,21,22,23,24,25,26,27,28,29,30,31. Among these phases, higher-order topological (HOT) insulators, which feature anomalous boundary states, have drawn particular attention recently32,33,34,35,36,37,38,39,40,41,42,43,44,45. In contrast to conventional bulk-boundary correspondence, i.e., a d-dimensional bulk topology corresponds to (dâââ1)-dimensional boundary states, the boundary states of HOT insulators are further restricted by crystalline symmetries and exhibit boundary states in a lower dimension, e.g., corner or hinge states. Recently, it has been found that HOT semimetals (e.g., Dirac, Weyl, and nodal-line) can be realized by nontrivially stacking higher-order insulators46,47,48,49,50,51,52,53,54,55,56,57,58. Specifically, higher-order Weyl semimetals were discovered by stacking two-dimensional (2D) HOT insulators, with broken time-reversal and/or inversion symmetry, in a third dimension50,51,52,53,54,55. Note that in this process, for electronic systems, the Kramers degeneracy of the original HOT insulators, if any, has to be broken before stacking50. In this regard, the higher-order Weyl semimetal phase depends on the 2D sub higher-order phases. The phase transition between these phases results in the higher-order Weyl points, which are typically of conventional type with chargeâ±â1.
Though topological gapless phases characterized by nontrivial band degeneracies, and topological gapped phases featuring anomalous boundary states, can be stabilized by the same crystalline symmetry, so far, their interplay remains to be explored.
In this work, we develop a theory to generate topological phases of matter, by intertwining existing ones protected by a crystalline symmetry. We discover intertwined Weyl phases, which feature the exotic interplay of unconventional Weyl fermions and higher-order topology. In the bulk, the two phases are independently stabilized by the crystalline symmetry: the Weyl fermions have a monopole charge larger than 1 by the symmetry, and a HOT phase was further superposed by the symmetry. The intertwining results in a specific distribution of the two phases in the bulk, i.e., the higher-order topology exists in the region outside pairs of unconventional Weyl points with opposite charges. Due to their independence, the two phases are separately tunable. The combination of different unconventional Weyl semimetallic and HOT phases results in distinct intertwined Weyl phases. Note that this mechanism is different from the higher-order Weyl semimetal which is determined by the higher-order topology as discussed above.
The intertwined Weyl phases are characterized by the drastic change of Fermi-arc topology on 2D surfaces upon rotating surface termination. This phenomenon comes from the intertwining between the two constituent phases. Even though HOT phases feature boundary states on 1D hinges, they can intertwine with the unconventional Weyl phase on 2D surfaces of the system: Fermi arcs form due to the underlying unconventional Weyl phases, while the topology of these Fermi arcs is drastically changed by the HOT phase. Thus, we identify a prominent topological feature of the intertwined phase: the change of Fermi-arc topology against surface orientation in a periodic pattern. The period is determined by rotation symmetry, e.g., a 2n-fold rotation symmetry leads to a period of Ï/n. Specifically, we discuss intertwined double-Weyl phases, where double-Weyl fermions and higher-order topology are intertwined due to a fourfold rotation symmetry. The topological phase is characterized by the periodic change of Fermi-arc topology with period Ï/2. The periodic behavior can be perfectly explained by an effective boundary Hamiltonian, in accordance with our theory. Finally, we show that intertwined double-Weyl phases are realizable in cold-atom experiments, which can serve as a promising platform to implement our theory in experiment.
Results
Intertwining topological phases by crystalline symmetry
We begin by showing how to intertwine topological phases by crystalline symmetry. The essential physics can be captured by the following simple but generic Hamiltonian in 3D momentum space,
where k±â=âkxâ±âiky, ϱâ=â(Ï1â±âiÏ2)/2, and Ïi and Ïi (iâ=â1,â2,â3) are Pauli matrices for (pseudo)spin and orbital degrees of freedom, respectively. HWeyl, responsible for generating unconventional Weyl points, is superposed with mÎh, a term for introducing higher-order topology. m is a real model parameter, and Îh a 4âÃâ4 matrix. After being projected to the boundary, mÎh acts as a mass term that is constrained by the symmetry, as we will discuss in Eqs. (3) and (4). Each of the two parts is invariant under the 2n-fold rotation symmetry C2n about the z axis, with 2nââ¤â6 by lattice restriction (nâââ{2,â3}).
As shown in Fig. 1, the intertwining can be induced by a crystalline symmetry that superposes the two phases at the same time:
-
1.
The bulk Weyl points with higher monopole charge are stabilized by rotation symmetry. Unconventional Weyl points with chargeâ±n (nâ>â1) can be generated by breaking time-reversal and/or inversion symmetry in HWeyl. After the substitution of \(({k}_{x},{k}_{y})=(k\cos \theta ,k\sin \theta )\) (Fig. 2a), HWeyl can be rewritten as HWeyl(θ)â=âkne+inθÏ3Ïââ+âkneâinθÏ3Ï+â+âkzÏ3Ï3. As can be seen, these unconventional Weyl points are protected by rotation symmetry14
$${\hat{C}}_{2n}{H}_{{{\mbox{Weyl}}}}(\theta ){\hat{C}}_{2n}^{-1}={H}_{{{\mbox{Weyl}}}}({R}_{2n}\theta ).$$(2)Here, \({\hat{C}}_{2n}={\tau }_{3}{\sigma }_{3}\), and R2nθâ=âθâ+âÏ/n acts in momentum space.
-
2.
The higher-order topology is further enforced by rotation symmetry. The topological protection by symmetry can be explicitly demonstrated at the boundary. After projecting the Hamiltonian to the boundary subspace, we can get an effective Hamiltonian on the surface in the form of (See Supplementary Note 2 for derivation details.)
$${h}_{{{\mathrm{surface}}}}=\mathop{\sum }\limits_{i = 0}^{n}{a}_{i}{k}_{\parallel }^{i}{\sigma }_{3}+m(\theta ){\gamma }_{{{\mathrm{h}}}},$$(3)where ai is a real coefficient, which may depend on kz, θ is the surface orientation, and k⥠the momentum parallel to the surface (Fig. 2a). The higher-order term m(θ)γh, originated from mÎh, is still rotation-symmetric, \({\hat{c}}_{2n}m({{{{\theta }}}}){\gamma }_{{{\mbox{h}}}}{\hat{c}}_{2n}^{-1}=m({R}_{2n}{{{{\theta }}}}){\gamma }_{{{\mbox{h}}}}\), with the projected rotation operator \({\hat{c}}_{2n}={\sigma }_{3}\). Note that m(θ), which stems from the constant m term in Eq. (1), becomes θ-dependent after the projection by the orientation-dependent boundary wavefunctions. γh anticommutes with Ï3, e.g., Ï1 or Ï2, so that it acts as a mass term in Eq. (3). Thus, we obtain
$$m({{{{\theta }}}})=-m({R}_{2n}{{{{\theta }}}})=-m(\theta +\pi /n).$$(4)This relation means m(θ) must have a zero value between θ and θâ+âÏ/n, which leads to a gapless point in the surface spectrum. The 1D Hamiltonian (3) for a fixed kz is characterized by a \({\mathbb{Z}}\) topological invariant. Note that the same mechanism leads to HOT insulators protected by crystalline symmetry, where gapless points correspond to corner or hinge modes43,44. Here, as a manifestation of higher-order topology, the gapless points result in 1D hinge states, as shown by the green solid lines in Fig. 1.
The two topological phases in superposition no longer act individually, but intertwine with each other, as we discuss next. Thus, we refer to the resultant phase as âintertwined Weyl phaseâ, which is different from each individual phase.
We emphasize that in the intertwined Weyl phase, the unconventional Weyl points and the higher-order topology are independently stabilized by the symmetry. This is different from the higher-order Weyl semimetals, where the higher-order Weyl points depend on the HOT phases and are typically conventional ones of topological charge-151,52,53,54,55. Due to the independent nature of the two phases, they can be tuned separately to generate various intertwined Weyl phases. Specifically, different unconventional Weyl fermion phases can be incorporated to build with HOT phase under the same or different crystalline symmetries, resulting in distinct kinds of intertwined Weyl phases, as we will discuss later.
Intertwining and the change of Fermi-arc topology in a periodic pattern
In the bulk, the intertwining results in a specific distribution of the two phases. This is because the higher-order topology cannot survive in regions with non-trivial Chern number. To have vanishing Chern number, it is required that the unconventional Weyl points with opposite charges appear in pairs. In this way, there can exist regions with trivial Chern numbers to host higher-order topology. As shown in Fig. 2 by green arrows, the HOT phase exists outside pairs of unconventional Weyl points with opposite charges.
Inside pairs of unconventional Weyl points, e.g., the lower and upper two sets of Weyl points in Fig. 2b, there is no m(θ), the surface Hamiltonian (3) describes Fermi arcs like those in unconventional Weyl semimetals. On 2D surfaces, the existence of these Fermi arcs is isotropic, i.e., they do not depend on surface orientation θ. However, in HOT regions, the HOT term m(θ) exists. The existence of Fermi arcs in this region depends on m(θ) and becomes anisotropic. We find that due to the intertwining of the two parts of Fermi arcs, the Fermi-arc topology on 2D surfaces is drastically changed. We use m(θ) in Eq. (4) to explain the phenomenon in detail below.
First, the gapless point enforced by the HOT phase in Eq. (4) crucially affects the Fermi-arc topology. We use the spectrum of the surface Hamiltonian (3), \({E}_{{{\mbox{surface}}}}=\pm \sqrt{{(\mathop{\sum }\nolimits_{i = 0}^{n}{a}_{i}{k}_{\parallel }^{i})}^{2}+{m}^{2}({{{{\theta }}}})}\), for illustration. If there is no gapless point [m(θ)ââ â0] in the HOT phase region for a specific angle θ, no Fermi arc can go through this HOT phase region because Esurface is gapped, e.g., as shown by Fig. 2b and d. In contrast, if there is a gapless point [m(θ)â=â0] in HOT region, Fermi arcs can go through this HOT region because Esurface is gapless for some kâ¥, e.g., as shown by Fig. 2c. In this regard, the gapless point [m(θ)â=â0] determines whether the Fermi arcs can go through the higher-order phase region at an angle θ. Thus, strikingly, the Fermi-arc topology changes with θ.
Second, m(θ) in Eq. (4) is antiperiodic. Remarkably, it becomes periodic in the surface spectrum Esurface after taking the square, since the square of m(θ) obeys
which has a period of Ï/n.
The periodic function m2(θ) indicates that the topology of Fermi arcs not only changes with the surface orientation θ, but also in a periodic pattern. We emphasize that the change of Fermi-arc topology comes from the intertwining of unconventional Weyl fermions and higher-order topology. Thus, this phenomenon is absent in any individual phase alone. This is distinct from ref. 59, where the edge states in different samples of different geometry are determined by higher-order topology. Thus, it constitutes the characteristic feature of intertwined Weyl phases.
After figuring out how m(θ) influences the behavior of boundary states on the 2D surfaces, it is also necessary to discuss its impact on the boundary states on the 1D hinges. The hinge states, seen in the HOT phases, are determined by the condition m(θ)â=â0. For m(θ)ââ â0, no localized zero modes can be observed. This is in accordance with the behavior of the HOT phase. Thus, the hinge states of higher-order topology can still be observed in 1D, as shown by the green lines in Fig. 1.
Exemplification of intertwined double-Weyl and triple-Weyl phases
The intertwined double-Weyl phase is generated by the fourfold rotation symmetry C4, i.e., nâ=â2, as we show previously. The Weyl points carry a monopole charge ofâ±â2, which emits two Fermi arcs on the surface. According to Eq. (5), on 2D surfaces, the change of Fermi-arc topology is of period Ï/2. In Fig. 2bâd, we show the topology of Fermi arcs in three representative orientations of θâ=â0, θâ=âÏ/4, and θâ=âÏ/2 in a full period of [0,âÏ/2], respectively. Clearly, by rotating surfaces, the Fermi-arc topology changes from Fig. 2b (θâ=â0) to Fig. 2c (θâ=âÏ/4), and then returns to the initial topology in Fig. 2d (θâ=âÏ/2). On 1D hinges, the boundary states occur at m2(θ)â=â0, e.g., θâ=ânÏ/4 with (nâ=â1,â3,â5,â7). This is in accordance with the boundary states of HOT phase with C4 symmetry, as shown in Fig. 1a.
A similar story applies to the intertwined triple-Weyl phase which is enforced by sixfold rotation symmetry C6, i.e., nâ=â3. On 2D surfaces, each triple-Weyl point emits three Fermi arcs on the surface, and the topology of Fermi arcs changes periodically with a period of Ï/3, as shown by Fig. 2eâg in a full period [0,âÏ/3]. On 1D hinges, the boundary states occur at, e.g., θâ=ânÏ/6 with (nâ=â0,â2,â4,â6,â8,â10). They are in accordance with the boundary states of the HOT phase with C6 symmetry, as shown in Fig. 1b.
Intertwined double-Weyl phase
We now apply our theory to the intertwined double-Weyl semimetals, which can be realized in cold-atom experiments or other artificial systems. The model Hamiltonian reads
where \(M({{{\bf{k}}}})={M}_{0}-2t(\cos {k}_{x}+\cos {k}_{y}+\cos {k}_{z})\). The four double-Weyl points are located on the kz axis at kzâ=â±kw1 and kzâ=â±kw2, where \({k}_{{{\mbox{w}}}1(2)}=\arccos [(M-(+)\sqrt{{\epsilon }^{2}+{m}^{2}})/2t]\) with Mâ=âM0âââ4t. By series expansion around the Weyl points at (0,â0,â±kw1(2)), i.e., \(({k}_{x},{k}_{y},{k}_{z})\to (\delta k\cos \theta ,\delta k\sin \theta ,\delta {k}_{z})\), the form of the low energy model H(k)â=âAδke+2iθÏ3Ïââ+âAδkeâ2iθÏ3Ï+â+âtδkzÏ3Ï3â+âϵÏ0Ï3â+âmÏ1Ï1 is the same as Eq. (1). Double-Weyl points are generated by the ϵ term that breaks time-reversal symmetry. They possess monopole charge of ±260,61. mÏ1Ï1 corresponds to mÎh in Eq. (1), which is responsible for the higher-order topology. The whole system is protected by fourfold rotation symmetry \({\hat{C}}_{4}={\sigma }_{3}{\tau }_{3}\). The Fermi arcs shown in Fig. 2bâd are numerically calculated by using the model (6) with open boundary conditions.
In Fig. 3a, the Chern number in kxky-plane against kz is plotted. Here, without loss of generality, we assume kw2â<âkw1, for the chosen parameters given in the caption of Fig. 3. We can see that the Chern number changes by 2 when passing a double-Weyl point. This is because the Chern number on the surface that encloses the Weyl point equals to the monopole charge, as it does for conventional Weyl phases. This explains the non-trivial Chern number between the upper pair kzâââ(kw2,âkw1) and lower pair kzâââ(âkw1, âkw2) of Weyl points. However, the distribution of the Chern number in the intertwined phase is different from conventional ones. In order to host HOT phases, it is required to have a region where the Chern number vanishes. Here, the region is between the lower and upper pairs of Weyl points, i.e., \({k}_{z}\in \left(-{k}_{{{\mbox{w}}}2},{k}_{{{\mbox{w}}}2}\right)\), as shown by the HOT region (green arrow) in Fig. 3. In this region, the existence of boundary states depends on the orientation θ, as shown by Fig. 3c. In contrast, in the region with non-zero Chern number, i.e., inside the two sets (upper and lower) of Weyl points, the Fermi arcs exist regardless of the angle θ.
Effective boundary theory
To understand the periodic behavior of Fermi-arc topology, a boundary theory applicable to any surface orientation is constructive. We can achieve this goal by firstly deriving two boundary states for each θ in the absence of the higher-order term mÏ1Ï1 in Eq. (6). The spinor part of two boundary states takes the form of \({\psi }_{1}\propto ({e}^{-2i\theta },\sqrt{2}+1,0,0)\) approximately in the region of kzâââ(âkw2,âkw2), and \({\psi }_{2}\propto (0,0,{e}^{-2i\theta },\sqrt{2}+1)\) in the region of kzâââ(âkw1,âkw1). Note that the contribution from spatial part of the boundary states does not affect the main results, and is neglected for simplicity. By projecting the whole Hamiltonian into the subspace spanned by the boundary states, we can obtain the effective boundary Hamiltonian as (See Supplementary Note 1 and 2 for derivation details.)
up to a constant term \(-1/\sqrt{2}\epsilon {\sigma }_{0}\), and \({k}_{\,{{\mbox{c}}}\,}^{2}=2t\cos {k}_{z}-M\). Clearly, the boundary Hamiltonian is in the form of Eq. (3). The effective Hamiltonian is valid in the region where Ï1 and Ï2 overlap, i.e., \({k}_{z}\in (-{k}_{{{\mbox{w}}}2},{k}_{{{{\mbox{w}}}}_{2}})\) between the middle two Weyl points, where the Chern number is trivial. On a 1D hinge, the gapless points at \(m\cos (2\theta )=0\) result in hinge Fermi arcs, as shown by the green solid lines in Fig. 1.
On 2D surfaces, clearly, the periodic change of Fermi-arc topology, shown in Fig. 2bâd, is caused by the higher-order term \(m\cos (2\theta )\). The period is Ï/2, as determined by \({m}^{2}{\cos }^{2}(2\theta )\) in the spectrum, in accordance with the general theory of Eq. (5). Within a single period, the higher-order topology enforces the appearance of gapless point [Eq. (4)], which is located at \(m\cos (2\theta )=0\). The gapless point drastically changes the topology of Fermi arcs, because it determines whether the arcs can go through the region of kzâââ(âkw2,âkw2) between the middle two Weyl points or not. Figure 2bâd shows three representative orientations in a full period of θâââ[0,âÏ/2]. The Fermi arcs cannot go through the region between the two middle Weyl points at θâ=â0. After rotating to θâ=âÏ/4 at the gapless point, the Fermi arcs are allowed to go through. Finally, at θâ=âÏ/2, the Fermi arcs return to their initial topology, completing one period.
Guidelines for material search
Two key ingredients of the discussed intertwined Weyl phases are fourfold (sixfold) rotation symmetry and the unconventional Weyl points on the rotation axis. Both requirements can be fulfilled in tetragonal, cubic, or hexagonal space groups, where the desired symmetries emerge. There is a correspondence between the absolute value of the chirality of a crossing and the rotation eigenvalues of the involved bands62, such that a crossing between certain rotation eigenvalues leads to the required unconventional Weyl points. The total phase of eigenvalues accumulated by such crossings for any pair of bands is restricted by the periodicity of Brillouin zone. To be more specific, to obtain crossings of the required charge of ±2 (±3), the phase of the eigenvalues of fourfold (sixfold) rotation must change by Ï. By looking at Fig. 1 it is evident that intertwined Weyl points require an even number of such changes of the eigenvalue phase. Thus, the necessary crossings can only occur if the total accumulated phase, which is consistent with the periodicity of the Brillouin zone, is equal to \(0\,{{\mathrm{mod}}}\,\,2\pi\), like it is fulfilled in Eq. (6). This excludes for nonsymmorphic systems certain filling factors, where the k-dependence of rotation eigenvalues would, for example, require an odd number of crossings. An appropriate filling b for a 2n-fold rotation with a fractional translation of \(\frac{m}{2n}\) has to fulfill \(\frac{bm}{2n}\in {\mathbb{Z}}\).
Commonly, a set of four crossings like in Fig. 1 would be accidental. Nevertheless, enforced crossings can be found within fourfold degenerate points, e.g., in space groups 106 and 133 on the path M-A. These comprise unconventional Weyl points related by mirror symmetry63. Hereby, mirror pairs of enforced unconventional Weyl points coincide. If a weak time-reversal and mirror symmetry breaking is introduced to such a system, one may obtain Weyl points in the configuration discussed in Fig. 1a with no other enforced bands at the Fermi energy. A comprehensive material search is left for future works.
Cold-atom experimental realization
Owing to technical advances, cold atoms have been widely applied in quantum simulations of topological matter64,65,66,67, and now are also readily available for realizing our intertwined double-Weyl semimetal described by Eq. (6). Here we present the realization proposal using fermionic atoms. We choose two hyperfine states as the pseudo-spins ââ for the Ï degrees of freedom. For our model Hamiltonian, which requires two extra degrees of freedom, we consider the atomic gases loaded in a 2D bilayer optical lattice, and thus, the layer index λ is used to represent the Ï subspace. The setup is illustrated in Fig. 4. The energy offset of the hyperfine states is prepared as \({{{\Gamma }}}_{\lambda }(\phi )={M}_{0}+{(-1)}^{\lambda }\epsilon -2t\cos \phi\) that is manually controllable by the cyclical parameter Ï. It not only generates the on-site energy of the model Hamiltonian, but also introduces the additional parameterized space represented by Ï. In this way, the band physics of the Weyl Hamiltonian can be captured in the 3D parameterized space (kx,âky,âÏ)68. Here, Ï corresponds to kz in Eq. (6).
In order to engineer the intra-layer spin-flipped hopping, we use laser fields of three modes to couple the pseudo-spins. The spatial modulations of the field modes are prepared as \({M}_{1}({{{\bf{r}}}})=i{M}_{1}\sin ({k}_{{{\mbox{L}}}}x)\sin ({k}_{{{\mbox{L}}}}y)\), \({M}_{2}({{{\bf{r}}}})={M}_{2}\cos ({k}_{{{\mbox{L}}}}x)\cos (3{k}_{{{\mbox{L}}}}y)\cos ({k}_{{{\mbox{L}}}}z)\), and \(M^{\prime}_{2} ({{{\bf{r}}}})=-{M}_{2}\cos (3{k}_{{{\mbox{L}}}}x)\cos ({k}_{{{\mbox{L}}}}y)\cos ({k}_{{{\mbox{L}}}}z)\), where kLâ=âÏ/d and d denotes the lattice constant. Due to the odd parity of M1(r) in the xy-plane69, the on-site and nearest-neighbor (NN) couplings vanish, while the next-NN coupling dominates, resulting in \(\sin {k}_{x}\sin {k}_{y}{\sigma }_{1}\). Due to the crystalline symmetry, the combination of M2(r) and \(M^{\prime}_{2} ({{{\bf{r}}}})\) leads to the NN coupling \((\cos {k}_{y}-\cos {k}_{x}){\sigma }_{1}\). After making operator transformations, all the intra-layer terms host opposite signs for different layers. Furthermore, the HOT term Ï1Ï1 is naturally introduced by the inter-layer hopping. The details of the realization proposal are shown in Supplementary Note 3.
Discussion
Based on the theory to create topological phases by intertwining existing ones, we have discovered the intertwined Weyl phase. The intertwined Weyl phase is different from the individual unconventional Weyl semimetallic phase or HOT phase, and exhibits its distinct characteristic topological features. The intertwining results in the drastic change of Fermi-arc topology in a periodic pattern, which constitutes the characteristic feature of the intertwined phase. We have proposed a feasible cold-atom experiment to verify our theory and to realize the intertwined Weyl phases.
Our theory could serve as a guiding principle to generate topological phases based on existing ones. A direct application would be to investigate the intertwining between topological gapless phases with emergent particles other than Weyl fermions and topological crystalline phases, that are protected by the same symmetry28,29,30,70,71.
Finally, we note that Fermi-surface topology is crucial for electronic properties of material. The change of Fermi-arc topology in the intertwined Weyl phases offers a direction of tuning Fermi-surface topology. It would inspire further research on electrical, magnetic, thermodynamic, and transport properties, that are determined by Fermi-surface topology. For instance, a potential application is to investigate the quantum oscillations of Fermi-arc surface states72. These oscillations are periodic against the inverse of magnetic field 1/B. Their frequency F is determined by the Onsager relation Fâ=âΦ0/(2Ï2)As. Here Φ0â=âh/2e is the magnetic flux quantum, and As is the Fermi surface cross section. In intertwined Weyl phases, the cross section As depends on the surface orientation θ. As shown in Fig. 3c, by changing θ, one may interpolate between small Fermi surfaces (localized around kw1 and kw2) and extended Fermi surfaces (connecting kw1(2) to âkw1(2)), which must yield a substantial change in the observed quantum oscillation spectra. Thus, the intertwined Weyl phases are expected to have a significant change of quantum oscillations upon rotating surface termination.
Methods
Analytical derivation and symmetry analysis
The derivation of the analytical results makes use of the kââ âp methods and symmetry analysis. The numerical calculation is based on tight-binding model and Greenâs function method. Details of these derivations are given in the Supplementary Notes 1, 2, and 4. This includes the calculation of Fermi arcs in intertwined double-Weyl semimetals, and the analysis of symmetry for intertwined double-Weyl and triple-Weyl semimetals. The Supplementary Note 3 also contains a detailed discussion of the cold-atom realization of intertwined double-Weyl semimetals.
Data availability
The numerical data of the plots within this paper are available from the corresponding author upon reasonable request.
Code availability
The numerical codes that support the findings of this paper are available from the corresponding author upon reasonable request.
References
Chiu, C.-K., Teo, J. C. Y., Schnyder, A. P. & Ryu, S. Classification of topological quantum matter with symmetries. Rev. Mod. Phys. 88, 035005 (2016).
Armitage, N. P., Mele, E. J. & Vishwanath, A. Weyl and Dirac semimetals in three-dimensional solids. Rev. Mod. Phys. 90, 015001 (2018).
Yang, S. A., Pan, H. & Zhang, F. Dirac and Weyl superconductors in three dimensions. Phys. Rev. Lett. 113, 046401 (2014).
Wan, X., Turner, A. M., Vishwanath, A. & Savrasov, S. Y. Topological semimetal and Fermi-arc surface states in the electronic structure of pyrochlore iridates. Phys. Rev. B 83, 205101 (2011).
Xu, S.-Y. et al. Discovery of a Weyl fermion semimetal and topological Fermi arcs. Science 349, 613â617 (2015).
Lv, B. Q. et al. Experimental discovery of Weyl semimetal TaAs. Phys. Rev. X 5, 031013 (2015).
Meng, T. & Balents, L. Weyl superconductors. Phys. Rev. B 86, 054504 (2012).
Morali, N. et al. Fermi-arc diversity on surface terminations of the magnetic Weyl semimetal Co3Sn2S2. Science 365, 1286â1291 (2019).
Liu, D. F. et al. Magnetic Weyl semimetal phase in a kagomé crystal. Science 365, 1282â1285 (2019).
Bradlyn, B. et al. Beyond Dirac and Weyl fermions: Unconventional quasiparticles in conventional crystals. Science 353, aaf5037 (2016).
Sanchez, D. S. et al. Topological chiral crystals with helicoid-arc quantum states. Nature 567, 500â505 (2019).
Yang, Y. et al. Topological triply degenerate point with double Fermi arcs. Nat. Phys. 15, 645â649 (2019).
He, H. et al. Observation of quadratic Weyl points and double-helicoid arcs. Nat. Commun. 11, 1820 (2020).
Fang, C., Gilbert, M. J., Dai, X. & Bernevig, B. A. Multi-Weyl topological semimetals stabilized by point group symmetry. Phys. Rev. Lett. 108, 266802 (2012).
Chen, W.-J., Xiao, M. & Chan, C. T. Photonic crystals possessing multiple Weyl points and the experimental observation of robust surface states. Nat. Commun. 7, 13038 (2016).
Huang, S.-M. et al. New type of Weyl semimetal with quadratic double Weyl fermions. Proc. Natl Acad. Sci. USA 113, 1180â1185 (2016).
Vaidya, S. et al. Observation of a charge-2 photonic Weyl point in the infrared. Phys. Rev. Lett. 125, 253902 (2020).
Yang, Y. et al. Ideal unconventional Weyl point in a chiral photonic metamaterial. Phys. Rev. Lett. 125, 143001 (2020).
Dantas, R. M. A., Peña Benitez, F., Roy, B. & Surówka, P. Non-abelian anomalies in multi-Weyl semimetals. Phys. Rev. Res. 2, 013007 (2020).
Fu, L. Topological crystalline insulators. Phys. Rev. Lett. 106, 106802 (2011).
Hsieh, T. H. et al. Topological crystalline insulators in the SnTe material class. Nat. Commun. 3, 982 (2012).
Tanaka, Y. et al. Experimental realization of a topological crystalline insulator in SnTe. Nat. Phys. 8, 800â803 (2012).
Dziawa, P. et al. Topological crystalline insulator states in PbSnSe. Nat. Mater. 11, 1023â1027 (2012).
Ando, Y. & Fu, L. Topological crystalline insulators and topological superconductors: from concepts to materials. Annu. Rev. Condens. Matter Phys. 6, 361â381 (2015).
Fang, C. & Fu, L. New classes of three-dimensional topological crystalline insulators: nonsymmorphic and magnetic. Phys. Rev. B 91, 161105 (2015).
Kruthoff, J. et al. Topological classification of crystalline insulators through band structure combinatorics. Phys. Rev. X 7, 041069 (2017).
Khalaf, E., Po, H. C., Vishwanath, A. & Watanabe, H. Symmetry indicators and anomalous surface states of topological crystalline insulators. Phys. Rev. X 8, 031070 (2018).
Po, H. C., Vishwanath, A. & Watanabe, H. Symmetry-based indicators of band topology in the 230 space groups. Nat. Commun. 8, 50 (2017).
Tang, F., Po, H. C., Vishwanath, A. & Wan, X. Comprehensive search for topological materials using symmetry indicators. Nature 566, 486â489 (2019).
Zhang, T. et al. Catalogue of topological electronic materials. Nature 566, 475â479 (2019).
Vergniory, M. G. et al. A complete catalogue of high-quality topological materials. Nature 566, 480â485 (2019).
Benalcazar, W. A., Bernevig, B. A. & Hughes, T. L. Quantized electric multipole insulators. Science 357, 61â66 (2017).
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. et al. Reflection-symmetric second-order topological insulators and superconductors. Phys. Rev. Lett. 119, 246401 (2017).
Peterson, C. W. et al. A fractional corner anomaly reveals higher-order topology. Science 368, 1114â1118 (2020).
Schindler, F. et al. Higher-order topological insulators. Sci. Adv. 4, eaat0346 (2018).
Schindler, F. et al. Higher-order topology in bismuth. Nat. Phys. 14, 918â924 (2018).
Serra-Garcia, M. et al. Observation of a phononic quadrupole topological insulator. Nature 555, 342â345 (2018).
Ezawa, M. Higher-order topological insulators and semimetals on the breathing kagome and pyrochlore lattices. Phys. Rev. Lett. 120, 026801 (2018).
Wang, Z. et al. 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).
Chen, R. et al. Higher-order topological insulators in quasicrystals. Phys. Rev. Lett. 124, 036803 (2020).
Zhang, S.-B., Calzona, A. & Trauzettel, B. All-electrically tunable networks of Majorana bound states. Phys. Rev. B 102, 100503 (2020).
Trifunovic, L. & Brouwer, P. W. Higher-order bulk-boundary correspondence for topological crystalline phases. Phys. Rev. X 9, 011012 (2019).
Khalaf, E. Higher-order topological insulators and superconductors protected by inversion symmetry. Phys. Rev. B 97, 205136 (2018).
Zhang, S.-B. et al. Topological and holonomic quantum computation based on second-order topological superconductors. Phys. Rev. Res. 2, 043025 (2020).
Lin, M. & Hughes, T. L. Topological quadrupolar semimetals. Phys. Rev. B 98, 241103 (2018).
Wang, K. et al. Boundary criticality of \({{{\mathcal{PT}}}}\)-invariant topology and second-order nodal-line semimetals. Phys. Rev. Lett. 125, 126403 (2020).
Ghorashi, S. A. A., Hu, X., Hughes, T. L. & Rossi, E. Second-order Dirac superconductors and magnetic field induced Majorana hinge modes. Phys. Rev. B 100, 020509 (2019).
Tiwari, A., Jahin, A. & Wang, Y. Chiral Dirac superconductors: Second-order and boundary-obstructed topology. Phys. Rev. Res. 2, 043300 (2020).
Roy, B. Antiunitary symmetry protected higher-order topological phases. Phys. Rev. Res. 1, 032048 (2019).
Wang, H.-X. et al. Higher-order Weyl semimetals. Phys. Rev. Lett. 125, 146401 (2020).
Ghorashi, S. A. A., Li, T. & Hughes, T. L. Higher-order Weyl semimetals. Phys. Rev. Lett. 125, 266804 (2020).
Rui, W. B. et al. Higher-order Weyl superconductors with anisotropic Weyl-point connectivity. Phys. Rev. B 103, 184510 (2021).
Wei, Q. et al. Higher-order topological semimetal in acoustic crystals. Nat. Mater. 20, 812â817 (2021).
Luo, L. et al. Observation of a phononic higher-order Weyl semimetal. Nat. Mater. 20, 794â799 (2021).
Zhang, R.-X., Hsu, Y.-T. & Das Sarma, S. Higher-order topological Dirac superconductors. Phys. Rev. B 102, 094503 (2020).
CÄlugÄru, D., JuriÄiÄ, V. & Roy, B. Higher-order topological phases: a general principle of construction. Phys. Rev. B 99, 041301 (2019).
Szabó, A. L. & Roy, B. Dirty higher-order Dirac semimetal: quantum criticality and bulk-boundary correspondence. Phys. Rev. Res. 2, 043197 (2020).
Zhang, X. et al. Second-order topology and multidimensional topological transitions in sonic crystals. Nat. Phys. 15, 582â588 (2019).
Bera, S., Sau, J. D. & Roy, B. Dirty Weyl semimetals: stability, phase transition, and quantum criticality. Phys. Rev. B 93, 201302 (2016).
Dantas, R. M. A., Peña Benitez, F., Roy, B. & Surówka, P. Non-abelian anomalies in multi-Weyl semimetals. Phys. Rev. Res. 2, 013007 (2020).
Tsirkin, S. S., Souza, I. & Vanderbilt, D. Composite Weyl nodes stabilized by screw symmetry with and without time-reversal invariance. Phys. Rev. B 96, 045102 (2017).
Hirschmann, M. M. et al. Symmetry-enforced band crossings in tetragonal materials: Dirac and Weyl degeneracies on points, lines, and planes. Phys. Rev. Mater. 5, 054202 (2021).
Zhang, D.-W. et al. Topological quantum matter with cold atoms. Adv. Phys. 67, 253â402 (2018).
Wu, Z. et al. Realization of two-dimensional spin-orbit coupling for Bose-Einstein condensates. Science 354, 83â88 (2016).
Song, B. et al. Observation of nodal-line semimetal with ultracold fermions in an optical lattice. Nat. Phys. 15, 911â916 (2019).
Wang, Z.-Y. et al. Realization of an ideal Weyl semimetal band in a quantum gas with 3D spin-orbit coupling. Science 372, 271â276 (2021).
Zhang, D.-W., Zhu, S.-L. & Wang, Z. D. Simulating and exploring Weyl semimetal physics with cold atoms in a two-dimensional optical lattice. Phys. Rev. A 92, 013632 (2015).
Zheng, Z. et al. Chiral magnetic effect in three-dimensional optical lattices. Phys. Rev. Res. 1, 033102 (2019).
Hasan, M. Z. et al. Weyl, Dirac and high-fold chiral fermions in topological quantum matter. Nat. Rev. Mater. 6, 784â803 (2021).
Yu, Z.-M. et al. Encyclopedia of emergent particles in three-dimensional crystals. Sci. Bull. https://doi.org/10.1016/j.scib.2021.10.023 (2021).
Potter, A. C., Kimchi, I. & Vishwanath, A. Quantum oscillations from surface Fermi arcs in Weyl and Dirac semimetals. Nat. Commun. 5, 5161 (2014).
Acknowledgements
This work was supported by the Key-Area Research and Development Program of GuangDong Province (Grant No. 2019B030330001), the NSFC/RGC JRS grant (No.: N_HKU774/21), the CRFs (Nos. C6005-17G and C6009-20G) and GRF (No. 17300220) of Hong Kong. The authors at HKU also thank support from Guangdong-Hong Kong Joint Laboratory of Quantum Matter.
Author information
Authors and Affiliations
Contributions
W.B.R., C.J.W., and Z.D.W. conceived the project. W.B.R., M.M.H., and S.B.Z. did the theoretical calculations. W.B.R. performed the numerical simulations. Z.Z. proposed the cold-atom realization, and M.M.H. analyzed the candidate materials. All authors authored, commented, and corrected the manuscript.
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 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
Rui, W.B., Zheng, Z., Hirschmann, M.M. et al. Intertwined Weyl phases emergent from higher-order topology and unconventional Weyl fermions via crystalline symmetry. npj Quantum Mater. 7, 15 (2022). https://doi.org/10.1038/s41535-022-00422-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41535-022-00422-0
This article is cited by
-
Orbital degree of freedom induced multiple sets of second-order topological states in two-dimensional breathing Kagome crystals
npj Quantum Materials (2023)
-
Square-root higher-order Weyl semimetals
Nature Communications (2022)