Temp

Pick four out of six problems below. Students who 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 {\sqrt {3+4\cos \left({\frac {\sqrt {3}}{2}}k_{y}a\right)\cos \left({\frac {3}{2}}k_{x}a\right)+2\cos({\sqrt {3}}k_{y}a)}}$

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:

File:Brick lattice.png

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 3$ 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

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} =B(-y,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}$.

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 Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle U(x)=0} otherwise, find transmission as a function of the incidence angle.

(b) For a step-like potential Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 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 Failed to parse (SVG with PNG fallback (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://en.wikipedia.org/api/rest_v1/":): {\displaystyle 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]