1. Simonka's ODE (Eq. 3.31, §3.4.1):
dCact/dt = −κ · (Cact − Ctot / (1 + Ctot/Css))
2. Let C∞ = Ctot / (1 + Ctot/Css) — the steady-state target. The ODE becomes:
dCact/dt = −κ · (Cact − C∞) → standard 1st-order linear ODE
3. With initial condition Cact(0) = 0, the analytical solution is:
Cact(t) = C∞ · (1 − e−κt) = Ctot / (1 + Ctot/Css) · (1 − e−κt)
4. Css(T) and κ(T) are Arrhenius: Z·exp(−Ea/kBT). Parameters from Table 3.5.
Css = solid solubility — the max dopant concentration the lattice can electrically activate at temperature T.
Ctot vs Cpeak — The model uses Ctot(x), the local concentration at depth x. Our calculator evaluates at x = Rp where Ctot = Cpeak (the worst-case point). At shallower/deeper depths where Ctot < Cpeak, activation ratio is higher.
Red dot = your current settings. Left chart sweeps T. Right chart sweeps t.