semic

Revamp lectures 13 and 14.

src/13_semiconductors.md

@@ -46,15 +46,15 @@ _(based on chapters 17–18 of the book)_
 
 Up until this point, we focused on calculating and understanding the band structures.
 However, the dispersion of a band is only part of the story.
-An empty band is not gonna lead to any interesting physical properties no matter how sophisticated it is.
+An empty band is not going to lead to any interesting physical properties no matter how sophisticated it is.
 Therefore, it is also important *how* bands are filled by the particles.
 By carefully controlling the distribution of particles in the bands, we are able to engineer material properties that we require.
-Without a doubt, the greatest example is *semiconductors* - the bedrock of all electronics. 
+Without a doubt, the greatest example is *semiconductors*—the bedrock of modern electronics. 
 In this lecture, we shall grasp the basics of semiconductors by learning how to treat bands at different levels of filling.
 
 ## Review of band structure properties
 
-Before proceeding further, it is vital to remind ourselves some of the concepts in band structures. 
+Before proceeding further, let us remind ourselves of important band structure properties.
 
 * Group velocity $v=\hbar^{-1}\partial E(k)/\partial k$. 
 Descibes how quickly electrons move within the lattice. 
@@ -66,68 +66,69 @@ The quantity is vital in order to calculate any bulk property of the material su
 
 In order to check that everything makes sense, we apply the concepts to the free electron model:
 
-$$H = \hbar^2 k^2/2m$$
+$$H = \frac{\hbar^2 k^2}{2m}$$
 
 The velocity is $\hbar^{-1}\partial E(k)/\partial k = \hbar k / m \equiv p/m$.  
 The effective mass is $m_{eff} = \hbar^2\left(d^2 E(k)/dk^2\right)^{-1} = m$.
 
-So in the simplest case the definitions match the usual expressions.
+So in this simplest case the definitions match the usual expressions.
 
 ## Filled vs empty bands
 
 We distinguish three different band filling types: filled, empty and partially filled. 
-Let us start with the most extreme cases - filled and empty bands.
-We treat these two cases together, because a completely filled band is very similar to a completely empty band.
-For example, both filled and empty bands lead to zero current:
+Despite being two opposite extreme extreme cases, filled and empty bands are very similar.
+For example, both filled and empty bands carry no electric current:
 $$
 \begin{align}
 j = 2e \frac{1}{2\pi} \int_{-\pi/a}^{\pi/a} v(k) dk = 2e \frac{1}{2\pi \hbar} \int_{-\pi/a}^{\pi/a} \frac{dE}{dk} \times dk = \\
 2e \frac{1}{2\pi \hbar} [E(-\pi/a) - E(\pi/a)] = 0.
 \end{align}
 $$
-An empty band has no electrons and thus no current.
 On the other hand, a filled band has an equal number of electrons going forwards and backwards which thus cancel and lead to zero current. 
-Similar results apply to many other physical quantities such as heat capacity and magnetisation. 
+Similar results apply to many other physical quantities such as heat capacity and magnetization. 
 Therefore, filled and empty bands do not affect most physical properties and can be disregarded.
-As a result, rather than to consider thousands of bands that a material contains, we neglect most of them and just focus on partially filled bands around Fermi level.
+As a result, rather than to consider hundreds of bands that a material contains, we neglect most of them and just focus on the handful of partially filled bands around Fermi level.
 
 ## From electrons to holes
 
-In order to understand partial filling, we begin with a simple analogy.
-Let's say we have 100 chairs: 99 are occupied and 1 is empty.
-To keep track which chair is occupied/empty, we could write down the occupation of each and every chair.
-However, that would require a lot of unnecessary writing!
-Instead, we only need to write down which chair is empty, because we know the rest wil be occupied. 
-The same philosophy is applied to band filling. 
+Because completely filled or completely empty bands have simple properties, we may search for a convenient way to describe a band that only has a few electrons missing or extra.
+While keeping track of a few electrons has no tricks, even a few electrons missing from a band seem to require considering all the other electrons in a band.
+
+A more efficient approach to describing a nearly filled band is motivated by the following analogy.
+Let us say we have 100 boxes: 99 are occupied and 1 is empty.
+To keep track which box is occupied/empty, we could write down the numbers of all 99 occupied boxes.
+If, on the other hand, we only keep track which single box is empty, we solve the problem with a lot less book-keeping.
+
+The same approach applies to band filling.
 Instead of describing a lot of electrons that are present in an almost filled band, we focus on those that are absent.
 The absence of an electron is called a **hole**: a state of a completely filled band with one particle missing.
 
 ![](figures/holes.svg)
 
 In this schematic we can either say that 8×2 electron states are occupied (the system has 8×2 electrons counting spin), or 10×2 hole states are occupied.
-A useful analogy to remember: glass half-full or glass half-empty.
 Electron and hole pictures correspond to two different, but equivalent ways of describing the occupation of a band.
-In practice, we choose to deal with electrons whenever a band is almost empty and holes when a band is almost full.
+Naturally, dealing with electrons is more convenient whenever a band is almost empty and with holes when a band is almost full.
 
 ## Properties of holes
 
-Let us compare the properties of electrons and holes.
+Let us compare the properties of an electron and a hole obtained by removing that electron.
+Since removing an electron reduces the total energy of the system, the hole's energy is opposite to that of an electron $E_h = -E$.
 The probability for an electron state to be occupied in equilibrium is given by $f(E)$:
 
 $$f(E) = \frac{1}{e^{(E-E_F)/kT} + 1}.$$
 
-On the other hand, since a hole is a missing electron, the probability for a hole state to be occupied is
+Since a hole is a missing electron, the probability for a hole state to be occupied is
 
 $$f_h(E) = 1 - f(E) = 1 - \frac{1}{e^{(E-E_F)/kT} + 1} = \frac{1}{e^{(-E+E_F)/kT} + 1},$$
 
-therefore for holes both energy $E_h= -E$ and $E_{F,h} = -E_F$.
+which is the Fermi distribution of particles with energy $E_h= -E$ and $E_{F,h} = -E_F$.
 
 The **momentum** $p_h$ of a hole should give the correct total momentum of a partially filled band if one sums momenta of all holes.
 Therefore $p_h = -\hbar k$, where $k$ is the wave vector of the electron.
 
-Similarly, the total **charge** should be the same regardless of whether we count electrons or holes, so holes have a positive charge $+e$ (electrons have $-e$).
+Similarly, the total **charge** should be the same regardless of whether we count electrons or holes, so holes have a positive charge $+e$ (electrons having $-e$).
 
-On the other hand, the velocity of a hole is **the same**:
+On the other hand, hole's velocity is **the same** as that of an electron:
 $$\frac{dE_h}{dp_h} = \frac{-dE}{-d\hbar k} = \frac{dE}{dp}.$$
 
 Finally, we derive the hole effective mass from the equations of motion:
@@ -138,42 +139,56 @@ Comparing with
 
 $$m_e \frac{d v}{d t} = -e (E + v\times B),$$
 
-we get $m_h = -m_e$.
+we get $m_h = -m_e$ (we could also obtain this by differentiating the hole's velocity).
 
 ## Semiconductors: materials with two bands.
 
-In semiconductors the Fermi level is between two bands.
-The unoccupied band is the **conduction band**, the occupied one is the **valence band**. 
-In the conduction band the **charge carriers** (particles carrying electric current) are electrons, in the valence band they are holes.
+Semiconductors are materials with all bands either nearly occupied or almost empty.
+Unlike in insulators, however, the band gap in semiconductors is sufficiently small, for it to be possible to create a few electrons in the lowest unoccupied band or the highest filled band.
+Because in the unoccupied band the **charge carriers** (particles carrying electric current) are electrons, it is called **conduction band**, while in the almost occupied **valence band** the charge carriers are holes.
 
-We can control the position of the Fermi level (or create additional excitations) making semiconductors conduct when needed.
+!!! note "Holes in semiconductors"
 
-Only the bottom of the conduction band has electrons and the top of the valence band has holes because the temperature is always smaller than the size of the band gap.
+    When introducting holes, we discussed holes obtained by removing *any* electron.
+    From this point on, similar to most other discussions of semiconductors in the literature, we will only speak of holes in valence band and electrons in conduction band.
 
+The occupation of the two bands is dictated by the Fermi distribution.
+Furthermore, the Fermi level of a semiconductor lies between the conduction and the valence bands, and the band gap $E_G \gg k_B T$ in most materials.
+As a result, only the bottom of the conduction band has electrons and the top of the valence band has holes.
 Therefore we can approximate the dispersion relation of both bands as parabolic.
 
 ![](figures/semiconductor.svg)
 
 Or in other words
 
-$$E_e = E_c + \frac{\hbar^2k^2}{2m_e}$$
-$$E_h = E_{v,h} + \frac{\hbar^2k^2}{2m_h} = -E_{v} + \frac{\hbar^2k^2}{2m_h},$$
-
-with the corresponding density of states
-
-$$ g(E) = (2m_e)^{3/2}\sqrt{E-E_c}/2\pi^2\hbar^3$$
-$$ g_h(E_h) = (2m_h)^{3/2}\sqrt{E_h+E_v}/2\pi^2\hbar^3.$$
+$$
+\begin{align}
+E_e &= E_c + \frac{\hbar^2k^2}{2m_e},\\
+E_h &= E_{v,h} + \frac{\hbar^2k^2}{2m_h} = -E_{v} + \frac{\hbar^2k^2}{2m_h}.
+\end{align}
+$$
 
 Here $E_c$ is the energy of an electron at the bottom of the conduction band and $E_v$ is the energy of an electron at the top of the valence band.
-
 Observe that because we are describing particles in the valence band as holes, $m_h > 0$ and $E_h > -E_v$.
 
-??? question "a photon gives a single electron enough energy to move from the valence band to the conduction band. How many particles does this process create?"
+The corresponding density of states of the two types of particles is
+
+$$
+\begin{align}
+g(E) &= (2m_e)^{3/2}\sqrt{E-E_c}/2\pi^2\hbar^3,\\
+g_h(E_h) &= (2m_h)^{3/2}\sqrt{E_h+E_v}/2\pi^2\hbar^3.
+\end{align}
+$$
+
+??? question "A photon gives a single electron enough energy to move from the valence band to the conduction band. How many particles does this process create?"
     Two: one electron and one hole.
 
 ## Semiconductor density of states and Fermi level
 
-### Part 1: Intrinsic semiconductor
+### Intrinsic semiconductor
+
+Our next task is to figure out *how many* electrons and holes there are, and for that we need to find where the Fermi level $E_F$ is located.
+Let us plot the density of states, the Fermi distribution function, and the density of particles at each energy in the same plot:
 
 ```python
 E = np.linspace(-3, 3, 1000)
@@ -197,6 +212,9 @@ ax.legend()
 draw_classic_axes(ax, xlabeloffset=.2)
 ```
 
+We know that by itself, the semiconductor should have no charge, and therefore the total numbers of electrons and holes must be equal.
+Since increasing the Fermi level increases the number of electrons and reduces the number of holes, we will use the charge neutrality condition to determine where the Fermi level is situated.
+
 **The key algorithm of describing the state of a semiconductor:**
 
 1. Compute the density of states of all types of particles.
@@ -205,29 +223,36 @@ draw_classic_axes(ax, xlabeloffset=.2)
 4. Apply approximations to simplify the equations (this is important!).
 5. Find $E_F$ and concentrations of electrons and holes
 
-Applying the algorithm:
+Applying the first two steps of the algorithm:
 
 $$n_h = \int_{-E_v}^\infty f_h(E_h) g_h(E_h) dE_h = \int_{-E_v}^\infty\frac{(2m_h)^{3/2}}{2\pi^2\hbar^3}\sqrt{E_h+E_v}\frac{1}{e^{(E_h+E_F)/kT}+1}dE_h$$
 
-$$n_e = \int_{E_c}^\infty f(E)g_e(E)dE = \int_{E_c}^\infty\frac{(2m_e)^{3/2}}{2\pi^2\hbar^3} \sqrt{E-E_c}\frac{1}{e^{(E-E_F)/kT}+1}dE$$
-
-We need to solve $n_e = n_h$
+$$n_e = \int_{E_c}^\infty f(E)g_e(E)dE = \int_{E_c}^\infty\frac{(2m_e)^{3/2}}{2\pi^2\hbar^3} \sqrt{E-E_c}\frac{1}{e^{(E-E_F)/kT}+1}dE.$$
 
-Simplification:
-Fermi level is far from both bands $E_F-E_v \gg kT$ and $E_c - E_F \gg kT$
+Note that whenever calculating the hole dependent quantities, we replace all the relevant physical quantities with their hole equivalents.
+Since the hole energy is opposite $E_h = -E$, we replace the Fermi energy $E_F \to -E_F$ and the bottom of the valance band by $E_v \to -E_v$ in the integration limits.
 
-Therefore Fermi-Dirac distribution is approximately similar to Boltzmann distribution.
+In the third step, we need to solve the equation under charge balance $n_e = n_h$.
+The equation is not a pleasant one and cannot be solved analytically unless an approximation is made.
+Therefore, the fourth step assumes that the Fermi level is far from both bands $E_F-E_v \gg kT$ and $E_c - E_F \gg kT$.
+As a result, the Fermi-Dirac distribution is approximately similar to Boltzmann distribution:
 
-$$f(E)_{e/h} \approx e^{-(E_{e/h}\pm E_F)/kT}$$
+$$
+f(E)_{e/h} \approx \exp\left[-(E_{e/h}\pm E_F)/kT\right].
+$$
 
-Now we can calculate $n_e$ and $n_h$:
+Now we can move to the last step and calculate $n_e$ and $n_h$:
 
-$$n_h \approx \frac{(2m_h)^{3/2}}{2\pi^2\hbar^3}e^{-E_F/kT} \int_{-E_v}^\infty\sqrt{E_h+E_v}e^{-E_h/kT}dE_h =
-N_V e^{E_v-E_F/kT},$$
+$$
+n_h \approx \frac{(2m_h)^{3/2}}{2\pi^2\hbar^3}e^{-E_F/kT} \int_{-E_v}^\infty\sqrt{E_h+E_v}e^{-E_h/kT}dE_h =
+N_V e^{E_v-E_F/kT},
+$$
 
 where we used $\int_0^\infty \sqrt{x}e^{-x}dx=\sqrt{\pi}/2$ and we defined 
 
-$$N_V = 2\left(\frac{2\pi m_h kT}{h^2}\right)^{3/2}$$
+$$
+N_V = 2\left(\frac{2\pi m_h kT}{h^2}\right)^{3/2}.
+$$
 
 We see that holes are exponentially activated into the valence band. 
 
@@ -236,16 +261,22 @@ We see that holes are exponentially activated into the valence band.
 
 Similarly for electrons:
 
-$$n_e = N_C e^{-(E_c - E_F)/kT},\quad N_C = 2\left(\frac{2\pi m_e kT}{h^2}\right)^{3/2}$$
+$$
+n_e = N_C e^{-(E_c - E_F)/kT},\quad N_C = 2\left(\frac{2\pi m_e kT}{h^2}\right)^{3/2}.
+$$
 
 
 Combining everything together:
 
-$$n_h \approx N_V e^{E_v-E_F/kT} = N_C e^{-(E_c-E_F)/kT} \approx n_e$$
+$$
+n_h \approx N_V e^{E_v-E_F/kT} = N_C e^{-(E_c-E_F)/kT} \approx n_e.
+$$
 
 Solving for $E_F$:
 
-$$E_F = \frac{E_c + E_v}{2} - \frac{3}{4}kT\log(m_e/m_h)$$
+$$
+E_F = \frac{E_c + E_v}{2} - \frac{3}{4}kT\ln(m_e/m_h).
+$$
 
 An extra observation: regardless of where $E_F$ is located, $n_e n_h = N_C N_V e^{-E_g/kT} \equiv n_i^2$, where $E_g=E_c-E_v$ is the band gap of the semiconductor.
 
@@ -255,6 +286,25 @@ $n_i$ is the **intrinsic carrier concentration**, and for a pristine semiconduct
 > $$n_e n_h = n_i^2$$
 > is the **law of mass action**. The name is borrowed from chemistry, and describes the equilibrium concentration of two reagents in a reaction $A+B \leftrightarrow AB$. Here electrons and hole constantly split and recombine.
 
+### Conduction
+
+Earlier, we deduced that empty and filled bands provide no current. 
+We finish the analysis by considering partially filled bands of an intrinsic (pristine) semiconductor.
+To calculate the current, we utilize the Drude model and sum the electron and hole contributions:
+
+$$j = -n_e e v_e + n_h e v_h $$
+
+$$ -m_e v_e /\tau_e = -eE;\quad -m_h v_h /\tau_h = eE.$$
+
+We see that despite opposite velocity signs for electrons and holes, they carry electric current in the same direction:
+
+$$ \sigma \equiv \frac{j}{E} = \left(\frac{n_e e^2 \tau_e}{m_e}+\frac{n_h e^2 \tau_h}{m_h}\right) = n_e e \mu_e + n_h e \mu_h.$$
+
+We know that for intrinsic semiconductors, the hole/electron densities are $n_e = n_h = n_i \propto e^{-E_G/kT}$.
+Therefore, it is possible to measure the band gap of an intrinsic semiconductor by looking at the temperature dependant conductivity $E_G \approx d \ln \sigma / d [kT]^{-1}$.
+Additional information can be obtained using Hall effect.
+However Hall effect is much more complex in semiconductors since only the current in the direction perpendicular to the applied electric field must vanish.
+This, however only means that the electron current is opposite of the hole current in that direction, not that the electrons and holes move parallel to the applied current.
 
 ## Exercises