** Graphene** Vol.3 No.2(2014), Article ID:45409,4 pages DOI:10.4236/graphene.2014.32004

Opening of a Gap in Graphene Due to Supercell Potential: Group Theory Point of View

E. Kogan, M. Kaveh

Department of Physics, Bar-Ilan University, Ramat-Gan, Israel

Email: Eugene.Kogan@biu.ac.il, Moshe.Kaveh@biu.ac.il

Copyright © 2014 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 25 February 2014; revised 22 March 2014; accepted 17 April 2014

ABSTRACT

We analyze in the framework of the space group theory the change of the dispersion law in grapenein and the vicinity of the (former) Dirac points due to application of supercell potential with the space priodicity and the same point symmetry as graphene.

**Keywords:**Graphene, Group Theory, Space Groups

1. Introduction

Graphene is a two-dimensional crystal of carbon atoms, which form a honeycomb lattice with the point symmetry described by the group. The first Brillouin zone (BZ) has a hexagonal form, and the conduction band touches the valence band in six BZ corners which form two non-equivalent triads of BZ corners, and. One of the routes toward tailoring the electronic properties of graphene is through the adsorption of metals [1] [2] . Recently, several types of adatoms were used to dope graphene in attempts to tailor properties of graphenebased devices [3] -[7] . The gap opening in the high symmetry points for the hexagonal lattices due to interaction with the interface was considered in [8] . In Ref. [9] , it was shown that, using ab initio density functional ory calculations, the adsorption of an alkali-metal submonolayer ongraphene occupying every third hexagon of the honeycomb lattice in a commensurate arrangement induces an energy gap in the spectrum of graphene. We decided to analyze this opening of the gap in the framework of the group theory. In our previous publications [10] -[12] we summed up the classification of the energy bands in graphene on the basis of the point group analysis. The only fact from that analysis we need in the present work, is the fact that the little group at the point is and the bands and realize representation of the group. The same can be said about the point. This fact by itself means that the bands and touch each other at the points and, and the electron states in the vicinity of these points are described by massless Dirac equation [11] .

For the purpose of the present paper, we should put the abovementioned fact into the framework of the theory of the space group symmetry [13] [14] , from which we will need only a few basic ideas. According to the theory of the space group symmetry, the bands and should be considered at the points and (these two points can be considered as the stars of the wave vector, and designated together, thus realizing a 4-dimensional representation of the group. Due to the identity

(1)

any element of the group can be presented as an element of the group, or as a product of and such element. Representation of the point group realized at the point defines representation of the space group realized at. The matrix representing an element is a super-matrix

(2)

super-indices 1 and 2 referring to the points and respectively. The matrix representing an element is

(3)

We will not need the exact form of the non-diagonal matrix elements in Equation (3); what we need is the fact that the trace of the matrix is equal to zero. Naturally, when we consider dispersion in grapheme as it is, space group symmetry point of view adds very little in comparison to point group symmetry point of view, because the Hamiltonian, which has the symmetry, is block-diagonal.

Now we apply the group theory to analyze what happens at the points in grapheme with a perfectly commensurate superlattice potential (which appears either because of the substrate or because of the absorbed atoms), which has the same point symmetry as graphene. We consider explicitly a superlattice, known as the Kekule distortion of the honeycomb lattice [9] . In this case we may consider the Bril-louine zone (BZ) of the superlattice as the folding of the original [9] [15] . The folding leads to the identification of the corners of the original BZ (and) with the center of the new BZ. The Hamiltonian is no longer block diagonal and, because the points and are now identical, has the full symmetry. We thus observe a paradox situation: due to decrease of the translational symmetry the point symmetry of the Hamiltonian has increased.

Because of the symmetry of the Hamiltonian, we need to decompose representation realized by Matrices (2) and (3) with respect to the irreducible representations of the group.To obtain the decomposition, it is convenient to use equation

(4)

which shows how many times a given irreducible representation is contained in a reducible one [16] . In Equation (4) is the number of elements in the group, is the character of an operator in the irreducible representation and is the character of the operator in the representation being decomposed. Actually, even without using Equation (4), just by inspection of the two lowest line of Table 1 we obtain the decomposition

(5)

We see that due to supercell potential two degenerate Dirac points disappear. At the point we have two merging bands realizing representation and another two merging bands realizing representation. We

Table 1. Characters table for irreducible representations of point groups.

may expect that representation, as being more symmetrical, is realized at the top of the valence band, and the representation is realized at the bottom of the conduction band.

To get the form of the energy spectrum of the electrons in the vicinity of the point let us consider both the term [17] and the supercell potential as a perturbation. The effective Hamoltonian is

(6)

where

(7)

and reduces to two independent real constants

(8)

The specific form of the operator follows from the symmetry of the base functions realizing representations and. By shifting origin of the energy axis these two constants can be chosen as

(9)

Forming and solving the secular equation from these matrix elements, we obtain

(10)

the sign plus corresponding to an upper pair of bands, and the sign minus corresponding to a lower pair of bands.

To resolve between the branches in each pair we should take into account corrections to the operator. The first order in corrections is equal to zero because the symmetry group contains the center of inversion. To the second order in we have

(11)

where is an Hermitian tensor operator (symmetrical in the suffixes and). These include the corrections from the terms linear in in the Hamiltonian in the second-order perturbation theory and the corections from the terms quadratic in in the first-order perturbation theory [17] . Notice that is small relative to both (because we consider the states in the vicinity of the point), and with respect to.

The relations exist between the matrix elements of the operator because of the requirements of symmetry. As regards their transformation law under the symmetry operations, the wave functions which form the basis of the representation (1, 2 branches) can be taken in the form

(12)

and the wave functions which form the basis of the representation can be taken in the form

(13)

From this, we easily conclude that in the first case the matrix elements of the reduce to three independent real constants

(14)

The matrix elements of the operator are

(15)

In the second case the matrix elements of the also reduce to three independent real constants

(16)

The matrix elements of the operator are

(17)

Forming and solving the secular equation from these matrix elements, we obtain the branches of the spectrum

(18)

The sign of is determined by the fact whether we are dealing with the lower or upper bands. The formula for the branches of the spectrumcan is obtained similarly.

The folding of the BZ, together with the destruction of previously existing gapless Dirac points, leads to appearance of the new ones. In fact, the new BZ is still a hexagon, and the same symmetry arguments used for graphene can be used to explain appearance of the gapless Dirac points at the corners of the new BZ. However, these new Dirac points are situated deep below or high above the Fermi level, hence, less of them manifest themselves than Dirac points of unreconstructed graphene.

References

- Uchoa, B., Lin, C.-Y. and Castro Neto, A.H. (2008) Tailoring Graphene with Metals on Top. Physical Review B, 77, Article ID: 035420.
- Giovannetti, G., Khomyakov, P.A., Brocks, G.V., Karpan, M., van den Brink, J. and Kelly, P.J. (2008) Doping Graphene with Metal Contacts. Physical Review Letters, 101, Article ID: 026803.
- Elias, D.C., Nair, R.R., Mohiuddin, T.M.G., Morozov, S.V., Blake, P., Halsall, M.P., Ferrari, A.C., Boukhvalov, D.W., Katsnelson, M.I., Geim, A.K. and Novoselov, K.S. (2009) Control of Graphene’s Properties by Reversible Hydrogenation: Evidence for Graphene. Science, 323, 610-613. http://dx.doi.org/10.1126/science.1167130
- Chen, J.-H., Jang, C., Adam, S., Fuhrer, M.S., Williams, E.D. and Ishigami, M. (2008) Charged-Impurity Scattering in Graphene. Nature Physics, 4, 377-381. http://dx.doi.org/10.1038/nphys935
- Zhou, S.Y., Siegel, D.A., Fedorov, A.V. and Lanzara, A., (2008) Metal to Insulator Transition in Epitaxial Graphene Induced by Molecular Doping. Physical Review Letters, 101, Article ID: 086402. http://dx.doi.org/10.1103/PhysRevLett.101.086402
- Schedin, F., Geim, A.K., Morozov, S.V., Hill, E.W., Blake, P., Katsnelson, M.I. and Novoselov, K.S. (2013) Detection of Individual Gas Molecules Adsorbed on Graphene. Nature Materials, 6, 652-655. http://dx.doi.org/10.1038/nmat1967
- Lamoen, D. and Persson, B.N.J. (1998) Adsorption of Potassium and Oxygen on Graphite: A Theoretical Study. Journal of Chemical Physics, 108, 3332-3341. http://dx.doi.org/10.1063/1.475732
- Vasseur, G., Fagot-Revurat, Y., Kierren, B., Sicot, M. and Malterre, D. (2013) Effect of Symmetry Breaking on Electronic Band Structure: Gap Opening at the High Symmetry Points. Symmetry, 5, 344-354. http://dx.doi.org/10.3390/sym5040344
- Farjam, M. and Rafii-Tabar, H. (2009) Energy Gap Opening in Submonolayer Lithium on Graphene: Local Density Functional and Tight-Binding Calculations. Physical Review B, 79, Article ID: 045417.
- Kogan, E. and Nazarov, V.U. (2012) Symmetry Classification of Energy Bands in Graphene. Physical Review B, 85, Article ID: 115418. http://dx.doi.org/10.1103/PhysRevB.85.115418
- Kogan, E. (2013) Symmetry Classification of Energy Bands in Graphene and Silicene. Graphene, 2, 74-80. http://dx.doi.org/10.4236/graphene.2013.22011
- Kogan, E., Nazarov, V.U., Silkin, V.M. and Kaveh, M. (2013) Energy Bands in Graphene: How good Is the Tight- Binding Model? Physical Review B (Submitted for Publication).
- Bradley, C.J. and Cracknell, A.P. (1972) The Mathematical Theory of Symmetry in Solids. Clarendon Press, Oxford.
- Aroyo, M.I., Capillas, C., De la Flor, G., Kirov, A.K., Orobengoa, D., Perez-Mato, J. and Wondraschek, H. (2010) Representations of Cristallographic Groups. http://www.crystallography.fr/mathcryst/pdf/nancy2010/Aroyo_reps2010.pdf
- Cheianov, V.V., Falko, V.I., Syljyuasen, O. and Altshuler, B.L. (2009) Hidden Kekulé Ordering of Adatoms on Graphene. Solid State Communications, 149, 1499-1501. http://dx.doi.org/10.1016/j.ssc.2009.07.008
- Landau, L.D. and Lifshitz, E.M. (1991) Landau and Lissitz: Course of Theoretical Physics on Quantum Mechanics. Pergamon Press, Oxford.
- Lifshitz, E.M. and Pitaevskii, L.P. (1979) Landau and Lissitz: Course of Theoretical Physics on Statistical Physics. Pergamon Press, Oxford.