Homework Set 2

Problem 2: Tight-binding Hamiltonian of one-dimensional nanowire on the lattice with a basis

Consider a tight-binding model of a 1D nanowire:

${\hat {H}}=\sum _{m}\cos \left(2\pi m{\frac {5}{3}}\right)|m\rangle \langle m|+t\sum _{m}\left(|m\rangle \langle m+1|+|m+1\rangle \langle m|\right)$,

The integer $m$ is indexing sites at which the atoms are located. The distance between two sites defines the lattice spacing $a$, while the nearest neighbor hopping $t$ sets the unit of energy. The ket $|m\rangle$ is quantum state of an electron on atom $m$, so that $\langle x|m\rangle =\psi (x-m)$ is the corresponding wave function in coordinate representation (or single "orbital" per site) which decays fast away from the position of an atom $m$.

(a) What is the periodicity of the Hamiltonian? (That is, after how many sites the chain starts to repeat itself? The atoms on those sites define the unit cell of the wire whose periodic repetition in both direction generates the whole wire.)

(b) Use the Bloch theorem to reduce the eigenvalue problem of an infinite matrix $\mathbf {H}$, obtained by representing the Hamiltonian in the basis of orbitals $|m\rangle$, to diagonalization of a small matrix [whose size $n$ is equal to the periodicity of the Hamiltonian found in (a)].

(c) The $n\times n$ matrix in (b) will depended on the Bloch wave vector $k$. For each value of $k$, diagonalize this matrix and plot $n$ bands $E_{n}(k)$ where $k$ vector belongs to the first Brillouin zone (this task will have to be carried our numerically).

Problem 3: Density of states in one-dimensional nanowire with a single impurity

A nanowire consists of 500 atoms described by a 1D tight-binding Hamiltonian:

${\hat {H}}=\sum _{m}\varepsilon _{m}|m\rangle \langle m|+t\sum _{m}\left(|m\rangle \langle m+1|+|m+1\rangle \langle m|\right)$.

(a) Compute numerically the density (DOS) of states for this wire assuming periodic boundary conditions and $\varepsilon _{m}=0$. In numerical calculations use $t=1$ as the unit of energy. How does DOS change if you increase the number of atoms from 500 to 5000?

(b) Replacement of the original atom in the middle of the chain by an impurity atom can be modeled by using $\varepsilon _{250}=10t$. Compute (DOS) for this case and comment on differences between (a) and (b). What is the highest eigenenergy in (a) vs. (b)? What is the physical meaning of this energy in the case (b)?

(c) Use the Green function method to compare its result for DOS to (a) and (b) obtained using eigenvalue method for 500 atom wire.

Problem 4: Hofstadter butterfly of electrons on square tight-binding lattice in external magnetic field

The Hofstadter butterfly is the energy spectrum of an infinite square lattice plotted as a function of the magnetic field. In this problem we will examine similar spectra for finite lattices modeling arrays of quantum dots in a magnetic field using methods that consider the appropriate molecular orbitals and compare their spectra to the Hofstadter butterfly.

(a) Reproduce panels (a)-(f) in Fig. 5 of American Journal of Physics 72, 5 (2004) for small arrays up to $10\times 10$.

(b) If you increase your quantum dot array size to $100\times 100$, does your $E\ \mathrm {vs.} \ \alpha$ plot resemble Hofstadter butterfly on an infinite lattice plotted in Fig. 6 of the same reference as in (a).

Problem 5: Shubnikov-de Haas oscillations in two-dimensional electron gas

Shubnikov-de Haas oscillations in the longitudinal voltage.

The physical properties of a two-dimensional electron gas (2DEG), in semiconductor heterostructures or within 2D materials, can be measured by embedding it into a multiterminal Hall bridge, as shown in the Figure above. The energy spectrum of 2D electrons in magnetic field $\mathbf {B} =B{\hat {z}}$ orthogonal to such bridge consists of discrete Landau lavels:

$E_{n}=\hbar \omega \left(n+{\frac {1}{2}}\right),\ \omega =eB/m^{*}$,

which are well-separated for large enough B. Thus, the step-like density of states associated with 2DEG $g(E)={\frac {m^{*}}{\pi \hbar ^{2}}}\theta (E-E_{b})$ ($E_{b}$ is the lowest level in the triangular potential that confines 2DEG) is then broken up into a series of peaks spaced by $\hbar \omega$:

$g(E,B)={\frac {2eB}{h}}\sum _{n=0}^{\infty }\delta \left[E-E_{b}-\hbar \omega \left(n+{\frac {1}{2}}\right)\right]$. (1)

Although the spikes are ideally delta functions, scattering processes spread them out in energy. As one changes the magnetic field B, longitudinal voltage (from which one computes the longitudinal resisitivity $\rho _{xx}=V_{x}W/IL$) in the Figure above goes through one cycle of oscillations as the Fermi energy moves from the center of one Landau level to the center of the next Landau level. These oscillations are refered to as the Shubnikov-de Haas (SdH) oscillations.

(a) Devise a procedure that will allow you to extract electron density $n$ in 2DEG from the values of magnetic field B at which $V_{x}$ (or the corresponding $\rho _{xx}$) reaches maximum in the course of SdH oscillations.

(b) What is $n$ for the 2DEG in the Figure above?

HINT: From Eq. (1) we see that the number of occupied Landau levels is $n/(2eB/h)$, if electron density $n$ is already known.

Problem 6: Energy-momentum dispersion in graphene

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 7: Density of state of graphene flake

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 Rev. Mod. Phys. 81, 109 (2009). [PDF]?