# Topological Properties of Electrons in Honeycomb Lattice with Kekulé Hopping Textures

###### Abstract

Honeycomb lattice can support electronic states exhibiting Dirac energy dispersion, with graphene as the icon. We propose to derive nontrivial topology by grouping six neighboring sites of honeycomb lattice into hexagons and enhancing the inter-hexagon hopping energies over the intra-hexagon ones. We reveal that this manipulation opens a gap in the energy dispersion and drives the system into a topological state. The nontrivial topology is characterized by the index associated with a pseudo time-reversal symmetry emerging from the symmetry of the Kekulé hopping texture, where the angular momentum of orbitals accommodated on the hexagonal ”artificial atoms” behaves as the pseudospin. The size of topological gap is proportional to the hopping-integral difference, which can be larger than typical spin-orbit couplings by orders of magnitude and potentially renders topological electronic transports available at high temperatures.

###### pacs:

03.65.Vf,73.43.-f,73.23.-b## I Introduction

Honeycomb lattice can host electrons with Dirac-like linear dispersion due to its crystal symmetry Wallace1947 (), and interests in questing for systems with honeycomb lattice structure flourished since the discovery of graphene produced by the scotch-tape technique Novoselov2004 (); GrapheneBook (); Geim2009 (). The Dirac dispersion and the associated chiral property of electronic wave functions accommodated on honeycomb lattice make it an ideal platform for exploring topological states QHE (); TKNN () without external magnetic field. It was shown first that a quantum anomalous Hall effect (QAHE) can be realized when complex hopping integrals among next-nearest-neighboring sites of honeycomb lattice are taken into account Haldane1988 (). Later on it was revealed that the intrinsic spin-orbit coupling (SOC) in honeycomb lattice can provide this complex hopping integrals, which drives spinful electrons into a topological state with preserved time-reversal (TR) symmetry, known as quantum spin Hall effect (QSHE) KaneRMP (); ZhangRMP (); Kane2005a (); Kane2005b (). Quite a number of activities have been devoted towards realizing topological states in electron systems on honeycomb lattice, such as QAHE by straining GeimNPhys2010 (); Gomes2012 (), twisting Hunt2013 (); He2013 () and decorating graphene QNiuPRB2010 (), and QAHE with spin-polarized edge currents in terms of the antiferromagnetic exchange field and staggered electric potential LiangNJP2013 (). Honeycomb lattice can also be tuned to support topological states in photonic crystals Haldane2008 (); Wu2015 () and cold atoms Tarruell2012 ().

In the present work, we explore possible topological properties in honeycomb lattice by introducing a Kekulé texture in hopping energy between nearest-neighboring (NN) sites. We take a hexagonal primitive unit cell and view the honeycomb lattice as a triangle lattice of hexagons [see the dashed red line in Fig. 1(a)]. When the inter-hexagon hopping is tuned to be larger than the intra-hexagon one , a topological gap is opened at the point accompanied by a band inversion between orbitals with opposite spatial parities accommodated on hexagons [see Fig. 1(b)]. A pseudo-TR symmetry associated with a pseudospin degree of freedom and Kramers doubling in the emergent orbitals are revealed based on point group symmetry, which generates the topology. For experimental implementations, we discuss that, along with many other possibilities, the molecular graphene realized by placing carbon monoxides (CO) periodically on Cu [111] surface Gomes2012 () is a very promising platform to realize the present idea, where the Kekulé texture can be controlled by adding extra CO molecules.

The paper is organized as follows. In Sec. II, we start with a symmetric Kekulé hopping texture among NN sites in honeycomb lattice, where a pseudo-TR symmetry and pseudospin emerge. In Sec. III, we derive an effective low-energy Hamiltonian of the system and reveal the topology. The topological properties of the system is demonstrated by the edge state and associated Hall and longitudinal conductances in sections IV and V. We discuss a promising experimental platform for realizing the topological state in Sec. VI. Concluding remarks are presented in Sec. VII.

## Ii Kekulé hopping texture and emergent orbitals

We start from a spinless tight-binding Hamiltonian on honeycomb lattice

(1) |

with the annihilation operator of electron at atomic site with on-site energy satisfying anti-commutation relations, and running over NN sites inside and between hexagonal unit cells with hopping energies and respectively [see Fig. 1(a)]. The orbitals are considered to be the simplest one without any internal structure, same as the electron of graphene. Below we are going to detune the hopping energy while keeping constant, and elucidate possible changes in the electronic state. In this case, the pristine honeycomb lattice of individual atomic sites is better to be considered as a triangular lattice of hexagons, with the latter characterized by symmetry.

Let us start with the Hamiltonian within a single hexagonal unit cell

(2) |

where [see Fig. 1(a)]. The eigenstates of Hamiltonian are given by

(3) |

with eigenenergies and respectively, up to normalization factors. As shown in Fig. 1(b), the emergent orbitals accommodated on the hexagonal “artificial atom” take the shapes similar to the conventional , , and atomic orbitals in solids.

A pseudo-TR symmetry operator can be composed in the present system with symmetry: with the complex conjugate operator and , where is the Pauli matrix. It can be checked straightforwardly that corresponds to rotation for orbitals and rotation for orbitals given in Eq. (3), which yields in the space formed by the and orbitals Wu2015 (). Therefore, the pseudo-TR symmetry satisfies the relation , which is same as that for fermionic particles even though the spin degree of freedom of electron has not been considered here. This indicates that electrons acquire a new pseudospin degree of freedom in the present system as far as the low-energy physics is concerned.

Explicitly the wave functions carrying pseudospins are given by the emergent orbitals with eigen angular momentum

(4) |

Distinguished from the intrinsic spin, the pseudospin is directly related to the chiral current density on the hexagon. For a lattice model, the current density between two sites is given by . The current distributions evaluated using wave functions in Eq. (4) for the pseudospin-up and -down states are shown in Figs. 2(a) and (b) with anticlockwisely and clockwisely circulating currents. By considering the hexagonal artificial atoms composed by six sites in honeycomb lattice, one harvests states with angular momentum merely from simple orbitals, such as electrons in graphene. The pseudo-TR symmetry is supported by the crystal symmetry, sharing the same underlying physics with the topological crystalline insulator FuPRL2011 (). However, crystal-symmetry-protected topological insulators addressed so far need strong SOCs to achieve band inversions Hsieh2012 (); Dziawa2012 (); Xu2012 (), which is different from the present approach as revealed below.

## Iii Topological phase transition

We calculate the energy dispersion of Eq. (1) for several typical values of (hereafter the on-site energy is put as without losing generality). As shown in Fig. 2, there are two two-fold degeneracies at the point corresponding to the two two-dimensional (2D) irreducible representations of point group. Projecting the wave functions for onto the orbitals given in Fig. 1(b), it is found that the topmost two valance bands show the character of orbitals whereas the lowest two conduction bands behave like orbitals [see Fig. 2(c)], with the order in energy same as those listed in Eq. (3). For , the and bands become degenerate at the point and double Dirac cones appear [see Fig. 2(d)], which are equivalent to the ones at and points in the unfolded Brillouin zone of honeycomb lattice with the rhombic unit cell of two sites. When increases further from , a band gap reopens at the point. As shown in Fig. 2(e) for , the valence (conduction) bands are now occupied by () orbitals around the point, opposite to the order away from the point, and to that before gap closing. Therefore, a band inversion between and orbitals takes place at the point when the inter-hexagon hopping energy is increased across the topological transition point , namely the pristine honeycomb lattice.

We can characterize the topological property of the gap-opening transition shown in Fig. 2 by a low-energy effective Hamiltonian around the point. Since the bands near Fermi level are predominated by and orbitals, it is sufficient to downfold the six-dimensional Hamiltonian associated with the tight-binding model (1) into the four-dimensional subspace . The second term in Eq. (1) is then simply given by

(5) |

in the four-dimensional subspace. Contributions from the third term in Eq. (1) to the effective Hamiltonian should be evaluated perturbatively. First, we list the inter-hexagon hoppings in terms of 66 matrices and with

on the basis of . Following the standard procedures Sakurai (), they can be projected to the subspace spanned by

(6) |

With Fourier transformations of matrices in Eqs. (5) and (6), one obtains the effective low-energy Hamiltonian on the basis in the momentum space. Expanding around the point up to the lowest-orders of , one arrives at

(7) |

with

(8) |

where , , , is a zero matrix, and is the lattice constant of the triangular lattice. For , the Hamiltonians and in Eq. (8) are the same as the well-known one for honeycomb lattice around the and points NetoRMP2009 (), where the quadratic terms of momentum in the diagonal parts can be neglected.

For , however, the quadratic terms are crucially important since they induce a band inversion Bernevig2006 (), resulting in the orbital hybridization in the band structures denoted by the rainbow colors in Fig. 2(e). Associated with a skyrmion in the momentum space for the orbital distributions in the individual pseudospin channels, a topological state appears characterized by the topological invariant Kane2005a (); Kane2005b (); Wu2015 (); FuPRB2007 (). It is clear that for there is no band inversion taking place and thus the band gap is trivial as shown in Fig. 2(c).

It is worthy noticing that, comparing with the Kane-Mele model for the honeycomb lattice Kane2005a (); Kane2005b (), the mass term in Eq. (8) can be considered as an effective SOC associated with the pseudospin, namely . For , a moderate Kakulé texture in hopping energy, the effective SOC is approximately 3000 times larger than the real SOC in magnitude in graphene where 0.1meV and eV. The huge effective SOC is due to its pure electronic character as compared with the intrinsic SOC originated from the relativistic effect. This is one of the fantastic aspects of the present approach, which renders a topological gap corresponding to temperature of thousands of Kelvin.

## Iv Topological edge states and associated conductances

### iv.1 Topological edge states

In order to check the edge state in the present system, we consider a ribbon of hexagonal unit cells of with its two edges cladded by hexagonal unit cells of . As can be seen in Fig. 3(a), additional states appear in the bulk gap as indicated by the red solid curves carrying double degeneracy. Plotting the spatial distribution of the corresponding wave functions, we find that the in-gap states are localized at the two interfaces between topological and trivial regions [see Fig. 3(b)]. As displayed in Fig. 3(c) [(d)], there is an excess upward (downward) edge current in the pseudospin-up (-down) channel associated with the state indicated by the red (green) dot in Fig. 3(a).

### iv.2 Hall and longitudinal conductances

At the interface between topological and trivial regimes, the crystal symmetry is reduced from to , which breaks the pseudo-TR symmetry in contrast to the real TR symmetry. As the result, a mini gap of [unnoticeable in the scale of Fig. 3(a)] opens in the otherwise gapless edge state at the point due to the coupling between two pseudospin channels. In order to quantitatively check possible backscatterings caused by this mini gap, we perform calculations on the longitudinal and Hall conductances based on a Hall bar system as sketched in Fig. 4(a). It is clear that the current injected from the left electrode divides itself into two branches according to the pseudospin states, namely pseudospin up (down) electrons can flow only in the upper (lower) edge of the Hall bar. By matching wave functions at the interfaces between the six semi-infinite electrodes and the topological scattering region AndoPRB1991 (); GrothNJP2014 (), one can evaluate the transmissions of plane waves scattered among all the six leads, and then the longitudinal and Hall conductances, and respectively, by the Landauer-Büttiker formalism Landauer1999 (), where and are the longitudinal and transverse resistances respectively. Similar to the case of QSHE with magnetic impurities TkachovPRL2010 (), in the present system the values of conductivity deviate from the quantized ones when Fermi level falls in the mini gap of as shown in Fig. 4(b). It is noticed, however, that both and heal quickly after several periods of oscillations that come from interferences between the two pseudospin channels. It is emphasized that almost perfectly quantized conductances and Bernevig2006 (); Konig2007 () are achieved for Fermi level beyond up to the bulk gap edge at , where the edge states with almost perfect linear dispersions hardly feel the existence of the mini gap and essentially no appreciable backscattering exists.

Now we discuss the hopping integral dependence of the longitudinal conductance. The size of scattering region is same as in Fig. 4(a) and fixed for all cases. As displayed in Fig. 5(a), saturates to the quantized value as expected for a topological state for all the cases with and in the topological region (whereas and in the trivial region correspondingly) when Fermi level is set away from the mini gaps, accompanied by oscillations due to interferences between two pseudospin channels.

Here we show the sample size dependence of the longitudinal conductance. We fix inter-hexagon hopping integrals at and in the topological and trivial regions respectively. As displayed in Fig. 5(b), saturates in all cases to the quantized value when Fermi level is shifted away from the mini gap. The topological edge transport remain unchanged when the topological sample becomes large.

## V Real spin and QSHE

In addition to the pseudospin, the real spin degree of freedom also contributes to transport properties in realistic systems. In absence of the real SOC, the results presented in Fig. 4 remain exactly the same, with an additional double degeneracy due to the two real spin and thus .

An intrinsic SOC is induced when next-nearest-neighbor hoppings in honeycomb lattice are taken into account Kane2005a (); Kane2005b (). The low-energy Hamiltonian around the point in Eq. (7) is then modified to

(9) |

with

where and stand for real spin-up and -down states respectively. Therefore, for up real spin SOC enhances (suppresses) the topological gap in the pseudospin up (down) channel presuming [see the left and central panels of Fig. 6(a)]. As far as , the system remains the topological state associated with the pseudospin, where electrons with up pseudospin and down pseudospin counter propagate at edges, both carrying on up and down real spins. The longitudinal conductance saturates to as displayed in Fig. 6(b).

When SOC is increased to , the pseudospin down (up) channel with real spin up (down) is driven into a semi-metallic state with zero band gap and the Dirac dispersion appears at the point. When SOC goes beyond , this Dirac dispersion opens a gap accompanied by a topological phase transition. The system now takes a QSHE state where at edges electrons with up real spin and pseudospin propagate oppositely to electrons with down real spin and pseudospin. Evaluating the longitudinal conductance, one finds that is quantized exactly to (see Fig. 6(b)), and as shown in the right panel of Fig. 6(a) there is no mini gap in the edge states, as protected by real time-reversal symmetry Kane2005a (); Kane2005b ().

## Vi Discussions

Much effort has been devoted towards realizing the Dirac-like energy dispersion in artificial honeycomb lattices Polini2013 (), ranging from optical lattices Tarruell2012 (); Wunsch2008 () to 2D electron gases modulated by periodic potentials Gomes2012 (); GibertiniPRB2009 (); SGLouie2009 (). All these systems provide promising platforms for realizing topological properties by detuning effective hopping energy among NN sites either by modulating muffin-tin potentials or bond lengthes periodically. To be specific, here we focus on how to achieve a topological state on the Cu [111] surface modulated by triangular gates of carbon monoxide (CO) molecules Gomes2012 (). When extra CO molecules are placed at specific positions over the pristine molecular graphene, the bonds of the hexagons surrounding them are elongated since the CO clusters enhance local repulsive potentials and push electrons away from them, which reduces the corresponding electron hopping energies SGLouie2009 (). It is extremely interesting from the present point of view that Kekulé hopping textures have already been achieved in experiments Gomes2012 (). We propose to place extra CO atoms in the pattern displayed in Fig. 7(a), where the intra-hexagon hopping energy (green thin bonds) surrounding the CO clusters is reduced and the inter-hexagon hopping energy (red thick bonds) is enhanced relatively. According to the above discussions, the system displayed in Fig. 7(a) with should take a topological state. The Kekulé hopping texture in Fig. 7(b), dual to that shown in Fig. 7(a), was realized in experiments Gomes2012 (), where the system takes a topologically trivial state because .

The underlying idea of the present scheme for achieving the topological state is to create artificial orbitals carrying on opposite orbital angular momenta and parities with respect to spatial inversion symmetry, and to induce a band inversion between them by introducing a Kekulé hopping texture on honeycomb lattice. In the sense that it does not require SOC, the present state may be understood as a quantum orbital Hall effect. The topological properties can also be extended to photonic crystals Wu2015 (), cold atoms, and even phonon systems where sound waves form band due to periodic configurations elastic materials.

## Vii Conclusion

We propose to derive topological properties by modulating electron hopping energy between nearest-neighbor sites of honeycomb lattice. Due to the Kekulé hopping texture with symmetry, atomic-like orbitals emerge, which carry a pseudospin degree of freedom characterizing a pseudo time-reversal symmetry and Kramers doubling. We reveal that the effective spin-orbit coupling associated with the pseudospin degree of freedom can be larger than the intrinsic one by orders of magnitude because of the pure electronic origin. The present work offers a new possibility for achieving topological properties and related novel quantum properties and functionalities at high temperatures.

###### Acknowledgements.

The authors thank Y.-Y. Wang, L. You, Z.-F. Xu, T. Taniguchi, K. Watanabe, L. Jiang and S. Mizuno for stimulating discussions. This work was supported by the WPI Initiative on Materials Nanoarchitectonics, Ministry of Education, Culture, Sports, Science and Technology of Japan.## References

- (1) P. R. Wallace, Phys. Rev. 71, 622 (1947).
- (2) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 666 (2004).
- (3) M. I. Katsnelson, Graphene Carbon in Two dimensions (Cambridge University Press, 2012).
- (4) A. K. Geim, Science 324, 1530 (2009).
- (5) K. v. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. Lett. 45, 494 (1980).
- (6) D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. d. Nijs, Phys. Rev. Lett. 49, 405 (1982).
- (7) F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988).
- (8) M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010).
- (9) X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011).
- (10) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 (2005).
- (11) C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 (2005).
- (12) F. Guinea, M. I. Katsnelson, and A. K. Geim, Nature Phys. 6, 30 (2010).
- (13) K. K. Gomes, W. Mar, W. Ko, F. Guinea, and H. C. Manoharan, Nature 483, 306 (2012).
- (14) B. Hunt et al., Science 340, 1427 (2013).
- (15) W. Yan, W.-Y. He, Z.-D. Chu, M. Liu, L. Meng, R.-F. Dou, Y. Zhang, Z. Liu, J.-C. Nie, and L. He, Nat. Commun. 4, 2159 (2013).
- (16) Z. Qiao, S. Yang, W. Feng, W.-K. Tse, J. Ding, Y. Yao, J. Wang, and Q. Niu, Phys. Rev. B 82, 161414 (2010).
- (17) Q.-F. Liang, L.-H. Wu, and X. Hu, New J. Phys. 15, 063031 (2013).
- (18) F. D. M. Haldane, and S. Raghu, Phys. Rev. Lett. 100, 013904 (2008).
- (19) L.-H. Wu, and X. Hu, Phys. Rev. Lett. 114, 223901 (2015).
- (20) L. Tarruell, D. Greif, T. Uehlinger, G. Jotzu, and T. Esslinger, Nature 483, 302 (2012).
- (21) L. Fu, Phys. Rev. Lett. 106, 106802 (2011).
- (22) T. H. Hsieh, H. Lin, J. Liu, W. Duan, A. Bansil, and L. Fu, Nat. Commun. 3, 982 (2012).
- (23) P. Dziawa et al., Nature Mater. 11, 1023 (2012).
- (24) S.-Y. Xu et al., Nat. Commun. 3, 1192 (2012).
- (25) J. J. Sakurai, Modern Quantum Mechanics (Addison Wesley, 1985).
- (26) A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 (2009).
- (27) B. A. Bernevig, T. L. Hughes, and S.-C. Zhang, Science 314, 1757 (2006).
- (28) L. Fu, and C. L. Kane, Phys. Rev. B 76, 045302 (2007).
- (29) T. Ando, Phys. Rev. B 44, 8017 (1991).
- (30) C. W. Groth, M. Wimmer, A. R. Akhmerov, and X. Waintal, New J. Phys. 16, 063065 (2014).
- (31) Y. Imry, and R. Landauer, Rev. Mod. Phys. 71, S306 (1999).
- (32) G. Tkachov, and E. M. Hankiewicz, Phys. Rev. Lett. 104, 166803 (2010).
- (33) M. König, S. Wiedmann, C. Brüne, A. Roth, H. Buhmann, L. W. Molenkamp, X.-L. Qi, and S.-C. Zhang, Science 318, 766 (2007).
- (34) M. Polini, F. Guinea, M. Lewenstein, H. C. Manoharan, and V. Pellegrini, Nature Nanotech. 8, 625 (2013).
- (35) B. Wunsch, F. Guinea, and F. Sols, New J. Phys. 10, 103027 (2008).
- (36) M. Gibertini, A. Singha, V. Pellegrini, and M. Polini, Phys. Rev. B 79, 241406 (2009).
- (37) C.-H. Park, and S. G. Louie, Nano Lett. 9, 1793 (2009).