Semiconductor Manufacturing Process…

Home› Knowledge Base› Semiconductor Manufacturing Process: Poisson Statistics & Mathematical Modeling

Semiconductor Manufacturing Process: Poisson Statistics & Mathematical Modeling

1. Introduction: Why Poisson Statistics?

Semiconductor defects satisfy the classical Poisson conditions:

Rare events — Defects are sparse relative to the total chip area
Independence — Defect occurrences are approximately independent
Homogeneity — Within local regions, defect rates are constant
No simultaneity — At infinitesimal scales, simultaneous defects have zero probability

1.1 The Poisson Probability Mass Function

The probability of observing exactly $k$ defects:

P(X = k) = \frac{\lambda^k e^{-\lambda}}{k!}

where the expected number of defects is:

\lambda = D_0 \cdot A

Parameter definitions:

$D_0$ — Defect density (defects per unit area, typically defects/cm²)
$A$ — Chip area (cm²)
$\lambda$ — Mean number of defects per chip

1.2 Key Statistical Properties

Property	Formula
Mean	$E[X] = \lambda$
Variance	$\text{Var}(X) = \lambda$
Variance-to-Mean Ratio	$\frac{\text{Var}(X)}{E[X]} = 1$

Note: The equality of mean and variance (equidispersion) is a signature property of the Poisson distribution. Real semiconductor data often shows overdispersion (variance > mean), motivating compound models.

2. Fundamental Yield Equation

2.1 The Seeds Model (Simple Poisson)

A chip is functional if and only if it has zero killer defects. Under Poisson assumptions:

\boxed{Y = P(X = 0) = e^{-D_0 A}}

Derivation:

P(X = 0) = \frac{\lambda^0 e^{-\lambda}}{0!} = e^{-\lambda} = e^{-D_0 A}

2.2 Limitations of Simple Poisson

Assumes uniform defect density across the wafer (unrealistic)
Does not account for clustering of defects
Consistently underestimates yield for large chips
Ignores wafer-to-wafer and lot-to-lot variation

3. Compound Poisson Models

3.1 The Negative Binomial Approach

Model the defect density $D_0$ as a random variable with Gamma distribution:

D_0 \sim \text{Gamma}\left(\alpha, \frac{\alpha}{\bar{D}}\right)

Gamma probability density function:

f(D_0) = \frac{(\alpha/\bar{D})^\alpha}{\Gamma(\alpha)} D_0^{\alpha-1} e^{-\alpha D_0/\bar{D}}

where:

$\bar{D}$ — Mean defect density
$\alpha$ — Clustering parameter (shape parameter)

3.2 Resulting Yield Model

When defect density is Gamma-distributed, the defect count follows a Negative Binomial distribution, yielding:

\boxed{Y = \left(1 + \frac{D_0 A}{\alpha}\right)^{-\alpha}}

3.3 Physical Interpretation of Clustering Parameter $\alpha$

$\alpha$ Value	Physical Interpretation
$\alpha \to \infty$	Uniform defects — recovers simple Poisson model
$\alpha \approx 1-5$	Typical semiconductor clustering
$\alpha \to 0$	Extreme clustering — defects occur in tight groups

3.4 Overdispersion

The variance-to-mean ratio for the Negative Binomial:

\frac{\text{Var}(X)}{E[X]} = 1 + \frac{\bar{D}A}{\alpha} > 1

This overdispersion (ratio > 1) matches empirical observations in semiconductor manufacturing.

4. Classical Yield Models

4.1 Comparison Table

Model	Yield Formula	Assumed Density Distribution
Seeds (Poisson)	$Y = e^{-D_0 A}$	Delta function (uniform)
Murphy	$Y = \left(\frac{1 - e^{-D_0 A}}{D_0 A}\right)^2$	Triangular
Negative Binomial	$Y = \left(1 + \frac{D_0 A}{\alpha}\right)^{-\alpha}$	Gamma
Moore	$Y = e^{-\sqrt{D_0 A}}$	Empirical
Bose-Einstein	$Y = \frac{1}{1 + D_0 A}$	Exponential

4.2 Murphy's Yield Model

Assumes triangular distribution of defect densities:

Y_{\text{Murphy}} = \left(\frac{1 - e^{-D_0 A}}{D_0 A}\right)^2

Taylor expansion for small $D_0 A$:

Y_{\text{Murphy}} \approx 1 - \frac{(D_0 A)^2}{12} + O((D_0 A)^4)

4.3 Limiting Behavior

As $D_0 A \to 0$ (low defect density):

\lim_{D_0 A \to 0} Y = 1 \quad \text{(all models)}

As $D_0 A \to \infty$ (high defect density):

\lim_{D_0 A \to \infty} Y = 0 \quad \text{(all models)}

5. Critical Area Analysis

5.1 Definition

Not all chip area is equally vulnerable. Critical area $A_c$ is the region where a defect of size $d$ causes circuit failure.

A_c(d) = \int_{\text{layout}} \mathbf{1}\left[\text{defect at } (x,y) \text{ with size } d \text{ causes failure}\right] \, dx \, dy

5.2 Critical Area for Shorts

For two parallel conductors with:

Length: $L$
Spacing: $S$

A_c^{\text{short}}(d) = \begin{cases} 2L(d - S) & \text{if } d > S \\ 0 & \text{if } d \leq S \end{cases}

5.3 Critical Area for Opens

For a conductor with:

Width: $W$
Length: $L$

A_c^{\text{open}}(d) = \begin{cases} L(d - W) & \text{if } d > W \\ 0 & \text{if } d \leq W \end{cases}

5.4 Total Critical Area

Integrate over the defect size distribution $f(d)$:

A_c = \int_0^\infty A_c(d) \cdot f(d) \, dd

5.5 Defect Size Distribution

Typically modeled as power-law:

f(d) = C \cdot d^{-p} \quad \text{for } d \geq d_{\min}

Typical values:

Exponent: $p \approx 2-4$
Normalization constant: $C = (p-1) \cdot d_{\min}^{p-1}$

Alternative: Log-normal distribution (common for particle contamination):

f(d) = \frac{1}{d \sigma \sqrt{2\pi}} \exp\left(-\frac{(\ln d - \mu)^2}{2\sigma^2}\right)

6. Multi-Layer Yield Modeling

6.1 Modern IC Structure

Modern integrated circuits have 10-15+ metal layers. Each layer $i$ has:

Defect density: $D_i$
Critical area: $A_{c,i}$
Clustering parameter: $\alpha_i$ (for Negative Binomial)

6.2 Poisson Multi-Layer Yield

Y_{\text{total}} = \prod_{i=1}^{n} Y_i = \prod_{i=1}^{n} e^{-D_i A_{c,i}}

Simplified form:

\boxed{Y_{\text{total}} = \exp\left(-\sum_{i=1}^{n} D_i A_{c,i}\right)}

6.3 Negative Binomial Multi-Layer Yield

\boxed{Y_{\text{total}} = \prod_{i=1}^{n} \left(1 + \frac{D_i A_{c,i}}{\alpha_i}\right)^{-\alpha_i}}

6.4 Log-Yield Decomposition

Taking logarithms for analysis:

\ln Y_{\text{total}} = -\sum_{i=1}^{n} D_i A_{c,i} \quad \text{(Poisson)}

\ln Y_{\text{total}} = -\sum_{i=1}^{n} \alpha_i \ln\left(1 + \frac{D_i A_{c,i}}{\alpha_i}\right) \quad \text{(Negative Binomial)}

7. Spatial Point Process Formulation

7.1 Inhomogeneous Poisson Process

Intensity function $\lambda(x, y)$ varies spatially across the wafer:

P(k \text{ defects in region } R) = \frac{\Lambda(R)^k e^{-\Lambda(R)}}{k!}

where the integrated intensity is:

\Lambda(R) = \iint_R \lambda(x,y) \, dx \, dy

7.2 Cox Process (Doubly Stochastic)

The intensity $\lambda(x,y)$ is itself a random field:

\lambda(x,y) = \exp\left(\mu + Z(x,y)\right)

where:

$\mu$ — Baseline log-intensity
$Z(x,y)$ — Gaussian random field with spatial correlation function $\rho(h)$

Correlation structure:

\text{Cov}(Z(x_1, y_1), Z(x_2, y_2)) = \sigma^2 \rho(h)

where $h = \sqrt{(x_2-x_1)^2 + (y_2-y_1)^2}$

7.3 Neyman Type A (Cluster Process)

Models defects occurring in clusters:

1. Cluster centers: Poisson process with intensity $\lambda_c$ 2. Defects per cluster: Poisson with mean $\mu$ 3. Defect positions: Scattered around cluster center (e.g., isotropic Gaussian)

Probability generating function:

G(s) = \exp\left[\lambda_c A \left(e^{\mu(s-1)} - 1\right)\right]

Mean and variance:

E[N] = \lambda_c A \mu

\text{Var}(N) = \lambda_c A \mu (1 + \mu)

8. Statistical Estimation Methods

8.1 Maximum Likelihood Estimation

8.1.1 Data Structure

Given:

$n$ chips with areas $A_1, A_2, \ldots, A_n$
Binary outcomes $y_i \in \{0, 1\}$ (pass/fail)

8.1.2 Likelihood Function

\mathcal{L}(D_0, \alpha) = \prod_{i=1}^n Y_i^{y_i} (1 - Y_i)^{1-y_i}

where $Y_i = \left(1 + \frac{D_0 A_i}{\alpha}\right)^{-\alpha}$

8.1.3 Log-Likelihood

\ell(D_0, \alpha) = \sum_{i=1}^n \left[y_i \ln Y_i + (1-y_i) \ln(1-Y_i)\right]

8.1.4 Score Equations

\frac{\partial \ell}{\partial D_0} = 0, \quad \frac{\partial \ell}{\partial \alpha} = 0

Note: Requires numerical optimization (Newton-Raphson, BFGS, or EM algorithm).

8.2 Bayesian Estimation

8.2.1 Prior Distribution

D_0 \sim \text{Gamma}(a, b)

\pi(D_0) = \frac{b^a}{\Gamma(a)} D_0^{a-1} e^{-b D_0}

8.2.2 Posterior Distribution

Given defect count $k$ on area $A$:

D_0 \mid k \sim \text{Gamma}(a + k, b + A)

Posterior mean:

\hat{D}_0 = \frac{a + k}{b + A}

Posterior variance:

\text{Var}(D_0 \mid k) = \frac{a + k}{(b + A)^2}

8.2.3 Sequential Updating

Bayesian framework enables sequential learning:

\text{Prior}_n \xrightarrow{\text{data } k_n} \text{Posterior}_n = \text{Prior}_{n+1}

9. Statistical Process Control

9.1 c-Chart (Defect Counts)

For constant inspection area:

Center line: $\bar{c}$ (average defect count)
Upper Control Limit (UCL): $\bar{c} + 3\sqrt{\bar{c}}$
Lower Control Limit (LCL): $\max(0, \bar{c} - 3\sqrt{\bar{c}})$

9.2 u-Chart (Defects per Unit Area)

For variable inspection area $n_i$:

u_i = \frac{c_i}{n_i}

Center line: $\bar{u}$
Control limits: $\bar{u} \pm 3\sqrt{\frac{\bar{u}}{n_i}}$

9.3 Overdispersion-Adjusted Charts

For clustered defects (Negative Binomial), inflate the variance:

\text{UCL} = \bar{c} + 3\sqrt{\bar{c}\left(1 + \frac{\bar{c}}{\alpha}\right)}

\text{LCL} = \max\left(0, \bar{c} - 3\sqrt{\bar{c}\left(1 + \frac{\bar{c}}{\alpha}\right)}\right)

9.4 CUSUM Chart

Cumulative sum for detecting small persistent shifts:

C_t^+ = \max(0, C_{t-1}^+ + (x_t - \mu_0 - K))

C_t^- = \max(0, C_{t-1}^- - (x_t - \mu_0 + K))

where:

$K$ — Slack value (typically $0.5\sigma$)
Signal when $C_t^+$ or $C_t^-$ exceeds threshold $H$

10. EUV Lithography Stochastic Effects

10.1 Photon Shot Noise

At extreme ultraviolet wavelength (13.5 nm), photon shot noise becomes critical.

Number of photons absorbed in resist volume $V$:

N \sim \text{Poisson}(\Phi \cdot \sigma \cdot V)

where:

$\Phi$ — Photon fluence (photons/area)
$\sigma$ — Absorption cross-section
$V$ — Resist volume

10.2 Line Edge Roughness (LER)

Stochastic photon absorption causes spatial variation in resist exposure:

\sigma_{\text{LER}} \propto \frac{1}{\sqrt{\Phi \cdot V}}

Critical Design Rule:

\text{LER}_{3\sigma} < 0.1 \times \text{CD}

where CD = Critical Dimension (feature size)

10.3 Stochastic Printing Failures

Probability of insufficient photons in a critical volume:

P(\text{failure}) = P(N < N_{\text{threshold}}) = \sum_{k=0}^{N_{\text{threshold}}-1} \frac{\lambda^k e^{-\lambda}}{k!}

where $\lambda = \Phi \sigma V$

11. Reliability and Latent Defects

11.1 Defect Classification

Not all defects cause immediate failure:

Killer defects: Cause immediate functional failure
Latent defects: May cause reliability failures over time

\lambda_{\text{total}} = \lambda_{\text{killer}} + \lambda_{\text{latent}}

11.2 Yield vs. Reliability

Initial Yield:

Y = e^{-\lambda_{\text{killer}} \cdot A}

Reliability Function:

R(t) = e^{-\lambda_{\text{latent}} \cdot A \cdot H(t)}

where $H(t)$ is the cumulative hazard function for latent defect activation.

11.3 Weibull Activation Model

H(t) = \left(\frac{t}{\eta}\right)^\beta

Parameters:

$\eta$ — Scale parameter (characteristic life)
$\beta$ — Shape parameter

$\beta < 1$: Decreasing failure rate (infant mortality)

$\beta = 1$: Constant failure rate (exponential)

$\beta > 1$: Increasing failure rate (wear-out)

12. Complete Mathematical Framework

12.1 Hierarchical Model Structure

<svg viewBox="0 0 760 470" xmlns="http://www.w3.org/2000/svg" font-family="-apple-system,BlinkMacSystemFont,Segoe UI,Roboto,sans-serif">
  <rect width="760" height="470" fill="#0d1117" rx="8"/>
  <text x="380" y="28" fill="#e6edf3" font-size="21" font-weight="700" text-anchor="middle">Poisson Yield Modeling &amp; Defect Distribution</text>
  <text x="380" y="48" fill="#8b98a5" font-size="12" text-anchor="middle">Statistical Defect Clustering, Murphy/Seeds Models, and Wafer Die Yield Prediction</text>

  <!-- Left: Poisson Distribution Model -->
  <g transform="translate(30, 70)">
    <rect width="220" height="350" fill="#161b22" stroke="#30363d" stroke-width="1.5" rx="6"/>
    <text x="110" y="24" fill="#79c0ff" font-size="13" font-weight="600" text-anchor="middle">1. Pure Poisson Model</text>

    <rect x="20" y="45" width="180" height="120" fill="#0d1117" stroke="#30363d" rx="4"/>
    <!-- Uniform random defects -->
    <rect x="35" y="55" width="150" height="100" fill="#161b22" stroke="#30363d"/>
    <circle cx="50" cy="70" r="3" fill="#f85149"/>
    <circle cx="110" cy="85" r="3" fill="#f85149"/>
    <circle cx="160" cy="120" r="3" fill="#f85149"/>
    <circle cx="80" cy="140" r="3" fill="#f85149"/>
    <circle cx="130" cy="65" r="3" fill="#f85149"/>

    <text x="110" y="185" fill="#e6edf3" font-size="11" text-anchor="middle">Random Uniform Defects</text>

    <rect x="20" y="200" width="180" height="135" fill="#0d1117" stroke="#30363d" rx="4"/>
    <text x="110" y="220" fill="#79c0ff" font-size="12" font-weight="600" text-anchor="middle">Y = e^(-A · D0)</text>
    <text x="110" y="245" fill="#8b98a5" font-size="10" text-anchor="middle">A = Die Area (cm²)</text>
    <text x="110" y="265" fill="#8b98a5" font-size="10" text-anchor="middle">D0 = Defect Density (/cm²)</text>
    <text x="110" y="285" fill="#8b98a5" font-size="10" text-anchor="middle">Assumes zero defect clustering</text>
    <text x="110" y="310" fill="#f85149" font-size="10" font-weight="600" text-anchor="middle">Pessimistic for Large Dies</text>
  </g>

  <!-- Middle: Negative Binomial (Clustered) -->
  <g transform="translate(270, 70)">
    <rect width="220" height="350" fill="#161b22" stroke="#30363d" stroke-width="1.5" rx="6"/>
    <text x="110" y="24" fill="#79c0ff" font-size="13" font-weight="600" text-anchor="middle">2. Negative Binomial</text>

    <rect x="20" y="45" width="180" height="120" fill="#0d1117" stroke="#30363d" rx="4"/>
    <!-- Clustered defects -->
    <rect x="35" y="55" width="150" height="100" fill="#161b22" stroke="#30363d"/>
    <circle cx="130" cy="110" r="3" fill="#f85149"/>
    <circle cx="135" cy="115" r="3" fill="#f85149"/>
    <circle cx="128" cy="122" r="3" fill="#f85149"/>
    <circle cx="140" cy="108" r="3" fill="#f85149"/>

    <text x="110" y="185" fill="#e6edf3" font-size="11" text-anchor="middle">Spatial Defect Clustering</text>

    <rect x="20" y="200" width="180" height="135" fill="#0d1117" stroke="#30363d" rx="4"/>
    <text x="110" y="220" fill="#d2a8ff" font-size="11" font-weight="600" text-anchor="middle">Y = (1 + A·D0 / α)^(-α)</text>
    <text x="110" y="245" fill="#8b98a5" font-size="10" text-anchor="middle">α = Cluster Parameter</text>
    <text x="110" y="265" fill="#8b98a5" font-size="10" text-anchor="middle">α → ∞ : Pure Poisson</text>
    <text x="110" y="285" fill="#8b98a5" font-size="10" text-anchor="middle">α → 1 : High Clustering</text>
    <text x="110" y="310" fill="#3fb950" font-size="10" font-weight="600" text-anchor="middle">Accurate Industrial Model</text>
  </g>

  <!-- Right: Murphy / Seeds Empirical Models -->
  <g transform="translate(510, 70)">
    <rect width="220" height="350" fill="#161b22" stroke="#30363d" stroke-width="1.5" rx="6"/>
    <text x="110" y="24" fill="#79c0ff" font-size="13" font-weight="600" text-anchor="middle">3. Empirical Models</text>

    <rect x="20" y="45" width="180" height="120" fill="#0d1117" stroke="#30363d" rx="4"/>
    <!-- Murphy & Seeds Equations -->
    <text x="110" y="75" fill="#58a6ff" font-size="11" font-weight="600" text-anchor="middle">Murphy Model:</text>
    <text x="110" y="95" fill="#e6edf3" font-size="10" text-anchor="middle">Y = ((1 - e^(-A D0)) / (A D0))²</text>
    <text x="110" y="125" fill="#58a6ff" font-size="11" font-weight="600" text-anchor="middle">Seeds Model:</text>
    <text x="110" y="145" fill="#e6edf3" font-size="10" text-anchor="middle">Y = 1 / (1 + A D0)</text>

    <text x="110" y="185" fill="#e6edf3" font-size="11" text-anchor="middle">Yield vs Die Area Curve</text>

    <rect x="20" y="200" width="180" height="135" fill="#0d1117" stroke="#30363d" rx="4"/>
    <text x="110" y="220" fill="#f0883e" font-size="11" font-weight="600" text-anchor="middle">Chiplet Advantage</text>
    <text x="110" y="240" fill="#8b98a5" font-size="10" text-anchor="middle">Splitting 800mm² Monolithic</text>
    <text x="110" y="260" fill="#8b98a5" font-size="10" text-anchor="middle">into 4 × 200mm² Chiplets</text>
    <text x="110" y="280" fill="#3fb950" font-size="11" font-weight="600" text-anchor="middle">Yield Increases from 35% → 82%</text>
    <text x="110" y="310" fill="#8b98a5" font-size="10" text-anchor="middle">KGD (Known Good Die) Test</text>
  </g>

  <text x="380" y="452" fill="#6b7684" font-size="11" text-anchor="middle">Wafer Yield Economics Engine Driving Reticle Limit Breakdown &amp; Modular Chiplet Partitioning</text>
</svg>

12.2 Summary of Key Equations

Concept	Equation
Poisson PMF	$P(X=k) = \frac{\lambda^k e^{-\lambda}}{k!}$
Simple Yield	$Y = e^{-D_0 A}$
Negative Binomial Yield	$Y = \left(1 + \frac{D_0 A}{\alpha}\right)^{-\alpha}$
Multi-Layer Yield	$Y = \prod_i \left(1 + \frac{D_i A_{c,i}}{\alpha_i}\right)^{-\alpha_i}$
Critical Area (shorts)	$A_c^{\text{short}}(d) = 2L(d-S)$ for $d > S$
Defect Size Distribution	$f(d) \propto d^{-p}$, $p \approx 2-4$
Bayesian Posterior	$D_0 \mid k \sim \text{Gamma}(a+k, b+A)$
Control Limits	$\bar{c} \pm 3\sqrt{\bar{c}(1 + \bar{c}/\alpha)}$
LER Scaling	$\sigma_{\text{LER}} \propto (\Phi V)^{-1/2}$

12.3 Typical Parameter Values

Parameter	Typical Range	Units
Defect density $D_0$	0.01 - 1.0	defects/cm²
Clustering parameter $\alpha$	0.5 - 5	dimensionless
Defect size exponent $p$	2 - 4	dimensionless
Chip area $A$	1 - 800	mm²

Poisson statisticsdefect distributionyield modelingcritical areaclustering

Explore 500+ Semiconductor & AI Topics

From EUV lithography to CUDA optimization — search the full knowledge base or chat with our AI assistant.

🔍 Search Topics 💬 Ask CFSGPT 📚 Browse All