8 Experiment 7 (Exploratory): Mass, Momentum, and the Dispersion Relation

Python: dispersion_poc.py

Experiments 1–6 establish the B-matrix formalism entirely within a massless, scale-free regime: B= 1√ --
  2 for any single hop, independent of chain length L or hop distance d (§4). That scale invariance is itself informative: a rest mass sets a scale, so a theorem that is exact and scale-free by construction has no room for mass to enter as an internal parameter. Momentum, via the hop phase k = π∕(2d), is already present in the formalism; energy and mass are not, because nothing in Experiments 1–6 plays the role of a generator of time evolution or of a second internal degree of freedom at a single site. This experiment asks whether the existing B-matrix machinery survives, unmodified, once such a degree of freedom is added, and whether the relativistic dispersion relation

E2 = (pc)2 + (m0c2 )2

can be recovered as a consequence rather than assumed.

Construction

We introduce a second internal band at a single site (rather than a second spatial site) and adopt the standard two-band lattice (Dirac-on-a-lattice / SSH-type) Hamiltonian

H (k ) = m σz + tsin (k )σx,

where m is an on-site band-splitting term and tsin(k) is a momentum-dependent inter-band coupling. This Hamiltonian is diagonalised numerically (not assumed analytically) at each k on a dense grid, for several values of m∕t. The lower-band eigenvector v(k) is mapped onto the existing Experiment 6 convention ψ(𝜃) = cos𝜃eA + isin𝜃eB by setting sin2𝜃 equal to the numerically computed Born-rule population |vB(k)|2 on the second band. The unmodified Experiment 6 formula B = ρ diag(ρ) is then applied to this ψ(𝜃).

Results

Four claims were checked numerically for m ∈{0,0.1,0.5,1,2} (with t = 1) over a dense grid k (π∕2,π∕2):

  1. Eigenvalues. The numerically diagonalised eigenvalues agree with E(k) = ∘ -------------
  m2 + (tsin k)2 to machine precision (< 1015) at every m and every k tested.
  2. Formalism survives unmodified. The boson norm computed directly from the 2 × 2 matrix B = ρ diag(ρ) agrees with the closed-form Experiment 6 result sin(2𝜃)√ --
  2 to machine precision (< 3 × 1016). No change to the boson-residual formalism was required to accommodate the reinterpretation of 𝜃 as a momentum/mass mixing angle rather than a spatial-hop angle.
  3. Mass–momentum identity. Working the eigenvector algebra of H(k) through by hand gives

                             2        |-------------------|
sin2𝜃 (1− sin2𝜃) = (tsin-k)-  = ⇒   |∥B (k)∥ = |√tsin-k|. |
                  4E (k)2         ------------2-E(k)--|

    This closed form agrees with the numerically computed boson norm to machine precision (< 3 × 1013) at every m and k tested. At m = 0 it reduces exactly to the universal B= 1√ --
  2 of Experiments 3–6, recovering the massless case as a special case rather than contradicting it. As m grows relative to tsink, B(k)shrinks monotonically toward zero: a more massive state produces a smaller compression residual at fixed momentum.

  4. Continuum limit. Identifying tsink pc and m m0c2 in the small-k regime, E(k)2 agrees with (pc)2 + (m0c2)2 to a relative error that shrinks as k 0 (e.g. 7 × 104 over |k|< 0.05 at m = 0, tightening to 5 × 107 at m = 2), consistent with the expected O(k2) lattice correction to the continuum dispersion relation.

Status of this result

This experiment shows that the B-matrix formalism of Papers VI–VII is compatible with a mass term and reproduces the relativistic dispersion relation as a continuum limit, once a second internal band and a Hamiltonian are introduced. It does not show that this construction is forced by the framework’s own first principles. Specifically:

Consistent with the three-tier accounting adopted in Paper XI, this result is reported as a strongly suggested correspondence, not an analytically exact theorem in the sense of Experiments 1–6: it is numerically exact given its assumptions, but its assumptions are not yet derived.