1 Introduction
Topological quantum materials (TQMs), which include topological insulators (TIs) and topological semimetals (TMs), have received perennial interest in condensed matter physics and material science over the past decades [
1–
9]. TIs, as the first class of TQMs, have been recognized to potentially exhibit numerous applications in various fields, owing to their intriguing topological electronic properties, which include quantum computation and high-temperature superconductivity. Although TIs have only been studied for approximately fifteen years [
3], many three-dimensional (3D) and two-dimensional (2D) TIs have already been experimentally confirmed and even more have been theoretically proposed [
10–
12]; however, desirable 3D TIs for application and especially those that are compatible with the current mainstream semiconductor industry are still rare. From the application point of view, a desirable TI would be nontoxic, easily manipulated and composed of widely used elements (e.g., Si, Ge, Sn). Unfortunately, the strong spin–orbit coupling (SOC) required by TIs generally limits the choice of materials to compounds containing heavy elements (e.g., the Bi
2Se
3 compounds [
13], ternary Heusler compounds [
14,
15]), complex heterostructures (e.g., a HgTe quantum well [
16]), or poisonous elements (e.g., mercury). Therefore, it is necessary to discover more desirable TIs for the purpose of a more effective integration with conventional semiconductor devices.
Considering these reasons, the allotropes of group-IV elements (e.g., C, Si, Ge, Sn, and Pb) seem to be good alternatives in the search for desirable TIs for at least two reasons: first, they can be efficiently integrated with various current semiconductor devices and modulated by standard semiconductor techniques; second, their numerous metastable polymorphs are very beneficial to find potentially promising TIs. In fact, the study of TQMs containing group-IV elements started almost simultaneously when graphene was successfully discovered, as it is well known that graphene is the prototype of TIs or TMs depending on whether SOC is turned on or off [
17,
18]. Because of the requirement of strong SOC in TIs, C and Si are not suitable for use in finding TIs, owing to their intrinsic weak SOC; for instance, the SOC-induced band gap in graphene and silicene is only 0.02 meV and 1.9 meV, respectively [
19]. Toxic Pb is also not a good choice despite it having the strongest SOC among the group-IV elements. Therefore, the most promising candidates are Ge and Sn.
However, the investigations on TIs composed of systems of Ge and Sn are far from complete. Only the graphene analogues of Ge and Sn (i.e., germanene and stanene) are 2D TIs and have successfully been grown on different substrates [
20]. As for the 3D structures, 3D Ge with a diamond cubic structure is an indirect band gap semiconductor at ambient conditions [
21]. Although many 3D Ge polymorphs have also been proposed using first-principles calculations or realized in experiments [
22–
27], they are usually trivial semiconductors and none of them is a TI. The diamond cubic structure of Sn, known as
α-Sn, is a low-temperature phase of Sn (stable below 13.2 °C). In contrast to its much more common metallic allotrope
β-Sn,
α-Sn is a zero-bandgap semiconductor with band inversion. Indeed,
α-Sn has received much attention as it can be engineered into different topological states, that is, a TI, Dirac semimetal and Weyl semimetal in the presence of external strains, magnetic fields and circularly polarized light [
28–
32]. Nevertheless, to the best of our knowledge, strained
α-Sn is the only elemental 3D TI composed from group-IV elements. Given the advantages and incompleteness of TIs in group-IV elements, identifying new group-IV elemental 3D TIs that are stable at room temperature is of great importance for both fundamental research and potential practical applications.
In this paper, based on first-principles calculations, we report new metastable 3D allotropes of Ge and Sn with topological electronic properties, named T5-Ge and T5-Sn. Their stabilities are thoroughly investigated through lattice dynamics and ab initio molecular dynamics (AIMD). Without considering SOC, T5-Ge and T5-Sn are highly anisotropic 3D Dirac semimetals; correspondingly, they are TIs with a negative bandgap of−21 meV and indirect bandgap of 31 meV if the SOC is switched on. Their topological properties are confirmed by calculation of the topological indices and topological surface states. To the best of our knowledge, T5-Ge/Sn are the first 3D allotropes of Ge/Sn that possess topological states in their equilibrium states at room temperature. Moreover, considering that group-IV elements usually share the same atomic structure, their C and Si isostructures are also studied. It is observed that T5-C and T5-Si are an indirect bandgap semiconductor and a direct bandgap semiconductor, respectively. Our findings not only propose the two 3D structures with topological states, but most significantly, we propose that the TIs in group-IV elements are nontoxic, easily manufactured and more compatible with the current semiconductor industry and thus present a new direction in the search for desirable TIs.
2 Results and discussion
As shown in Figs. 1(a–c), T5-Ge and T5-Sn share the same tetragonal crystal system with space group P-4M2 (No. 115) and contain five atoms in the unit cell. One of the atoms occupies the vertex of the cuboid, which corresponds to the 1a (0,0,0) Wyckoff position, and the other four atoms distributed at two equivalent lateral facets are related by S4 rotoinversion symmetry along the tetragonal axis (defined as the c axis), which represents the 4j (0, 0.2776, 0.3803) Wyckoff site, and forms a buckled square-octagon lattice in the a–b plane [see Fig. 1(a)]. The corresponding first Brillouin Zone (BZ) with high symmetry points is shown in Fig. 1(d). After full relaxation of the structures, the optimized lattice constants are a =5.63/6.49 Å and c = 5.23/5.99 Å for T5-Ge/Sn. According to symmetry, there are three kinds of bonds in T5-Ge/Sn; the first is between Ge/Sn 1 and Ge/Sn 2 (d1 = 2.53/2.91 Å), the second is between two Ge/Sn 2 atoms within one facet (d2 = 2.50/2.89 Å), and the last bond is that between two nearby Ge/Sn 2 atoms belonging to two neighboring facets (d3 = 2.54/2.93 Å), which are presented in Table 1 along with other crystalline information.
Stability examination is crucial to prove the existence of the theoretically proposed phase. Here, the dynamic, thermodynamic and mechanical stabilities of T5-Ge/Sn are thoroughly assessed through state-of-theart theoretical techniques. Their phonon spectra show there are no soft modes throughout the entire BZ, which demonstrates that they are dynamically stable [see Figs. 2(a) and (c)]. The thermodynamic stability of T5-Ge/Sn is verified through the AIMD simulations at 300 K [see Figs. 2(b) and (d)] in the framework of the NVT ensemble. One can observe that their structures are well maintained after 10 ps, with a time step of 1 fs. For the tetragonal crystal phase, the independent elastic constants are
C11 = 88 GPa/47 GPa,
C33 = 88 GPa/44 GPa,
C44 =7 GPa/3 GPa,
C66 = 4 GPa/2 GPa,
C12 = 3 GPa/5 GPa,
C13 = 12 GPa/12 GPa for T5-Ge/T5-Sn, respectively (see Table 1). Obviously, these results meet the mechanical stability criteria for tetragonal systems [
33]:
C11 >0,
C33 >0,
C44>0,
C66 >0, (
C11−
C12) >0, (
C11 +
C33−2
C13) >0, [2(
C11 +
C12) +
C33+4
C13] >0, which confirms that they are also mechanically stable. We also present their formation energies in comparison with those of diamond structures in Table S1, which shows that they are metastable group-IV elemental phases.
Noticing the similarity of the topological electronic properties between T5-Ge and T5-Sn, as well as the more pronounced SOC effect of T5-Sn, here, we mainly focus on the topological electronic properties of T5-Sn, and those of T5-Ge will be briefly discussed later. Without a special statement, the electronic properties below are calculated with the HSE06 functional. The calculated electronic band structure of T5-Sn without including SOC around the Fermi level is shown in Fig. 3(a). One observes that the valence band and conduction band cross each other to form a point along Z-A. Interestingly, one Dirac band is nearly linear but the other is parabolic-like around the Dirac point, which implies that the Dirac cone is highly anisotropic [see the inset of Fig. 3(a)]. Because we do not include the degree of spin, each band is doubly degenerate, the crossed point should be quadruple degenerate, thus T5-Sn is a 3D Dirac semimetal. Given the nontrivial topology of a 3D Dirac semimetal, the projected 2D Dirac points and Fermi arcs are expected to be observed on some surfaces of T5-Sn. As shown in Figs. 3(b) and (c), a pair of Dirac points (colored in blue) can be clearly observed at the electronic band structure of the (001) surface along [see Fig. 3(b)] and the Fermi arcs across the BZ boundary can also be clearly seen between them [see Fig. 3(c)]. Furthermore, there are no extraneous bands that pass through the Fermi level [see Fig. 3(a)], except the Dirac bands, and no extraneous trivial surface states are entangled with the nontrivial ones [see Fig. 3(b) and (c)], which is very beneficial for detecting the Dirac point and its topological nontrivial surface states in experiments.
It is known that the SOC strength of Sn is sufficient to induce observable physical properties at ambient conditions; for instance, the Dirac-point gap by SOC in stanene is approximately 101 meV [
19]. Similarly, the Dirac-point gap of T5-Sn is 31 meV after including SOC [see Fig. 3(d)]. Using the topological quantum chemistry method proposed by Bradlyn
et al.. [
34], we verified that T5-Sn is a TI with topological index Z
2 =1, which belongs to the strong topological class 1. The appearance of the topologically guaranteed surface states is one of the most significant consequences of TIs. To further unveil the topological nature of T5-Sn, we calculate its electronic states on the (001) surface, from which, the Dirac cone formed at the
point can be clearly observed [see Fig. 3(e)]. Fixed at the Fermi level, its isoenergy surface of the projected surface states on the 2D BZ of the (001) surface distinctly shows that only the two topological surface bands exist around the Fermi level, which is very helpful for experimentally finding the topological surface states.
Strain engineering has been a ubiquitous paradigm to tailor the electronic band structure and harness the associated new or enhanced fundamental properties in TQMs, just as
α-Sn demonstrates different topological states under compressive or tensile strain [
30]. From the stress-strain relationship, we find that T5-Sn can sustain at least 2% uniaxial strain along the
a and
c directions [see Fig. 4(a)]. The two kinds of strains have a different impact on their topological electronic properties. Without including SOC, the Dirac points are ruined after applying uniaxial strain along the
a direction while they remain intact along the
c direction [see Figs. 4(b–e)]. The S
4 rotoinversion symmetry would be destroyed when the strain is imposed along the
a direction, whereas the symmetry would remain unchanged along the
c direction. Therefore, the Dirac semimetal of T5-Sn is protected by S
4 rotoinversion symmetry. In addition, we calculate the topological surface states of T5-Sn under 2% compressive and tensile strains along the
c direction, which confirms that T5-Sn is still a TI under strain along the
c direction [see Figs. 4(f) and (g), their phonon spectra are shown in Fig. S1].
Despite the topological electronic properties of T5-Ge being similar to those of T5-Sn, they still demonstrate a few small differences. For example, without including SOC, T5-Ge is an even more notable anisotropic Dirac semimetal, of which one band is nearly flat but the other is linear [see Fig. 5(a)]; and more interestingly, a negative band gap of 21 meV appears after switching SOC on [see Fig. 5(b)]. As a 3D Dirac semimetal, we also calculate its surface states and Fermi arc to confirm the topology [see Fig. 5(c)], and the topological surface states of T5-Ge as a TI are given as well [see Fig. 5(d)].
Before concluding, we would like to provide a few remarks about our work. Group-IV elements have an s2p2 valence electronic configuration, which usually leads to common crystalline structures and similar chemical features, but can also lead to significant differences. Here, we also studied the T5 isostructures of C/Si and found that they are metastable phases, which is confirmed by the calculation of their phonon spectra and independent elastic constants [see Figs. 6(a, b) and Table 1]. In contrast, T5-C is a normal semiconductor with an HSE (PBE) indirect bandgap of 2.83 eV (1.82 eV), while T5-Si is a semiconductor with a direct bandgap of 1.58 eV (0.99 eV), as shown in Figs. 6(c, d). The calculated optical absorption of T5-Si as well as the air mass 1.5 solar spectral irradiance (see Fig. S2) indicate the potential of T5-Si as an absorber of sunlight in a photovoltaic device.
In addition, we would like to point out that the possibility of alloying different group-IV elements to optimize the material properties has been experimentally considered for many years, such as SiGe or GeSn alloys for better optoelectronic performance [
35,
36]. Therefore, the implication of alloys of group-IV elements with T5 structures on the electronic, optical and topological properties is very worthy of being studied in the future.
Finally, we briefly discuss the TQMs in group-IV elements (e.g., C, Si, Ge and Sn). The most distinct features for group-IV elements as TQMs are they are nontoxic, easily manipulated and compatible with current semiconductor devices. The increasing strength of SOC from C to Sn propels them to fit different kinds of TQMs. Given the strong SOC requirement of TIs, Ge, and Sn are more suitable than C and Si to induce topological insulating states, as discussed in this work; inversely, the negligible SOC in C and Si makes them better candidates as TMs because the SOC-induced band inversion mechanism in TIs is not mandatory for TMs [
37,
38]. In fact, many TMs, including nodal-line semimetals or nodal-surface semimetals, have been reported in carbon or silicon allotropes [
39–
45]. More recently, even 2D higher-order TIs have been claimed in 2D carbon materials (i.e., in graphdiyne) [
46,
47], which is beyond the paradigm of SOC-related topological states discussed here. Thus, we are deeply convinced that group-IV elements can spawn numerous desirable and novel TQMs.
3 Conclusions
In summary, through first-principles calculations, we reveal a new 3D metastable crystalline form with a tetragonal space group (named as T5), which can be adopted by all group-IV elements. Among them, T5-C and T5-Si are indirect and direct bandgap semiconductors, while T5-Ge and T5-Sn are 3D Dirac semimetals or elemental TIs depending on whether the SOC is switch off or on. T5-Ge is the first 3D Ge allotrope possessing topological states and T5-Sn is the first 3D Sn allotrope that has topological states at ambient temperature and can sustain at least a 2% external uniaxial strain along the c axis. Meanwhile, the topological electronic properties are intrinsic to T5-Sn, compared with those of α-Sn that should be induced by external methods. In conclusion, we expect that our work of discovering TIs in nontoxic, easily manipulated and widely used group-IV elements can bring a fresh perspective to the search for TQMs.
4 Computational methods
We performed our first-principles calculations within the density functional theory formalism as implemented in the Quantum ESPRESSO (QE) code [
48,
49]. The electron-ion interaction and the exchange-correlation functional between the valence electrons are described, respectively, by the projector augmented wave (PAW) and the generalized gradient approximation (GGA) parametrized by the Perdew–Burke–Ernzerhof (PBE) approach [
50]. The kinetic energy cutoffs of 200 eV, 200 eV, 300 eV, and 500 eV were employed for C, Si, Ge, and Sn, respectively. The atomic positions were optimized until the residual forces converged at 10−2 eV/Å. The BZ was sampled with 8×8×8 k-point meshes in the scheme of Monkhorst–Pack [
51]. Hybrid functional calculations of HSE06 [
52] were mainly used to obtain the electronic properties. Phonon spectra were calculated using Phonopy [
53] interfaced with QE. The topological surface properties were computed with WannierTools [
54], of which the tight-binding Hamiltonian was generated by the maximally localized Wannier functions [
55].