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

Brief solution

2. Reasoning and answer.

\[\boxed{E=\hbar^2k^2/(2m_0)=p^2/(2m_0)}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. Substitute \(e^{-jkx}\) in the free Schrödinger equation to get \(\boxed{E=\hbar^2k^2/(2m_0)=p^2/(2m_0)}\).

Detailed step 2. 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

Brief solution

2. Reasoning and answer.

\[p=N_ve^{-(E_F-E_v)/kT}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. 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\).

Detailed step 2. 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

Brief solution

1. Method. The working uses algebraic rearrangement and dimensional checks.

2. Reasoning and answer.

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

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.

Detailed step 1. 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)\).

Detailed step 2. 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

Brief solution

1. Method. The working uses algebraic rearrangement and dimensional checks.

2. Reasoning and answer.

\[R=r[(n_0+\Delta n)(p_0+\Delta n)-n_i^2]\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. \(p_0=n_i^2/n_0\); low-injection lifetime is \(\tau=[r(n_0+p_0)]^{-1}\).

Detailed step 2. 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,

Detailed step 3. 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

Brief solution

2. Reasoning and answer.

\[\boxed{E_q=\hbar^2q^2\pi^2/(2md^2)}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. 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)}\).

Detailed step 2. A finite well has lower energies,

Detailed step 3. evanescent tails,

Detailed step 4. 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

Brief solution

1. Method. The working uses algebraic rearrangement and dimensional checks.

2. Reasoning and answer.

\[\boxed{E_{Fc}-E_{Fv}>h\nu}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. Thermal equilibrium has one Fermi level and detailed balance makes absorption larger.

Detailed step 2. In quasi-equilibrium,

Detailed step 3. 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

Show detailed stepsHide detailed steps

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.

Detailed step 1. Differentiate the equilibrium direct-gap form \(\alpha\propto\sqrt{h\nu-E_g}/\nu\); its maximum occurs at \(h\nu=2E_g\).

Detailed step 2. 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

Brief solution

1. Method. The working uses algebraic rearrangement and dimensional checks.

2. Reasoning and answer.

\[r_D=a_0\epsilon_r/(m^*/m_0)\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. Use \(E_D=13.606(m^*/m_0)/\epsilon_r^2\) eV and \(r_D=a_0\epsilon_r/(m^*/m_0)\).

Detailed step 2. Results (energy,

Detailed step 3. radius) are Si (88.1 meV, 0.664 nm),

Detailed step 4. GaAs (5.63 meV, 9.82 nm),

Detailed step 5. GaN (69.7 meV, 1.65 nm),

Detailed step 6. and polyacetylene (1.51 eV, 0.159 nm).

Detailed step 7. The last two radii approach lattice scale,

Detailed step 8. 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

Brief solution

2. Reasoning and answer.

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

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.

Detailed step 1. Setting \(n=p\) gives \(\boxed{E_i=(E_c+E_v)/2+(3/4)kT\ln(m_h/m_e)}\).

Detailed step 2. With nondegenerate doping,

Detailed step 3. 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

Brief solution

2. Reasoning and answer.

\[\boxed{\Delta n(t)=\Delta n_0/[1+r\Delta n_0(t-t_0)]}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. After the source turns off, \(d\Delta n/dt=-r(\Delta n)^2\).

Detailed step 2. Separation gives \(\boxed{\Delta n(t)=\Delta n_0/[1+r\Delta n_0(t-t_0)]}\)—a reciprocal power law,

Detailed step 3. 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

Brief solution

1. Method. The working uses algebraic rearrangement and dimensional checks.

2. Reasoning and answer.

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

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.

Detailed step 1. 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.

Detailed step 2. 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

Brief solution

1. Method. The working uses algebraic rearrangement and dimensional checks.

2. Reasoning and answer.

\[\boxed{d=4\sqrt{2\hbar^2/(mV_0)}}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. Thirty-percent Al raises total gap by \(30(12.47)=374.1\) meV;

Detailed step 2. 60% gives electron barrier \(\boxed{V_0=224.5\ \mathrm{meV}}\).

Detailed step 3. 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

Brief solution

2. Reasoning and answer.

\[T_2=1\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. At 300 K plot the lifetime Lorentzian width \(1/(\pi T_2)\) beside the \(kT/h\) widths of occupation and joint DOS.

Detailed step 2. For \(T_2=1\) ps the lineshape is much narrower,

Detailed step 3. 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

Brief solution

2. Reasoning and answer.

\[\boxed{h\nu_p=E_g+kT/2}\]
Show detailed stepsHide detailed steps

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.

Detailed step 1. Maximize \(\sqrt{h\nu-E_g}e^{-h\nu/kT}\) to obtain \(\boxed{h\nu_p=E_g+kT/2}\).

Detailed step 2. Substitution gives the printed closed peak rate.

Detailed step 3. 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

Brief solution

2. Reasoning and answer.

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

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.

Detailed step 1. Set \(x=h\nu-E_g\) and use the supplied gamma integral to obtain the printed \((kT)^{3/2}e^{-E_g/kT}\) rate.

Detailed step 2. The peak-times-width estimate has the same scaling.

Detailed step 3. 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]}\);

Detailed step 4. GaAs evaluation is of order \(10^{-10}\ \mathrm{cm^3/s}\),

Detailed step 5. 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.