# Landau levels in deformed bilayer graphene at low magnetic fields

###### Abstract

We review the effect of uniaxial strain on the low-energy electronic dispersion and Landau level structure of bilayer graphene. Based on the tight-binding approach, we derive a strain-induced term in the low-energy Hamiltonian and show how strain affects the low-energy electronic band structure. Depending on the magnitude and direction of applied strain, we identify three regimes of qualitatively different electronic dispersions. We also show that in a weak magnetic field, sufficient strain results in the filling factor being the most stable in the quantum Hall effect measurement, instead of in unperturbed bilayer at a weak magnetic field. To mention, in one of the strain regimes, the activation gap at is, down to very low fields, weakly dependent on the strength of the magnetic field.

## I Introduction

Electrons in monolayer graphene are chiral quasiparticles with Berry phase , Dirac-like linear dispersion and Landau quantization in magnetic field leading to a peculiar spectrum with a Landau level at ”zero energy” (the Fermi point separating conduction and valence bands in a neutral structure) geim_science_2009 (); novoselov_nature_2005 (); zhang_nature_2005 (); abergel_advphys_2010 (). Electrons in bilayer graphene also exhibit an exotic behaviour. They can be viewed as ‘massive chiral fermions’, that is, quasiparticles with a parabolic spectrum but also Berry phase of geim_science_2009 (); abergel_advphys_2010 (); mccann_prl_2006 (); novoselov_natphys_2006 (), and their high-magnetic-field spectrum features an eight-fold degenerate Landau level at zero energy. It has been noticed that the topology of the constant-energy lines in the dispersion bands in bilayer graphene may change, from single-connected, almost circular at higher energies, meV, into several disconnected at lower energies meV: a transformation known in the physics of metals as the Lifshitz transition lifshitz_jetp_1960 (), which occurs when the Fermi level passes the saddle point in the spectrum. Theoretically, details of this transformation in the spectrum and the energy at which topology of the energy band changes are determined by the interplay between the nearest-neighbour and skew interlayer hoppings of electrons. For energies , the single-particle electron spectrum consists of four Dirac cones with a linear dispersion, three with a Berry phase and one with . In this article, we show that the topology of the low-energy single-particle dispersion of bilayer graphene critically depends on mechanical deformations of the crystal. Strain not only influences the energy at which the topology of the dispersion changes but more critically, determines the number of Dirac cones in the low-energy part of BLG spectrum, from four in the unperturbed bilayer, down to two (both with the Berry phase equal to ) in a strongly strained crystal. To mention, no such dramatic changes occur in strained monolayer graphene where deformations of the crystalline lattice induce merely a shift in the position of the Dirac point in the momentum plane manes_prb_2007 (); pereira_prb_2009 (). We also track the strain-induced spectral changes in the bilayer down to the evolution of the Landau levels for electrons, which enables us to predict what features in the quantum Hall effect in bilayer graphene crystals would persist down to the lowest magnetic fields and lowest carrier densities.

## Ii Tight-binding model for strained bilayer graphene

Bilayer graphene consists of two honeycomb layers of carbon atoms arranged according to Bernal stacking bernal_porsa_1924 () (see Fig. 1(a)) with four inequivalent carbon sites in the unit cell, , in the bottom layer and , in the top layer, and with positioned over . Hexagonal symmetry in the plane leads to three nearest neighbours at the same distance surrounding each site, their position described with the help of vectors , and , where is the carbon-carbon distance.

For the unperturbed bilayer, the electronic spectrum is gapless, with the conduction and valence bands touching in the corners of the hexagonal Brillouin zone, denoted as the K points at the energy conventionally used as the reference level, which also coincides with the Fermi level in undoped bilayer graphene. Two additional bands are split from the neutrality point by energy (see Fig. 1(b)). At very low energies, meV, the isoenergetic line undergoes a topological transition as it splits from a single-connected line into four disconnected parts, each resulting from a separate Dirac cone, as shown in the top of Fig. 1(b) .

In the suspended graphene devices du_natnano_2008 (); bolotin_ssc_2008 (); feldman_natphys_2009 (); weitz_science_2010 (); martin_prl_2010 () deformations can develop in the atomically thin carbon membranes due to the motion of contacts as the sample is cooled down. Understanding the influence of deformations of bilayer graphene in particular on its electronic properties is, now, becoming a pressing matter, in order to distinguish strain effects in suspended devices from the spectral changes expected in the recently predicted low-temperature phase transitions induced by the electron-electron interactions in undoped bilayer castro_prl_2008 (); zhang_prb_2010 (); vafek_prb_2010 (); lemonik_prb_2010 (); nandkishore_prl_2010 (); kotov_arxiv_2010 () which are being searched for in suspended bilayer graphene devices feldman_natphys_2009 (); weitz_science_2010 (); martin_prl_2010 (). Below, we characterise strain using angle between the principal axis of strain tensor and coordinates in Fig. 1, and its eigenvalues, and . Recently investigated interlayer shear deformations son_arxiv_2010 () can be included by allowing a shift of the top layer with respect to the bottom one, by . As in monolayers, stretching of the lattice changes its symmetry and makes the intralayer hops - (-) dependent on the direction of the hop, from a single value, in the Slonczewski-Weiss-McClure parametrisation slonczewski_physrev_1958 (); dresselhaus_advphys_1981 (), to three different values for the hops along bonds , and , , where , or , and takes into account unilateral deformation while the second term, , accounts for the ‘hydrostatic’ rescaling of the lattice period and leads only to a change in the value of Dirac velocity, . Here, quantifies the change of the intralayer - hopping upon the change of the distance of carbon atoms on the plane. To describe the effect of strain in bilayers, we need to account for additional changes in the interlayer hoppings. The direct - interlayer coupling, , may be changed by shear, , which has no bearing on the topology of electron bands at low energies. More importantly, strain and shear make the value of the skew interlayer coupling dependent on the direction of the hop, leading, again, to three distinct couplings , where , for the three hops indicated in Fig. 1 with black dashed lines.

The variation of the hopping elements in the tight-binding model by strain and the lower symmetry of the bilayer lattice result in modification of the dispersion relation for electrons at low energies. The latter is most convenient to capture using an effective Hamiltonian, which describes in details the pair of low-energy bands near the Dirac point in bilayer graphene corresponding to electron states located predominantly on the sublattices and (those coupled by the skew hop , but not by the direct hop ). The calculation leading to such effective Hamiltonian includes two steps. One consists in the Schrieffer-Wolff transformation schrieffer_physrev_1966 () projecting the four-band model onto the effective two-band model mccann_prl_2006 (). The other takes into account that cancellation between the contributions of the and hops in directions , and for the coupling of plane wave Bloch states on and sublattices in the corners and of the Brillouin zone is no more exact. Note that lattice deformations also distort the hexagonal shape of the Brillouin zone, what is automatically taken into account in our analysis. As a result, the hopping integrals between sublattices are modified from to for the intralayer - hops (where and ), and from to for the interlayer skew hops (where and with denoting the anticlockwise angle between negative direction of the axis and vector ). Then, the resulting two-band model can be described using the Hamiltonian

where and . Following the suggestion manes_prb_2007 () that in monolayers the effect of homogeneous strain is equivalent to a constant vector potential, which can be eliminated by a gauge transformation equivalent to a shift of the valley centre from the Brillouin zone corners [ () or ()], we introduce here valley momentum , and arrive at the effective Hamiltonian for strained bilayer graphene:

(1) |

where

## Iii The influence of the electron-electron interaction on the strain-induced term in the low-energy Hamiltonian

In this section, we investigate the influence of electron-electron interactions on the low-energy Hamiltonian of (strained) bilayer graphene using the renormalisation group (RG) approach. The derivation of the RG equations is based on the analysis of the leading logarithmic divergences in the perturbation theory and their subsequent resummation. It follows the line of Ref. lemonik_prb_2010 (), though in this paper we will not consider the renormalizations of the marginal short range interaction giving rise to the instability towards nematic state lemonik_prb_2010 (). As usual, the bookkeeping of the perturbation theory is greatly simplified by using the diagrammatic technique.

We pictorially represent the Hamiltonian (1) as

(2) |

In the energy region where the RG flow occurs, two last terms in Eq. (2) can be considered as the perturbation so the Green function has the form

(3) |

The Coulomb interaction denoted by

(4) |

is divergent at so all the terms in approximation should be summed which gives

(5) |

The explicit calculation of the polarization operator yields

(6) |

where is the total degeneracy of the electronic spectrum. All the calculation is controlled in approximation. Using Eqs. (5) and (6), we obtain

(7) |

The last equation is written under assumption that ( being the screening radius), which is always true in the region of applicability of the Hamiltonian (1) of the main text.

To derive the RG equations we, first, calculate the correction to the self-energy. On each step of the RG the integration over along the interaction line is restricted by . Its energy derivative gives the correction to the quasiparticle weight

(8) |

which can be always eliminated by the rescaling of the Fermionic operators and thus does not affect by itself any gauge invariant quantity. We will perform the rescaling in such a way that the quasiparticle weight is set to unity, so that all the corrections have the meaning of the observable spectrum (note in passing that all the corrections to the observables are logarithmic whereas is divergent more severely). As a bonus, the scalar vertex is also not renormalized due to the gauge invariance

The correction to the effective mass is found from

(9) |

The other corrections are

(10) |

Calculating the integrals in Eqs. (9)-(10), substituting , and , we obtain:

(11) |

where and is the running energy scale. The electron-electron interaction enhances the strain-induced term stronger than other parameters, and at energies , where the influence of strain plays a dominant role in determining the electron spectrum, we substitute , in the Hamiltonian (1).

## Iv Low-energy electronic band structure of strained bilayer graphene

The resulting low-energy electronic dispersions of the conduction band [the band structure described by Hamiltonian in equation (1) is always electron-hole symmetric] of strained bilayer graphene in the vicinity of the valley are displayed for representative values of in Fig. 2. The dispersion for the case of the unperturbed bilayer, that is for , is shown in Fig. 2(a). It is quasi-parabolic at energies , with responsible for the trigonal warping of isoenergetic lines. At , the isoenergetic line undergoes a Lifshitz transition and splits into four Dirac cones mccann_prl_2006 (). Remaining dispersions in the top row depict the bottom of the conduction band for , . As the value of is increased (Fig. 2(b)-(c)), the two side cones positioned off the axis and the central cone move closer, as shown in the graph for . Those three cones collide for , and for only two Dirac cones remain, as shown in Fig. 2(c) for . If is negative (graphs (d) and (e) in Fig. 2), the central cone and the side cone positioned on the axis approach each other and collide for , creating a local, quasi-parabolic minimum. This minimum persists for some range of the strain, although it lifts off the plane, as shown in the graph 2(d) for . Eventually, for , the minimum merges with a saddle point and with further decrease of , again only two Dirac cones remain in the spectrum, as in Fig. 2(e) for . A contrasting situation of , , is presented in the graphs in Fig. 2(f)-(h) (due to symmetry, dispersions for opposite values of are mirror reflections of each other with respect to the plane and hence only situation of is described here). Again, starting from the unperturbed system and increasing the magnitude of leads to two of the cones moving closer to each other (Fig. 2(f) for ) and colliding with a creation of a local, quasi-parabolic minimum which lifts off the plane (graph (g) in Fig. 2 for ) and eventually disappears, leaving only two Dirac cones (graph (h) in Fig. 2 for ). We identify three qualitatively different regimes in the space: (i) spectrum contains four Dirac cones, (ii) spectrum contains two Dirac cones and a local minimum, created by a collision of two of the cones and which in general does not touch the plane, and (iii) spectrum contains only two Dirac cones. The extent of those regimes is shown in Fig. 3.

Strain-induced deformation of the low-energy electronic dispersion gives rise to distinctive features in the low-energy density of states in each of the strain regimes. In particular, existence of a local parabolic minimum results in a sharp step in the density of states, while saddle points lead to van Hove singularities. Note that for , the symmetry between the side cones is broken and the saddle points are no longer all found at the same energy , leading to more than one van Hove singularity and more than one Lifshitz transition. Densities of states representative of all the regimes are shown in Fig. 4 for , , and , that is points that lie on the line. Indeed, for only one van Hove singularity is present, whereas three can be traced for . A sharp step, corresponding to the contribution of a quasi-parabolic part of the band structure to the density of states, is visible for . It is followed by a van Hove singularity, as for this dispersion the saddle point lies not far above the position of the local minimum. Finally, for , the sharp step is washed out (as the local minimum merged with the saddle point) and only one van Hove singularity remains.

## V Landau levels of strained bilayer graphene

Presence of strain in bilayer graphene may also lead to new features in the low-energy Landau level (LL) spectrum. To obtain Landau level spectra for bilayer graphene in a perpendicular external magnetic field, we use the two-band Hamiltonian (1), where the electron momentum, , incorporates the electromagnetic vector potential . For simplicity, we use the Landau gauge, . The operators and fulfil the role of raising and lowering operators when acting on the magnetic oscillator functions , where the magnetic length . If we neglect the second and third terms in the effective Hamiltonian (1), we can describe the low-energy Landau levels as a sequence of levels mccann_prl_2006 (); abergel_prb_2007 (); ando_jpsj_2007 (); koshino_prb_2008 (); nakamura_prb_2008 (); mucha-kruczynski_jpcm_2009 ()

(12) |

with energies , where , and distinguishes between the conduction () and the valence () band. To incorporate the and the strain-induced terms into this picture, we follow the approach employed before to describe the effect of the coupling on the LL structure of graphite nakao_jpsj_1976 (). We describe the new eigenstates, each of them a linear combination of an infinite number of functions , with vectors which th and th component create a minimal subspace required to describe the th LL in the absence of . In the corresponding infinite Hamiltonian matrix,

(13) |

where

(14) | ||||

(15) | ||||

(16) |

the first term in the low-energy Hamiltonian results in decoupled diagonal blocks, the second term results in some off-diagonal couplings between functions and , while the strain-induced third term leads to off-diagonal couplings between functions and footnote (). We then truncate the infinite matrix, restricting the calculation to a given LLs and diagonalise the resulting Hamiltonian numerically. The number of the basis vectors required in the calculation in order to describe correctly the low-energy LL structure increases with decreasing magnetic field (reflecting growing importance of the terms at weaker fields). In this work, matrices of the dimension 1200 have been used to produce the spectra, which are shown in the top row of Fig. 5 for , , and .

The additional degeneracy of the th LL in bilayer graphene at weak fields is connected to the number of Dirac cones in the electronic spectrum. At magnetic fields such that the inverse of the magnetic length is smaller than the distance in the reciprocal space between a pair of cones, each of them provides additional four-fold degeneracy to the LL (just like in monolayer graphene novoselov_nature_2005 (); zhang_nature_2005 ()). At stronger fields, contributions of separate cones can no longer be resolved, leading to an 8-fold degenerate th LL mccann_prl_2006 (); novoselov_natphys_2006 (). Hence, for the strain within the region shown in yellow in Fig. 3, the zeroth LL becomes 16-fold degenerate at very weak but nonzero fields as three LLs, 8-, 4- and 4-fold degenerate, merge together, reflecting the existence of four Dirac cones in the spectrum. However, for regions of greater strain, shown in Fig 3 with dark and light shading, the zeroth LL is 8-fold degenerate at all fields due to only two cones in the electronic spectrum. As the size of the gap between two LLs decreases with the LL index, the filling factor determining the biggest activation gap at low fields is closely connected to the additional degeneracy of the zeroth LL. Graphs showing the magnitude of activation gaps for filling factors , , and as a function of magnetic field for the values of representing all regions from Fig. 3, are shown in the bottom row in Fig. 5. For the region of weak strain, here portrayed with , the biggest activation gap occurs at , whereas for the two other regions (, and ), filling factor is the most stable. For the region of intermediate strain (), where a local parabolic minimum exists in the electronic dispersion, the activation energy between the th and the first LL experiences a very unusual intermediate saturation, indicating that one of the LLs gets stuck in that minimum.

## Vi Acknowledgements

This project has been funded by the EPSRC grants EP/G041954 and Science & Innovation Award EP/G035954. This review is based upon the lecture course ”Electronic properties of graphene” taught by V. Fal’ko at the 5th Windsor Summer School, , in Windsor, United Kingdom windsor ().

## References

- (1) A.K. Geim, Science 324, 1530 (2009).
- (2) K.S. Novoselov et al, Nature 438, 197 (2005).
- (3) Y. Zhang, Y.-W. Tan, H.L. Stormer, and P. Kim, Nature 438, 201 (2005).
- (4) D.S.L. Abergel, V. Apalkov, J. Berashevich, K. Ziegler, and T. Chakraborty, Adv. Phys. 59, 261 (2010).
- (5) E. McCann and V.I. Fal’ko, Phys. Rev. Lett. 96, 086805 (2006).
- (6) K.S. Novosleov et al, Nature Phys. 2, 177 (2006).
- (7) I.M. Lifshitz, Sov. Phys. JETP 11, 1130 (1960).
- (8) J.L. Manes, Phys. Rev. B 76, 045430 (2007).
- (9) V.M. Pereira, A.H. Castro Neto, and N.M.R. Peres, Phys. Rev. B 80, 045401 (2009).
- (10) J.D. Bernal, Proc. R. Soc. A 106, 749 (1924).
- (11) X. Du, I. Skachko, A. Barker, and E.Y. Andrei, Nature Nanotech. 3, 491 (2008).
- (12) K.I. Bolotin et al, Sol. St. Commun. 146, 351 (2008).
- (13) B.E. Feldman, J. Martin, and A. Yacoby, Nature Phys. 5, 889 (2009).
- (14) R.T. Weitz, M.T. Allen, B.E. Feldman, and A. Yacoby, Science 330, 812 (2010).
- (15) J. Martin, B.E Feldman, R.T. Weitz, M.T. Allen, and A. Yacoby, Phys. Rev. Lett. , (2010).
- (16) E.V. Castro, N.M.R. Peres, T. Stauber, N.P. Silva, Phys. Rev. Lett. 100, 186803 (2008).
- (17) F. Zhang, H. Min, M. Polini, and A.H. MacDonald, Phys. Rev. B 81, 041402 (2010).
- (18) O. Vafek and K. Yang, Phys. Rev. B 81, 041401 (2010).
- (19) Y. Lemonik, I.L. Aleiner, C. Toke, and V.I. Fal’ko Phys. Rev. B 82, 201408 (2010).
- (20) R. Nandkishore and L. Levitov, Phys. Rev. Lett. 104, 156803 (2010).
- (21) V.N. Kotov, B. Uchoa, V.M. Pereira, A.H. Castro Neto, and F. Guinea, arXiv.:10123484.
- (22) Y.-W. Son, S.-M. Choi, Y. P. Hong, S. Woo, S.-H. Jhi, arXiv:1012.0643 (2010).
- (23) J.C. Slonczewski and P.R. Weiss, Phys. Rev. 109, 272 (1958).
- (24) M.S. Dresselhaus and G. Dresselhaus, Adv. Phys. 30, 139 (1981).
- (25) J.R. Schrieffer and P.A. Wolff, Phys. Rev. 149, 491 (1966).
- (26) D.S.L. Abergel and V.I. Fal’ko, Phys. Rev. B 75, 155430 (2007).
- (27) T. Ando, J. Phys. Soc. Jpn. 76, 104711 (2007).
- (28) M. Koshino and T. Ando, Phys. Rev. B 77, 115313 (2008).
- (29) M. Nakamura, L. Hirasawa, and K.-I. Imura, Phys. Rev. B 78, 033403 (2008).
- (30) M. Mucha-Kruczynski, D.S.L. Abergel, E. McCann, and V.I. Fal’ko, J. Phys. Condens. Mat. 21, 344206 (2009).
- (31) K. Nakao, J. Phys. Soc. Jap. 40, 761 (1976).
- (32) Note, that two rows and columns in the above matrix, identically equal to zero, give rise to two solutions at zero energy, which correspond to unphysical eigenstates (zero vectors) and should not be confused with the zero energy Landau levels.
- (33) V.I. Fal’ko, Lecture notes on ”Electronic properties of graphene”, 5th Windsor Summer School, , Windsor, United Kingdom.