# Resonant modes in strain-induced graphene superlattices

###### Abstract

We study tunneling across a strain-induced superlattice in graphene. In studying the effect of applied strain on the low-lying Dirac-like spectrum, both a shift of the Dirac points in reciprocal space, and a deformation of the Dirac cones is explicitly considered. The latter corresponds to an anisotropic, possibly non-uniform, Fermi velocity. Along with the modes with unit transmission usually found across a single barrier, we analytically find additional resonant modes when considering a periodic structure of several strain-induced barriers. We also study the band-like spectrum of bound states, as a function of conserved energy and transverse momentum. Such a strain-induced superlattice may thus effectively work as a mode filter for transport in graphene.

###### pacs:

81.05.ue, 72.80.Vp, 85.30.MnGraphene is a single layer of carbon atoms in the hybridization state, arranged according to a honeycomb lattice Novoselov et al. (2004, 2005). Transport properties in graphene are largely determined by its reduced dimensionality, which characterizes its remarkable electronic properties Castro Neto et al. (2009); Abergel et al. (2010). These include low-energy quasiparticles with a Dirac-like spectrum and a linearly vanishing density of states (DOS) at the Fermi level. Evidence of such an unconventional behaviour is to be found in several electronic properties, such as Klein tunneling Pereira et al. (2006); Barbier et al. (2009); Barbier et al. (2010a, b); Peres (2009), the optical conductivity Wang et al. (2008); Stauber et al. (2008); Pellegrino et al. (2010a), and the plasmon dispersion relation Hwang and Das Sarma (2007); Abedinpour et al. (2011); Pellegrino et al. (2010b, 2011a). These have been predicted to depend quite generally on applied strain Pellegrino et al. (2011b), following the earlier suggestion that suitably deformed graphene sheets could be engineered into nanodevices with the desired electron properties Pereira and Castro Neto (2009). For instance, it has been recently demonstrated that the electrical properties of epitaxial graphene on SiC strongly depend on the local strain induced in graphene by the substrate Low et al. (2012). One thus expects that a suitable pattern of periodically repeating stripes, with alternating values of strain, i.e. a strain-induced superlattice, may produce coherent effects on single particle transport, depending on the energy and momentum of the incident electrons. Here, we therefore study the possible occurrence of resonant states within a strain-induced superlattice in graphene.

We consider quasiparticle transmission across identical barriers, each of width , the inter-barrier separation being also , such that (Fig. 1). Let denote the coordinate orthogonal to the barriers, forming an angle with the graphene zig-zag direction. Thus, (resp., ) will refer to a superlattice oriented along the zig-zag (resp., armchair) direction. Such a superlattice is usually obtained via a step-wise varying gate potential , with within each barrier [, ], and between two neighboring barriers [, ]. Here, we will additionally consider a nonuniform profile of uniaxial strain applied along the direction, with strain modulus alternating between the values inside and outside a barrier, as above. Such a dependence approximates a smooth periodic strain wave with period , as a train of sharp steps.

Within each barrier, strain is described by the tensor , where is the strain modulus, is Poisson’s ratio for graphene Farjam and Rafii-Tabar (2009), and , with () denoting the Pauli matrices. Within a quite general tight-binding approach Castro Neto et al. (2009), strain then enters the electronic properties through the dependence of the hopping parameters on the lattice vectors Pereira et al. (2009). Expanding such a tight-binding Hamiltonian to linear order in the strain modulus, one finds that the low-lying spectrum can still be described by a Dirac-like Hamiltonian, but now (i) applied strain shifts the location of the Dirac points in reciprocal space with respect to at the vertices of the first Brillouin zone, and (ii) it induces a deformation of the Dirac cones, which can be accounted in terms of an anisotropic Fermi velocity . Specifically, one finds for the Hamiltonian under applied strain

(1) |

where measures the wave vector displacement from the shifted Dirac points , () take into account of the strain-induced deformation of the Fermi velocity, with , , is the unitary matrix performing a rotation mapping the zig zag direction onto the the direction of applied strain, is related to the logarithmic derivative of the nearest-neighbor hopping parameter with respect to the lattice parameter at zero strain, and (cf. Ref. Pellegrino et al., 2011c).

Since the strain superlattice is uniform along the coordinate orthogonal to the direction of applied strain, say , stationary eigenmodes will be characterized by constant energy and transverse wave vector . The stationary Dirac equation associated to Eq. (1), with appropriate matching conditions for the quasiparticle spinor due to the continuity of its associated current density at the barriers’ edges, can then be equivalently recast using the transfer matrix formalism Titov (2007); Bruus and Flensberg (2004). Following Ref. Pellegrino et al., 2011b, for the transfer matrix across the first, say, barrier in Fig. 1, one finds , where , and is a unimodular matrix, . Specifically, one obtains

(2a) | |||||

(2b) | |||||

(2c) | |||||

where is always real, whereas can be real or purely imaginary, depending and . More compactly, one also finds

(3) |

In Eqs. (2) and (3), we have employed the definitions , , and . Making use of the Chebyshev identity for the th power of a unimodular matrix Yeh et al. (1977), for the evolution matrix across identical barriers, one finds Pellegrino et al. (2011c)

(4) |

where . Finally, the transmission can be related to the evolution matrix as

(5) |

We are now in the position to discern whether an electronic mode is characterized by an oscillating or evanescent behavior far from the barrier superlattice. To this aim, we preliminarly observe that, depending on and , one has a propagating (resp., evanescent) wave for (resp., ), where the subscript refers to the region between two consecutive barriers [, ], and the subscript refers to the region within a barrier [, ] (Fig. 1).

Fig. 2 shows the single electron transmission across a single barrier, Eq. (5) with , as a function of the transverse wave vector and energy , each scaled by the characteristic energy

(6) |

where . Here and in the following, strain is applied along the armchair direction, , and we set , , and . In Fig. 2, cyan dashed lines delimit the two (deformed) Dirac cones defined by (left cone) and (right cone), corresponding to regions I+III and II in Fig. 1), respectively. One finds that is defined within the left cone and is exponentially vanishing within the intersection between both cones. This corresponds to having propagating modes in all the three regions. In this case, resonant modes, i.e. propagating modes with unit transmission, are characterized by the condition for stationary waves

(7) |

where , and is an integer.

Fig. 3 shows the single electron transmission across a superlattice composed of five identical barriers, Eq. (5) with . Again, nonzero values of the transmission are to be found within the intersection of the Dirac cones corresponding to the region inside a barrier and between two consecutive barriers. However, at variance of the case , because of multiple scatterings, a nonzero transmission is also possible beyond the cone . This corresponds to having evanescent modes within the barriers. Such a phenomenon is analogous to what happens to photons propagating across a 1D photonic crystal with alternative layers of a left-handed and a right-handed material (1D metamaterial) Wu et al. (2003). As for resonant modes, , besides the ones given by Eq. (7) regardless of , additional resonant modes are given by the condition

(8) |

where is defined by Eq. (2b), and . The latter condition implies that these resonant modes have globally propagating behavior. Comparing Figs. 2 and 3, one finds that, in the domain within both Dirac cones, in addition to the resonant modes given by Eq. (7), in the case there exist new resonant modes given by Eq. (8). It should also be noted that in the domain within the left cone but outside the second the resonant modes, which are only given by Eq. (8), are characterized by quite narrow lines in the transmission plots.

Outside the left Dirac cone, it is still possible to find bound states Pereira and Castro Neto (2009); Ramezani Masir et al. (2009); Barbier et al. (2010c). Within the transfer matrix method, these are given by the condition Bliokh et al. (2010) . For one finds evanescent modes outside the barriers, and therefore also far from the superlattice structure. In the case , one finds several such confined modes within the second cone (Fig. 2, solid lines outside the left cone), whereof only one such mode survives in the region outside both cones. The latter is the surface mode analyzed in Ref. Pereira and Castro Neto, 2009. In the case (Fig. 3, solid lines outside the right cone), one finds that to each bound mode in the single barrier case there correspond exactly bound states outside either cones. This is remindful of electron bands in solids, where the overlap of periodically arranged atomic orbitals give rise to a band of states.

In conclusion, we have found that a strain-induced superlattice in graphene can accomodate additional resonant quasiparticle states, analytically characterized by Eq. (8), besides the ones usually found across a single barrier, given by Eq. (7). One finds that applied strain modifies the kinetic part of the quasiparticle Hamiltonian, which preserves its Dirac-like character, but around shifted and deformed Dirac cones. This can be described in terms of a coordinate-dependent, periodic, profile of the Fermi velocity, which produces coherent effects on the quasiparticle transmission. Specifically, we find resonant modes with globally propagating behavior far from the superlattice, for conserved energy and transverse momentum within the intersection of the two deformed Dirac cones corresponding to the two alternating strained regions. Other modes are exponentially suppressed, and we also discuss the spectrum of bound states, which arrange themselves as ‘bands’, depending on the overall number of barriers making up the superlattice. We thus surmise that a strain-induced superlattice in graphene can be used as a filter for the resonant modes here discussed.

## References

- Novoselov et al. (2004) 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).
- Novoselov et al. (2005) K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature 438, 197 (2005).
- Castro Neto et al. (2009) A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 000109 (2009).
- Abergel et al. (2010) D. S. L. Abergel, V. Apalkov, J. Berashevich, K. Ziegler, and T. Chakraborty, Adv. Phys. 59, 261 (2010).
- Pereira et al. (2006) J. M. Pereira, V. Mlinar, F. M. Peeters, and P. Vasilopoulos, Phys. Rev. B 74, 045424 (2006).
- Barbier et al. (2009) M. Barbier, P. Vasilopoulos, and F. M. Peeters, Phys. Rev. B 80, 205415 (2009).
- Barbier et al. (2010a) M. Barbier, P. Vasilopoulos, and F. M. Peeters, Phys. Rev. B 81, 075438 (2010a).
- Barbier et al. (2010b) M. Barbier, P. Vasilopoulos, and F. M. Peeters, Phys. Rev. B 82, 235408 (2010b).
- Peres (2009) N. M. R. Peres, J. Phys.: Cond. Matter 21, 323201 (2009).
- Wang et al. (2008) F. Wang, Y. Zhang, C. Tian, C. Girit, A. Zettl, M. Crommie, and Y. R. Shen, Science 320, 206 (2008).
- Stauber et al. (2008) T. Stauber, N. M. R. Peres, and A. K. Geim, Phys. Rev. B 78, 085432 (2008).
- Pellegrino et al. (2010a) F. M. D. Pellegrino, G. G. N. Angilella, and R. Pucci, Phys. Rev. B 81, 035411 (2010a).
- Hwang and Das Sarma (2007) E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 (2007).
- Abedinpour et al. (2011) S. H. Abedinpour, G. Vignale, A. Principi, M. Polini, W.-K. Tse, and A. H. MacDonald, Phys. Rev. B 84, 045429 (2011), URL http://link.aps.org/doi/10.1103/PhysRevB.84.045429.
- Pellegrino et al. (2010b) F. M. D. Pellegrino, G. G. N. Angilella, and R. Pucci, Phys. Rev. B 82, 115434 (2010b).
- Pellegrino et al. (2011a) F. M. D. Pellegrino, G. G. N. Angilella, and R. Pucci, High Press. Res. 31, 98 (2011a).
- Pellegrino et al. (2011b) F. M. D. Pellegrino, G. G. N. Angilella, and R. Pucci, Phys. Rev. B 84, 195407 (2011b).
- Pereira and Castro Neto (2009) V. M. Pereira and A. H. Castro Neto, Phys. Rev. Lett. 103, 046801 (2009).
- Low et al. (2012) T. Low, V. Perebeinos, J. Tersoff, and P. Avouris, Phys. Rev. Lett. 108, 096601 (2012), URL http://link.aps.org/doi/10.1103/PhysRevLett.108.096601.
- Farjam and Rafii-Tabar (2009) M. Farjam and H. Rafii-Tabar, Phys. Rev. B 80, 167401 (2009).
- Pereira et al. (2009) V. M. Pereira, A. H. Castro Neto, and N. M. R. Peres, Phys. Rev. B 80, 045401 (2009).
- Pellegrino et al. (2011c) F. M. D. Pellegrino, G. G. N. Angilella, and R. Pucci, Phys. Rev. B 84, 195404 (2011c).
- Titov (2007) M. Titov, Eur. Phys. Lett. 79, 17004 (2007).
- Bruus and Flensberg (2004) H. Bruus and K. Flensberg, Many-Body Quantum Theory in Condensed Matter Physics: An Introduction (Oxford University Press, Oxford, 2004).
- Yeh et al. (1977) P. Yeh, A. Yariv, and C. Hong, J. Opt. Soc. Am. 67, 423 (1977).
- Wu et al. (2003) L. Wu, S. He, and L. Shen, Phys. Rev. B 67, 235103 (2003).
- Ramezani Masir et al. (2009) M. Ramezani Masir, P. Vasilopoulos, and F. M. Peeters, New J. Phys. 11, 095009 (2009).
- Barbier et al. (2010c) M. Barbier, P. Vasilopoulos, and F. M. Peeters, Phil. Trans. R. Soc. A 368, 5499 (2010c).
- Bliokh et al. (2010) Y. P. Bliokh, V. Freilikher, and F. Nori, Phys. Rev. B 81, 075410 (2010), URL http://link.aps.org/doi/10.1103/PhysRevB.81.075410.