Topology on a new facet of bismuth.

Bismuth-based materials have been instrumental in the development of topological physics, even though bulk bismuth itself has been long thought to be topologically trivial. A recent study has, however, shown that bismuth is in fact a higher-order topological insulator featuring one-dimensional (1D) topological hinge states protected by threefold rotational and inversion symmetries. In this paper, we uncover another hidden facet of the band topology of bismuth by showing that bismuth is also a first-order topological crystalline insulator protected by a twofold rotational symmetry. As a result, its [Formula: see text] surface exhibits a pair of gapless Dirac surface states. Remarkably, these surface Dirac cones are "unpinned" in the sense that they are not restricted to locate at specific k points in the [Formula: see text] surface Brillouin zone. These unpinned 2D Dirac surface states could be probed directly via various spectroscopic techniques. Our analysis also reveals the presence of a distinct, previously uncharacterized set of 1D topological hinge states protected by the twofold rotational symmetry. Our study thus provides a comprehensive understanding of the topological band structure of bismuth.

Bismuth is well known for its peculiar physical properties. It was long considered to be the stable element with the highest atomic mass, but relatively recent experiments have shown that bismuth is in fact weakly radioactive (1). It is a semimetal with a vanishingly small carrier density ( cm−3) but an exceptionally high electron mobility (s−1⋅V−1) (2–5). As a result, the ultraquantum regime is reached in bismuth at a magnetic field as small as 9 T, beyond which a number of correlated electron states have been observed (4–7). Because of bismuth’s large spin–orbit coupling, bismuth-based materials have also played a fundamental role in topological physics (8, 9). A bismuth–antimony alloy (Bi1−Sbx) was the first experimental realization of a 3D topological insulator (TI) (10). The family supports the prototypical TI state with a single surface Dirac cone (11–14). Doped gives rise to the quantum anomalous Hall effect (15) and unconventional nematic superconductivity (16, 17). Bi is a 3D Dirac semimetal (18, 19). Despite being a crucially important chemical component in many topological materials, pure bismuth has long been thought to be topologically trivial. However, distinct from an atomic band insulator, bismuth does have an even number of band inversions (10, 20). The existence of these band inversions implies that bismuth cannot be smoothly connected to an atomic band insulator without undergoing closing of its band gap. This suggests that the band topology of bismuth must be far from being trivial (8, 9, 20–22).

Theoretical advances on topological crystalline insulators (TCIs) (21) have greatly expanded the topological classification of band insulators beyond the TIs (8, 9). After the success of theoretical prediction and experimental realization of the mirror-symmetry–protected TCIs (20–35), two novel TCI phases, rotational-symmetry–protected TCIs (36) and higher-order TIs (37–43), have been theoretically proposed recently. Rotational-symmetry–protected TCIs are predicted to harbor “unpinned” Dirac surface states in that the surface normal to an -fold rotational axis () can host Dirac cones whose Dirac points appear at generic points in the surface Brillouin zone. In contrast, higher-order TIs are generated through the consideration of higher-order bulk-boundary correspondence. For instance, a 3D second-order TI supports 1D topological hinge states in a rotational symmetry-preserving rod. It is important to note that the first-order and higher-order topologies are not mutually exclusive. In particular, a TCI could support both topological 2D surface states and 1D hinge states. Thus, strictly speaking, a pure higher-order TI should refer to the TCIs that do not show the presence of topological surface states. Beyond these new TCI phases, another important theoretical advance is the development of methods to systematically diagnose topological invariants in terms of the symmetry eigenvalues of the electronic states (44–51). In this connection, Song et al. (47) and Khalaf et al. (48) found that, when certain additional symmetry is present, topological invariants of TCIs protected by symmetry can be inferred from the -related symmetry eigenvalues of the energy bands. Such proposals of symmetry indicators and topological quantum chemistry have facilitated first-principles studies of new topological materials (52–58).

Building upon these theoretical advances, a recent work (41) showed that pure bismuth hosts a second-order band topology that is protected by threefold rotational and inversion symmetries. As a result, it supports 1D topological hinge states in a crystal whose shape preserves the threefold rotational and inversion symmetries. Here we show that the band topology of bismuth is even richer in that it also hides a first-order TCI state, which is protected by its twofold rotational symmetries, resulting in a pair of unpinned topological Dirac surface states on its surface. We also reveal the presence of a distinct, previously uncharacterized set of 1D topological hinge states protected by twofold rotational symmetry.

Results

First-principles calculations were carried out using the VASP package (59, 60). An ultrasoft pseudopotential and the generalized-gradient approximation were applied in the self-consistency process. The experimental crystal structure of Bi was used (61) in the calculations. A Wannier-basis–based tight-binding model, where the and orbitals of Bi were included, was obtained via the Wannier90 package (62).

Bismuth crystalizes in a rhombohedral structure (61). Its space group and point group are R-3m () and , respectively. Each layer in this structure forms a buckled honeycomb lattice. The three principle lattice vectors () are tilted with respect to the Cartesian directions (Fig. 1). In this construction, the out-of-plane direction is [111], and the , , and directions lie in plane. Key point-group symmetries of bismuth are space inversion ; out-of-plane threefold rotation axis ; in-plane twofold rotational axes , , and ; and mirror planes , , and . The bulk Brillouin zone (BZ) has the shape of a truncated octahedron (Fig. 1). The time-reversal invariant momentum (TRIM) points include one , one , three , and three symmetry points. Our band structure (Fig. 1) confirms that bismuth is a compensated semimetal with a continuous band gap. To unfold the nature of band inversions in bismuth, we calculated the parity eigenvalues at the TRIM points, and with the exception of the point, the lowest three valence bands at all of the TRIM points were found to have one state with positive parity and two states with negative parity. In contrast, at , all three valence bands have positive parity eigenvalues, so that bismuth has two band inversions at the point.

Fig. 1.

(A) Crystal structure of bulk bismuth. The principle lattice vectors , , and are expressed in the Cartesian coordinates. (B) The bulk Brillouin zone. The principle reciprocal space vectors (, , and ) and high-symmetry points are noted. (C) Band structure of bismuth including spin–orbit coupling. The red and green circles represent positive and negative parity eigenvalues, respectively.

We turn next to discuss the symmetry-based indicators. As shown in refs. 47 and 48, crystals in the space group are characterized by four symmetry indicators, three invariants and one invariant (). By enumerating the symmetry eigenvalues of the electronic states at high-symmetry points, we obtain . Table 1 shows the corresponding topological states and the associated topological invariants. Importantly, the specific symmetry indicators () point to two possible topological states, one being a purely rotational-symmetry–protected TCI with , while the other is a purely mirror-symmetry–protected TCI with . To uniquely determine the topological state of bismuth, we have further calculated the mirror Chern number and found for bismuth; see for details. Our analysis thus reveals that bismuth is a purely rotational-symmetry–protected TCI with . Note that there are three twofold rotational axes (, , and ) involved here, which are related by the out-of-plane threefold rotation . We elaborate on the nontrivial topology protected by . The other two twofold rotational axes are associated with the same nontrivial topology, so that the associated surfaces can also be expected to host a pair of unpinned Dirac surface states. Moreover, the symmetry indicators point to a nontrivial topological invariant associated with space-inversion symmetry (). This guarantees 1D topological hinge states protected by inversion symmetry, which is consistent with previous findings based on the additional constraint (41). We emphasize that the high-order TI predicted in ref. 41 does not support first-order surface states (36, 47, 48) because there are no topological invariants associated with alone.

Table 1.

Two possible topological states for the symmetry indicator (47)

(ν_0;ν_1ν_2ν_3)	n_2_[11¯0]	n_M_[11¯0]	n_I
(0;000)	1	0	1
(0;000)	0	2	1

are the invariants for 3D TIs. is a invariant for the twofold rotational symmetry . corresponds to a rotational-symmetry–protected TCI, which features two unpinned Dirac surface states on the surface. is the mirror Chern number for mirror plane and is a invariant. features N Dirac surface states. Finally, is the classified inversion () protected TCI index. features the hinge state in a 3D finite geometry preserving .

To obtain a deeper understanding of , we provide an intuitive physical picture that connects the band inversion property to the rotational-symmetry–protected topology along the lines of the physics discussed in ref. 41. Recall that a single band inversion at can lead to a strong TI state. Therefore, the double band inversions in bismuth can be viewed as representing two sets of strong TI band structures, where each strong TI contributes a surface Dirac cone at (Fig. 2). The nontrivial topology protected by the twofold rotation then yields the double Dirac cones, which remain gapless on the surface. To further connect the above intuitive picture to the symmetry indicators, we provide a more detailed analysis of the representations at the TRIM points. For this purpose, it is helpful to rearrange the Fu–Kane formula (63) (Eq. ) aswhere is the strong TI (STI) index, and is the number of occupied bands with negative parity eigenvalue. For space group 166, only TRIMs are invariant under one of the three axes. However, states at the other pairs of and can be linearly combined into states invariant under . Thus, the occupied bands can be separated into two subspaces according to their respective eigenvalues (), which are and (Eq. ). This rearrangement may lead to new topology beyond the TI state, but it can be made only when the system is trivial (). From Eq. , we see that plays a role similar to that of within each of the subspaces. can be directly obtained from the symmetry indicator (37, 46–48), usingEq. shows that corresponds to . In other words, the band structure in each subspace can be viewed as a strong TI. The hybridization between the two strong TI surface states in the presence of the rotational symmetry thus leads to the pair of unpinned Dirac surface states on the surface.

Fig. 2.

(A and B) The projected () surface Brillouin zone on which several points, paths, and the predicted positions of the unpinned surface DPs are detailed. (C) The schematic shows the unpinned surface DPs on the surface preserving . (D) Surface band spectrum along the high-symmetry path, . Inset shows that the surface states are gapped along this path. (E and F) Surface band spectrum along the paths containing the DP (see the trajectories in B). The spectrum in F shows that the surface Dirac fermion is type II: The Dirac cone is strongly tilted so that the velocities of the two bands that cross have the same sign. (G and H) WCCs of the loop integrated around () at (G) and (H) with respect to , where A and B indicate two WCC sheets that connect only at generic ().

With the preceding discussion in mind, we computed the surface band structure throughout the surface BZ (Fig. 2). As expected, we found two sets of surface Dirac cones, whose Dirac points are located at () of () axes indicated in Fig. 1. These are generic points, which are related only by the twofold rotational symmetry . The calculated surface states along the paths passing through the Dirac points (DPs) (Fig. 2 ) directly show the presence of gapless Dirac surface states. Interestingly, these unpinned Dirac fermions are of type II (64) in that the velocities of the two surface bands that cross possess the same sign. By contrast, along the path that does not go through the DP, the surface states are gapped (Fig. 2). To further confirm the nontrivial nature of these states, we studied the Wannier charge centers (WCCs) (65–67) on the surface, where is the th Wannier state in the unit cell of = 0 and is the position operator (68). Two dominant WCC sheets, marked by A and B in Fig. 2 , are seen to be disjointed in the whole surface Brillouin zone, except at (Fig. 2). These results unambiguously demonstrate the presence of unpinned Dirac surface states associated with the twofold rotational-symmetry–protected topology.

We now demonstrate the 1D topological hinge states protected by . For this purpose, we constructed two types of bismuth rods (Fig. 3 ), which are periodic along the axis but finite sized within the plane that is perpendicular to . In the first rod, the side surfaces are the (111) and surfaces and the cross-section is rectangular (Fig. 3). In the second rod, the side surfaces are (001) and (111) and the cross-section is a parallelogram (Fig. 3). Both rods are invariant under the twofold rotation . The calculated band structures for rods 1 and 2 are shown in Fig. 3 , from which we can clearly identify the existence of 1D helical edge states lying inside the bulk band gap. We also investigated the real-space distribution of the wave function of these 1D helical edge states. Fig. 3 shows that these helical states are localized on the edges shared by adjacent side surfaces, further confirming that these are topological hinge states.

Fig. 3.

(A and B) Two types of 1D bismuth rods that preserve the twofold rotation symmetry . These rods are periodic along the axis but finite sized within the plane that is perpendicular to . In the first rod (A), the side surfaces are the (111) and surfaces and the cross-section is a rectangle. (C and D) Energy dispersions of rods 1 and 2. The topological hinge states are marked by red lines. (E and F) Real-space distributions of the hinge states of rods 1 and 2.

Finally, we show that the rotational-symmetry–protected topology sheds light on the large Rashba surface states observed on the (111) and (001) surfaces (2, 3). The double bulk band inversion at also leads to a pair of surface Dirac cones on the (111) and (001) surfaces. However, as (111) and (001) surfaces do not respect , these Dirac cones are gapped, which naturally gives rise to the Rashba surface bands (Fig. 4). Note that since the Rashba surface states are nontopological, they can be created or removed by adjusting the surface potential (69). Our analysis gives insight into the origin of Rashba bands on bismuth’s (111) and (001) surfaces, which have been observed experimentally (2, 3, 69), in terms of the underlying symmetries and the double band inversions in the band topology of bismuth.

Fig. 4.

(A and B) The projected (111) and (001) surface Brillouin zones (SBZs) on which some of the high-symmetry k points are indicated. (C) Schematic illustration of the gapped surface states due to the double band inversions without the protection as discussed in the text. (D and E) Band spectrum for (111) and (001) surfaces, respectively. The green dotted lines show the Rashba bands due to the double band inversions.

In summary, we have investigated topological properties of the bulk and surface band structures of bismuth. We show that bismuth is a TCI with multiple nontrivial topological invariants. In the first order, bismuth features unpinned Dirac surface states on the , , and surfaces that are protected by their twofold rotational symmetries. These unpinned Dirac surface states would be amenable to direct imaging via photoemission spectroscopy. In the second order, bismuth features two sets of independent topological helical hinge states. The first set is protected by the twofold rotational symmetries whereas the second set is protected by space inversion symmetry. Our study thus provides a comprehensive picture of the rich band topology of bismuth, which is arguably the most important element involved in the field of topological materials.

Note Added in Proof.

After we completed our study, we became aware of refs. 54 and 55, which also propose bismuth as a potential TCI. Our analysis, however, uniquely identifies bismuth as a rotational-symmetry–protected TCI.