Intuition What this page is for
The parent note gave you two formulas:
the exact hole count p = 2 1 ( N A + N A 2 + 4 n i 2 ) and the shortcut p ≈ N A .
A shortcut is only trustworthy if you know exactly when it fails . This page walks through every situation a doping problem can hand you — heavy doping, light doping, no doping, freeze-out, temperature changes, compensation, and a real device question — so you never meet a case you haven't already seen solved.
Before we start, let us re-anchor every symbol so no notation is used before it is earned.
Definition The symbols, in plain words
N A = number of acceptor (boron) atoms per cubic centimetre we deliberately added. Units: cm − 3 (things counted inside one cubic centimetre).
N D = number of donor (phosphorus) atoms per cm³ . A donor is the mirror of an acceptor — it gives away a spare electron (see N-type doping with donor atoms (phosphorus) ). For pure p-type there are none, so N D = 0 ; it only matters in the compensation case (Cell H).
n i = intrinsic carrier concentration : how many free electrons (and equally, holes) pure undoped silicon makes on its own from heat alone. In silicon at room temperature n i ≈ 1.5 × 1 0 10 cm − 3 . See Intrinsic Semiconductors .
p = how many free holes per cm³ actually exist after doping.
n = how many free electrons per cm³ actually exist after doping.
q = the charge of one electron in magnitude, q = 1.6 × 1 0 − 19 coulombs (C). We will only need q once, in Cell G , where it converts a carrier count into an actual electric current — until then every equation is pure counting and q does not appear.
μ p = hole mobility , units cm 2 / ( V⋅s ) — how fast (drift speed, cm/s) a hole moves per unit of electric push (field, V/cm). Bigger μ p = faster holes. In silicon μ p ≈ 480 cm 2 / ( V⋅s ) .
μ n = electron mobility , same units cm 2 / ( V⋅s ) . In silicon electrons are actually faster than holes (μ n ≈ 1350 cm 2 / ( V⋅s ) vs μ p ≈ 480 ), which we will need in Cell G when we argue why minority electrons still carry negligible current.
Dimensional check for later: in σ = q p μ p the units multiply as C × cm − 3 × cm 2 / ( V⋅s ) = C / ( V⋅s⋅cm ) = A / ( V⋅cm ) = ( Ω ⋅cm ) − 1 , using C/s = A and A/V = Ω − 1 . So conductivity comes out in ( Ω ⋅cm ) − 1 automatically — good to keep in mind before we ever plug numbers.
The two laws we lean on the whole way (both from the parent note):
Every doping problem is really one of these cells . The whole point of a "deep dive" is to fill every box, not just the easy one.
Cell
Regime
What controls p ?
Which formula?
Example
A
Heavy doping N A ≫ n i
the dopant
p ≈ N A
Ex 1
B
Comparable N A ∼ n i
dopant and heat
full quadratic
Ex 2
C
Zero doping N A = 0
pure heat
p = n i
Ex 3
D
Trace doping N A ≪ n i
mostly heat, tiny tilt
full quadratic + series expansion
Ex 4
E
Temperature raised (hot)
n i grows, may overtake N A
recompute n i , then quadratic
Ex 5
F
Temperature lowered (freeze-out)
only some acceptors ionized
N A − < N A
Ex 6
G
Real device: resistance/current
conductivity from p
σ = q p μ p
Ex 7
H
Compensation, net p-type (N A > N D )
N A − N D
net-doping neutrality
Ex 8
I
Compensation, net n-type (N D > N A )
N D − N A
net-doping neutrality
Ex 9
J
Extreme doping (degenerate)
statistics break down
formulas above invalid
note in Cell J
The figure below is the map we will keep returning to. Figure 1 plots, on log–log axes, the hole count p (purple) and electron count n (coral) as we sweep N A (horizontal axis, cm − 3 ) from far below n i to far above it; both vertical and horizontal axes are carrier concentrations in cm − 3 . The green dashed line marks the n i level, and the shaded bands + text labels tell you exactly which cell and which example each region belongs to.
Figure 1 — Carrier concentrations p and n (both cm⁻³, vertical) versus acceptor doping N A (cm⁻³, horizontal), log–log. Mint band = Cells C/D (below n i ); butter band = Cell B bend (N A ∼ n i ); lavender band = Cell A (N A ≫ n i ).
Worked example Example 1 — Cell A
Silicon doped with N A = 1 0 16 cm − 3 boron, n i = 1.5 × 1 0 10 cm − 3 . Find p and n . (This is the far-right end of Figure 1, where the purple line sits flat on N A .)
Forecast: Will p be close to 1 0 16 , or noticeably less? Guess before reading.
Compare sizes: N A / n i = 1 0 16 / ( 1.5 × 1 0 10 ) ≈ 6.7 × 1 0 5 .
Why this step? The whole decision "shortcut vs quadratic" rests on whether N A ≫ n i . Here it is 600,000× bigger — deep in Cell A .
So the N A 2 + 4 n i 2 term is essentially N A 2 = N A , giving p ≈ N A = 1 0 16 cm − 3 .
Why this step? 4 n i 2 is ∼ 1 0 − 11 of N A 2 — invisible under the square root.
n = n i 2 / p = ( 1.5 × 1 0 10 ) 2 /1 0 16 = 2.25 × 1 0 20 /1 0 16 = 2.25 × 1 0 4 cm − 3 .
Why this step? Once p is fixed, mass-action pins n automatically.
Verify: Plug into neutrality: n + N A = 2.25 × 1 0 4 + 1 0 16 ≈ 1 0 16 = p . ✓ The electron term is a rounding speck. Units: cm − 3 throughout. ✓
Worked example Example 2 — Cell B
N A = 2 × 1 0 10 cm − 3 , so N A is barely above n i = 1.5 × 1 0 10 . Find p and n . (This is the bend in Figure 1, where purple and coral start to separate.)
Forecast: The shortcut says p = 2 × 1 0 10 . Do you trust it here? By how much might it be off?
Do not drop n i . Compute the discriminant: N A 2 + 4 n i 2 = ( 2 × 1 0 10 ) 2 + 4 ( 1.5 × 1 0 10 ) 2 = 4 × 1 0 20 + 9 × 1 0 20 = 1.3 × 1 0 21 .
Why this step? 4 n i 2 is now bigger than N A 2 — the term we would have thrown away dominates, so we keep it.
1.3 × 1 0 21 = 3.606 × 1 0 10 .
Why this step? The exact formula needs the square root of that discriminant; taking it turns the "sum of squares" back into a quantity with the same units (cm − 3 ) as N A , so it can be added to N A in the next step.
p = ( 2 × 1 0 10 + 3.606 × 1 0 10 ) /2 = 2.803 × 1 0 10 cm − 3 .
Why this step? This is literally the quadratic-formula answer p = 2 1 ( N A + … ) — we average N A with the root we just found to get the one physical (positive) hole count.
n = n i 2 / p = 2.25 × 1 0 20 /2.803 × 1 0 10 = 8.03 × 1 0 9 cm − 3 .
Why this step? Mass-action fixes n once p is known.
Verify: neutrality check n + N A = 8.03 × 1 0 9 + 2 × 1 0 10 = 2.803 × 1 0 10 = p . ✓
The shortcut p ≈ N A = 2 × 1 0 10 would be 28% too low — real lesson of Cell B.
Worked example Example 3 — Cell C
Set N A = 0 . What do the formulas give, and is it consistent with Intrinsic Semiconductors ? (This is the far-left of Figure 1, where purple and coral merge onto the green n i line.)
Forecast: With no boron at all, should p equal n i , zero, or something in between?
Put N A = 0 into the exact formula: p = 2 1 ( 0 + 0 + 4 n i 2 ) = 2 1 ( 2 n i ) = n i .
Why this step? This is the sanity boundary — a good formula must reduce to the undoped answer when doping vanishes.
n = n i 2 / p = n i 2 / n i = n i .
Why this step? Pure silicon must have equal electrons and holes.
Verify: p = n = n i = 1.5 × 1 0 10 cm − 3 . Neutrality: p = n + 0 . ✓ The p-type formula gracefully becomes the intrinsic result — no discontinuity.
Before crunching numbers, let us do the promised analytical expansion so we understand the shape of the answer when doping is tiny.
Worked example Example 4 — Cell D
N A = 3 × 1 0 9 cm − 3 , five times below n i = 1.5 × 1 0 10 . Find p and n both ways (exact and the expansion above), and show the material is barely p-type. (Figure 1: just right of the far-left merge — purple is a hair above coral.)
Forecast: Is p closer to N A (tiny) or to n i (big)? Which carrier wins?
Exact route — discriminant: N A 2 + 4 n i 2 = ( 3 × 1 0 9 ) 2 + 4 ( 1.5 × 1 0 10 ) 2 = 9 × 1 0 18 + 9 × 1 0 20 = 9.09 × 1 0 20 .
Why this step? N A 2 is 1% of 4 n i 2 — the crystal's own heat dominates, so we must keep the n i term.
9.09 × 1 0 20 = 3.015 × 1 0 10 , so p = ( 3 × 1 0 9 + 3.015 × 1 0 10 ) /2 = 1.657 × 1 0 10 cm − 3 .
Why this step? Direct quadratic-formula answer p = 2 1 ( N A + … ) .
Expansion route — plug into p ≈ n i + N A /2 + N A 2 / ( 8 n i ) : = 1.5 × 1 0 10 + 1.5 × 1 0 9 + ( 9 × 1 0 18 ) / ( 1.2 × 1 0 11 ) = 1.5 × 1 0 10 + 0.15 × 1 0 10 + 7.5 × 1 0 7 = 1.6575 × 1 0 10 cm − 3 .
Why this step? This shows the series and the exact formula agree to 4 figures — confirming the expansion is trustworthy in Cell D and giving physical insight (the N A /2 tilt).
n = n i 2 / p = 2.25 × 1 0 20 /1.657 × 1 0 10 = 1.357 × 1 0 10 cm − 3 .
Why this step? Mass-action.
Verify: n + N A = 1.357 × 1 0 10 + 3 × 1 0 9 = 1.657 × 1 0 10 = p . ✓ Exact and expansion match to < 0.1% . ✓
Notice p / n = 1.22 — holes only slightly outnumber electrons. As N A → 0 , p / n → 1 (back to intrinsic, Cell C). The doping merely tilts the balance.
The intrinsic concentration is not a constant — it rises steeply with temperature because heat rips more bonds. If you heat p-type silicon enough, n i can catch up to N A , dragging you from Cell A back toward Cell B. This is why power devices lose their doping character when they overheat. Figure 2 plots n i climbing (coral) toward the fixed N A line (purple) against temperature — the yellow crossing marker is exactly the temperature used in Example 5.
Figure 2 — Temperature T (K, horizontal) vs concentration (cm⁻³, vertical log axis). Coral = n i rising with heat; lavender dashed = fixed N A = 1 0 15 cm − 3 ; butter dot = the Example 5 operating point where n i has grown to half of N A .
Worked example Example 5 — Cell E
A wafer with N A = 1 0 15 cm − 3 is heated until n i rises to 5 × 1 0 14 cm − 3 . Find p and n , and compare to the room-temperature shortcut p = N A .
Forecast: At this temperature n i is half of N A . Will p still be near N A ?
Discriminant: N A 2 + 4 n i 2 = ( 1 0 15 ) 2 + 4 ( 5 × 1 0 14 ) 2 = 1 × 1 0 30 + 1 × 1 0 30 = 2 × 1 0 30 .
Why this step? Now 4 n i 2 equals N A 2 — we are firmly in the "must use quadratic" zone (Cell B territory, reached by heating instead of by light doping).
2 × 1 0 30 = 1.4142 × 1 0 15 .
Why this step? Take the root of the discriminant so it carries units cm − 3 and can be added to N A .
p = ( 1 0 15 + 1.4142 × 1 0 15 ) /2 = 1.2071 × 1 0 15 cm − 3 .
Why this step? Quadratic-formula answer p = 2 1 ( N A + … ) .
n = n i 2 / p = 2.5 × 1 0 29 /1.2071 × 1 0 15 = 2.071 × 1 0 14 cm − 3 .
Why this step? Mass-action fixes n from the p we just computed; note n is now only ~6× below p , showing the material is drifting toward intrinsic.
Verify: n + N A = 2.071 × 1 0 14 + 1 0 15 = 1.2071 × 1 0 15 = p . ✓
The shortcut p = N A = 1 0 15 is now 21% low . Heat pushes any doped crystal toward intrinsic.
At room temperature boron's acceptor level sits just above the valence band, so essentially all acceptors have grabbed an electron and are ionized: N A − ≈ N A . Cool the crystal and there isn't enough thermal energy to complete that hop — some boron atoms sit un-ionized , holding onto their empty bond without releasing a hole. Then N A − < N A and p < N A . This is why the parent note warns "p = N A breaks at low T."
Worked example Example 6 — Cell F
N A = 1 0 16 cm − 3 , but at a cold temperature only 40% of acceptors are ionized, and n i has dropped to a negligible 1 × 1 0 6 cm − 3 . Find p and n .
Forecast: Will p be 1 0 16 , or the fraction 0.4 × 1 0 16 , or something else?
Ionized acceptors: N A − = 0.40 × 1 0 16 = 4 × 1 0 15 cm − 3 .
Why this step? Only ionized boron has released a hole; neutral boron contributes none. Charge neutrality uses N A − , not N A .
Since n i = 1 0 6 ≪ N A − , we are deep in Cell A relative to the ionized count: p ≈ N A − = 4 × 1 0 15 cm − 3 .
Why this step? The quadratic with N A − still applies; n i is utterly tiny so the shortcut is fine.
n = n i 2 / p = ( 1 0 6 ) 2 / ( 4 × 1 0 15 ) = 1 0 12 / ( 4 × 1 0 15 ) = 2.5 × 1 0 − 4 cm − 3 .
Why this step? Mass-action fixes n from p ; the answer being far below 1 per cm³ just tells us electrons are essentially absent at this cold temperature.
Verify: p = 4 × 1 0 15 , well below the full N A = 1 0 16 — freeze-out has hidden 60% of the holes. As temperature rises back to room level, N A − → N A and p climbs back to 1 0 16 . ✓
This is the cell where the charge q (defined in the symbols block) finally does work: it turns a carrier count (cm − 3 ) into an electric current . Recall the dimensional check from the symbols block — σ = q p μ p comes out in ( Ω ⋅cm ) − 1 .
Worked example Example 7 — Cell G
A bar of p-type silicon: length L = 2 mm = 0.2 cm , cross-section A = 0.01 cm 2 , doping N A = 1 0 16 cm − 3 (so p ≈ N A from Cell A), hole mobility μ p = 480 cm 2 / ( V⋅s ) , q = 1.6 × 1 0 − 19 C . Find the conductivity σ and the resistance R .
Forecast: Which carrier carries the current — holes or the minority electrons? Guess the order of magnitude of R .
Conductivity from holes: σ = q p μ p = ( 1.6 × 1 0 − 19 ) ( 1 0 16 ) ( 480 ) = 0.768 ( Ω ⋅cm ) − 1 .
Why this step? Current in p-type is holes drifting (see Drift and Diffusion Currents ). Minority electrons (n ≈ 2.25 × 1 0 4 ) would add q n μ n = ( 1.6 × 1 0 − 19 ) ( 2.25 × 1 0 4 ) ( 1350 ) ∼ 5 × 1 0 − 12 — about 1 0 12 × smaller, so we drop it even though μ n > μ p .
Resistivity ρ = 1/ σ = 1/0.768 = 1.302 Ω ⋅cm .
Why this step? Resistivity is just the reciprocal of conductivity.
Resistance R = ρ L / A = 1.302 × 0.2/0.01 = 26.04 Ω .
Why this step? Standard geometry law: longer bar = more resistance, fatter bar = less.
Verify: Units: ( Ω cm ) ⋅ ( cm ) / ( cm 2 ) = Ω . ✓ Recompute σ backward: q p μ p = 1.6 × 1 0 − 19 ⋅ 1 0 16 ⋅ 480 = 0.768 . ✓
Sometimes a wafer already has some phosphorus donors (N-type doping with donor atoms (phosphorus) ) when we add boron. The donors' free electrons cancel an equal number of holes. Only the excess dopant survives. Neutrality is the full p + N D = n + N A , so the effective doping is N A − N D .
Worked example Example 8 — Cell H
A crystal has N A = 8 × 1 0 15 cm − 3 boron and N D = 3 × 1 0 15 cm − 3 phosphorus, n i = 1.5 × 1 0 10 . Is it p-type or n-type? Find p and n .
Forecast: Which dopant wins? Will the leftover be big enough for the shortcut?
Net acceptor doping: N A − N D = 8 × 1 0 15 − 3 × 1 0 15 = 5 × 1 0 15 cm − 3 (positive → still p-type ).
Why this step? Each donor electron fills one acceptor hole; only the surplus of holes remains mobile.
Since 5 × 1 0 15 ≫ n i , Cell A shortcut applies to the net : p ≈ N A − N D = 5 × 1 0 15 cm − 3 .
Why this step? The net doping plays exactly the role N A did before.
n = n i 2 / p = ( 1.5 × 1 0 10 ) 2 / ( 5 × 1 0 15 ) = 2.25 × 1 0 20 / ( 5 × 1 0 15 ) = 4.5 × 1 0 4 cm − 3 .
Why this step? Mass-action fixes n from p ; because p is large, n is heavily suppressed — holes are still the majority.
Verify: neutrality p + N D = n + N A ⇒ 5 × 1 0 15 + 3 × 1 0 15 = 8 × 1 0 15 , and n + N A = 4.5 × 1 0 4 + 8 × 1 0 15 ≈ 8 × 1 0 15 . ✓ Both sides = 8 × 1 0 15 .
If the donors outnumber the acceptors, the surplus is now electrons, and the crystal is n-type despite containing boron. The same neutrality equation handles it — you just find the sign of N A − N D flips, so you solve for n (majority) instead of p .
Worked example Example 9 — Cell I
A crystal has N A = 2 × 1 0 15 cm − 3 boron and N D = 7 × 1 0 15 cm − 3 phosphorus, n i = 1.5 × 1 0 10 . Find the majority carrier, n and p .
Forecast: With more donors than acceptors, which carrier now dominates — and by roughly how much?
Net donor doping: N D − N A = 7 × 1 0 15 − 2 × 1 0 15 = 5 × 1 0 15 cm − 3 (positive donor surplus → n-type ).
Why this step? Each acceptor hole is filled by a donor electron; the leftover 5 × 1 0 15 electrons are the mobile majority.
Since 5 × 1 0 15 ≫ n i , the n-type shortcut applies: n ≈ N D − N A = 5 × 1 0 15 cm − 3 .
Why this step? This is the mirror of Cell A with electrons as majority — see N-type doping with donor atoms (phosphorus) .
p = n i 2 / n = ( 1.5 × 1 0 10 ) 2 / ( 5 × 1 0 15 ) = 2.25 × 1 0 20 / ( 5 × 1 0 15 ) = 4.5 × 1 0 4 cm − 3 .
Why this step? Mass-action fixes the minority holes from the majority n .
Verify: neutrality p + N D = n + N A ⇒ 4.5 × 1 0 4 + 7 × 1 0 15 ≈ 7 × 1 0 15 and n + N A = 5 × 1 0 15 + 2 × 1 0 15 = 7 × 1 0 15 . ✓ Both sides = 7 × 1 0 15 . Adding boron did not make it p-type — the donor surplus wins.
Everything above uses the non-degenerate mass-action law n p = n i 2 . That law quietly assumes carriers are dilute enough to obey Boltzmann statistics . When doping gets very high, that assumption fails and we need a quantitative flag for when .
Definition The quantitative degeneracy criterion
Let N V = the effective density of states in the valence band (silicon: N V ≈ 1.04 × 1 0 19 cm − 3 at room temperature) — roughly "how many hole-seats the valence band edge offers." The non-degenerate approximation is safe only while
p ≲ 0.1 N V ⟺ E V − E F ≳ 3 k B T
where E F is the Fermi level, E V the valence band edge, k B Boltzmann's constant and T temperature (this energy language is developed in Fermi Level in Doped Semiconductors ). Once p approaches N V — i.e. N A ≳ 1 0 19 cm − 3 in silicon — the Fermi level enters the valence band, the crystal is degenerate , and it behaves almost like a metal.
Worked example Example 10 — Cell J (recognise, don't misuse)
A wafer is doped N A = 5 × 1 0 19 cm − 3 . Someone plugs it into p ≈ N A then n = n i 2 / p . Why is this wrong , and what should be used instead?
Forecast: Is p really 5 × 1 0 19 , larger, or smaller than that?
Test the criterion: N A / N V = 5 × 1 0 19 /1.04 × 1 0 19 ≈ 4.8 > 0.1 .
Why this step? This is the single check that decides "legal vs illegal formula." Here N A is nearly 5× the band's capacity — deeply degenerate.
So n p = n i 2 is invalid ; you cannot get p from the simple formula. Instead p must be found from the Fermi–Dirac integral
p = N V F 1/2 ( k B T E V − E F )
where F 1/2 is the Fermi–Dirac integral of order 2 1 (it replaces the simple exponential of Boltzmann statistics). Not every acceptor even releases a free hole, so p < N A here.
Why this step? At these densities the holes fill a finite number of valence-band seats; Fermi–Dirac counting (from Mass-Action Law and Carrier Statistics and Fermi Level in Doped Semiconductors ) is required.
Verify (sanity, not a plug-in): the warning sign N A > 0.1 N V is what flags degeneracy. Numerically 5 × 1 0 19 > 0.1 × 1.04 × 1 0 19 = 1.04 × 1 0 18 . ✓ So the honest answer is "degenerate — non-degenerate formulas do not apply; use p = N V F 1/2 ( … ) ."
Recall Which cell am I in? (decision flow)
Answer these in order.
First: is N D present? ::: If yes, replace N A with N A − N D (Cell H if positive → p-type, Cell I if negative → n-type).
Next: is the net doping near 0.1 N V ≈ 1 0 18 –1 0 19 cm − 3 or above? ::: Then it is degenerate (Cell J) — the n p = n i 2 formulas do not apply; use p = N V F 1/2 ( … ) .
Then: how does the net doping compare to n i at this temperature? ::: ≫ n i → shortcut p ≈ N A (Cell A). Comparable or below → full quadratic or its series expansion (Cells B, C, D).
Cold crystal? ::: Use ionized N A − , not full N A (Cell F).
Hot crystal? ::: Recompute the larger n i first, then decide (Cell E).
When is the shortcut p ≈ N A safe? Only when N A ≫ n i and all acceptors ionized and below degeneracy (N A ≪ N V ) — Cell A.
In Example 2 (N A = 2 × 1 0 10 ) why is p = 2.8 × 1 0 10 not 2 × 1 0 10 ? Because n i ∼ N A , so 4 n i 2 dominates the discriminant and cannot be dropped, pushing p about 28% above the naive N A .
What does the exact formula give when N A = 0 ? p = n = n i — it reduces to the intrinsic result.
Series expansion of p for N A ≪ n i ? p ≈ n i + N A /2 + N A 2 / ( 8 n i ) — leading intrinsic term, then a half-doping tilt.
Why does heating p-type silicon reduce its "p-ness"? n i rises steeply, catches up to N A , and the carrier ratio p / n falls toward 1.
What is freeze-out? At low temperature only a fraction of acceptors ionize, so N A − < N A and p < N A .
For a compensated crystal, what sets the majority carrier? The net dopant N A − N D : positive → p-type (Cell H), negative → n-type (Cell I).
Quantitative test for when n p = n i 2 breaks down? Degeneracy when p ≳ 0.1 N V (i.e. N A ≳ 1 0 19 cm − 3 in Si); then use the Fermi–Dirac integral p = N V F 1/2 (( E V − E F ) / k B T ) (Cell J).
Conductivity of p-type silicon formula, units, and why holes only? σ = q p μ p in ( Ω ⋅cm ) − 1 ; minority electrons contribute q n μ n which is ~1 0 12 × smaller even though μ n > μ p .