Silicon Processing, Powder Preparation and Sintering
Segregation and zone refining
When a melt freezes, the solid generally holds a different solute concentration from the liquid. The equilibrium segregation coefficient is k = C_s/C_l, the ratio of solute concentration in the solid to that in the liquid at the interface. For a solute that lowers the melting point k < 1: the solute is rejected into the liquid. For one that raises it k > 1. Most impurities in silicon have k ≪ 1, which is why freezing purifies it. Normal freezing of a bar of initial solute concentration C₀, with complete mixing in the liquid and no diffusion in the solid, gives at fraction solidified f: C_s = k C₀ (1 − f)^(k−1). The same law governs a Czochralski melt: for k < 1 the melt grows richer as the crystal is pulled, so the dopant concentration in the crystal rises along its length.
Zone refining (Pfann) passes a narrow molten zone of length l along a bar. At the freezing end the solid takes k C_l and at the melting end the zone picks up solute from the bar, so for k < 1 solute is swept to the far end of the bar. For the first pass, with starting concentration C₀ uniform, the solid concentration at distance x from the start (valid until the last zone length) is C(x) = C₀ [1 − (1 − k) exp(−k x/l)]. At x = 0 this gives k C₀ (the first solid is purest); it climbs towards C₀ as x grows. Many passes, with a new zone started at the beginning each time, push the purified region further and finally yield an ultra-pure bar. Worked: k = 0.1 and x/l = 5 give C/C₀ = 1 − 0.9 e^(−0.5) = 1 − 0.9 × 0.6065 = 0.454.
Single-crystal growth: Czochralski and float zone
In the Czochralski (CZ) method polysilicon and the desired dopant are melted in a fused-silica crucible (silicon melts at 1414 °C) under inert argon. A small seed crystal of the chosen orientation is dipped into the melt and pulled upward while seed and crucible rotate in opposite senses. The melt freezes on the seed and copies its orientation. A narrow neck is grown first, so that dislocations from thermal shock grow out; then the diameter is widened to the shoulder and held constant in the body by adjusting the pull rate and the heater power. The ingot is then ground, sliced into wafers, lapped, etched and polished.
The quartz crucible dissolves slightly, so CZ silicon contains oxygen (of order 10¹⁸ atoms/cm³), which is not always harmful: under heat treatment it can precipitate and trap metal impurities (gettering). In the float-zone (FZ) method a polysilicon rod is held vertically; a radio-frequency induction coil melts a short zone, touching a seed at the bottom, and the zone is moved up the rod. The melt is held by surface tension, so there is no crucible and no oxygen pick-up; zone refining happens during growth. FZ silicon therefore has higher purity and resistivity, and suits power devices and detectors, but the diameter is limited and the cost is higher. CZ supplies most wafers, because it gives large diameters, is cheaper and allows uniform doping from the melt.
| Feature | Czochralski | Float zone |
|---|---|---|
| Container | Fused-silica crucible | None; zone held by surface tension |
| Oxygen content | High | Very low |
| Diameter and cost | Large, lower cost | Smaller, higher cost |
| Dopant uniformity along ingot | Varies with k, as in normal freezing | Also affected by segregation; doping is by gas phase or by neutron irradiation (transmutation), which gives very uniform n-type material |
Silicon oxidation: the Deal–Grove model
Silicon oxidises to SiO₂ in dry O₂ (Si + O₂ → SiO₂) or in steam (Si + 2H₂O → SiO₂ + 2H₂), typically at 900 to 1200 °C. The oxide is the gate insulator, the mask against dopants and the isolation layer. Oxidation happens at the Si–SiO₂ interface: the oxidant must diffuse through the existing oxide, so the oxide thickens more slowly as it grows. The Deal–Grove model combines oxidant diffusion through the oxide with the interface reaction and gives, for oxide thickness x after time t, x² + A x = B (t + τ), where τ accounts for an initial oxide. Solving, x = (A/2)[√(1 + 4B(t + τ)/A²) − 1].
- Thin oxide, short time (t + τ ≪ A²/4B): linear regime. The interface reaction is the bottleneck and x ≈ (B/A)(t + τ). B/A is the linear rate constant.
- Thick oxide, long time (t ≫ A²/4B): parabolic regime. Diffusion through the oxide controls and x ≈ √(B t). B is the parabolic rate constant; doubling the time raises x only by √2.
Worked example (A = 0.20 µm, B = 0.050 µm²/h, τ = 0, t = 4 h; constants assumed for illustration): x² + 0.20x − 0.20 = 0, so x = [−0.20 + √(0.04 + 0.80)]/2 = (−0.20 + 0.9165)/2 = 0.358 µm. The parabolic law alone would give √(0.050 × 4) = 0.447 µm, an over-estimate of 25 %, because the oxide is not yet thick enough for the diffusion limit. Silicon is consumed: since SiO₂ occupies about 2.27 times the volume of the silicon it replaces, a thickness x of oxide uses up 0.44 x of silicon, here 0.158 µm. Wet oxidation is faster than dry because H₂O is far more soluble in SiO₂ than O₂ (larger B); dry oxide is denser and is chosen for thin gate oxides.
Doping and the photolithographic process
Doping by diffusion uses D = D₀ exp(−E_a/kT), so depth is set by √(Dt) and rises steeply with temperature. With a constant surface concentration C_s (predeposition from a gas or glass source) the profile is C(x, t) = C_s erfc[x/(2√(Dt))]. The junction lies where C equals the background concentration C_b of the opposite type, so x_j = 2z√(Dt) with erfc(z) = C_b/C_s. With a fixed dose Q per unit area and no further supply (drive-in) the profile is Gaussian, C = [Q/√(πDt)] exp(−x²/4Dt). Worked: C_b/C_s = 0.01 gives z = 1.82 (erfc(1.82) ≈ 0.01); for D = 2 × 10⁻¹³ cm²/s and t = 1800 s, √(Dt) = 1.90 × 10⁻⁵ cm, so x_j = 2 × 1.82 × 0.190 µm = 0.69 µm. Diffusion is isotropic: dopant spreads sideways under the mask edge by a comparable distance.
Ion implantation accelerates dopant ions (keV to MeV) into the wafer. The dose is Q = It/(qA) (ions per unit area) and the profile is approximately Gaussian, centred at the projected range R_p with straggle ΔR_p, with peak concentration C_max = Q/(√(2π) ΔR_p). The peak lies below the surface, depth is set by ion energy and dose by beam current and time, and the process is at low temperature, so a thin oxide or resist masks it. It damages the lattice, so an anneal restores crystallinity and puts the dopant on substitutional sites. Diffusion needs no vacuum and is cheap but is isotropic and hard to control for shallow junctions; implantation gives precise dose and depth control.
Photolithography transfers a pattern into the wafer:
(1) grow or deposit the film (for example SiO₂);
(2) clean and spin-coat photoresist, then soft-bake;
(3) align the mask and expose with ultraviolet light;
(4) develop: a positive resist becomes soluble where exposed, so the pattern is the same as on the mask; a negative resist crosslinks where exposed, so the unexposed part dissolves and the pattern is reversed;
(5) hard-bake;
(6) etch the unprotected film, by wet chemical etching (isotropic, undercuts the resist) or by dry plasma etching (reactive-ion etching, anisotropic);
(7) strip the resist. The smallest printable feature scales as k₁λ/NA, so shorter wavelength and larger numerical aperture improve resolution.
Preparation of powders and compaction
| Route | Principle | Powder produced |
|---|---|---|
| Gas atomisation | A melt stream is broken by a high-pressure inert-gas jet | Spherical, clean, fine; lower oxygen |
| Water atomisation | Same, with water jets | Irregular, cheaper, more oxide; compacts well |
| Mechanical milling | Ball impacts fracture brittle particles; ductile ones weld and fracture repeatedly | Fine, angular; mechanical alloying and nanostructures |
| Chemical routes | Reduction of oxides by H₂ or C; electrolysis; decomposition of carbonyls; precipitation or sol–gel | Sponge or fine powders; high purity, nanoscale possible |
A powder is described by size distribution, particle shape, apparent density and flow. Pure powder is compacted in a hard die (typically a few hundred MPa) to a green compact, using a lubricant, which cuts die-wall friction. Pressure raises density by particle rearrangement, then plastic deformation; because friction makes the density non-uniform, pressing from both ends and cold isostatic pressing (fluid pressure from all sides) give a more uniform green density. The green compact is weak, held together by cold welding and mechanical interlocking; the strength comes from sintering.
Sintering: driving force, stages, mechanisms and variables
Sintering heats the compact, typically to 0.7 to 0.9 of the absolute melting temperature, below melting, so that particles bond and the compact strengthens and usually densifies. The driving force is the reduction of total surface (solid–vapour) free energy: powder has a very large surface area, and replacing surface by grain boundary, which has lower energy, lowers the free energy. The matter transport is by atomic diffusion. In the initial stage necks grow at the contacts (the neck-to-particle radius ratio x/a rises to about 0.3), with little shrinkage. In the intermediate stage pores form a connected network of channels and most densification occurs. In the final stage pores are isolated spheres; they close slowly, and grain growth, which can trap pores inside grains, competes with densification.
Mechanisms differ in whether they densify. Surface diffusion, lattice diffusion from the surface and evaporation–condensation move matter from the particle surface to the neck: the neck grows but the particle centres do not approach, so there is no shrinkage. Grain-boundary diffusion, lattice diffusion from the boundary and plastic flow take matter from the contact region, so the centres approach and the compact shrinks. For a pair of spheres the neck radius follows x^n ∝ t (with a size factor), so x ∝ t^(1/n); n = 5 for lattice diffusion, 6 for grain-boundary diffusion, 7 for surface diffusion, 3 for evaporation–condensation and 2 for viscous flow. Neck growth therefore slows with time, and the rate depends on temperature through the diffusivity (Arrhenius), so a modest rise in temperature shortens the time greatly.
Particle size: finer powder has a larger surface area per unit mass and shorter diffusion paths, so it sinters faster and at lower temperature; for the same neck ratio Herring’s scaling law gives a time ∝ a³ for lattice diffusion and ∝ a⁴ for surface and grain-boundary diffusion (a is the particle radius). Shrinkage: for isotropic shrinkage with no mass loss, the relative density ρ/ρ_g = (1 − ΔL/L₀)⁻³ times the green relative density ρ_g. Worked: a compact of 60 % relative density that shrinks linearly by 5 % reaches 0.60/0.95³ = 0.70, that is 70 %. Entrapped gas and grain growth limit the final density; sintering aids, liquid phases, pressure (hot pressing) and fine powders raise it.
Key takeaways
- MG-Si comes from SiO₂ + 2C → Si + 2CO at about 1800 °C (98 to 99 %); EG-Si via the Siemens route: SiHCl₃ is made with HCl at about 300 °C, distilled, and reduced by H₂ on hot silicon rods at about 1100 °C.
- k = C_s/C_l. Normal freezing: C_s = kC₀(1 − f)^(k−1). First-pass zone refining: C = C₀[1 − (1 − k)e^(−kx/l)]; for k < 1 impurity goes to the end where the zone finishes, and zone refining fails when k is near 1.
- CZ pulls a crystal from a melt in a silica crucible (oxygen present, large diameters, cheap); FZ uses an RF zone with no crucible (very pure, smaller, costlier).
- Deal–Grove: x² + Ax = B(t + τ); linear x ≈ (B/A)(t + τ) when thin, parabolic x ≈ √(Bt) when thick; oxide consumes 0.44 x of silicon; wet is faster than dry.
- Diffusion doping gives erfc (constant source) or Gaussian (fixed dose) profiles with x_j = 2z√(Dt); implantation gives a Gaussian peak C_max = Q/(√(2π)ΔR_p) below the surface and needs an anneal. Lithography is coat, expose, develop, etch, strip.
- Powders: atomisation, milling, chemical routes; green compacts are weak. Sintering lowers surface energy by diffusion; neck x ∝ t^(1/n); only grain-boundary routes shrink the compact; finer powder and higher temperature speed it.
Practice questions (18)
Attempt each one before opening the answer. Every explanation names the tempting wrong option as well as the right one, because that is where marks are lost.
Which sequence correctly describes the production of semiconductor-grade polysilicon by the Siemens route from metallurgical-grade silicon?
Show answer
Answer: A — React with HCl to give liquid SiHCl₃, purify it by fractional distillation, then reduce it with H₂ on heated silicon rods
The Siemens route turns silicon into a volatile liquid compound, SiHCl₃ (made with HCl near 300 °C), distils it repeatedly to remove the chlorides of Fe, B and P, and deposits silicon by hydrogen reduction on rods at about 1100 °C. Carbothermic reduction alone gives only 98 to 99 % Si, which is the tempting but insufficient route. Oxidation followed by stripping would remove silicon itself, not the impurities.A silicon bar contains an impurity with segregation coefficient k = 0.05 (C_s/C_l at the interface). A molten zone is passed along the bar several times in one direction. After many passes, the impurity is most concentrated:
Show answer
Answer: C — at the end of the bar where the zone finishes, because the liquid retains the solute that the freezing solid rejects
For k < 1 the freezing solid takes only a fraction k of the liquid concentration, so the liquid zone becomes richer as it travels and deposits its load at the end where it freezes last. The first solid is the purest, since it takes k C₀. For k > 1 the result would be reversed.Why does float-zone (FZ) silicon contain much less oxygen than Czochralski (CZ) silicon?
Show answer
Answer: B — The molten zone touches no crucible, whereas in CZ the fused-silica crucible dissolves slightly into the melt
In CZ growth the silica crucible is the main oxygen source. In FZ the zone is held by surface tension, so there is no crucible and no oxygen pick-up. FZ is not a rapid-cooling process, and oxygen in silicon has k of about 1.25, so zone passes do not sweep it out, and the coil heats only a short zone near the melting point, not above the boiling point of silicon.Silicon is oxidised so that the oxide is already thick and growth follows the parabolic Deal–Grove law, x ≈ √(Bt). The oxide is 0.30 µm thick after 3 h. Taking τ = 0 and no change in temperature, the thickness after a total of 12 h is:
Show answer
Answer: D — 0.60 µm, because the thickness varies as the square root of time and the time has increased by a factor of four
In the parabolic regime x² = Bt, so x₂/x₁ = √(t₂/t₁) = √4 = 2, giving 0.60 µm. The linear-in-time options, 0.90 and 1.20 µm, describe the thin-oxide linear regime, which does not apply here. Quadrupling the time doubles, not quadruples, the thickness.Which statement correctly compares doping by ion implantation with doping by thermal diffusion?
Show answer
Answer: C — Implantation allows the dose to be set by beam current and time and the depth by ion energy, at low wafer temperature
Dose is Q = It/(qA) and depth follows the ion energy (projected range). The wafer need not be hot, so a thin resist or oxide can mask it. The profile peaks at R_p below the surface, and the damaged lattice must be annealed. Thermal diffusion is isotropic and does spread laterally.In photolithography, a positive photoresist is one in which:
Show answer
Answer: B — the exposed regions become soluble in the developer and are removed, so the pattern on the wafer is the same as the mask
A positive resist is made soluble by exposure, so developing removes the exposed areas and the pattern matches the mask. A negative resist behaves in the opposite way, crosslinking where exposed. The word positive has no connection with charge or with the use of ions.What is the thermodynamic driving force for the sintering of a loose compact of metal powder?
Show answer
Answer: A — The reduction of the total surface free energy as solid–vapour surface is replaced by lower-energy grain boundary
Powder has a very large surface area, and replacing surface by grain boundary lowers the free energy; atomic diffusion at 0.7 to 0.9 of the absolute melting temperature supplies the transport. Sintering is a solid-state process, so melting at contacts is not the mechanism. Stored strain relief is a separate effect, and not every powder has a polymorphic transformation.Two identical compacts of the same powder are sintered. In one, the neck between particles grows and strength rises, but the density and the dimensions do not change at all. Which mechanism is mainly operating?
Show answer
Answer: D — Surface diffusion, because it moves atoms from the particle surface to the neck without bringing the centres closer
Matter that is taken from the free surface and placed on the neck leaves the particle centres where they were, so the neck grows without shrinkage. Surface diffusion is one such non-densifying path (evaporation–condensation is another). Grain-boundary diffusion and lattice diffusion from the boundary both take matter from the contact region and so densify, and plastic flow also densifies.The sintering time needed to reach a given neck-to-particle radius ratio by lattice diffusion scales with particle radius a as a³ (Herring’s scaling law, for particles of the same shape and material at the same temperature). A powder of 20 µm particles needs 16 h. The time for a powder of 10 µm particles is:
Show answer
Answer: C — 2 h, because the time scales as the cube of the particle size
With t ∝ a³, halving a gives t = 16 × (1/2)³ = 2 h. The linear and square options use the wrong exponent, and the longer-time option reverses the trend: finer powder has a larger driving force per unit mass and shorter diffusion distances, so it sinters faster.Compared with water-atomised metal powder, gas-atomised powder typically has:
Show answer
Answer: C — a spherical shape with lower oxygen content, because the melt is broken up by an inert-gas jet
An inert-gas jet breaks the liquid stream into droplets that become spherical before freezing and that meet little oxygen. The irregular, oxide-rich shape belongs to water atomisation, the flake shape to milling of ductile metals and the sponge structure to solid-state oxide reduction.Which of the following statements about segregation and crystal growth of silicon are correct? Take k = C_s/C_l.
Show answer
Answer: A — For an impurity with k < 1, the first part of a normally frozen bar to solidify is purer than the starting material; B — In a Czochralski growth with a dopant of k < 1, the dopant concentration in the crystal rises from seed end to tail end; D — A neck is grown at the start of a Czochralski pull so that dislocations produced by thermal shock grow out of the crystal
The first solid takes k C₀ < C₀, so it is purer; with k < 1 the melt is enriched as the crystal grows, so the crystal dopant level rises along its length (C_s = kC₀(1 − f)^(k−1)); and the thin neck lets dislocations escape through its sides. Zone refining relies on rejection of solute, so it fails when k is near 1, which is the reason the statement about k close to 1 is wrong.Which of the following statements about thermal oxidation of silicon are correct?
Show answer
Answer: A — New oxide forms at the Si–SiO₂ interface, because the oxidant diffuses through the oxide already grown; C — An oxide of thickness x consumes about 0.44 x of silicon, because the oxide occupies roughly 2.27 times the volume of the silicon replaced
Growth occurs at the inner interface, and the volume ratio gives 1/2.27 = 0.44 for the silicon consumed. The regimes are the other way round: thin oxide is reaction-limited (linear) and thick oxide diffusion-limited (parabolic). Wet oxidation is faster, not slower, because H₂O is much more soluble in SiO₂ than O₂ is.Which of the following statements about sintering of metal powder compacts are correct?
Show answer
Answer: A — Raising the sintering temperature speeds every diffusion mechanism, because the diffusivity follows an Arrhenius law; B — In the final stage the pores are isolated, and grain growth that sweeps past them can leave them trapped inside grains
Diffusion is thermally activated, so a hotter compact sinters faster; and grain growth can detach grain boundaries from pores, leaving the pores isolated in the grains and hard to remove. Conventional sintering stays below the melting temperature. Surface diffusion does not shrink the compact at all, so it cannot shrink it faster than grain-boundary diffusion.A bar of silicon of uniform impurity concentration C₀ is zone refined in a single pass with a molten zone of length l. For the first pass, the solid concentration at distance x from the start is C(x) = C₀[1 − (1 − k) exp(−kx/l)]. The impurity has k = 0.2. At x = 4l the ratio C/C₀ is ____ (answer to two decimal places).
Numerical answer — type the value.
Show answer
Answer: 0.64
kx/l = 0.2 × 4 = 0.8. exp(−0.8) = 0.4493. C/C₀ = 1 − 0.8 × 0.4493 = 1 − 0.3595 = 0.6405, which is 0.64. Check: at x = 0 the formula gives k = 0.2 and at large x it tends to 1, so 0.64 lies between, as it should.Silicon is oxidised at a fixed temperature for which the Deal–Grove equation x² + Ax = B(t + τ) has A = 0.10 µm, B = 0.020 µm²/h and τ = 0. The oxide is grown from bare silicon for 5 h. Solve the quadratic exactly (do not use the parabolic approximation). The oxide thickness is ____ nm (answer as a whole number).
Numerical answer — type the value.
Show answer
Answer: 270
x² + 0.10x − 0.10 = 0 (since B t = 0.020 × 5 = 0.10). x = [−0.10 + √(0.01 + 0.40)]/2 = (−0.10 + 0.6403)/2 = 0.2702 µm = 270 nm. The parabolic approximation √(Bt) = 0.316 µm = 316 nm would overestimate, because the oxide is not yet far beyond A/2.Boron is predeposited into n-type silicon from a source that holds the surface concentration constant, so that C(x, t) = C_s erfc[x/(2√(Dt))]. The substrate concentration is C_b = 0.05 C_s, and erfc(z) = 0.05 at z = 1.39 (use this value). Take D = 4 × 10⁻¹⁴ cm²/s and t = 3600 s. The junction depth is ____ µm (answer to two decimal places).
Numerical answer — type the value.
Show answer
Answer: 0.33
The junction is where C = C_b, i.e. erfc(z) = 0.05, so x_j = 2z√(Dt). Dt = 4 × 10⁻¹⁴ × 3600 = 1.44 × 10⁻¹⁰ cm², so √(Dt) = 1.2 × 10⁻⁵ cm = 0.12 µm. x_j = 2 × 1.39 × 0.12 = 0.3336 µm, which is 0.33 µm. Unit slip to avoid: 1 cm = 10⁴ µm.Two copper spheres sinter by lattice diffusion, for which the neck radius x follows x⁵ ∝ t at fixed temperature and particle size. The neck radius is 3 µm after 2 h. Neglect any change of the geometry other than this scaling. The time at which the neck radius reaches 6 µm is ____ h (answer as a whole number).
Numerical answer — type the value.
Show answer
Answer: 64
x⁵ ∝ t gives t₂/t₁ = (x₂/x₁)⁵ = 2⁵ = 32. So t₂ = 2 × 32 = 64 h. Doubling the neck radius needs a 32-fold longer time, which shows how strongly neck growth slows. Assuming x ∝ t would wrongly give 4 h.A powder compact has a green relative density of 55 % of the theoretical density. On sintering it shrinks isotropically, with a linear shrinkage ΔL/L₀ = 8 %, and no mass is lost. The relative density after sintering, in per cent, is ____ (answer to one decimal place).
Numerical answer — type the value.
Show answer
Answer: 70.6
Mass is conserved and volume falls by the cube of the linear factor: V/V₀ = (1 − 0.08)³ = 0.7787. Density therefore rises by 1/0.7787 = 1.2842, so the relative density is 55 × 1.2842 = 70.6 %. Using 8 % as a volume shrinkage would give 59.8 %, which is wrong.