Strain-induced quantum spin Hall effect in methyl-substituted germanane GeCH3.

Quantum spin Hall (QSH) insulators exhibit a bulk insulting gap and metallic edge states characterized by nontrivial topology. We investigated the electronic structure of an isolated layer of methyl substituted germanane GeCH3 by density functional calculations (DFT), and its dynamic stability by phonon dispersion calculations. Our results show that an isolated GeCH3 layer has no dynamic instability, and is a QSH insulator under reasonable strain. This QSH insulator has a large enough band gap (up to 108 meV) at 12% strain. The advantageous features of this QSH insulator for practical room-temperature applications are discussed.

Recently, two-dimensional (2D) topological insulators (TIs), namely, QSH insulators1234567, have attracted much attention. These insulators possess a bulk insulting gap but metallic edge states. The edge states exhibit a linear dispersion near the Fermi energy giving rise to a Dirac Fermion behavior, which can lead to an efficient transport of charge carriers. The Dirac Fermions are protected by time-reversal symmetry and hence are robust against nonmagnetic defects, disorder and surface contamination4567. The QSH effect was first proposed by Kane and Mele for graphene, considering the opening of a band gap at the Dirac point by spin-orbit coupling (SOC)1. However, this proposal was shown to be unrealistic89, because the strength of SOC in graphene is extremely small. Subsequently, Beneviget et al. proposed that the quantum well semiconductor HgTe-CdTe would be a QSH insulator5, which was soon confirmed experimentally by observing ballistic edge channels3. Nevertheless, these quantum-well-based QSH insulators have a very small bulk band gap so that the QSH effect is observed only at low temperatures (below 10 K)5. Besides, the HgTe-CdTe quantum wells as well as their analogues1011 have other serious limitations such as toxicity, difficulty of preparation, and incompatibility with current silicon-based electronic technology.

In searching for practical QSH insulators, two key ingredients are essential: (1) the materials must have a large bulk gap to ensure the realization of a QSH effect at room temperature, and (2) they must be compatible with the current silicon-based electronic technology. The latter condition narrows down our attention to the group IV honeycomb lattices. Silicene can be a nontrivial QSH insulator in its native structure, and can be observed in an experimentally accessible low temperature12. However, up to now, no free-standing silicene has been prepared experimentally, and its QSH effect can be observed only at low temperature13. Closely related to silicene in terms of compatibility and structure are germanene and stanene, which have been proposed as QSH insulators with a relative large band gap induced by spin-orbit coupling (SOC)12131415. However, they have not yet been realized experimentally. As a stable free-standing species that is experimentally accessible, only the germanane, GeH, can be well produced1617. When strained, germanane can display a QSH effect18. However, it begins to become amorphous upon annealing at 75°C and becomes completely amorphous above 175°C16, which is undesirable because most transistors should operate up to hundreds of degrees °C. Recently, a methyl-substituted germanane, GeCH3, has been synthesized19 (Fig. 1a, b). In contrast to the case of germanane GeH, it is found that GeCH3 has considerably enhanced thermal stability19. Therefore, it is worthwhile to probe whether GeCH3 has a QSH effect.

Figure 1

(a) Top and (b) side views of a GeCH3 monolayer. The red, blue and green balls denote Ge, C and H atoms, respectively. (c) Phonon band dispersion relations calculated for a GeCH3 monolayer. (d) Energy gap calculated for a GeCH3 monolayer as a function of strain.

In this work, we study the electronic states of a single layer of GeCH3 by density functional calculations to show that reasonable strain could induce the s-p band inversion in the GeCH3 layer and makes it a QSH insulator, and that SOC plays an important role in opening a local band gap at the crossing points near the Fermi level. In the band-inverted GeCH3 layer under 12% strain, the band gap is 108 meV, so that this QSH insulator can be useful in an experimentally accessible temperature regime.

Results

Stability of GeCH3 monolayer

GeCH3 crystallizes in a hexagonal structure in which GeCH3 layers stack along the c axis19. GeCH3 is thermally stable up to 250°C, so it is quite promising as a building block for future nanodevices15. Due to its large interlayer separation the electronic interaction between adjacent GeCH3 layers is extremely weak19, contrary to the case of graphite, which has a substantial interlayer overlap between the C 2pz orbitals, and MoS2 at the Γ point, which has a substantial interlayer overlap between the S 3pz orbitals20. The interlayer binding energy is calculated to be 5.7 meV/Å2 for GeCH3, which is even smaller than that of graphite (9.3 meV/Å2). Such weak coupling between GeCH3 layers provides a natural cleaving plane. Thus, bulk GeCH3 can be easily exfoliated to produce stable GeCH3 monolayers. The stability of GeCH3 monolayer is further confirmed by the phonon dispersion curves, presented in Fig. 1c, which shows that all branches of the phonon dispersion curves have positive frequencies, and no imaginary phonon modes are found. This confirms the stability of a GeCH3 monolayer.

Electronic structures of GeCH3 monolayer

As shown in Fig. 1a, b, each sp3-hybridized Ge atom is bonded to three Ge atoms and one CH3 group such that the buckled sheet of Ge atoms is sandwiched between two sheets of CH3 groups. Each GeCH3 monolayer has a hexagonal unit cell, and there are two Ge atoms and two CH3 groups per unit cell. In the optimized structure, the lattice constant is 3.954 Å with Ge-Ge = 2.415 Å and Ge-C = 1.972 Å.

The band structures calculated for a GeCH3 monolayer without and with SOC are plotted in Fig. 2. Without SOC, the band structure has a direct band gap (1.07 eV) at the Γ point (Fig. 2a). The Ge-4s and Ge-4p states are located near Fermi level, with the Ge-4p states lying lower than the Ge-4s states. The valence band maximum (VBM) is doubly degenerate at Γ because they are derived from the Ge-4px,y states. When SOC is included, the degenerate level is split by a gap of 193 meV (Fig. 2e, f). Except for this, no significant change is caused by SOC for the band structure of a GeCH3 monolayer; the conduction band minimum (CBM) at Γ is still occupied by the Ge-4s states while the two valence bands near the Fermi level at Γ point are still dominated by Ge-4p states. Thus, no inverted band order is observed in this 2D lattice, showing that a GeCH3 monolayer is a trivial insulator.

Figure 2

The band structures calculated without SOC for an isolated GeCH3 monolayer structures with (a) 0%, (b) 4%, (c) 8% and (d) 12% strain. The band structures calculated with SOC for an isolated GeCH3 monolayer structures with (e) 0% and (g) 12% strain. (f) and (h) show zoomed-in views of (e) and (g), respectively. The black and olive lines denote the Ge-4s and Ge-4p orbitals, respectively. The horizontal dashed lines indicate the Fermi level.

Strain effects

Mechanical strain can induce the s-p band inversion and hence lead to a nontrivial topological state212223. For example, under external strain, the layered compound InSe becomes a three dimensional TI22. In addition, strain can be readily exerted on nanosheets. Thus, one might consider if a GeCH3 monolayer becomes a QSH insulator under mechanical strain. We probe this question by introducing 0–12% strain into a GeCH3 monolayer maintaining all crystal symmetries and the overall honeycomb-like structure. Here we define the strain as ε = Δa/a, where a and a = Δa + a are the lattice constants of the unstrained and strained cell, respectively. The GeCH3 lattice is expanded uniformly along both zigzag and armchair directions in the ab plane and is relaxed along the c axis. To clarify the stability of strained GeCH3, we also calculate phonon dispersion curves of GeCH3 under different strain. The corresponding results are shown in Fig. S1. We find that, the strain can reach to 12% without destroying the feature that imaginary phonon mode is absent in the phonon dispersion curves. This suggests that strained systems maintain their stability within the strain range of 12%.

The evolution of the band gap of a GeCH3 monolayer as a function of the applied mechanical strain, presented in Fig. 1d, shows that with increasing the mechanical strain, the CBM is driven continuously to shift downward to the Fermi level thus significantly decreasing the energy gap. Throughout the changes under strain, it has a direct band gap with the VBM and CBM both located at the Γ point. The band structures of a GeCH3 monolayer under the strain of 4%, 8% and 12% are presented in Fig. 2b–d. From Fig. 2, we can see that, with 12% strain, the GeCH3 monolayer becomes a semimetal since the VBM and CBM is degenerate at Γ point, and the s-p band inversion occurs around Fermi level; for the other two cases with strain below ε = 12%, the lowest-conduction band is represented by the Ge-4s states, and the highest two valence bands by the Ge-4p states. For the GeCH3 monolayer with ε = 12%, the lowest-conduction band and the highest-valence band are represented by the Ge-4p orbitals, but the next-highest-valence band by the Ge-4s states. This establishes the existence of band inversion for the GeCH3 monolayer with ε = 12%. The inversion of band orders is a strong indication for the formation of nontrivial topological phases. The s-p band order inversion in our case is caused predominantly by strain, but the SOC is still indispensable because it lifts the band degeneracy at Γ, as shown in Fig. 2h. With SOC, a finite energy gap of 161 meV opens at Γ (an indirect gap of 108 meV) while the band order stays unchanged. This large SOC-induced nontrivial band gap reflects the fact that the states near the Fermi level are made up of σ orbitals rather than π orbitals182425. The SOC-induced nontrivial energy gap is significantly greater than kBT at room temperature ( = 26 meV), so the QSH effect can be measured at room temperature. Thus GeCH3 might have useful applications in fields such as phase-change random access memory, quantum computing and spintronics.

The topological phase transition in a GeCH3 monolayer is further confirmed by calculating the Z2 topological invariant before and after the band inversion. Nontrivial topological insulators can be distinguished from normal insulators by nontrivial Z2 topological invariant associated with the “bulk” electronic structure: Z2 = 1 characterizes a nontrivial band topology, but Z2 = 0 a trivial band topology. All the structures of a GeCH3 monolayer considered here possess inversion symmetry. According to the method developed by Fu and Kane26, the existence of inversion symmetry allows us to calculate the Z2 topological invariant directly from knowledge about the parity of the Bloch wave functions for all the filled bands at all the time-reversal-invariant momenta (TRIM). For a GeCH3 monolayer, there are four TRIM points, namely, one Γ point and three M points in the Brillouin zone. The Z2 invariant is determined by the quantity , where is the parity eigenvalue of the occupied band at the TRIM points ki. The Z2 invariant is then given by the product of all the δi's, . For the GeCH3 monolayer without strain or with strain below 8%, they have trivial band topology with the topological invariant Z2 = 0. When the strain is 12%, the calculated Z2 topological invariants is 1, which confirms that the strain induces a topological phase transition.

The 2D nontrivial insulating states in the GeCH3 monolayer with ε = 12% should support an odd number of topologically protected Dirac-like edge states connecting the conduction and valence band edges at certain k-points. To see the nontrivial topological features of the GeCH3 monolayer with ε = 12% explicitly, we performed calculations of its edge-state band structure by constructing an zigzag nanoribbon structure with symmetric edges, for which we saturate the dangling bonds of all the edge atoms by hydrogen atoms. The width of the nanoribbon structure adopted here is 8.2 nm, which is enough to separate the interaction between its two edges. The corresponding band structure of the nanoribbon plotted in Fig. 3a shows that the nontrivial topological gapless edge states connecting the band edge states and form a single Dirac cone at the Γ point. For this nanoribbon structure, the two edges are the same, leading to two energetically degenerate Dirac cones located at the opposite edges. At a given edge, two counter-propagating edge states exhibit opposite spin-polarizations, a manifestation of the one-dimensional nature of the system. Thus, the surface-state calculation is consistent with the above parity analysis, and confirms conclusively the topologically non-trivial nature of this material.

Figure 3

(a) Edge states of a GeCH3 monolayer calculated by using a nanoribbon model (see the text). Schematic diagrams showing the relative positions at Γ of the Ge-4s and Ge-4p states of a GeCH3 monolayer near the Fermi level for the case of (b) 0% and (c) 12% strain in the absence of SOC. The horizontal dashed lines indicate the Fermi level.

Let us now consider how strain influences the band gap and the s-p band order. There are two factors affecting the band order; one is SOC, and the other is the strain. The latter factor has been discussed in many previous studies1518222327. The previous showed that the band inversion in a GeCH3 monolayer can be driven by mechanical strain instead of SOC. Thus the physics of band inversion can be understood without considering SOC, while the effect of SOC is to create an energy gap afterwards. As discussed above, the valence and conduction bands away from the Γ point are well separated, therefore, one needs to focus on the bands near the Γ point. The relevant bands are dominated by the Ge-4s and Ge-4px,y states. Due to the chemical bonding between adjacent Ge atoms, these states are split into bonding and antibonding states, which we denote as and , with the subscript ± representing the parity (Fig. 3b, c). For the GeCH3 monolayer without strain (Fig. 3b), the bands near the Fermi level are contributed by the and orbitals, with the orbital being above the orbital. For the GeCH3 monolayer with strain, increasing the lattice constant weakens the interaction between the Ge atoms. This decreases the splitting between the bonding and antibonding states, which lowers the state and raises state. Consequently, depending on the strength of strain, the energy gap of GeCH3 monolayer can be continuously tuned. At 12% strain, the level shifts above the level, which implies the normal band order is inverted under this condition, as depicted in Fig. 3c. Then the SOC then can lift the band degeneracy and create an energy gap at the Γ point, signifying that a GeCH3 monolayer with ε = 12% is a QSH insulator.

Discussion

Our density functional calculations show that a topological phase transition can be induced in a recently synthesized GeCH3 monolayer by strain, and our phonon dispersion analysis shows that this layer is stable. The analysis of the Z2 topological invariant confirms the topological insulating nature of the insulating state of the GeCH3 monolayer with 12% strain. This QSH insulator has a large enough band gap that exceeds the thermal energy at room-temperature, making them suitable for room-temperature spintronics applications. It should be mentioned that there are many different ways to apply strain. For example, the strain can be realized by exerting mechanical forces to suspended samples, while very strong strain can also be introduced in thin films via the lattice mismatch with the substrate. For deposition onto substrates, one main concern is whether interface effect on the topological properties can be ignored. To this end, lots of efforts are still needed. However, as the first step for its realization, it is necessary to investigate the topological nature of free-standing GeCH3 as discussed in the present work. In actual conditions, there are many other factors that may affect the topological properties of GeCH3, so a more comprehensive study concerning these problems should be carried out in further works.

It is of interest to comment on the advantages of this GeCH3 QSH insulator over other QSH insulators. First, the GeCH3 monolayer is made up of the commonly used elements and is thermally stable up to 250°C, offering compatibility with current silicon-based electronic technology. Second, on the experimental side, GeCH3 has been synthesized whereas other 2D group IV lattice-based QSH insulators with comparable nontrivial gap (such as silicene, germanene, stanene, silicane, and stanane) have not yet been produced in the free-standing form. Third, the CH3 terminated surface can be resilient against oxidation, whereas the -H terminated surface (such as germanane and silicane) are extremely air sensitive. Finally, strain can switch between a normal insulating edge state and an almost linear dispersion for the metallicg edge state. Therefore, the GeCH3 monolayer allows one to construct a potential QSH field-effect transistor.

The peculiar judgment of the nontrivial topological nature is the existence of band inversion. Standard DFT calculations tend to underestimate the gap, and this could give an enhanced and not reliable estimate of the strain for band inversion. To overcome this drawback, we repeat the simulations using the Heyd–Scuseria–Ernzerhof (HSE) functional28, and the corresponding results are plotted in Fig. S2. We find that the band inversion occurs at 13% strain based on HSE, which increases only a few percent with respect to the LDA case. From this we can argue that standard DFT is still a good approximation to study the band inversion feature in GeCH3 monolayer.

Methods

Our density functional calculations for an isolated GeCH3 layer were performed using the projector augmented wave (PAW) method2930 implemented in the Vienna ab initio Simulation Package (VASP)3132 with the cutoff energy of 450 eV for the plane-wave basis set and the local density approximation (LDA)33 to describe the exchange-correlation functional. Repeated image interactions were made negligible by including a vacuum layer of 18 Å in the simulations. The Monkhorst–Pack set34 of 17 × 17 × 1 k-points was used to sample the Brillouin zone for both geometry optimizations and total energy calculations. The lattice constants and the atom coordinates were fully optimized until the forces on each ion are reduced below 0.02 eV/Å. The optimization with the inclusion of van der Waals interactions are also tested, and the results show the main conclusions obtained here are not affected (see Fig. S3). Complementary calculations were performed to determine the phonon dispersion relations using the finite displacement method with the density functional perturbation theory (DFPT) as implemented in the CASTEP code3536.

Author Contributions

Y.D. and B.B.H. designed the research. Y.D.M. and W.W. performed the density functional theory calculations. Y.D.M., Y.D. and M.-H.W. prepared the manuscript. All authors contributed to the scientific discussions.