Level 3 — ProductionDoping & PN Junctions

Doping & PN Junctions

45 minutes60 marksprintable — key stays hidden on paper

Difficulty: Level 3 (from-scratch derivations, explain-out-loud, code-from-memory) Time limit: 45 minutes Total marks: 60

Use ...... / ...... for all mathematics. Constants where needed: kB=1.381×1023 J/Kk_B = 1.381\times10^{-23}\ \mathrm{J/K}, q=1.602×1019 Cq = 1.602\times10^{-19}\ \mathrm{C}, thermal voltage VT=kBT/qV_T = k_BT/q.


Question 1 — Doping & band diagram (explain-out-loud) [10 marks]

(a) Explain from scratch, atom by atom, what happens when a phosphorus atom substitutes for a silicon atom in the lattice. State the group of P and Si, the number of valence electrons involved, and why a nearly-free electron results. [4]

(b) Do the same for boron in silicon: explain the creation of a mobile hole. [3]

(c) Draw (describe in words with energy values relative to EcE_c and EvE_v) where the donor level EDE_D and acceptor level EAE_A sit in the band gap, and explain why these levels are "shallow" and easily ionised at room temperature. [3]


Question 2 — Built-in potential derivation from scratch [12 marks]

(a) Starting from the requirement that the total current is zero at equilibrium, derive the built-in potential of a step PN junction: Vbi=VTln ⁣(NANDni2)V_{bi} = V_T \ln\!\left(\frac{N_A N_D}{n_i^2}\right) State clearly each physical principle you invoke (drift–diffusion balance, mass-action law). [7]

(b) Compute VbiV_{bi} at T=300 KT = 300\ \mathrm{K} for silicon with NA=1017 cm3N_A = 10^{17}\ \mathrm{cm^{-3}}, ND=1016 cm3N_D = 10^{16}\ \mathrm{cm^{-3}}, ni=1.5×1010 cm3n_i = 1.5\times10^{10}\ \mathrm{cm^{-3}}. Give the numeric answer in volts. [5]


Question 3 — Depletion region & space charge [10 marks]

(a) Explain why a depletion region forms and why it is "depleted" of mobile carriers. Describe the sign of the space charge on each side. [3]

(b) Using charge neutrality, derive the relation between depletion widths xnx_n, xpx_p and the dopings, and hence show which side of an asymmetric junction (NANDN_A \gg N_D) the depletion region extends into. [4]

(c) State the depletion width formula and compute the total width WW for the junction of Q2(b) with VbiV_{bi} from your answer, using εs=1.04×1012 F/cm\varepsilon_s = 1.04\times10^{-12}\ \mathrm{F/cm}. Formula: W=2εsq(1NA+1ND)VbiW = \sqrt{\frac{2\varepsilon_s}{q}\left(\frac{1}{N_A}+\frac{1}{N_D}\right)V_{bi}} [3]


Question 4 — Shockley equation & I–V (code-from-memory) [12 marks]

(a) State the Shockley diode equation and define every symbol. [3]

(b) A diode has IS=1×1014 AI_S = 1\times10^{-14}\ \mathrm{A} at T=300 KT=300\ \mathrm{K}, ideality factor n=1n=1. Compute the current at V=0.6 VV = 0.6\ \mathrm{V}. [4]

(c) Write a short Python function (from memory, using numpy) diode_current(V, Is, n, T) that returns the Shockley current for an array of voltages, and describe what the I–V curve looks like in forward vs reverse. [5]


Question 5 — Reverse behaviour & breakdown [10 marks]

(a) Explain reverse saturation current ISI_S: what carriers cause it, and why is it nearly independent of reverse voltage but strongly temperature dependent? [4]

(b) Contrast avalanche and Zener breakdown: mechanism, doping dependence, and the sign of the temperature coefficient of breakdown voltage for each. [6]


Question 6 — Junction capacitance [6 marks]

(a) Name the two capacitance contributions of a diode and state which dominates in reverse bias vs forward bias. [3]

(b) The depletion capacitance is Cj=εsAWC_j = \dfrac{\varepsilon_s A}{W}. For A=104 cm2A = 10^{-4}\ \mathrm{cm^2} and WW from Q3(c), compute CjC_j in pF. [3]


Answer keyMark scheme & solutions

Q1 [10]

(a) [4] P is Group V (5 valence e⁻), Si Group IV (4 valence e⁻) [1]. P forms 4 covalent bonds with neighbouring Si atoms using 4 of its electrons [1]; the 5th electron is not needed for bonding [1] and is weakly bound (donor binding energy ~0.045 eV), so at room temperature it is thermally ionised into the conduction band as a nearly-free electron → n-type [1].

(b) [3] B is Group III (3 valence e⁻) [1]. It forms only 3 bonds, leaving one bond incomplete — an acceptor state [1]. An electron from a neighbouring bond hops in to complete it, leaving behind a mobile hole; B becomes a fixed negative ion → p-type [1].

(c) [3] Donor level EDE_D lies just below EcE_c (≈0.045 eV below) [1]; acceptor level EAE_A lies just above EvE_v [1]. "Shallow" because the ionisation energy (\simtens of meV) is comparable to kBT26k_BT\approx26 meV at 300 K, so nearly all dopants are ionised at room temperature [1].


Q2 [12]

(a) [7]

  • At equilibrium, net electron current = 0: drift balances diffusion [1]: Jn=qμnnE+qDndndx=0J_n = q\mu_n n E + qD_n\frac{dn}{dx} = 0
  • Use Einstein relation Dn/μn=VTD_n/\mu_n = V_T and E=dψ/dxE=-d\psi/dx [1]: dψdx=VT1ndndx\frac{d\psi}{dx} = V_T\frac{1}{n}\frac{dn}{dx}
  • Integrate across junction from p-side to n-side [1]: Vbi=ψnψp=VTlnnnnpV_{bi} = \psi_n-\psi_p = V_T\ln\frac{n_n}{n_p}
  • On n-side nnNDn_n \approx N_D; on p-side (minority) np=ni2/NAn_p = n_i^2/N_A by mass-action law [2]:
  • Substitute [2]: Vbi=VTlnNDNAni2V_{bi} = V_T\ln\frac{N_D N_A}{n_i^2}

(b) [5] VT=kT/q=0.02585V_T = kT/q = 0.02585 V [1]. Argument =10171016(1.5×1010)2=10332.25×1020=4.44×1012= \dfrac{10^{17}\cdot10^{16}}{(1.5\times10^{10})^2} = \dfrac{10^{33}}{2.25\times10^{20}} = 4.44\times10^{12} [2]. ln(4.44×1012)=29.12\ln(4.44\times10^{12}) = 29.12 [1]. Vbi=0.02585×29.12=0.753 VV_{bi} = 0.02585\times29.12 = 0.753\ \mathrm{V} [1].


Q3 [10]

(a) [3] Near the metallurgical junction, electrons diffuse from n→p and holes from p→n, recombining and leaving behind uncovered fixed ionised dopants [1]. This region has no mobile carriers → depleted [1]. n-side: positive space charge (ionised donors); p-side: negative space charge (ionised acceptors) [1].

(b) [4] Overall charge neutrality requires total positive = total negative charge [1]: qNDxn=qNAxpNDxn=NAxpq N_D x_n = q N_A x_p \Rightarrow N_D x_n = N_A x_p [2] If NANDN_A \gg N_D then xnxpx_n \gg x_p: the depletion region extends mostly into the lightly-doped n-side [1].

(c) [3] (1NA+1ND)=1017+1016=1.1×1016 cm3\left(\frac1{N_A}+\frac1{N_D}\right) = 10^{-17}+10^{-16} = 1.1\times10^{-16}\ \mathrm{cm^3} [1]. W=2(1.04×1012)1.602×1019(1.1×1016)(0.753)W = \sqrt{\frac{2(1.04\times10^{-12})}{1.602\times10^{-19}}(1.1\times10^{-16})(0.753)} Inner =1.298×107×1.1×1016×0.753=1.075×109= 1.298\times10^7 \times 1.1\times10^{-16}\times0.753 = 1.075\times10^{-9} [1]. W=1.075×109=3.28×105 cm=0.328 μmW = \sqrt{1.075\times10^{-9}} = 3.28\times10^{-5}\ \mathrm{cm} = 0.328\ \mu\mathrm{m} [1].


Q4 [12]

(a) [3] I=IS ⁣(eV/(nVT)1)I = I_S\!\left(e^{V/(nV_T)}-1\right). ISI_S = reverse saturation current, VV = applied voltage, nn = ideality factor, VT=kBT/qV_T = k_BT/q thermal voltage. [3] (1 for eqn, 2 for symbols)

(b) [4] VT=0.02585V_T=0.02585 V. V/VT=0.6/0.02585=23.21V/V_T = 0.6/0.02585 = 23.21 [1]. e23.21=1.20×1010e^{23.21}=1.20\times10^{10} [1]. I=1014(1.20×10101)1.20×104 A=0.12 mAI = 10^{-14}(1.20\times10^{10}-1) \approx 1.20\times10^{-4}\ \mathrm{A} = 0.12\ \mathrm{mA} [2].

(c) [5]

import numpy as np
def diode_current(V, Is, n, T):
    kB = 1.381e-23; q = 1.602e-19
    VT = kB*T/q
    return Is*(np.exp(V/(n*VT)) - 1.0)

[3] for correct function. Forward: exponential rise, current ~0 below ~0.6 V then steep knee [1]. Reverse: current saturates at IS-I_S (tiny, ~constant) until breakdown [1].


Q5 [10]

(a) [4] ISI_S is carried by minority carriers (thermally generated electrons in p-region, holes in n-region) swept across the junction by the field [2]. Nearly voltage-independent because it's limited by the rate of thermal generation, not the field, once carriers are collected [1]. Strongly T-dependent because generation ni2eEg/kT\propto n_i^2 \propto e^{-E_g/kT} — roughly doubles every ~10 °C [1].

(b) [6]

Avalanche Zener
Mechanism Impact ionisation: high-field carriers knock out further e-h pairs (multiplication) [1] Direct band-to-band tunnelling across thin junction [1]
Doping Lighter doping, wider depletion, higher VBRV_{BR} (>~6 V) [1] Heavy doping, very thin junction, low VBRV_{BR} (<~5 V) [1]
Temp. coeff. Positive (VBRV_{BR} rises with T — phonon scattering shortens mean free path) [1] Negative (VBRV_{BR} falls with T — band gap shrinks) [1]

Q6 [6]

(a) [3] Depletion (junction) capacitance and diffusion capacitance [1]. Depletion capacitance dominates in reverse bias [1]; diffusion capacitance dominates in forward bias (stored minority charge) [1].

(b) [3] Cj=εsAW=1.04×1012×1043.28×105C_j = \dfrac{\varepsilon_s A}{W} = \dfrac{1.04\times10^{-12}\times10^{-4}}{3.28\times10^{-5}} [1] =1.04×10163.28×105=3.17×1012 F= \dfrac{1.04\times10^{-16}}{3.28\times10^{-5}} = 3.17\times10^{-12}\ \mathrm{F} [1] =3.17 pF= 3.17\ \mathrm{pF} [1].


[
  {"claim":"V_bi = 0.753 V for NA=1e17,ND=1e16,ni=1.5e10 at 300K",
   "code":"kB=1.381e-23;q=1.602e-19;T=300;VT=kB*T/q;NA=1e17;ND=1e16;ni=1.5e10;Vbi=VT*log(NA*ND/ni**2);result=abs(float(Vbi)-0.753)<0.01"},
  {"claim":"Depletion width W approx 3.28e-5 cm",
   "code":"eps=1.04e-12;q=1.602e-19;NA=1e17;ND=1e16;Vbi=0.753;W=sqrt(2*eps/q*(1/NA+1/ND)*Vbi);result=abs(float(W)-3.28e-5)<0.1e-5"},
  {"claim":"Diode current at 0.6V, Is=1e-14, n=1, T=300 approx 1.2e-4 A",
   "code":"kB=1.381e-23;q=1.602e-19;T=300;VT=kB*T/q;Is=1e-14;I=Is*(exp(0.6/VT)-1);result=abs(float(I)-1.2e-4)/1.2e-4<0.05"},
  {"claim":"Junction capacitance approx 3.17 pF",
   "code":"eps=1.04e-12;A=1e-4;W=3.28e-5;Cj=eps*A/W;result=abs(float(Cj)*1e12-3.17)<0.1"}
]