Weibull Distribution Mathematics in Semiconductor Manufacturing

A comprehensive guide to the mathematical foundations and applications of Weibull distribution in semiconductor reliability engineering.

1. Fundamental Weibull Mathematics

1.1 The Core Equations

Two-parameter Weibull Probability Density Function (PDF):

Cumulative Distribution Function (CDF) — probability of failure by time $t$:

Reliability (Survival) Function:

Parameter Definitions:

• $t \geq 0$ — random variable (typically time or stress cycles)
• $\beta > 0$ — shape parameter (Weibull slope/modulus)
• $\eta > 0$ — scale parameter (characteristic life, where $F(\eta) = 0.632$)

1.2 Three-Parameter Weibull

Adding a location parameter $\gamma$ (threshold/minimum life):

1.3 The Hazard Function (Instantaneous Failure Rate)

Physical Interpretation of Shape Parameter $\beta$:

This directly models the semiconductor bathtub curve.

2. Semiconductor-Specific Applications

2.1 Time-Dependent Dielectric Breakdown (TDDB)

Gate oxide breakdown follows Weibull statistics. The area scaling law derives from weakest-link theory:

Where:

• $A_1$ — reference test area
• $A_2$ — target device area
• $\eta_1$ — characteristic life at area $A_1$
• $\eta_2$ — predicted characteristic life at area $A_2$

Typical $\beta$ values for oxide breakdown:

• Intrinsic breakdown: $\beta \approx 10$–$30$ (tight distribution)
• Extrinsic/defect-related: $\beta \approx 1$–$5$ (broader distribution)

2.2 Electromigration

Metal interconnect failure combines Black's equation with Weibull statistics:

Where:

• $MTF$ — median time to failure
• $j$ — current density ($A/cm^2$)
• $n$ — current density exponent (typically 1–2)
• $E_a$ — activation energy (eV)
• $k_B$ — Boltzmann constant ($8.617 \times 10^{-5}$ eV/K)
• $T$ — absolute temperature (K)

Typical $\beta$ values: 2–4 (wear-out behavior)

2.3 Hot Carrier Injection (HCI)

Degradation follows power-law kinetics:

Where:

• $\Delta V_{th}$ — threshold voltage shift
• $t$ — stress time
• $n$ — time exponent (typically 0.3–0.5)

2.4 Negative Bias Temperature Instability (NBTI)

For PMOS transistors:

3. Statistical Analysis Methods

3.1 Weibull Probability Plotting

Linearization transformation — take double logarithm of CDF:

Plotting $\ln[-\ln(1-F)]$ vs $\ln(t)$:

• Slope = $\beta$
• Intercept at $F = 0.632$ gives $t = \eta$

Bernard's Median Rank Approximation for ranking data:

Where:

• $r$ — rank of the $r$-th ordered failure
• $n$ — total sample size

3.2 Maximum Likelihood Estimation (MLE)

Log-likelihood function for $n$ samples with $r$ failures and $(n-r)$ censored units:

MLE Estimator for $\eta$:

MLE Equation for $\beta$ (solve numerically):

4. Accelerated Life Testing Mathematics

4.1 Acceleration Factors

Arrhenius Model (Thermal Acceleration):

Exponential Voltage Acceleration:

Power-Law Voltage Acceleration:

Life Extrapolation:

4.2 Combined Stress Models (Eyring)

Where:

• $RH$ — relative humidity
• $m$ — humidity exponent
• Additional stress factors can be included

5. Competing Failure Modes

5.1 Series (Competing Risks) Model

Device fails when the first mechanism fails:

Combined CDF:

5.2 Mixture Model

Different subpopulations with different failure characteristics:

Where:

• $p_i$ — proportion in subpopulation $i$
• $\sum_{i=1}^{k} p_i = 1$
• $F_i(t)$ — CDF for subpopulation $i$

PDF for mixture:

6. Key Derived Quantities

6.1 Moments of the Weibull Distribution

$k$-th Raw Moment:

Mean (MTTF — Mean Time To Failure):

Variance:

Standard Deviation:

6.2 Percentile Lives (B$X$ Life)

Time by which $X\%$ have failed:

Common Percentile Lives:

6.3 Characteristic Life Significance

At $t = \eta$:

This means 63.2% of units have failed by the characteristic life, regardless of $\beta$.

7. Confidence Bounds

7.1 Fisher Information Matrix Approach

Information Matrix:

Asymptotic Variance-Covariance Matrix:

Fisher Matrix Elements:

Where $\gamma_E \approx 0.5772$ is the Euler-Mascheroni constant.

7.2 Likelihood Ratio Bounds (Preferred for Small Samples)

Approximate $(1-\alpha)$ Confidence Interval:

8. Order Statistics

8.1 Expected Value of Order Statistics

For $n$ samples, the expected value of the $r$-th order statistic:

8.2 Plotting Positions

Bernard's Approximation (recommended):

Hazen's Approximation:

Mean Rank:

9. Practical Example: Gate Oxide Qualification

9.1 Test Setup

• Sample size: 50 oxide capacitors
• Stress conditions: 125°C, 1.2× nominal voltage
• Test duration: 1000 hours
• Failures: 8 units at times: 156, 289, 412, 523, 678, 734, 891, 967 hours
• Censored: 42 units still running at 1000h

9.2 Analysis Steps

Step 1: Calculate Median Ranks

Step 2: MLE Results

Step 3: Calculate Acceleration Factor

Given: $E_a = 0.7$ eV, voltage exponent $n = 40$

Step 4: Extrapolate to Use Conditions

Step 5: Calculate B0.1 Life

10. Key Equations

10.1 Quick Reference Table

10.2 Why Weibull Works for Semiconductors

1. Physical meaning of $\beta$ — directly indicates failure mechanism type 2. Area/volume scaling — derives from extreme value theory (weakest-link) 3. Censored data handling — essential since most test units don't fail 4. Acceleration compatibility — seamlessly integrates with physics-based models 5. Competing risks framework — models complex multi-mechanism devices

Gamma Function Values

Common values of $\Gamma(1 + 1/\beta)$ for mean life calculations:

Common Activation Energies