Zero modes, energy gap, and edge states of anisotropic honeycomb lattice in a magnetic field
We present systematic study of zero modes and gaps by introducing effects of anisotropy of hopping integrals for a tight-binding model on the honeycomb lattice in a magnetic field. The condition for the existence of zero modes is analytically derived. From the condition, it is found that a tiny anisotropy for graphene is sufficient to open a gap around zero energy in a magnetic field. This gap behaves as a non-perturbative and exponential form as a function of the magnetic field. The non-analytic behavior with respect to the magnetic field can be understood as tunneling effects between energy levels around two Dirac zero modes appearing in the honeycomb lattice, and an explicit form of the gap around zero energy is obtained by the WKB method near the merging point of these Dirac zero modes. Effects of the anisotropy for the honeycomb lattices with boundaries are also studied. The condition for the existence of zero energy edge states in a magnetic field is analytically derived. On the basis of the condition, it is recognized that anisotropy of the hopping integrals induces abrupt changes of the number of zero energy edge states, which depend on the shapes of the edges sensitively.
pacs:71.70.Di, 73.43.-f, 81.05.Uw
Recent experiments on graphene(1); (2); (3); (4); (5) have led to renewed interest in physical properties of electrons on the honeycomb lattice. Despite its simple structure, the honeycomb lattice provides non-trivial physical phenomena which can not be observed in the ordinary square lattice. Among them, much attention has been paid to its peculiar dispersion. In the absence of a magnetic field, the honeycomb lattice has zero modes at the corners K and K’ of the Brillouin zone. By treating these zero modes as Dirac fermions, the unconventional quantization of the Hall conductance observed for graphene was explained(7); (6); (8), although the full proper theoretical treatment of the Hall conductance on the honeycomb lattice was made very recently(9). Moreover, when the system has a boundary, there are edge modes localized on the boundary. The existence of the edge modes depends on a choice of the boundary, and for zigzag and bearded edges there occur large density of states localized on these edges at the Fermi energy due to their flat band structures(10); (11); (12); (13); (14); (15); (16); (17); (18); (19); (20). In addition, it is suggested that the edge modes induce charge accumulation on these edges(21); (22).
In this paper, we study properties of these zero modes in the presence of anisotropy of the hopping integrals in the honeycomb lattice. Recently, the anisotropy of the hopping integrals was introduced by replacing one of the hopping integrals with a general value (23); (24); (25); (26); (9); (28); (27) in order to investigate the unconventional quantum Hall effects on graphene (see Fig.1). (For , we have the anisotropic honeycomb lattice.) In Ref.(9), by using topological arguments, an algebraic expression of the quantum Hall conductance was obtained for almost all gaps including subband gaps, and it was shown that the unconventional quantization of the Hall conductance in a weak magnetic field is realized for weak (), while only the conventional quantization is obtained for strong (). Furthermore, for the graphene case (), the unconventional quantization was found to persist up to the Van Hove singularity(9); (19).
The anisotropy of the hopping parameters is also known to change the peculiar dispersion mentioned above(24); (25); (26); (27). However, in the absence of a magnetic field, its influence is restrictive: Although a gap opens around zero energy for , there remain two zero modes in the Brillouin zone for (24). Therefore, a large anisotropy is needed to change the zero mode structure. As well as the zero mode structure, that of zero energy edge states was shown to change by a large anisotropy(28).
In this paper, it will be shown the situation is drastically changed in the presence of a magnetic field. We analytically derive the condition for the existence of zero modes in a magnetic flux ( and are mutually prime integers), and from the condition it is found that, in the limit of , zero modes exist only for , but a gap around zero energy opens for . In other words, a small anisotropy () is sufficient to open a gap in the presence of a weak magnetic field.
For , the gap around zero energy in a weak magnetic field behaves as a non-perturbative and exponential form as a function of . It will be shown that this behavior is naturally explained in terms of the spontaneous breaking of supersymmetry (29); (30). In particular, an explicit form of the gap around zero energy for is obtained by the WKB method. At , the gap around zero energy in a weak magnetic field is found to make a transition from an exponential (non-perturbative) to a power-law (perturbative) behavior as a function of , and for , the energy bands in a weak magnetic field show linear dependence on .
We will also show that the structure of edge states in the presence of a magnetic field is different from that in the absence of a magnetic field. The condition for the existence of zero energy edge states in a magnetic field is analytically derived, and it is found that the anisotropy of the hopping integrals induces abrupt changes of the number of zero energy edge states, which also sensitively depend on shapes of the edges.
The organization of this paper is as follows. In Sec.II, we present our model. The condition for the existence of zero modes in a magnetic field is analytically derived in Sec.III, both from the secular equation and from the normalizability condition of states with zero energy. On the basis of the condition for the existence of zero modes, the energy spectrum near zero energy in a weak magnetic field is systematically examined in Sec.IV. In Sec.V, zero energy edge states are analyzed, where crucial roles of the anisotropy of the hopping integrals are recognized again. Finally, we summarize our results and discuss possible experimental realization of anisotropy of the hopping integrals in Sec.VI.
Ii Tight-binding model on the honeycomb lattice in a magnetic field
Let us consider the tight-binding model on the honeycomb lattice with nearest-neighbor hopping in a magnetic field as shown in Fig.1. By denoting wave functions on two sublattices of the honeycomb lattice as and , respectively, the tight-binding model is given by
where a magnetic flux through the unit hexagon is given by . Here we have introduced anisotropy of the hopping integrals: The hopping integrals of the horizontal bonds are , and those for the other bonds are 1. For simplicity, we neglect the spin degrees of freedom in the following.
Iii The condition for the existence of zero modes
For the isotropic case (), it was found that zero modes exist for all (rational) values of (31). We now derive the condition for the existence of zero modes in the anisotropic case.
From the Bloch’s theorem, the wave functions are written as
where is given by
The eigenenergies are given by
From (6) with , we find two zero modes at
for . By expanding and around and in (7),
is given by
For , the two Dirac zero modes merge into a confluent point
and for , we have a gap around zero energy.
iii.1 Derivation from a secular equation
Now we consider . We suppose that is a rational number, ( and are mutually prime integers). Since Eq.(1) has translational symmetry in the -direction, the wave functions are written as
and (1) becomes
By the gauge transformation and , (13) is rewritten as
where with . Since the spectrum is found to be invariant under the transformation , we can restrict the range of to without loss of generality. Moreover, in (14), (,) and (,) obey the same equation, thus from the Bloch’s theorem, we have
where satisfies . Therefore (14) reduces to the eigenequation of a matrix. In the secular equation of this, all non-constant terms containing less than factors of should cancel out each other since the eigenvalue has periodicity with respect to . From this property, it is found that the secular equation is written as the following form:
where is a th-order polynomial of with , and it is independent of . The secular determinant for (14) with determines as
Here we have used
which is derived from .
the range of is and there exist two independent ’s with . From the secular equation (16), we have two modes at these ’s. On the other hand, if satisfies
we have and there is no with . We have a gap around zero energy in this case.
iii.2 Derivation from the normalizability condition of states
In Sec. III.1, we derived the condition for the existence of zero modes from the secular equation. Here, we re-derive it from the normalizability condition of states with zero energy.
Let us first consider (14) with ,
where . Taking the absolute values of the both sides in (22), we obtain
For , (23) gives
From (24), we see that diverges for , and diverges for . Thus these states are not normalizable, and no relevant zero modes exist. We have a gap around in this case. On the other hand, for , (23) gives
at . Thus there exist two zero modes for . These results coincide with those of Sec. III.1.
Iv Spectrum near zero energy in a weak magnetic field
In this section, we examine the spectrum near zero energy in a weak magnetic field. Although some numerical study was presented in Ref.(23), we perform detailed analytical study here. On the basis of the condition for the existence of zero modes obtained in the previous section, we consider the following four cases separately:
, where the condition (19) is always satisfied and we have zero modes for all rational values of .
, where zero modes disappear and a gap around opens for .
, where one zero mode exists for .
, where no zero modes exist.
We show two examples of the energy bands as a function of in Fig.2. For , we have zero modes for all rational values of . As shown in the following, the energy bands in a weak magnetic field are well described by the continuum approximation.
where is given by
Here, and are those near and , respectively. Since and commute each other, we can replace with a c-number . Then putting , we obtain
Therefore, we have
around , and
around , where is the -component of the Pauli matrix.
and that around is given by
In a similar manner, the energy levels for around and are given by
For , a gap around opens for . This implies that in a weak magnetic field (), a gap around opens by a tiny distortion of graphene, . Since we do not have states, the expressions (34) and (37) need to be modified. We show two examples of energy bands as a function of for in Fig.5.
Let us focus on the gap around . In order to see how it behaves, we plot the natural logarithm of the gap around as a function of in Fig.6. From this, we find that it behaves as
The values of are obtained from Fig.6 as and for and , respectively.
The non-analytic behavior (38) can be understood as breaking of supersymmetry (29); (30); (32) in our model. The operator transforms to and vice versa, which is seen from (28) and (29). By identifying with generators of supersymmetry, in (30) can be considered as sypersymmetric Hamiltonians, ( is “boson”, and is “fermion”). Due to the supersymmetry, there is no perturbative (or power-law) correction with respect to for the states. However, tunneling effects break the supersymmetry spontaneously and the non-perturbative correction (38) appears as a gap around .
When two Dirac zero modes at are close to each other in the momentum space, namely, , the gap around can be estimated by the WKB method. For , the two Dirac zero modes at are located at (7) with
Then in (5) is given by
where is given by
In a similar manner as Sec.IV.1, we replace with a c-number and put . Then we obtain
where is given by
Identifying with a generator of supersymmetry, we have the supersymmetric Hamiltonian , which satisfies
By using the following variable ,
(46) can be rewritten as
Therefore, the potential terms for and are given by
respectively (see Fig.7).
and around , (48) becomes
Therefore, the energy levels for around and are given by
respectively, and those for around and are given by
Now take into account tunneling effects between the energy levels around and . The tunneling effects can be estimated by the WKB method presented in Appendix D of Ref.(32). Here we consider only the equation for since the equation for gives the same result. The solution of (52) which vanishes for is given by
For near zero energy (), we have
thus the solution of (61) for is obtained as
This implies that the gap around is given by
and the exponent in (38) is given by
At , the two Dirac zero modes at merge into the confluent point (11). As a consequence, a gap around in a weak magnetic field makes a transition from an exponential (non-perturbative) to a power-law (perturbative) behavior as a function of .
For , we obtain the exponent as , , and for the lowest, the second lowest, and the third lowest bands in , respectively. Thus for we have a power-law behavior .
For , we do not have zero modes but have a gap around . In Fig.10, we show two examples of the energy bands as a function of for .
Let us study behavior of energy bands in a weak magnetic field by the continuum approximation. For , we expand and around ,
then (6) gives
Thus for small and (), (71) is written as
which is quadratic in both and . In the continuum approximation, we put and with , and . Then we have
where we have neglected higher order corrections of . The energy bands in the weak magnetic field limit