Holstein → CRETIN · an explorable explanation docs · team · refs · other gifts

Holstein → CRETIN

Eight layers of radiation trapping. From a single atomic resonance to a CRETIN-class 3D NLTE solver. Then a LightCell 819 nm Na extraction calculator falls out.

You can read Molisch & Oehry's 443-page Radiation Trapping in Atomic Vapours (Oxford 1998) end to end. Most people don't. The mathematics is not the obstacle — it's that nobody has ever made it manipulable. Bret Victor would tear into this book in a paragraph. Howard Scott already did the production-scale version (CRETIN, JQSRT 71:689, 2001) but the binary is export-controlled and the code is hard to teach with.

So we built it — interactively, in eight layers, with the constants from Molisch & Oehry Appendix D and the architectural shape of CRETIN. Every parameter is a slider. Every equation has a worked numerical case. Each layer compounds: Layer N reuses the objects from Layers 1..N-1.

The terminal layer is the engineering case study that motivated this whole effort: LightCell Energy's sodium-vapour 819 nm extraction problem. You'll be able to drag bath temperature, optical depths, cavity-Q, and recuperator efficiency and watch the design surface live.

The eight layers

Layer 1
The Holstein equation
Lineshapes (Doppler / Lorentz / Voigt), 1D modes for slab / cyl / sph, trapping factors $g_j$ from M&O fitting equations. Wing-escape physics interactive.
shipped
Layer 2
Higher-order modes & the escape factor's accidental success
3-term vs 10-term modal expansion. Why narrow-laser excitation breaks single-mode fits. Why the escape factor works in cylindrical lamps only by coincidence.
shipped
Layer 3
CFR vs PFR
Redistribution functions $R_I, R_{II}, R_{III}, R_V$. JW approximation breakdown at $k_0L > 500$. Live decision flowchart.
shipped
Layer 4
Multi-level atoms, hfs, isotope doping
Hg-196 doping optimum at 2.6 %, reproducing the GTE 1985 result that increased fluorescent-lamp efficiency by >1 %.
shipped
Layer 5
Saturation & three-level depletion
Ground-state depletion, Bennett holes, the three-level metastable depletion regime that matters for LightCell.
shipped
Layer 6
NLTE + CRETIN architecture
Three-phase continuum / line / spectral. Auer–Heasley linearization + ALI. Tokamak-divertor Lyman-α suppression by order of magnitude.
shipped
Layer 7
2D / 3D solver
Live 2D Monte Carlo photon-tracing. Source-function heatmap, escape angular distribution, $p_i$ histogram. Scott 2001 §5.1 case discussion.
shipped
Layer 8
LightCell 819 nm extraction · TERMINAL FEAT
Eight sliders. Live $\eta_{819}$ readout. Boltzmann ladder + loss-channel breakdown + sensitivity-to-$T_h$ curve. LTE-bath framing (operator-corrected).
shipped

Open questions surfaced along the way

Team: Feynman (theorist-lead) · Bret Victor (UX-lead) · Alan Kay (architecture) · Howard Scott (CRETIN computational) · Andreas Molisch + Bernhard Oehry (canon) · N.N. Bezuglov (Russian analytical extension) · Holstein / Biberman / Milne (historical anchors)

Source material: M&O 1998 (OCR'd via Gemini 3.1 Pro — see the bench essay) · 6 CRETIN PDFs (Scott 2001, Sequoia-Tillack-Scott 2006, Langer-Scott-Marinak-Landen 2002, Scott 2011/2017, fusion-plasma Ch II) · 3 deep-research runs (Claude, Gemini, ChatGPT — cross-referenced for trust)

Built by an LLM and a human in a few sessions, ~2026-05-26.