aidoesscience
aidoessciencefindings › Energy quantization
ValidatingOracle-validated

Energy quantization

Why is energy quantized? When a quantum particle is trapped, can it have any energy — and what decides which energies are allowed?

Measured by the lab
6.0247e-20
Known value
6.0247e-20
Relative error
4.05e-6

Units: J — E₁ = h²/(8 m_e L²) = 0.376030 eV for an electron in a 1 nm infinite square well (Schrödinger 1926; Griffiths §2.2)

▶ Run this simulationRead how it works

The finding

Energy quantization: no-oracle → validated — the ladder weighed by a clock, not an eigensolver: with NOTHING coded but i ψ_t = −½ψ_xx + Vψ (Crank–Nicolson, exactly unitary) and a wavepacket, the energy levels appear as spectral LINES of the autocorrelation C(t) = ⟨ψ(0)|ψ(t)⟩, and an electron in a 1 nm hard-wall box returns E₁ = 0.376032 ± 1.5e-6 eV vs the exact h²/(8m_eL²) = 0.376030 eV (rel 4.1e-6); the n² law emerges from a FREE fit (p = 1.999996 ± 2.5e-6), E₁ ∝ L⁻² to a free-fit q = −2.000000, and the SAME machinery on a harmonic trap returns the even ladder a = ω AND the zero-point b = ω/2 from one free linear fit — Planck's universal equal-step ladder dies at 393,786 SE on the box, and his original no-zero-point spectrum dies at 16,129 SE on the oscillator; between the box's lines the spectrum is EMPTY (mid-gap power 1.3e-4 of the weakest line) — the classical continuum has nowhere to live

Method

Dynamical, blind spectral recovery — deliberately independent of the module's method (the module diagonalizes the discretized Hamiltonian by Sturm-sequence bisection; the oracle never solves an eigenproblem). The generator integrates the time-dependent Schrödinger equation i ψ_t = −½ψ_xx + V(x)ψ (ℏ = m = 1) by Crank–Nicolson on a Dirichlet grid (complex Thomas solve with a precomputed constant LU — exactly unitary), launches a Gaussian wavepacket, and records the autocorrelation C(t) = ⟨ψ(0)|ψ(t)⟩ = Σ|c_n|²e^{−iE_n t} with 1% complex measurement noise. Energy levels are found BLIND: Hann-windowed 4x-zero-padded FFT peak detection (log-parabola interpolated), then iterative refinement — per line, subtract the complex least-squares model of every other line (Gram via exact geometric sums), demodulate, and fit phase vs time across 16 Hann-weighted segments with incremental unwrap (segment increments < π; no per-sample unwrap, no cumulative drift). The Crank–Nicolson phase map is deconvolved exactly per line (E = (2/dt)·tan(ω dt/2) — a property of the scheme, not of the law). Levels are indexed by RANK ORDER of energy (node count), never by the law; the n² exponent, 1/L² scaling, even harmonic spacing and ½ zero-point all come from FREE fits. No eigensolver, no n²π²/2L², no ω(n+½) anywhere in the generator or estimator. 9 gates in scripts/qwell-derisk.mjs (~32 s); known values in scripts/oracles/qwell.reference.json, loaded ONLY to score; tamper ⇒ exit 1. ?world=qwell.

Measurements, controls & cross-checks

Recovered se

2.4000e-25

Worst seed rel error

1.9000e-5

Box spectrum

Exponent free fit
1.999996
Exponent se
2.5400e-6
Worst n2 ratio rel
2.4700e-5
Note
free fit of ln E_n vs ln n over the 5 lowest lines, indexed by rank only — the n² law is in the data, not the code; 6 seeds of packet position/width/momentum jitter + 1% noise on C(t)

L sweep

Exponent
-2
Note
free power fit over L ∈ {0.7, 1, 1.3, 1.6} returns q = −2.000000 — halve the box, quadruple the confinement energy

Harmonic ladder

Spacing
0.999949
Spacing se
1.6600e-5
Zero point
0.50004
Zero point se
3.1000e-5
Omega
1
Note
one FREE linear fit E_k = a·k + b over the 6 lowest lines returns BOTH the even spacing a = ω and the zero-point intercept b = ω/2 — the particle is never at rest; ω-sweep {0.5, 1, 2}: spacing tracks trap stiffness with free log-log slope 1.000000

Discreteness

Midgap power rel
0.000131
Note
the strongest spectral power in the central half of the E₁–E₂ gap is 1.3e-4 of the weakest gated line — the spectrum is lines, not a band

Rival equal step ladder

Rival exponent
1
Measured exponent
1.999996
Z
393786
E5 over E1 measured
24.9998
E5 over E1 rival
5
Note
Planck's 1900 quantum applied as a universal rule (E_n = n·ε for every trap) predicts the box levels spread linearly; the box refuses at ~394,000 SE — the SHAPE of the trap sets the quantization law, which is what Schrödinger's equation delivers and the old quantum theory could not

Rival no zero point

Rival intercept
0
Measured intercept
0.50004
Z
16129
Note
Planck's original ladder started at zero; the free fit's intercept sits at ω/2, 16,129 SE away — the zero-point energy Mulliken's 1924 BO band spectra proved real

Module certificate

Summary
The derisk EXECUTES the shipped QuantumWellModule.ts (sha256-pinned, 35 exact strip pairs + 2 counted bulk passes, run via new Function under Babylon/DOM recorder stubs) and certifies a fully STATIC world in the blackbody/hydrogen frozen-proof pattern — gates I–M, zero module edits (the display was honest as shipped).
Exec bit exact
Executed init measurement === op-order statics replica BIT-FOR-BIT: all 5 box eigenvalues (E1 4.934790963858765), all 7 oscillator eigenvalues (E0 0.4999898950180409, gap 0.9998585140771867), every sample of all 12 eigenfunctions (5×600 + 7×1000), the grids, potentials and both Emax caps; all 3 thin-instance buffers (potentials+hard walls, box levels+psi, HO levels+psi) bit-exact vs replica geometry — the ladders are drawn at exactly the measured eigenvalues; executed _countBelow/_eigenvalue/_thomas/_eigenvector/_solve match the replica at fresh probes including a tilted well never displayed.
Frozen proof
600 engine render() calls: the module declares NO fixedUpdate, ZERO DOM writes after init at every call, all 3 buffers frozen (same reference, setBuffer called once), HUD text and chart SVG equal to exact template replicas and to a fresh executed _buildChart(); style/HUD/chart sha256 cert-pinned. Statics-by-design: the init measurement IS the whole display.
Earned
TRIPLE rounds-onto-label — every shown number equals its theory-label string at display precision while the underlying double differs: shown ratios '1.00, 4.00, 9.00, 16.00, 25.00' (ratio5 double 24.99863, delta -1.37e-3), shown gap '1.000' (double 0.99985851, delta -1.41e-4), shown E0 '0.50' (double 0.49998990, delta -1.01e-5). Fresh-physics probes through the EXECUTED solver: E1(L=2)*4/E1(L=1) deviates 1.7e-12 (the FD box spectrum scales EXACTLY as 1/L^2 at fixed N — a closed identity of the tridiagonal, so the probe floor is float rounding); omega=2 trap returns gap/2 to 8.1e-5 and zero-point/1 to 4.0e-5; node theorem holds — psi_k has exactly k sign changes for all 12 displayed states (rank IS node count).
Priced
Display bias priced with NO free constants, two techniques: (1) BOX EXACT — the shipped matrix has closed spectrum lambda_n = (1-cos(n*pi/(N+1)))/dx^2; worst |E_n - lambda_n| = 2.88e-11 against a 1.6e-10 derived float floor (bisection bracket 1e-11 + eps*||T||), so the entire on-screen ratio bias (worst 1.37e-3 absolute, invisible at 2 dp) is 100% the priced FD discretization gap (residual 1.45e-10 < 1.05e-8 propagated floor); (2) HO RICHARDSON — two-level defect ceilings derived in-run from the executed solver at N=1000/2001/4003 (exact dx halving): worst level bias 8.59e-4 <= its ceiling, gap bias 1.41e-4 <= 1.45e-4, zero-point bias 1.01e-5 <= 1.01e-5 (the |ext1-ext2| term provably covers the alternating h^4 tail). Wall truncation evidenced: worst edge |psi|/max 2.15e-13. The TWO law-fed elements — the display-axis caps Emax_box = (5.6)^2*pi^2/2 and Emax_HO = 7.6 — are disclosed, bit-verified against their expressions, and proven geometry-only (never read in the HUD/chart numeric path); the eigensolver block is pure linear algebra (no pi, no closed form, no toFixed; Sturm pivot + node-count bisection + inverse-iteration landmarks asserted present); no answer digits in the scanned code, planted-violation self-test green.
Tampers
Three surgical hand-tampers, restored by hand: known_value*1.02 → ONLY gate B fails, exit 1, recovered unchanged; module_sha256 first-hex flip → ONLY gate I fails; QWELL_ULP=1 (post-strip 1-ulp on one solver line, sha still pinned) → bit-exact gates I/J/L fail while PRICED gate M stays green (box FD gap rises to exactly the predicted eps*inv ≈ 9e-11, still under the derived floor) — and gate K stays green because a 1-ulp shift truly changes NO displayed string (2-3 dp), which is precisely why the doubles-level pins in J/L exist.
Gates
13/13 (A-H science + I-M certification) in ~39 s

What it reduces to

The two canonical bound-state spectra of introductory quantum mechanics: the infinite square well E_n = n²h²/(8mL²) (levels spread as n², scale as 1/L²) and the harmonic oscillator E_n = ℏω(n+½) (even ladder with a zero-point) — Schrödinger 1926, Griffiths §2.2–2.3. It VALIDATES, not derives: the generator is bare unitary time evolution of the Schrödinger equation plus a trap shape; the closed forms appear only in the reference file, loaded to score. Non-circularity is checked five ways in-script: the recovery is dynamical (no eigensolver — a genuinely different method from the module's Sturm diagonalization), levels are indexed by rank order rather than by the law, the n²/L⁻²/even-spacing/half-quantum structure emerges from four FREE fits (gates C, D, E, E′), the scoring self-test flips the exit code with the recovered value unchanged (gate H), and the decisive controls are PHYSICAL and historically real — Planck's universal equal-step ladder and his original no-zero-point spectrum are both fitted by the same machinery and rejected at 393,786 and 16,129 SE. This closes the lab's bound-state gap in the quantum arc: blackbody/photoelectric/franckhertz established the quantum of ACTION from counting and currents; qwell establishes that a TRAP turns it into a discrete energy ladder whose spacing law is set by the trap's shape — the fact hydrogen (already validated) inherits as E_n ∝ −1/n².

Module systematics

CERTIFIED (refiner honest-module run, zero module edits): the derisk now EXECUTES the shipped module — see result.module_certificate. The on-screen box ratios, oscillator gap and zero-point are MEASURED from the module's own Sturm-sequence diagonalization (init === replica bit-for-bit; display frozen over 600 engine calls). The old estimate that discretization bias reaches ~5.7e-5 relative on the box's 5th level is now MEASURED by execution: 5.5e-5 relative (1.37e-3 absolute on E5/E1 = 25) — the estimate was honest to within 4%, and the bias is 100% priced by the tridiagonal's exact closed spectrum. The oscillator ladder carries -1.41e-4 (gap) and -1.01e-5 (zero-point) FD bias, bounded by in-run Richardson defect ceilings. The law ANNOTATIONS next to measured numbers ('→ the n² law', '= ½ω') and the two law-fed display-axis caps (Emax rulers) are disclosed and proven geometry-only; the module/oracle method independence (stationary Sturm eigensolve vs dynamical Crank-Nicolson spectroscopy) is unchanged and remains the oracle's non-circularity backbone.

Notes

Zero gate retries — prototyped in the scratchpad first, which caught two estimator traps at design time: (1) plain demodulated-mean amplitude estimates leak badly at the oscillator's gap ω (2/(ΔE·T) ≈ 5%), poisoning subtraction and collapsing adjacent lines — fixed with exact complex least-squares amplitudes (Gram matrix via geometric sums) before every subtraction; (2) per-sample phase unwrap of a noisy residual random-walks into 2π slips that drag the slope fit across a whole level gap — fixed by fitting phase across 16 Hann-weighted SEGMENT means (increments guaranteed < π), the segment-level analogue of doppler's incremental cycle indexing. Estimator-checklist compliance: no division near zero, no naive argmax (lines are log-parabola-interpolated FFT peaks refined by phase-slope fits), untracked ghost lines removed by a 1%-amplitude cut, and the recovery is noise-limited at the Cramér–Rao bound (measured E₁ SE matches the CRB for 1% complex noise on 8000 samples over T = 40). Tolerances set at ≈6–25× the measured statistics (gate table in the reference note), five orders tighter than the 100% / half-quantum rival gaps. Hand tamper (known ×1.079) ⇒ exit 1 with recovered unchanged, restored by hand; in-script scoring self-test gate H green. Runtime ~32 s, 9/9 gates. CERTIFICATION RUN (refiner, zero-edit): gates I-M appended to the derisk; QWELL_CAL calibration dump; three tamper self-tests surgical (known→B only; sha→I only; post-strip ulp→I/J/L with M green and K string-insensitive as designed). New cert phenomena: (1) the earned L-probe has a FLOAT-floor tolerance (1e-9) because the FD box spectrum obeys 1/L^2 EXACTLY at fixed N — the discretization cancels in the ratio; (2) first world priced by TWO techniques at once — exact closed FD spectrum (box) + in-run two-level Richardson defect ceilings (oscillator); (3) triple rounds-onto-label in one HUD.

Confidence & reproduction

Confidence
high
Validation
derisk-pass
Re-run the check
npm run derisk -- qwell (scripts/qwell-derisk.mjs)
Oracle
scripts/oracles/qwell.reference.json

Sources

E. Schrödinger, 'Quantisierung als Eigenwertproblem (Erste Mitteilung)', Annalen der Physik 79, 361–376 (1926); D. J. Griffiths & D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., §2.2–2.3; M. Planck, Verh. Dtsch. Phys. Ges. 2, 237 (1900) — the equal-step ladder (rival 1) without a zero-point (rival 2); R. S. Mulliken, 'The band spectrum of boron monoxide', Nature 114, 349–350 (1924) — zero-point energy is real.

One finding from the lab's 104 catalogued results — each an experiment run end to end by an AI: a question, a method, measured data, a control, and a confidence.