Temp HW 3

Pick four out of six problems below. Students who try to solve all six problems will be given extra credit.

Problem 1

The energy momentum-dispersion of graphene valence and conduction bands is given by the formula derived in class:

$E(k_{x},k_{y})=\pm t{\sqrt {1+4\cos \left({\frac {\sqrt {3}}{2}}k_{x}a\right)\cos \left({\frac {k_{y}a}{2}}\right)+4\cos ^{2}\left({\frac {k_{y}a}{2}}\right)}}$

which assumes tight-binding Hamiltonian projected to a basis set of single $\pi$ orbital per carbon atom and hopping allowed only between the nearest neighbor carbon atoms. The honeycomb lattice spacing is $a$.

(a) Using MATLAB or Mathematica plot this function as a surface in 3D k-space within the range of $(k_{x},k_{y})$ vector values defined by the first Brillouin zone.

(b) Near one of the two inequivalent K points [the so-called Dirac points at which $E(k_{x},k_{y})\equiv 0$] in the corners of the first Brillouin zone, where the conduction (top surface) and valence (bottom surface) bands touch, the dispersion has a conical shape characteristics of massless relativistic particles (such as photons or neutrinos). Up to what energy away from the Dirac point $E=0$ (which is the Fermi energy of ideal undoped graphene) can the real dispersion shown in your figure be approximated by linear dispersion of the "Dirac cones"? You can subtract $E(k_{x},k_{y})$ surface around one of the K points from the conical surface $E(q_{x},q_{y})=\pm v_{F}\hbar {\sqrt {q_{x}^{2}+q_{y}^{2}}}$ (corrections are $O[(q/K)]^{2}$; here $\mathbf {q} =\mathbf {K} +\mathbf {k}$ and the Fermi velocity is $v_{F}={\frac {\sqrt {3}}{2}}ta/\hbar \approx 10^{6}m/s$) and find at which energy $\pm E_{0}$ their difference starts to be more than 5%.

Problem 2

The honeycomb lattice of graphene is topologically equivalent to a brick lattice shown in the Figure below:

Honeycomb lattice of graphene (top left), its equivalent brick lattice (top right) and submatrices of the matrix representation of the tight-binding Hamiltonian defined on either of these two lattices.

Thus, the matrix representation of the Hamiltonian defined on this lattice can be obtained in the same fashion as for the square lattice (used in Problem 4. of Homework Set 2), i.e., by loading its submatrices $\mathbf {H} _{0}$ describing graphene supercells coupled to adjacent supercells via the hopping submatrices $\mathbf {H} _{1}$. The carbon atoms within the supercells as well as the structure of these submatrices are shown in the Figure. Here it is assumed that electron can hop between only the nearest neighbor carbon atoms with probability $t$.

(a) Using the honeycomb (or equivalent brick) lattice of $59\times 30$ carbon atoms (for example, the honeycomb lattice in the top inset of the Figure above is $5\times 4$ in this terminology), construct the full Hamiltonian matrix $\mathbf {H}$ of finite size graphene sheet and compute its density of states (DOS) via the Green function method used in Homework Set 2 in the energy range $[-3.5t,3.5t]$.

(b) Introduce periodic boundary into $\mathbf {H}$ (this requires to input two additional submatrices $\mathbf {H} _{1}$ into $\mathbf {H}$, as well as to introduce extra hoppings $t$ into $\mathbf {H} _{0}$) and compute its DOS via the same method as in (a).

(c) Does any of your numerically obtained DOS for finite-size graphene sheet resembles analytical DOS shown in Fig. 5 (page 114) of A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys. 81, 109 (2009). [PDF]?

HINT: Test if your Hamiltonian is properly set up with the MATLAB function visual_graphene_H.m posted on the Computing page of PHYS 824 wiki. The function will plot real space position of the carbon atoms and the hoppings between them, based on the matrix content and lattice size you supply as an input.

Problem 3

Massive Dirac particles in graphene: Consider a tight-binding model on a honeycomb lattice with on-site potential1 different on the sublattices A and B, $V_{A}=-V_{B}=\lambda$.

a) Develop a general solution of this problem in terms of plane wave states, expressing the spectrum of excitations in terms of the nearest neighbor hopping amplitude $t$ and the asymmetry parameter $\Delta ={\frac {1}{2}}\left(V_{A}-V_{B}\right)$.

b) Linearize the solution found in part a) near the points $K$ and $K^{\prime }$ treating the asymmetry as a small perturbation, and find how the massless Dirac picture (valid for $V_{A}=V_{B}=0$) of low energy states of graphene is altered. That is, show that $\Delta$ generates a finite mass of Dirac quasiparticles whose Hamiltonian is given by:

${\hat {H}}=v_{F}{\hat {\mathbf {\sigma } }}\cdot {\hat {\mathbf {p} }}+\lambda {\hat {\sigma }}_{z}$,

where $\lambda =V_{A}=-V_{B}$.

Problem 4

Bilayer graphene: The graphene bilayers are expected to play an important role in nanoelectronic device fabrication since one can manipulate their electronic structure and band gaps (see, e.g., Nature Materials 7, 151 (2007)] Suppose that two graphene layers are stacked atop one another according to the "Bernal stacking" (the stacking fashion of graphite). This means that the B sublattice sites of the upper layer (layer 1) are directly above the A sublattice sites of the lower layer (layer 2). The A sublattice sites of the upper layer and the B sublattice sites of the lower layer have no "partner" atoms below/above them, as shown in the Figure below.

Lattice of bilayer graphene.

(a) Formulate and solve the tight-binding model for this system consisting of the usual hopping $t$ between nearest-neighbor sites in each layer, and an additional smaller hopping $t_{b}$ between the B1-A2 sites in different layers. Choose the zero of energy to equal that of an isolated atom. You should find 4 bands.

(b) Two of the bands found above touch at the Fermi energy $E_{F}=0$ at the Brillouin zone boundary points ±K. Find the effective mass around these points.

(c) Calculate the Berry phase acquired by an electron encircling one of the zone boundary points in each of the two touching bands. Use your result to argue that the bands cannot split if a small perturbation is applied provided inversion and time-reversal symmetry are preserved. Where is the inversion center for the bilayer? REFERENCE: The Berry phase in graphene and graphite multilayers.

(d) Check this conclusion in a simple case: Because the B1 and A2 sites have more atoms close by than the B2 and A1 sites, there will generally be some difference of the site energies for these orbitals. Add an energy $+U_{0}$ to the electrons on the B1, A2 sites and Failed to parse (syntax error): {\displaystyle −U_0 } to the electrons on the B2, A1 sites. Show that for small enough $U_{0}$, the two bands still touch.

REFERENCE: B. Partoens and F. M. Peeters, From graphene to graphite: Electronic structure around the K point, Phys. Rev. B 74, 075404 (2006). [PDF]

Problem 5

Chiral dynamics and Klein tunneling of low-energy quasiparticles in graphene: Consider motion of massless Dirac fermions in an external potential, as described by the Hamiltonian:

${\hat {H}}_{\mathrm {Weyl} }=v_{F}{\hat {\mathbf {\sigma } }}\cdot \mathbf {p} +U(x)$

(a) For a rectangular potential barrier, $U(x)=V_{0},\,0<x<D$ and $U(x)=0$ otherwise, find transmission as a function of the incidence angle.

(b) For a step-like potential $U(x)=U_{0}sign(x)$ find transmission as a function of the incidence angle. For simplicity, consider the case of particle at zero Fermi energy (zero doping).

(c) What happens at normal incidence $k_{y}=0$ in both (a) and (b)?

REFERENCE: M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Chiral tunnelling and the Klein paradox in graphene, Nature Physics 2, 620 (2006). [PDF]

Problem 6

Graphene in magnetic field: Find the energy spectrum of massless Dirac fermions (as the low-energy quasiparticles in graphene) in the external magnetic field $B=B{\hat {z}}$ which is orthogonal to the graphene xy-plane. The (Weyl) Hamiltonian in magnetic field is given by:

${\hat {H}}_{\mathrm {Weyl} }=v_{F}{\hat {\mathbf {\sigma } }}\cdot (\mathbf {p} -e\mathbf {A} )$

where you can employ the so-called Landau gauge $\mathbf {A} =(-By,0,0)$ to diagonalize this Hamiltonian. The vector of Pauli matrices is ${\hat {\mathbf {\sigma } }}=({\hat {\sigma }}_{x},{\hat {\sigma }}_{y},{\hat {\sigma }}_{z})$. HINT: ${\hat {H}}_{\mathrm {Pauli-Schrodinger} }=2m{\hat {H}}_{\rm {Weyl}}^{2}$.