Dimensional plasma fluxes and a material stack in; coupled transport, surface chemistry, self-consistent charging, and a moving three-dimensional profile out — fast enough to matter, accurate against real published experiments, and built to predict process conditions it has never seen.
One continuous calculation: energy–angle-resolved ion and neutral arrival → in-feature ballistic and diffuse transport → material-specific surface chemistry with finite, conserved surface state → mask erosion, passivation, and substrate etch → 3-D multi-material level-set motion, straight through pinch-off, enclosed-cavity formation, and reopening.
Same case, same hardware, same etch depth — petch against ViennaPS, the leading open-source profile simulator. Not projections: wall-clock recordings of both tools running.
Krüger 2024: SiO₂ trench etch in fluorocarbon plasma, 60 s process, measured by SEM, driven by the paper's own measured reactor fluxes and ion energy–angle distributions. Two interpretable physical closures calibrated; everything else is physics.
| Observable | Experiment | petch | Error |
|---|---|---|---|
| Minimum mask opening | 45.000 nm | 45.085 nm | +0.19% |
| Etch depth | 825.000 nm | 853.219 nm | +3.42% |
Electron shading charges deep features; the resulting in-feature field steers ions and produces notching, twisting, and profile distortion — the dominant failure modes of high-aspect-ratio dielectric etch. petch computes this self-consistently: field ↔ charged trajectories ↔ accumulated surface charge, converged until currents balance on every dielectric node, on the evolving profile.
ViennaPS is an excellent geometry engine — petch matches it exactly where it is strong, and carries the physics and guarantees it doesn't have.
| petch | ViennaPS | |
|---|---|---|
| Self-consistent charging (fields, deflection, notching) | Converged I⁺=I⁻ on the evolving 3-D profile | Not computed — structurally blind to charging effects |
| Surface chemistry | Finite, conserved multi-material surface inventories (mask film, passivation, polymer) with exact ledgers | Parametrized rate/coverage models |
| Reactor boundary | Ingests measured energy–angle-resolved fluxes (checksummed) — predicts new conditions from their physics | Parametrized source distributions, retuned per condition |
| Reproducibility & audit | Bit-identical reruns, conservation receipts, checksummed provenance, freeze/reveal validation gates | No determinism or audit contract |
| Speed on matched cases | ~50× (2-D trench), ~100× (3-D via), depth-matched, same hardware | Baseline |
| Geometric transport & topology | Parity — verified directly: ARDE transport agrees (0.727 vs 0.73 at AR 8.6) and pinch-off/reopen event times match ViennaLS within one grid cell. We benchmark against it because it's good. | |
Developed against published wafer experiments spanning cryo SF₆/O₂ deep-trench ARDE, fluorocarbon SiO₂ trenches, charging microtrenches and notching, and Bosch cycling — one engine, one set of physics contracts.
Exact material conservation (zero residual, enforced every step), watertight topology-complete surface meshes, float64-certified ray transport, deterministic seeded sampling — bit-identical on re-run. And the observables converge: refine the grid 2× and the answers stay.
Speed comparisons (sections 02, 06) are measured wall-clock results on matched cases from the engine's GPU-accelerated transport generation; the certified full-physics operator used in the validation campaign runs slower today and inherits that acceleration stack next. The Krüger base result is a two-observable calibrated reproduction; the blind-transfer campaign is what converts it into held-out validation. Charging notching is demonstrated causally with the experimentally observed trend; absolute notch depth is not yet calibrated against a wafer measurement. The de Boer ARDE result shown is from the engine generation that produced the speed demonstrations; it is being re-verified through the current certified operator as part of the ongoing validation program. Profile-edge texture in plots is grid-scale numerics, not a physical roughness prediction.