Chapter 16: Semiconductor Optics

Source: Saleh and Teich, Fundamentals of Photonics, second edition, Chapter 16.

In-text exercises

Exercise 16.1-1 — Free-electron dispersion

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.1-1, Free-electron dispersion

Figure 93 — Exercise 16.1-1: Free-electron dispersion. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses exponential, logarithmic, and phasor identities and algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Substitute \(e^{-jkx}\) in the free Schrödinger equation to get \(\boxed{E=\hbar^2k^2/(2m_0)=p^2/(2m_0)}\). Group velocity \(\hbar^{-1}dE/dk=\hbar k/m_0=p/m_0\) matches the particle velocity.

Step 4 — State the numbered result. The principal result obtained in the working is

(1)\[\boxed{E=\hbar^2k^2/(2m_0)=p^2/(2m_0)}\]

Step 5 — Check. Equation (1) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Exercise 16.1-2 — Boltzmann limit of Fermi occupation

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.1-2, Boltzmann limit of Fermi occupation

Figure 94 — Exercise 16.1-2: Boltzmann limit of Fermi occupation. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses integration identities, exponential, logarithmic, and phasor identities, and algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

When \(E-E_F\gg kT\), \(f\simeq e^{-(E-E_F)/kT}\); similarly \(1-f\simeq e^{-(E_F-E)/kT}\) below \(E_F\). Integrating these with the band densities of states gives the chapter’s \(n=N_ce^{-(E_c-E_F)/kT}\) and \(p=N_ve^{-(E_F-E_v)/kT}\).

Step 4 — State the numbered result. The principal result obtained in the working is

(2)\[p=N_ve^{-(E_F-E_v)/kT}\]

Step 5 — Check. Equation (2) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Differentiation of an antiderivative, or normalization of a definite integral, checks the integration step.

Exercise 16.1-3 — Quasi-Fermi levels

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.1-3, Quasi-Fermi levels

Figure 95 — Exercise 16.1-3: Quasi-Fermi levels. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

At zero temperature fill the 3-D density of states to \(k_F=(3\pi^2n)^{1/3}\): \(E_{Fc}=E_c+\hbar^2(3\pi^2n)^{2/3}/(2m_e)\) and \(E_{Fv}=E_v-\hbar^2(3\pi^2p)^{2/3}/(2m_h)\). In the Boltzmann regime invert the preceding exponentials instead.

Step 4 — State the numbered result. The principal result obtained in the working is

(3)\[E_{Fv}=E_v-\hbar^2(3\pi^2p)^{2/3}/(2m_h)\]

Step 5 — Check. Equation (3) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Exercise 16.1-4 — Injected GaAs carriers

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.1-4, Injected GaAs carriers

Figure 96 — Exercise 16.1-4: Injected GaAs carriers. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

\(p_0=n_i^2/n_0\); low-injection lifetime is \(\tau=[r(n_0+p_0)]^{-1}\). Steady excess solves \(R=r[(n_0+\Delta n)(p_0+\Delta n)-n_i^2]\); insert the stated values and select the positive quadratic root, then compute the two quasi-Fermi levels from Exercise 16.1-3.

Step 4 — State the numbered result. The principal result obtained in the working is

(4)\[R=r[(n_0+\Delta n)(p_0+\Delta n)-n_i^2]\]

Step 5 — Check. Equation (4) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Exercise 16.1-5 — Infinite quantum well

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.1-5, Infinite quantum well

Figure 97 — Exercise 16.1-5: Infinite quantum well. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses trigonometric and small-angle identities and algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Boundary conditions \(\psi(0)=\psi(d)=0\) select \(\psi_q=\sqrt{2/d}\sin(q\pi x/d)\) and \(\boxed{E_q=\hbar^2q^2\pi^2/(2md^2)}\). A finite well has lower energies, evanescent tails, and only finitely many bound roots of its tangent/cotangent equations.

Step 4 — State the numbered result. The principal result obtained in the working is

(5)\[\boxed{E_q=\hbar^2q^2\pi^2/(2md^2)}\]

Step 5 — Check. Equation (5) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. The zero-angle or paraxial limit supplies an independent sign and magnitude check whenever that limit is part of the model.

Exercise 16.2-1 — Semiconductor gain condition

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.2-1, Semiconductor gain condition

Figure 98 — Exercise 16.2-1: Semiconductor gain condition. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Thermal equilibrium has one Fermi level and detailed balance makes absorption larger. In quasi-equilibrium, emission exceeds absorption when \(\boxed{E_{Fc}-E_{Fv}>h\nu}\) for the same-k states—the Bernard–Duraffourg population-inversion condition.

Step 4 — State the numbered result. The principal result obtained in the working is

(6)\[\boxed{E_{Fc}-E_{Fv}>h\nu}\]

Step 5 — Check. Equation (6) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Exercise 16.2-2 — Peak direct absorption

Step 1 — Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Illustrated calculation map for Exercise 16.2-2, Peak direct absorption

Figure 99 — Exercise 16.2-2: Peak direct absorption. The diagram identifies the input quantities, physical operation, requested result, variable meanings, and an independent verification route. Every symbol in the variable strip is labeled on the model itself.

Step 2 — Mathematical formulas used. The working uses stationary-value condition, product, quotient, and chain rules, and algebraic rearrangement and dimensional checks.

Step 3 — Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Differentiate the equilibrium direct-gap form \(\alpha\propto\sqrt{h\nu-E_g}/\nu\); its maximum occurs at \(h\nu=2E_g\). Therefore \(\boxed{\lambda_p=hc/(2E_g)}\); for GaAs \(E_g=1.42\) eV, \(\boxed{\lambda_p=436.6\ \mathrm{nm}}\).

End-of-chapter problems

Step 4 — State the numbered result. The principal result obtained in the working is

(7)\[\boxed{\lambda_p=436.6\ \mathrm{nm}}\]

Step 5 — Check. Equation (7) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale. Repeat the substitution with unrounded intermediate values and retain the displayed units; the final unit must have the requested dimension.

Problem 16.1-6 — Hydrogenic donors

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Use \(E_D=13.606(m^*/m_0)/\epsilon_r^2\) eV and \(r_D=a_0\epsilon_r/(m^*/m_0)\). Results (energy, radius) are Si (88.1 meV, 0.664 nm), GaAs (5.63 meV, 9.82 nm), GaN (69.7 meV, 1.65 nm), and polyacetylene (1.51 eV, 0.159 nm). The last two radii approach lattice scale, where bulk dielectric/effective-mass theory is least credible.

Numbered result. The principal result obtained in the working is

(8)\[r_D=a_0\epsilon_r/(m^*/m_0)\]

Check. Equation (8) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale. Repeat the substitution with unrounded intermediate values and retain the displayed units; the final unit must have the requested dimension.

Problem 16.1-7 — Intrinsic and doped Fermi levels

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses exponential, logarithmic, and phasor identities and algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Setting \(n=p\) gives \(\boxed{E_i=(E_c+E_v)/2+(3/4)kT\ln(m_h/m_e)}\). With nondegenerate doping, charge neutrality gives \(E_F=E_i+kT\ln(n/n_i)\) for n type or \(E_F=E_i-kT\ln(p/n_i)\) for p type.

Numbered result. The principal result obtained in the working is

(9)\[\boxed{E_i=(E_c+E_v)/2+(3/4)kT\ln(m_h/m_e)}\]

Check. Equation (9) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Problem 16.1-8 — Strong-injection decay

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses exponential, logarithmic, and phasor identities and algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

After the source turns off, \(d\Delta n/dt=-r(\Delta n)^2\). Separation gives \(\boxed{\Delta n(t)=\Delta n_0/[1+r\Delta n_0(t-t_0)]}\)—a reciprocal power law, not an exponential.

Numbered result. The principal result obtained in the working is

(10)\[\boxed{\Delta n(t)=\Delta n_0/[1+r\Delta n_0(t-t_0)]}\]

Check. Equation (10) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Problem 16.1-9 — Alloy bowing

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

At each plotted composition solve \(\boxed{b=[xE_{AC}+(1-x)E_{BC}-E_g(x)]/[x(1-x)]}\) and average consistent points for every listed alloy. Large \(b\) means gap tuning is strongly nonlinear even while lattice constant follows Vegard; it changes which composition simultaneously provides a desired gap and substrate lattice match.

Numbered result. The principal result obtained in the working is

(11)\[\boxed{b=[xE_{AC}+(1-x)E_{BC}-E_g(x)]/[x(1-x)]}\]

Check. Equation (11) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Problem 16.1-10 — GaAs/AlGaAs well

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Thirty-percent Al raises total gap by \(30(12.47)=374.1\) meV; 60% gives electron barrier \(\boxed{V_0=224.5\ \mathrm{meV}}\). From \(\sqrt{mV_0d^2/(2\hbar^2)}=4\), \(\boxed{d=4\sqrt{2\hbar^2/(mV_0)}}\); solve the finite-well even/odd equations to place the levels inside the drawn conduction/valence offsets.

Numbered result. The principal result obtained in the working is

(12)\[\boxed{d=4\sqrt{2\hbar^2/(mV_0)}}\]

Check. Equation (12) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale. Repeat the substitution with unrounded intermediate values and retain the displayed units; the final unit must have the requested dimension.

Problem 16.2-3 — Delta-lineshape approximation

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses Fourier-transform and convolution identities and algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

At 300 K plot the lifetime Lorentzian width \(1/(\pi T_2)\) beside the \(kT/h\) widths of occupation and joint DOS. For \(T_2=1\) ps the lineshape is much narrower, validating replacement by a delta function for both emission and absorption; numerical convolution quantifies the small error.

Numbered result. The principal result obtained in the working is

(13)\[T_2=1\]

Check. Equation (13) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Transform dimensions and the expected even/odd or conjugate symmetry provide an independent check on signs and scale factors.

Problem 16.2-4 — Thermal spontaneous peak

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses exponential, logarithmic, and phasor identities and algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Maximize \(\sqrt{h\nu-E_g}e^{-h\nu/kT}\) to obtain \(\boxed{h\nu_p=E_g+kT/2}\). Substitution gives the printed closed peak rate. Nondegenerate doping shifts Fermi factors but cancels under thermal mass action until degeneracy; insert GaAs parameters in that formula for the requested numerical rate.

Numbered result. The principal result obtained in the working is

(14)\[\boxed{h\nu_p=E_g+kT/2}\]

Check. Equation (14) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Check that dimensions agree term by term, then test the simplest symmetry or limiting case for the expected sign and scale.

Problem 16.2-5 — Integrated radiative rate

Definitions and setup. Symbols are local to this item and follow the chapter convention. Each physical quantity and supplied numerical value is introduced at its first use below; angles are in radians unless a degree symbol is shown, and units are retained through numerical substitution.

Mathematical formulas used. The working uses integration identities, exponential, logarithmic, and phasor identities, and algebraic rearrangement and dimensional checks.

Worked derivation. The calculation is kept in symbolic form until the governing relation has been rearranged for the requested quantity.

Set \(x=h\nu-E_g\) and use the supplied gamma integral to obtain the printed \((kT)^{3/2}e^{-E_g/kT}\) rate. The peak-times-width estimate has the same scaling. Equating it to \(r_rn_i^2\) gives \(\boxed{r_r=\sqrt2\pi^{3/2}\hbar^3/ [(m_e+m_h)^{3/2}(kT)^{3/2}T_r]}\); GaAs evaluation is of order \(10^{-10}\ \mathrm{cm^3/s}\), consistent with the table.

Numbered result. The principal result obtained in the working is

(15)\[\boxed{r_r=\sqrt2\pi^{3/2}\hbar^3/ [(m_e+m_h)^{3/2}(kT)^{3/2}T_r]}\]

Check. Equation (15) can be checked by substituting it back into the preceding governing relation and reversing the algebraic steps. Differentiation of an antiderivative, or normalization of a definite integral, checks the integration step. Repeat the substitution with unrounded intermediate values and retain the displayed units; the final unit must have the requested dimension.