Mills-Randal: Mathematician review (Step 2)

Reviewer: PhD mathematician (checks derivations, notation, units, hidden assumptions). Date: 2026-10-09. Subject: R. L. Mills, The Grand Unified Theory of Classical Physics (GUT-CP), project copy = the "December 10, 2018 Edition" (1,876 pp, file name "2020 Ed Web"). The inventor has a commercial stake (BlackLight Power, now Brilliant Light Power, hydrino technology). Everything below about what the book says is the author's own account; the checks are mine.

Tags: [read] from the project documents or the deep reads (which I re-checked where stated), [calc] my own computation in the project venv (mpmath/sympy/numpy; scratch scripts, not deliverables), [web] with URL, [memory] with confidence.

Source gap and coverage (repeated from claims.md)

About 40 percent of the book's pages were deep-read. Digest only: Ch.3-4, Ch.9 (except Appendix VI), Ch.11 (the hydrogen molecule and the 4.478 eV bond), Ch.15 pdf 801-1082, Ch.17-21, 23, 24-30 (statistics, superconductivity, quantum Hall, Aharonov-Bohm, pair production, positronium), Ch.31, 33-35 and the first half of Ch.32. Claim 7 therefore rests on digests and [memory] alone, and so does any statement below about Ch.11. Displayed equations do not survive text extraction; every number the deep readers checked came from page images. I did not re-render page images; I re-computed the deep readers' numbers from their printed formulas and found them to reproduce (details per claim).

Verdicts at a glance

Corrections to Step 1 (facts.json, claims.md, deep reads)

Each item: document and passage, what is wrong or incomplete, the evidence.

1. claims.md, "Observations ... Hydrino evidence" and deep read pp.1214-1252 (EPR): the "average error 0.097 G" is the book's number, not the paper's. The Hagen and Mills 2022 paper (IJHE 47(56), 23751; Investigator's file) prints in its Table I eleven main-line positions at 9.82029 GHz with differences of 0.74, 0.30, 0.31, 0.23, 0.44, 0.30, 0.20, 0.09, 0.29, 0.58, 1.22 gauss and the summary "Average (st. dev.): 0.004 +/- 0.550", and says "The GUTCP fit to the X-band spectrum of H2(1/4) is semi-quantitative" [read: research/files/Hagen-Mills-2022..., Table I and "Remaining questions"]. The 0.004 is a signed mean; the mean absolute difference of the 11 lines is 0.43 G [calc]. The book's 0.097 G (Tables 16.49-16.50) is a different statistic over an unstated peak set. Neither should be quoted alone. The paper also says the fluxon (sub-line) separation "is not a constant" and that the number and intensities of fluxon lines "are presently not understood" [read].
2. claims.md, Source gap: the Hagen-Mills paper is now in hand (already recorded by the Investigator). The Bruker sentence has a second statement (Billerica, MA, "two different spectrometers") but both statements come from papers Mills co-authored and cite Mills's own preprint [read]. Addition, not a reversal.
3. Deep read (Ch.1), item "Eqs. 1.15-1.20: wavelengths conflict; vanishing depends on lambda_n = r_n": confirmed and extended. I confirm with sympy that a rotating spherical-harmonic charge wave rho = P_l^m(cos theta) cos(m(phi + omega t)) on a sphere of radius R satisfies the stated two-dimensional wave equation only if the wave speed is v = m omega R / sqrt(l(l+1)). For l = m = 1 that is omega R / sqrt 2; l = m = 2 gives 0.816 omega R; l = m = 3 gives 0.866 omega R [calc]. Mills sets v = omega R (his Eq. 1.35) and v = omega R is never equal to m omega R / sqrt(l(l+1)), because l(l+1) = m^2 has no solution with l >= 1 and |m| <= l. So the l >= 1 "orbitsphere" charge waves (the basis for excited states, Eq. 1.12, Eq. 48 of Appendix I) are not solutions of his own Eq. 1.56, whatever the radial factor. This is a stronger statement than the deep reader's and independent of the radial delta function.
4. Mills's 2008 reply to Rathke (Investigator's file) is half right. Mills notes that Rathke (New J. Phys. 7, 127, 2005) prints the wrong sign in front of the time-derivative term in his Eq. (9) [read: Millsian-Mills-response, p.90 and footnote 1]. I confirm: separating rho = f(r) A(theta, phi, t) gives (1/r^2) Delta A - (1/v^2) A_tt = 0, with a minus sign. Mills's conclusion that "his other comments ... are made moot" does not follow: with the correct sign the rotating-harmonic failure above remains, and Rathke's other point (a radial delta function is not a solution of the radial equation (r^2 f')' = 0) is Mills's own admission: the reply says the delta function "is not a radial solution of Eq. 4; rather, it defines a constant charge density on a spherical shell" [read]. That is a boundary-condition choice, not a derivation from the wave equation. Anyone citing Rathke should cite the sign slip and the surviving points together.
5. Deep read Ch.5, section 7a (Schroedinger: "Non-integer n gives solutions that blow up ... so n = 1 is the lowest eigenvalue"). Correct for the regular solutions, but incomplete in one mathematically important way, which the Investigator's arXiv files bring in. For the l = 0 radial equation the two solutions near the origin behave like r (regular) and like a constant (u = r psi, i.e. psi ~ 1/r, irregular). Both are square-integrable near r = 0 in three dimensions (the integral of r^2 * (1/r)^2 dr is finite), so the l = 0 radial Hamiltonian is "limit-circle" at the origin and is not essentially self-adjoint; it has a one-parameter family of self-adjoint extensions, and some extensions contain an extra, very deep bound state [memory, high on the limit-circle statement; this is how Naudts, Khelashvili and Nadareishvili argue]. Standard quantum mechanics selects the regular boundary condition. Dombey (arXiv physics/0608095) shows the deep Klein-Gordon state "disappears if the nuclear charge distribution is taken to have an arbitrarily small but non-zero radius" [read: abstract], and Khelashvili and Nadareishvili conclude for the Dirac equation that the extra solution "must be ignored" because the norm is not conserved and orthogonality fails [read: abstract]. So the deep reader's statement stands for physical (finite-size) protons, but "no such state exists" is a statement about a boundary-condition choice and should say so.
6. claims.md claim 6 and the deep read (Ch.32): "H0 = 78.5 km/s/Mpc ... compared with measurements". I reproduce that Mills's expression reduces to H = 8 pi^2 / T = 4 pi c^3 / (G m_U) independent of time, giving 78.27 km/s/Mpc for m_U = 2e54 kg and needing m_U = 2.32e54 kg for 67.4 [calc]. The number is therefore set by the one chosen quantity m_U, which the book says was chosen "for internal consistency" with H0 and other data [read: deep read, Ch.32 section 3]. Step 1 reported the comparison; it should also carry this dependence.
7. Deep read p.1218 (Ch.5 section 5) and claims.md: "Ar+ 27.62 vs 27.63". Trivial; skip.
8. Deep read (Ch.2) "Lamb shift 1057.09 vs 1057.845 MHz ... about 80 standard deviations" [memory, high] is now [web] checked for the uncertainty: the Lundeen and Pipkin value is 1057.845(9) MHz (PRL 46, 232, 1981), "the most precise direct determination" per a review; other determinations are 1057.842(12) and 1057.852(15) MHz; the mean of direct measurements is 1057.8439(72) MHz [web: https://arxiv.org/pdf/0809.2985]. Mills's 1057.09 MHz is about 0.75 MHz low, 80 times the quoted error. The "natural linewidth" argument (100 MHz) is about how sharp the line is, not how well its centre is known.
9. Deep read (Ch.2) "muonium uncertainty '12 ppm' should be ppb": confirmed [web: Liu et al., PRL 82, 711 (1999), 4 463 302 765(53) Hz, 12 ppb, via https://www.doi.org/10.1103/PHYSREVLETT.82.711 listing]. Mills's 4.46330328 GHz is 515 Hz, about 10 standard errors, above it.
10. Deep read Appendix III, open item (CERN II muon g-2). A search finds the Bailey et al. CERN result (Nuovo Cimento 9A, 1972) with a value consistent with the recalled 116616 x 10^-8 and agreement with theory to (240 +/- 270) ppm; and that the third CERN experiment (1974-79) is the one with electrostatic quadrupole focusing [web: https://archives.web.cern.ch/?p=479 and the 1972 paper listing]. I could not confirm from the sources found that the 1972 ring had no electric focusing; the gamma = 12 weak-focusing design suggests so [memory, medium]. Treat "a magnetic-focusing measurement already differs from the electron value" as probable, not verified. Either way the circularity of Appendix III Eq. 26 (below) does not depend on it.
11. Deep read pp.121-158, Section 2.3, "eleven-figure agreement". Confirmed and tightened (see claim 1b). No correction to the deep read; correction to the book's phrase only.
12. Digest corrections already listed in claims.md (digest .28, .110, .111, .121) are consistent with what I read; none contradicted.
13. Not a correction but a gap in Step 1: the digest-only chapters carry two quantitative statements I could check from the digests alone and which the claims.md list omits: Ch.25 gives Tc = 30.8 K "for conventional three-dimensional metallic superconductors" against 23.2 K for Nb3Ge, and three values (7 K, 22 K, 180 K) for ceramic oxides by dimension [read: digests p.1554]. Conventional elemental superconductors have Tc below about 10 K and BCS-type compounds top out near the A15 values [memory, high]; a single formula giving a spread of 7-180 K by choosing a dimension is a free discrete choice. Flagged under claim 7.

---

Claim 1. Classical electron, closed-form atoms

1a. Is the mathematics internally consistent?

What is consistent (checked): the Bohr-type force balance (m_e v^2/r = e^2/(4 pi eps0 r^2) with v = hbar/(m_e r)) gives r = a_0 and T = 13.606 eV; Box 1.1 numbers reproduce; the magnetostatic boundary-value solution (Eqs. 1.146-1.153) is correct for the current assumed; the Dirac 1s formula (Eq. 1.294) is exact [read; calc].

Where it is not:

• hbar is an input, not a result. Eq. 1.37 integrates a speed that was set by Eq. 1.35 so that m_e v r = hbar, and returns hbar [read: deep read Ch.1, section 3]. The hbar/2 of "spin 1/2" comes from a construction of great-circle currents whose normals are spread over one hemisphere; a full-sphere spread gives zero. Condition (4) states the target in advance [read].
• mu_B needs an extra postulate: a trapped resonant photon of angular momentum hbar precessing at pi/3, with the angle fixed by an asserted energy balance (footnote 8). The factor N = 3/2 (Eq. 1.140) turns the natural 2/3 mu_B of a rigid rotating shell into 1 [read; I re-derived: dipole moment of a uniformly charged shell turning at omega = hbar/(m_e r^2) is (2/3) mu_B].
• Units: factor 4 pi in the field, (4 pi)^2 in the energy. Eqs. 1.132-1.162 combine Gaussian-form fields with the SI mu_B. I confirm by direct integration that the SI field energy of a shell carrying moment m is mu_0 m^2 / (4 pi R^3), so Eq. 1.162 (4 pi mu_0 mu_B^2 / r^3) is 16 pi^2 = 157.9 times too large; at r = a_0 it is 0.0572 eV against 3.62e-4 eV [calc]. The E_mag of Eq. 7.46 and 8 pi mu_0 mu_B^2 / r^3 are used downstream (helium, hydride ion), so the error propagates unless the later chapters re-scale; I did not trace every use.
• Nonradiation is not proved for l >= 1. The l = 0 spin shell is static and nonradiating by standard electrodynamics (true, trivial). For l >= 1:
• the wave equation fails for the rotating harmonics unless v = m omega R/sqrt(l(l+1)) (Corrections, item 3);
• the Fourier "proof" uses a zero of one factor (sinc) to kill a convolution, which does not follow; the correct radial factor is j_l(k r_n), and j_1(2 pi) = -0.159, j_2(2 pi) = -0.076 are not zero at Mills's own k r_n = 2 pi [read: deep read App. I; I re-evaluated these Bessel values: j_1(2 pi) = sin(2 pi)/(2 pi)^2 - cos(2 pi)/(2 pi) = -0.1592, correct];
• for hydrogen n = 1 the surface speed is alpha c, so k r_1 = alpha and j_0 = sin(alpha)/alpha = 0.999991, not zero;
• the Poynting-vector version sets the electric multipole coefficient to zero by dropping a charge term that the book's own continuity equation makes nonzero; for l = m = 1 the charge density is a rotating dipole of size e r_n/3 which radiates P = p^2 omega^4/(6 pi eps0 c^3) ~ 5e-9 W at hydrogen scale [read; I re-did the arithmetic: p = 2.8e-30 C m, omega = 4.13e16 rad/s gives 5.1e-9 W, in agreement];
• the three definitions of wavelength (2 pi r_n = lambda_n; lambda_n = r_n; s = v T_n = r_n = lambda_n) cannot all hold [read].
• Self-force. The book's force balance omits the electrostatic self-repulsion of the shell, justified by postulates (zero thickness, "perfect conductor", "no remote action of any surface point on another"). Purcell's own derivation (reproduced in Appendix II) gives an outward pressure sigma^2/(2 eps0) independent of layer thickness. Including it halves the net inward Coulomb force, and the same balance then gives r = 2 a_0 instead of a_0 [calc: sympy, ratio exactly 2]. The ground-state radius of the whole theory therefore depends on the postulate that a surface charge does not push on itself. The field-energy route (Eq. 56 of App. II) silently includes the shell's self-energy e^2/(8 pi eps0 r) = 13.6 eV at a_0, which the force route excludes; the two agree only because that energy numerically equals the kinetic energy at a_0 [read: deep read App. II].

1b. Consistent with established experimental results?

I recomputed the quantities the deep readers tabulated.

Electron g factor. Mills's Eq. 1.228: g/2 = 1 + alpha/(2 pi) + (2/3) alpha^2 (alpha/(2 pi)) - (4/3)(alpha/(2 pi))^2. With CODATA 2018 alpha it gives 1.00115965246723; the measured value is 1.00115965218059 (Fan et al. 2023 [memory, high]). Difference 2.87e-10, about 2,000 times the experimental uncertainty (1.3e-13) [calc]. In powers of alpha/pi, Mills's second- and third-order coefficients are -1/3 = -0.3333 and pi^2/3 = 3.2899; QED's are -0.328479 and +1.181241 [calc for Mills's; QED values: memory, high]. The two errors, -2.6e-8 and +2.6e-8, nearly cancel, which is why the sum lands near experiment. The 1987 test used alpha known to 8e-9 relative in g/2, so the claim "eleven-figure agreement" was agreement inside a +/-7e-9 band [read: deep read; the propagation is a one-line derivative]. The first term (alpha/2 pi) is Schwinger's (1948). Several modelling choices set the coefficients (2 pi radius scaling, replacing Phi_0/area by the applied B, v = c in one term and v = alpha c in another, internal field only for the electric term, sign of the dissipation term) [read].

Fine structure 2P1/2-2P3/2. Mills's leading term is alpha^5 (2 pi)^2 sqrt(3/4) m_e c^2 / 8 = 10,927.01 MHz; the Dirac result is alpha^4 m_e c^2/32 = 10,949.28 MHz. Their ratio is 8 sqrt(3) pi^2 alpha = 0.99797 [calc]: the book's formula has a different power of alpha, and matches Dirac to 0.2 percent only because 8 sqrt(3) pi^2 alpha is numerically near 1. Since fine structure scales as (Z alpha)^4 in hydrogen-like ions, Mills's Z-scaling is not given in the pages read.

Lamb shift, muonic hydrogen, 21 cm, muonium. Arithmetic reproduces (deep reader's recomputation; I also confirmed the recoil term uses the 10.2 eV Lyman-alpha energy and that without it the result is 1.3 percent low). Residuals: H Lamb 0.75 MHz (80 sigma [web]), muonium 515 Hz (10 sigma [web]), 21 cm 1.07 kHz, 0.75 ppm. The geometric factor multiplying the stored magnetic energy is (1 + (2/3)^2 + alpha/4) for hydrogen and (1 + (1/3)^2 + alpha) for muonium; swapping them spoils each result by 0.42 percent [read: deep read]. The muonic-hydrogen result depends on lepton-counting coefficients (1, 2, -3) in (2 pi alpha)^2 + 2(2 pi alpha)^3 - 3(2 pi alpha)^4; changing the last by one unit moves the answer 94 GHz, more than the 75 GHz proton-radius discrepancy it claims to explain [read].

Relativistic ionization energies. Eq. 1.294, IP = m_e c^2 (1 - sqrt(1 - (Z alpha)^2)), is Dirac's point-nucleus 1s result [memory, high]. Table 1.5's 39 rows reproduce to every printed digit; the residual against "experiment" grows like Z^4 and is the Lamb shift plus nuclear size, which the theory does not compute. For U91+ the book's residual is 431 eV against a measured Lamb shift of about 460 eV [read; memory, high]. So Table 1.5 is agreement with Dirac theory, not independent support.

Multi-electron atoms (Ch.7, 10). Helium: the force balance gives r_1 = a_0 (1 - sqrt(3/4)/2) = 0.566987 a_0 [calc, exact]; the magnetic factor 1/Z is argued, not fitted (with factor 1 the radius would be 0.134 a_0 [calc]). The ionization energy formula for helium imports a bracket from the muonium chapter. Neutral atoms Al-Ca: errors 0.03-0.65 percent. Isoelectronic ions: up to 4.6 percent (Ar3+), pattern smooth in Z, with many "experimental" high-Z values being themselves interpolations (the table footnote says so) [read]. The per-atom integers A, B (Table 10.18, 10.21) are chosen "to give an energy minimum"; in the 2p series I accept the deep reader's finding that Mills's own energy-minimum rule would pick the configurations he rejects for neutrals [read; I did not repeat that test]. The hydride ion: reproduced at 0.754175 eV [calc], against Lykke's 0.754195(19) eV (about 1 sigma), but with an ordinary reduced mass the same formula gives 0.754104 eV [calc], 5 sigma off, and with infinite nuclear mass 0.753643 eV. The special reduced mass m_e m_p/(m_e/sqrt(3/4) + m_p) is therefore load-bearing at the 5-sigma level; it is introduced in Ch.7 and not derived in the pages read. Mainstream ab initio calculations also reproduce the hydride affinity [memory, medium-high].

1c. Has it been tested? By whom? Replicated?

Tests of the numbers: the deep readers' recomputations (arithmetic reproduces everywhere I checked) and Booker's chapter-by-chapter verification for the company (Investigator: he reports reproducing Tables 1.1-1.5, 7.1-7.2 "to a very high degree of accuracy"). I regard these as consistent: nobody finds that the book's equations fail to give its printed numbers. That is a test of arithmetic, not of the postulates, and Booker was also the author of a company power-validation report (interested party; compensation not stated in pages read). Tests of the physics: the comparison with experiment is made by the author; the CODATA g-2 value and the 2018-era alpha were never used by the book (it uses 1987-era constants). No outside group is known to me to have adopted the framework except Selke (2025, extension of the force-balance ionization energies to heavy ions, Investigator; not read).

1d. If true, what in mainstream physics would have to change?

Quantum electrodynamics would need to be an accidentally good numerical fit: the g-2 series is computed to five loops and agrees with atom-interferometry alpha at 10^-13 [memory, high], which Mills's three-term closed form misses by 2.9e-10. The Lamb shift, hyperfine and Dirac fine structure would be reproducible by a different mechanism with different alpha-powers. The whole of atomic spectroscopy (selection rules, the Rydberg-atom n^2 scaling of size, which Mills replaces by n a_H) would need re-derivation; the book's own Stark-effect energies reproduce the textbook first-order Stark levels exactly [read: deep read Ch.2], so not every claim differs in numbers.

Verdict, claim 1

• Nonradiation of l >= 1 states and the no-self-force step: Internally inconsistent (specific defects above).
• g factor, Lamb, fine structure, muonium as closed-form predictions: Tested, not supported (correct arithmetic; residuals 10 to 2,000 times the experimental uncertainty; choices of factors per case).
• Ionization energies of 2-20 electron atoms, "exact" solutions: Tested, not supported as exact and parameter-free; reasonable as a post-hoc semi-empirical scheme for neutrals.
• One-electron relativistic ionization energies: Established (as Dirac's result).

---

Claim 2. Hydrinos

2a. Is the mathematics internally consistent?

What checks: the energy ladder. With Z_eff = p + m the shell radius is a_H/(p+m), kinetic energy +13.6 (p+m)^2 eV, potential -27.2 (p+m)^2 eV, binding 13.6 (p+m)^2 eV; the released energy is [(p+m)^2 - p^2] x 13.6 eV, of which m x 27.2 eV goes to the catalyst and the rest [(p+m)^2 - p^2 - 2m] x 13.6 eV is the "continuum edge" [read; calc]. Examples: p = 1, m = 1: 40.8 eV released, edge 13.6 eV (91.2 nm); m = 2: 108.8 eV, edge 54.4 eV (22.8 nm); m = 3: 204.0 eV, edge 122.4 eV (10.13 nm) [calc]. Tables 5.1 and 5.2 add up. The H2(1/p) numbers reproduce (31.678 eV total, 4.478 eV bond with E(2H) = 27.20 eV; with 2 x 13.598 the bond is 4.481 eV [read; calc by the deep reader, simple]).

What does not:

• The force balance does not select the integer. Substituting r = a_H/Z_eff into hbar^2/(m_e r^3) = Z_eff e^2/(4 pi eps0 r^2) holds identically for any real Z_eff (the factors cancel), so Eq. 5.28 is satisfied for 1.5 or 2.7 as well as for 2 and 3. The integers come from an additional postulate: the "trapped photon" adds an integer multiple of the proton's central field (Eq. 5.27) [read: deep read Ch.5, 2b]. It is the same algebra as a hydrogen-like ion of nuclear charge Z, with the extra charge attributed to the electron's own photon field.
• The delta shell and the rotating harmonics are not solutions of the stated wave equation (Corrections 3-4; Rathke's radial argument).
• The "27.2 eV per unit" hole is justified by counting multipole degeneracy p^2 and assuming one multipole must be excited (Eqs. 5.44-5.46) [read].
• The lower limit is inconsistent within the book. The relativistic radius at p = 137 is (a_0/137) sqrt(1 - (137 alpha)^2) = 8.853e-15 m (Eq. 5.113) [calc, reproduces the book]; Eq. 6.14 (alpha a_0/(2 pi)) gives 6.15e-14 m; Eq. 6.13 gives 3.86e-13 m; the proton radius is 0.84e-15 m [memory, high]. The pages do not reconcile them [read: deep read]. And the 2025 Mills and Hagen paper states p <= 126 ("limited by the energy in the electric field of the proton ... whereafter the proton and electron decay") [read: Investigator, FeOOH paper p.1]; the book (p <= 137) and the paper disagree by 11 states, and the physical reason changed (speed reaching c versus field energy). Also note 137 here is the value at which a point nucleus Dirac 1s binding energy m c^2(1 - sqrt(1 - (Z alpha)^2)) becomes complex (Z alpha = 1 at Z = 137.036); Mills's p limit coincides with the Dirac point-nucleus singularity at Z = 1/alpha [calc/memory, high]. The radius at p = 126 is 1.65e-13 m [calc].
• Table 5.2 (catalysts) is arithmetic-consistent, but the criterion "approximately m x 27.2 eV" has no stated tolerance. The deviations of the 48 numeric entries from the nearest multiple of 27.2114 eV have median 0.58 eV and mean 0.86 eV, maximum 2.84 eV [calc, from the deep-read table]. A random enthalpy lies within 1 eV of a multiple of 27.2114 eV with probability 2/27.2 = 7.4 percent (within 3 eV: 22 percent), so with several choices of how many electrons to ionize for each element the table is what chance supplies; a handful of very close entries (He+ 0.005 eV, Ti 0.02 eV, Pb 0.04 eV) are the expected tail. The table cannot, by itself, test the resonance claim; a non-resonant mismatch of 0.3-3 eV cannot be absorbed by phonons in a gas (the book's own footnote says phonons absorb the mismatch in solids).

2b. Consistent with established experimental results?

Hydride ion 0.75418 eV and H2 -28.0 ppm NMR, vibration, bond energy: these are ordinary hydrogen quantities that standard theory also reproduces; they do not discriminate (Mills's hydrino claims are the p-scaled analogues, for which no measured value appears in Ch.5). The book's Rydberg-atom radius n a_H (instead of n^2 a_0) conflicts with Rydberg-atom physics [memory, high] (Physicist's ground). The hydrino hyperfine table (Table 6.1: 34.3 GHz for p = 2, l = 0, etc.) is internally consistent (each row E/h = f, f/c = wavenumber) and the predicted spontaneous rate for p = 2 is 4.05e-11 /s [calc by deep reader, reproduced] (lifetime ~780 years): extremely weak, with no observation reported.

2c. Tested? By whom? Replicated?

Record as presented to me, all with the stake of the author's company noted: every experimental citation for hydrino evidence in Ch.5 ([1]-[9], [28]) has Mills as an author; one is "submitted" [read: deep read]. The Investigator adds: NASA Lewis TM-107167 (Niedra et al., 1996: apparent excess heat up to 11 W on 60 W input, gains 1.06-1.68, with the authors saying recombination of dissolved H2/O2 is an ordinary explanation they could not exclude); Thermacore SBIR Phase I for the USAF (1993-94; 25 W; Lehigh ESCA "55 eV signature"; Mills consulted); the EarthTech 1997 negative replication (no excess within about 10 percent); the NIAC Phase I study (2002); company-hosted validation reports (Nansteel 2021: one test, local gain about 4.2; Booker 2020; Tse 2020, "hired to observe"; Glumac 2014 DSC); Hagen and Mills 2022 EPR; and the IIT Kanpur 2026 abstract announcing a replication programme with no results yet. Absence of evidence: thirty-five years of filings and reports, no sale, one non-company group beginning to replicate; the records that decide are the outside laboratories' own reports (AECL, Westinghouse, MIT Lincoln, INEL, Hokkaido, Brookhaven, Lehigh, Harvard CfA, Bruker) and the IIT Kanpur results (Investigator). I add no verdict on calorimetry (physicist/contrarian). My mathematical remarks on the EPR evidence:

• The prediction is a doublet separated by 3.9943 G (frequency independent), spin-orbit satellites at +/- m x 7.9885 G with m = 0.5, 1, 2, ..., a downfield-only magnetic correction, and sub-lines of 0.31/0.62/0.93 G with an "unbounded" integer j [read: deep read, Hagen-Mills]. A line set with a finest spacing s has mean random miss s/4: 0.078 G for s = 0.31, 0.155 G for 0.62 [calc]. So an "average error" of 0.097 G (book) on sub-line positions is at the level of chance if any sub-line combination may be assigned; the book's number needs the assignment table to weigh. The main-line comparison in the paper (11 lines, differences 0.1-1.2 G, mean absolute 0.43 G) is on a grid of 2-6 G spacing, where random placement gives 0.5-1.5 G, so it is better than chance but the authors themselves call it "semi-quantitative" and note irregular fluxon spacing.
• The paper offers an ordinary alternative that "approximately reproduces" the spectrum with two isotropic radicals (hyperfine doublets) and argues it is "highly unlikely" on chemical grounds [read]. That is a chemistry argument (Physicist's).
• Hagen is an outside EPR specialist at TU Delft (Investigator [web]); the sample comes from Mills's reactor; the interpretation and theory are co-published with Mills; competing-interest statement as printed.

2d. If true, what in mainstream physics would have to change?

The hydrogen ground state would have to be a metastable state of a Hamiltonian whose spectrum, as mathematics currently defines it, is bounded below by -13.6 eV. In standard quantum mechanics the hydrogen Hamiltonian with a physical (finite-size) proton has no such states; the only formal route is the singular-solution/self-adjoint-extension family (Naudts 2005, Khelashvili and Nadareishvili 2012-2015), which the authors themselves, or Dombey, show to be unphysical or non-normalizable in the Dirac case [read: abstracts]. A hydrino would also require stability of the hydrogen atom to be a function of the environment rather than a theorem (Lieb-Thirring-type stability bounds on the binding energy of order 20 eV exist for the non-relativistic isolated atom, per Rathke) [read: Rathke, citing a stability bound; memory, medium]. The radiation condition (Haus) is accepted by mainstream; the new content is the claim that a "trapped photon" adds integer central charge.

Verdict, claim 2

• The CP derivation of states below n = 1: Internally inconsistent (integers not selected by the force balance; wave-equation failures; self-force).
• Standard-QM version with a finite nucleus: Refuted by calculation (Dombey; Khelashvili and Nadareishvili) [read abstracts].
• Energy-release algebra: arithmetic consistent (no verdict needed).
• Laboratory evidence for hydrinos: Never properly tested on the record read: every discriminating claim (H2(1/p) rotation, p-scaled lines) rests on the inventor's or his company's reports, outside measurements exist (Hagen EPR; calorimetry by hired or consulting persons; NASA/EarthTech with ordinary explanations not excluded) but no neutral controlled test of the p-scaled signatures is on record. Not "dismissed without testing": the 2000 USPTO withdrawal and 2008 UK IPO refusal were examinations of patent applications (Investigator), not tests, and the UK IPO's "contrary to generally accepted physical law" is a legal characterization.

---

Claim 3. Hydrinos as dark matter and cosmic emission

3a. Mathematics

The numbers are internally consistent: 3481.6 eV = 256 x 13.6 eV = 16^2 x 13.6 eV; 3481.6 + 231.2 = 3712.8 eV (Eq. 5.79); 0.35625 nm printed against 0.3561 nm computed (0.04 percent); 122.4 eV for the 10.1 nm cutoff [read; calc]. The pairing "10.1, 22.8, 91.2 nm for m = 1, 2, 3 respectively" is reversed relative to Eq. 5.57 (lambda = 91.2/m^2 nm) [read: deep read]. The disproportionation power and rate examples (Eqs. 5.91, 5.104, 5.111) reproduce (8e21 /s; 100.7 kW; 144.1 kW) but only because the Foerster-type coupling factors (orientation factor, transition probability, overlap) are set to 1 by assertion, and the efficiency 1e-4 in Eq. 5.111 is chosen; 3e19 /m^3 hydrinos in Eq. 5.103 are assumed present at 0.28 Pa [read; calc]. Eq. 6.39 sets the radiative decay rate A of a hyperfine line equal to the collisional rate, which with A = 4e-11 /s would imply an absurdly low collision rate [read, observation by the deep reader].

3b. Consistent with experiment?

The 3.48 keV match is a single energy (a predicted continuum edge compared with a line); the line observed by Bulbul et al. was reported at about 3.55-3.57 keV and Boyarsky et al. at 3.52 keV, and later work did not confirm it [memory, medium]; this matter, and 511 keV and GRB origin, are the Physicist's. No predicted continuum shape or intensity is printed.

3c. Tested?

No spectrum, flux or line-shape comparison appears in the pages read; the cited independent papers (Bulbul, Boehm) are cited for facts, and the readings are Mills's.

3d. If true

Dark matter would be baryonic hydrogen in sub-ground states; the CMB-peak and BBN baryon budget would conflict with the standard fit [memory, high] (Physicist).

Verdict, claim 3: Never properly tested. Note the p <= 137 limit on which the 511 keV / gamma chain relies is withdrawn in the author's 2025 paper (p <= 126).

---

Claim 4. One method solves molecules and matter

4a. Mathematics

All checked arithmetic reproduces. The method: each bond is an H2-type prolate-spheroidal orbital; per-group energies are summed. The deep reads document where discrete choices enter: factors 0.75, 0.25 (charge splits), 0.91771 and 0.93172 (sp3 hybridization), "N factor 0.936127", the CH2 extra magnetic energy (0.14803 eV), double-bond multipliers, experimental vibrational energies in E_osc, experimental sub-radical bond energies in the CH2-CH4 and NH3 chains, experimental angles for boron bridges (85.4, 58.9 degrees), and "energy minimum" assigned per group [read: deep reads Ch.13, Ch.22]. Ionic groups reuse the neutral bond geometry by construction, so the comparison of the ion bond lengths is not independent [read].

4b. Consistent with experiment?

431 summary rows: mean absolute relative error 0.172 percent, median 0.056 percent, 317 within 0.1 percent [read: deep read; the deep reader's own tally]. Because total bond energies are sums of 25-300 eV, 0.17 percent is 0.04-0.5 eV; chemical accuracy (about 0.04 eV, 1 kcal/mol) is the usual benchmark [memory, medium]; for boron compounds (mean 0.3-0.75 percent, max 1.78 percent) the error exceeds the 0.87 eV bridge bond energy itself. Classes with 1-3 percent: cytosine, trimethylphosphine (2.3 percent, while triethylphosphine is -0.04 percent), D-ribose, 13 of 20 amino acids; 12 of the 13 and D-ribose have crystal-phase experimental values compared with gas-phase calculations. Nine rows in Table 15 have relative errors not matching (exp - calc)/exp [read]. Geometry: comparisons often against neutral basis-element bond lengths or computed (AM1) "experimental" values; water geometry compared against liquid-phase neutron data [read]. The book also reports no aggregate statistic itself.

4c. Tested?

Not blind. The book's comparisons are after the fact; the only independent verifications are Booker/Nansteel arithmetic checks (company-hosted). I know of no test in which the method predicted a molecular energy before the measurement (record read).

4d. If true

Quantum chemistry (Hartree-Fock, DFT, coupled cluster) would become a roundabout way to compute what closed-form algebra gives; the molecular orbital as a probability amplitude would be replaced by a charge surface. Quantum chemistry agreement is not a consequence of "no fitted parameters": group-additivity schemes (Benson) also reproduce thermochemistry to 1 kcal/mol with fitted values [memory, high].

Verdict, claim 4: strong form Never properly tested; weak form Established as reproducible arithmetic (a bookkeeping scheme). Not "dismissed without testing": no independent group is known to have run it blind.

---

Claim 5. Quantum mechanics is wrong or unnecessary

5a. Mathematics

• "The Schroedinger charge density must radiate" (Ch.40). Mills computes a transform of r exp(-r/a_0) (my sympy check: the closed form agrees with Mills's up to a factor 2) that is nonzero at the light cone, then convolves it with angular and time transforms with time dependence delta(omega - omega_n). The premise is that the Schroedinger charge density oscillates at omega_n. For a stationary state |psi|^2 is time independent; there is no oscillation, so no Fourier component at omega = c|k| exists [memory, high]. Also, r exp(-r/a_0) is not |psi|^2 of the 1s state (which is exp(-2r/a_0)). Any smooth non-shell radial profile gives a nonzero transform, so the "result" is generic. Refuted by calculation (accepted electrodynamics applied to a stationary density).
• Aspect/Bell. Eq. 42.23, R(phi)/R_0 = 1/4(...)(...) + 1/4(eps_M1-eps_m1)(eps_M2-eps_m2) F cos 2phi, is the formula Aspect, Grangier and Roger give as the QM prediction (deep reader, medium memory; not verified by me against PRL 47, 460). The mechanism described in words (a shared random linear polarization, independent Malus-law transmission) gives modulation 0.1089 cos 2 phi, exactly half the book's 0.2178, and satisfies the inequalities (delta = -0.096, S = -0.192) [read: deep reader's calculation, simple Malus average; plausibly correct]. The text's own statement of Bell's theorem uses deterministic outcome functions, contradicting "Bell's theorem does not apply to deterministic classical physics" [read]. Jackson's equations (17.70)/(17.71) are quoted with labels reversed [read].
• Muon g-2 (App. III). Eq. 26: g_mu/2 = 1 + (g_e/2)(1 - beta^2), with beta^2 = 2/g_mu taken from the measured g_mu. I recomputed: for any input g/2 between 1.00116 and 1.00117 the output differs from the input by -4e-10 to -1.2e-8 (for the Fermilab value, 1.00116592059 in, 1.00116591329 out, difference -7.3e-9) [calc]. The output simply repeats the input; the "agreement with Carey et al." tests nothing. The assumed E = beta B balance is not the BNL apparatus (electrostatic quadrupoles; field zero on the central orbit) [memory, medium].
• Heisenberg, Durr experiment. The chapter restates the orthogonal-internal-state account (the standard quantum account of Durr et al.) and then calls it classical [read].

5b. Experiments

Loophole-free Bell tests (Hensen 2015 Delft; Giustina 2015 Vienna; Shalm 2015 NIST; later) [memory, high] are not addressed; the newest Bell-test reference in the book is 2007 [read: deep-read text search]. Bell's inequality is violated under strict locality in those experiments, which the book's photon model (polariser cascade, per-pair conservation) does not reach.

5c. Tested? Replicated?

Mills reinterprets other groups' data; no new experiment. The 1982-2018 Bell tests are the tests.

5d. If true

Bell violations would have a local deterministic explanation reproducing the entire quantum correlation curve (it cannot, by Bell's theorem, for local hidden-variable models, unless the model is not local or has loopholes); atomic stability would need a classical proof (App. II postulates). Nothing in the book supplies one.

Verdict, claim 5: radiation claim Refuted by calculation; Aspect/Bell classical derivation Internally inconsistent (+ Tested, not supported by loophole-free tests); muon g-2 Internally inconsistent (circular).

---

Claim 6. Spacetime, gravity, cosmology, particles

6a. Mathematics (cosmology)

Reproduced (deep reader; I recomputed the central ones): T = 2 pi G m_U / c^3 = 3.113e19 s (9.86e11 yr; book 3.10e19, 9.83e11); aleph_0 = 1.866e28 m; aleph(t) = r_g + aleph_0 (1 - cos(2 pi t / T)); H0 = 8 pi^2/T = 78.27 km/s/Mpc, independent of t at fixed m_U [calc]. Defects: - mass divided by mass rate gives seconds; it becomes a length "by multiplying by c" by declaration (Eq. 32.150) [read]; - "km/s^2 = km/s/Mpc" (Eq. 32.155) and 1 Gyr = 3.1358e16 s, prefactor 1229 vs 12,290 (Eq. 32.188) [read]; - rho_min = 8.52e27 m (from T_CMB) is set "equal to" r_g = 2.96e27 m, a factor 2.9 apart [read]; - Mills's H divides by c t (10 Gyr), not by the model's radius aleph; with aleph the same function gives 2.46 km/s/Mpc, and the model's own deceleration parameter is -39.7 [read: deep reader's calculation]; - T_CMB = 3.08 K from the equations (13 percent above 2.7255 K) vs about 2.8 K read off Fig. 32.9; one m_U cannot give both T_CMB and H0 = 78.5 (T needs 2.56e54 kg, giving H0 = 61.0) [read]; - the CMB peak formula 77 sinc[pi(l - 197)/140] uK has positive lobes at l = 197, 541, 824 [read: calc by the deep reader], close to observed 220, 540, 810 [memory, high], but both 197 and 140 come from t = 10 Gyr; with 13.8 Gyr the lobes move to 143, 369, 555 [read]; the relative heights 1 : 0.5 : 0.3 are not those of the squared formula (1 : 0.016 : 0.005) without an unstated additive constant; Eq. 33.46 is an implicit equation for aleph(t) (aleph appears in its own argument) and no substitution into the wave equation is shown [read].

6b. Mathematics (particles)

I verified the pure-alpha ratios: m_mu/m_e = (alpha^-2/(2 pi))^(2/3) = 207.488 (0.35 percent high), corrected by (1 + pi alpha^2)/(1 + alpha/2) to 206.768280 against CODATA 206.7682830(46): miss -1.6e-8, inside the measurement uncertainty [calc; the measured value memory, high]. Look-elsewhere test: I enumerated 57,600 two-factor corrections of the form (1 +/- c alpha^j pi^b)^(+/-1) with c in {1/4, 1/3, 1/2, 2/3, 3/4, 1, 3/2, 2, 3, 4}, j in {1, 2}, b in {0, 1, 2}; 24,148 lie within 1 percent of the needed correction; exactly one family (Mills's two factors, with the variant (1 - pi alpha^2)^-1 that differs at order alpha^4) falls within 2.2e-8, and the chance expectation, from the local density of candidates, is about 0.02 [calc]. So within this hindsight-built family the muon match is not a trivial look-elsewhere artefact; but the family was constructed around Mills's formula, the base exponent (2/3) and the placement of 2 pi in the base were also free, and the tau ratios (miss 5.6e-5, uncertainty 5e-5) are an easy target. Absolute masses use "sec": m_e = (h alpha/(sqrt 2 sec c^2))^(1/2) m_P^(1/2) gives 9.09924e-31 kg with sec = 1 s (0.111 percent low) [calc]; the formula is equivalent to the statement that a time unit t_0 = h alpha m_P/(sqrt 2 c^2 m_e^2) = 0.9978 s equals 1.85e43 Planck times [calc]. In minutes the formula would give m_e smaller by sqrt 60. Fitting "sec = 0.9975 s" (Eq. 36.4) removes any test: with it fitted from m_e, later absolute masses carry one fitted constant, as the deep reader found. A fit within 0.25 percent of the SI second, a unit defined by a caesium transition, is either a numerical coincidence or a statement that the second is fundamental; no mechanism is given. - Nucleon: mu_p = (8 pi/9) mu_N = 2.79253 vs 2.792847 measured (-1.15e-4, about 4e5 times the uncertainty); the book prints the experimental value as 2.79268 [read]. mu_n = -1.91253 vs -1.91304 (-2.7e-4); the 1/5, 1/25 and factor 3 steps have no derivation I could follow [read]. m_N/m_e (Eq. 38.31) misses by 5.2e-6, about 10^4 sigma, after a fitted sec [read]. Z, W, Higgs: 91.16 (uses Mills's own muon mass; 90.97 with the real one), 80.51, 128.75 GeV against 91.188, 80.369, 125.20 [read; memory, medium-high]. Top: 196 - 2(4.58) = 187 GeV (188 from the printed formula) against 172.57 +/- 0.29 GeV today [memory, medium-high]. Deuteron binding 2.181 vs 2.2246 MeV (-2 percent), hidden when expressed as a mass [read]. X17 17.11 MeV (written after the 2016 and 2019 reports, so not a prediction; status unconfirmed) [read; memory, medium-low].

6c. Consistency with experiment and 6d. Mainstream change

Hubble constant 78.5 against 67.4 +/- 0.5 (Planck) and 73.0 +/- 1.0 (SH0ES) [memory, high]; B-mode r = 0.16 versus r < 0.036 (BICEP2's r = 0.20 was withdrawn as a primordial detection) [memory, high]; expansion age 10 Gyr vs 13.8 Gyr. The claimed 1995 prediction of acceleration cannot be checked here (the November 1995 edition is not in the project folder; record needed: that edition's Ch.22, and any earlier edition, to date H0 = 78.5 and the acceleration statement; the top-quark announcement was 2 March 1995 [memory, high]). The model's acceleration is positive from the start of the expansion with no deceleration epoch, whereas the book's own cited supernova paper title mentions a glimpse of "the epoch of deceleration" [read]. If true, general relativity's field equations would not hold as stated, and cosmological parameters would be set by one mass m_U = 2e54 kg.

Verdict, claim 6: cosmological predictions Tested, not supported; absolute lepton/nucleon masses with fitted "sec" Not testable as stated; muon/electron ratio Never properly tested (a notable match, unexplained, and no prospective prediction); nucleon moments, top quark, Higgs, deuteron Tested, not supported.

---

Claim 7. Collective phenomena are classical (digest only)

Digest-level observations [read: digests p.1549-1564]: the "band-pass filter" argument says a cubic array of magnetic dipoles on lattice nuclei passes all frequencies, so no energy is lost; Tc = 30.8 K for conventional metals vs 23.2 K for Nb3Ge, 7/22/180 K for ceramic oxides by dimension; the Hall resistance h/(n e^2) is argued from a flux-linkage and band-pass condition; the Aharonov-Bohm effect is treated as magnetoresistance of a current-carrying ring with a changing flux. Mathematically: a Fourier transform of a lattice of dipoles whose lattice factor is nonzero only on reciprocal-lattice vectors does not by itself exclude scattering, and the conclusion "no frequency is attenuated" needs the condition k = k_z to be the only one that matters, not shown in the digest; the free choice of dimension (1, 2, 3) changes the predicted Tc by a factor of 25; and the Aharonov-Bohm effect measured in vacuum electron interferometry (a shift of the interference fringes by a flux where B = 0 along the beams) [memory, high] cannot be a ring-resistance effect. I did not read these chapters in depth. Provisional verdict: Never properly tested, with the strong reservation that superconductivity, quantum Hall and AB are described by mainstream theory with experimental confirmation, and no outside test of Mills's version is known to me. The order from claims.md (Opus deep read of Ch.25-27) should be run before this verdict is relied on.

---

Investigator's findings I weighed (kept whatever the verdict)

Priority 1, inventor's work: a newer edition (February 2025, 1,888 pp) exists; changes between editions are evidence (the p limit 137 -> 126 is a concrete change); the 165-page reply to Rathke (read in part, Corrections 4); Hagen-Mills 2022/2024 EPR (Corrections 1). Priority 2, support by others: NASA TM-107167, Thermacore/USAF, NIAC, company-hosted validation reports, Booker's arithmetic verifications, IIT Kanpur 2026. Priority 3, related work: Naudts, Khelashvili-Nadareishvili, Giri (graphene), Phillips/Loureiro-Amorim (fast H atoms), Dombey, de Castro, Selke, EarthTech negative. Priority 4: USPTO withdrawal 2000, Federal Circuit 2002, UK IPO 2008, funding history, absence of products.

Where the mathematics bears on these: the related work shows the formal route to sub-ground states in standard theory exists only through singular solutions that disappear with a finite nucleus (Corrections 5); Booker's checks agree with ours that the printed numbers follow from the printed equations, but his phrase "to a very high degree of accuracy" does not address the alpha-dependence (2,000 sigma against modern alpha), the circular muon input, the factor 16 pi^2, or the nonradiation proof; the NASA result is a measured excess with an ordinary explanation not excluded by its authors, which is why I keep it separate from any statement about Mills's theory.

What would settle each point (calculations and records)

1. Claim 1: recompute Eq. 1.228 with an alpha measured independently of g-2 (done above: 2.9e-10 off); evaluate the Fourier transform of Eq. 48 (l = 1, m = 1) at k = omega/c numerically; solve the shell-with-self-force balance (r = 2 a_0).
2. Claim 2: a blind test of p-scaled signatures (H2(1/4) rotational lines at 1950 cm^-1 in Table 16.52, or H^-(1/p) NMR shifts) by a laboratory without a stake; EPR on a control sample made in the same reactor without hydrogen, with spin counting; the assignment table and peak list behind the 0.097 G figure; the IIT Kanpur results.
3. Claim 3: a predicted spectrum shape and flux for the 3.48 keV continuum edge, to be compared with Hitomi/XMM data.
4. Claim 4: a blind prediction, written down before measurement, of a molecule class not in the book.
5. Claim 6: the November 1995 edition (Ch.22); and a derivation of the correction factors 1/(1 + alpha/2) and (1 + pi alpha^2) that does not use the target value.

Open items

Deep reads still needed for claims 4 and 7 (Ch.11, Ch.15 pdf 801-1082, Ch.17-21, Ch.24-31, Ch.33-35, Ch.32 first half). The 2025 edition's changes. The ESI of the 2022 EPR paper. The Lykke/Pekeris comparison for the hydride was not re-run beyond the formula. I did not verify Aspect's printed QM numbers, Pearle (1977), or the 3.48/3.5 keV measurements (Physicist). Web checks done: Lamb shift uncertainty, muonium uncertainty (ppb), CERN muon g-2 series (partial); not done: Fermilab final a_mu, BICEP2 withdrawal, top mass (all [memory, medium-high], none changes a mathematical verdict).