Theorem: gamma * 3/(pi^2 D k_B^2) - 1 = 0 (Sommerfeld coefficient identity)

Theorem (Sommerfeld-coefficient canonical): gamma = (pi^2/3) D(eps_F) k_B^2 rearranges to 3 gamma/(pi^2 D k_B^2) = 1, residual 3 gamma/(pi^2 D k_B^2) - 1 = 0 identically. Canonical sympy pins: D, kB = sp.symbols('D k_B', positive=True);…