Calculation: The GWT Lagrangian
From lattice springs to the Standard Model action — every term derived, every coefficient computed.
§1 — The Microscopic Lattice Lagrangian
Every node of the elastic lattice has two kinds of degrees of freedom: displacement (where the node sits) and orientation (its internal yin–yang angle). The Lagrangian for a single lattice element is:
Lattice Element Lagrangian
ℒ = T − V
Kinetic Energy (T)
1
Translational kinetic energy
Ttrans = ½ η (∂u/∂t)²
η = 1.385 × 10−8 kg (inertial response per node)
u = displacement vector of the node from equilibrium
Ttrans = ½ η (∂u/∂t)²
η = 1.385 × 10−8 kg (inertial response per node)
u = displacement vector of the node from equilibrium
2
Internal rotation kinetic energy
Trot = ½ I (∂θ/∂t)²
I = moment of inertia of the yin–yang oscillator at each node
θ = orientation angle (the internal “spin” degree of freedom)
Trot = ½ I (∂θ/∂t)²
I = moment of inertia of the yin–yang oscillator at each node
θ = orientation angle (the internal “spin” degree of freedom)
Total kinetic energy:
T = ½ η (∂u/∂t)² + ½ I (∂θ/∂t)²
T = ½ η (∂u/∂t)² + ½ I (∂θ/∂t)²
Potential Energy (V)
3
Elastic stretching (bond distortion)
Velastic = ½ k (∂u/∂x)² a
k = 4.77 × 1078 N/m (lattice stiffness per bond)
a = 1.616 × 10−35 m (lattice spacing = Planck length)
The factor of a converts from strain to energy per node.
Velastic = ½ k (∂u/∂x)² a
k = 4.77 × 1078 N/m (lattice stiffness per bond)
a = 1.616 × 10−35 m (lattice spacing = Planck length)
The factor of a converts from strain to energy per node.
4
Orientation coupling (neighbour interaction)
Vorient = ½ κ (θi − θj)²
κ = orientation stiffness (torsional spring between adjacent nodes)
(θi − θj) = angular mismatch between neighbours
This term resists orientation gradients — it will become the gauge field.
Vorient = ½ κ (θi − θj)²
κ = orientation stiffness (torsional spring between adjacent nodes)
(θi − θj) = angular mismatch between neighbours
This term resists orientation gradients — it will become the gauge field.
5
Double-well self-energy (yin–yang bistability)
VDW(θ) = −½ k2 θ² + ¼ k4 θ4
Each node prefers one of two orientations: yin (+) or yang (−).
The double well has minima at θ = ±√(k2/k4).
This term will become the Higgs potential.
VDW(θ) = −½ k2 θ² + ¼ k4 θ4
Each node prefers one of two orientations: yin (+) or yang (−).
The double well has minima at θ = ±√(k2/k4).
This term will become the Higgs potential.
Complete microscopic Lagrangian per element:
ℒ = ½η(∂u/∂t)² + ½I(∂θ/∂t)² − ½k(∂u/∂x)²a − ½κ(θi − θj)² − VDW(θ)
ℒ = ½η(∂u/∂t)² + ½I(∂θ/∂t)² − ½k(∂u/∂x)²a − ½κ(θi − θj)² − VDW(θ)
In Full 3D
6
Promote to 3D vector field
u → u(x,y,z,t) = displacement vector (3 components)
θ → Θ(x,y,z,t) = orientation field (SU(2) doublet [4 real DOF] + U(1) phase [1 DOF] = 5 total internal DOF)
The full 3D lattice Lagrangian density sums over all three spatial directions:
ℒ3D = ½η |∂u/∂t|² + ½I |∂Θ/∂t|² − ½k a ∑α |∂u/∂xα|² − ½κ ∑α |∂Θ/∂xα|² − VDW(Θ)
u → u(x,y,z,t) = displacement vector (3 components)
θ → Θ(x,y,z,t) = orientation field (SU(2) doublet [4 real DOF] + U(1) phase [1 DOF] = 5 total internal DOF)
The full 3D lattice Lagrangian density sums over all three spatial directions:
ℒ3D = ½η |∂u/∂t|² + ½I |∂Θ/∂t|² − ½k a ∑α |∂u/∂xα|² − ½κ ∑α |∂Θ/∂xα|² − VDW(Θ)
Physical Meaning
This is the Lagrangian of a 3D elastic medium with internal structure. Every term is a spring — literal Hooke’s law. The displacement sector will produce gravity. The orientation sector will produce gauge fields, fermions, and the Higgs mechanism. There are no quantum postulates, no path integrals, no renormalization — just springs.
§2 — The Continuum Limit
At scales much larger than the lattice spacing a, we replace discrete differences with derivatives. This is the standard long-wavelength approximation used in condensed matter physics.
1
Discrete → Continuous
(uj − ui) / a → ∂u/∂x
(θj − θi) / a → ∂θ/∂x
∑nodes a³ → ∫ d³x
(uj − ui) / a → ∂u/∂x
(θj − θi) / a → ∂θ/∂x
∑nodes a³ → ∫ d³x
2
Two sectors emerge cleanly
Displacement sector (longitudinal + transverse compression): → Gravity
Orientation sector (internal angle gradients + double-well): → Gauge fields + Higgs
These sectors couple only through the metric (spacetime curvature = lattice compression).
Displacement sector (longitudinal + transverse compression): → Gravity
Orientation sector (internal angle gradients + double-well): → Gauge fields + Higgs
These sectors couple only through the metric (spacetime curvature = lattice compression).
3
Natural UV cutoff
The lattice spacing a = 1.616 × 10−35 m provides a hard ultraviolet cutoff.
Maximum wave number: qmax = π/a = 1.943 × 1035 m−1
Maximum energy: Emax = ℏc × qmax = 1.221 × 1019 GeV = EPlanck
No wavelength shorter than 2a can propagate on the lattice.
This is the Brillouin zone boundary — the same physics that prevents ultraviolet catastrophes in crystals.
Result: No infinities. Renormalization in QFT is an approximation to the finite lattice sum.
The lattice spacing a = 1.616 × 10−35 m provides a hard ultraviolet cutoff.
Maximum wave number: qmax = π/a = 1.943 × 1035 m−1
Maximum energy: Emax = ℏc × qmax = 1.221 × 1019 GeV = EPlanck
No wavelength shorter than 2a can propagate on the lattice.
This is the Brillouin zone boundary — the same physics that prevents ultraviolet catastrophes in crystals.
Result: No infinities. Renormalization in QFT is an approximation to the finite lattice sum.
The continuum Lagrangian density (both sectors):
ℒ = ℒgravity(gμν) + ℒgauge(Aμ) + ℒmatter(ψ) + ℒHiggs(φ)
ℒ = ℒgravity(gμν) + ℒgauge(Aμ) + ℒmatter(ψ) + ℒHiggs(φ)
Why No Infinities
In standard QFT, loop integrals diverge because momentum is integrated to infinity. On the lattice, the Brillouin zone boundary qmax = π/a truncates every integral. The “renormalization group” is a correct accounting tool, but the underlying physics is simply: the medium is discrete, so infinite momenta do not exist.
§3 — The Gravity Sector (Einstein–Hilbert)
Longitudinal compression of the lattice — nodes pushed closer together or pulled apart — curves spacetime. The elastic energy stored in this compression is proportional to the Ricci scalar R, which measures intrinsic curvature.
Gravity Lagrangian
ℒgravity = (ka/32) R = c4/(16πG) × R
Key Identity: ka/32 = c4/(16πG)
1
Compute c4/(16πG) from known constants
c = 2.998 × 108 m/s
G = 6.674 × 10−11 m³/(kg·s²)
c4 = (2.998 × 108)4 = 8.077 × 1033 m4/s4
16πG = 16 × 3.14159 × 6.674 × 10−11 = 3.353 × 10−9 m³/(kg·s²)
c4/(16πG) = 8.077 × 1033 / 3.353 × 10−9
c4/(16πG) = 2.409 × 1042 kg/(m·s²)
c = 2.998 × 108 m/s
G = 6.674 × 10−11 m³/(kg·s²)
c4 = (2.998 × 108)4 = 8.077 × 1033 m4/s4
16πG = 16 × 3.14159 × 6.674 × 10−11 = 3.353 × 10−9 m³/(kg·s²)
c4/(16πG) = 8.077 × 1033 / 3.353 × 10−9
c4/(16πG) = 2.409 × 1042 kg/(m·s²)
2
Compute ka/32 from lattice constants
k = 4.77 × 1078 N/m
a = 1.616 × 10−35 m
ka = 4.77 × 1078 × 1.616 × 10−35 = 7.708 × 1043 N
ka/32 = 7.708 × 1043 / 32
ka/32 = 2.409 × 1042 N = 2.409 × 1042 kg/(m·s²)
k = 4.77 × 1078 N/m
a = 1.616 × 10−35 m
ka = 4.77 × 1078 × 1.616 × 10−35 = 7.708 × 1043 N
ka/32 = 7.708 × 1043 / 32
ka/32 = 2.409 × 1042 N = 2.409 × 1042 kg/(m·s²)
GWT (lattice)
ka/32 = 2.409 × 1042
Standard GR
c4/(16πG) = 2.409 × 1042
EXACT
3
Algebraic proof (using G = 2c4/(πka))
From the master equation G = 2c4/(πka), solve for ka:
ka = 2c4/(πG)
Therefore: ka/32 = 2c4/(32πG) = c4/(16πG) ✓
The identity is exact, not approximate. The Einstein–Hilbert coefficient is literally half a lattice bond stiffness divided by 16.
From the master equation G = 2c4/(πka), solve for ka:
ka = 2c4/(πG)
Therefore: ka/32 = 2c4/(32πG) = c4/(16πG) ✓
The identity is exact, not approximate. The Einstein–Hilbert coefficient is literally half a lattice bond stiffness divided by 16.
Result: The Einstein–Hilbert action for gravity
Sgravity = ∫d4x √(−g) × (ka/32) R = ∫d4x √(−g) × c4/(16πG) × R
This IS general relativity. The Ricci scalar R measures lattice curvature (compression gradients). The coefficient ka/32 sets the stiffness of spacetime. Einstein’s gravity is Hooke’s law for the medium.
Sgravity = ∫d4x √(−g) × (ka/32) R = ∫d4x √(−g) × c4/(16πG) × R
This IS general relativity. The Ricci scalar R measures lattice curvature (compression gradients). The coefficient ka/32 sets the stiffness of spacetime. Einstein’s gravity is Hooke’s law for the medium.
§4 — The Electromagnetic Sector (Maxwell)
Oscillations of the yin–yang orientation angle θ propagate as transverse waves through the lattice. These are photons. The orientation gradient defines the electromagnetic field tensor Fμν.
EM Lagrangian
ℒEM = −(1/4μ0) Fμν Fμν
1
Orientation gradient → field tensor
Define a 4-potential from the orientation field: Aμ ∝ ∂θ/∂xμ
The field tensor is the antisymmetric curl:
Fμν = ∂μAν − ∂νAμ
This is automatically gauge-invariant: θ → θ + const leaves Fμν unchanged.
Gauge symmetry is not postulated — it is a geometric identity of the orientation field.
Define a 4-potential from the orientation field: Aμ ∝ ∂θ/∂xμ
The field tensor is the antisymmetric curl:
Fμν = ∂μAν − ∂νAμ
This is automatically gauge-invariant: θ → θ + const leaves Fμν unchanged.
Gauge symmetry is not postulated — it is a geometric identity of the orientation field.
2
The photon as a lattice wave
A propagating orientation oscillation has:
• Spin-1: because θ is a vector-like rotation (orientation has direction)
• Massless: because VDW does not affect the U(1) phase degree of freedom
• Speed c: the orientation waves travel at the same speed as displacement waves
• Transverse: orientation oscillations are perpendicular to propagation direction
These are precisely the properties of the photon. Light is a transverse orientation wave.
A propagating orientation oscillation has:
• Spin-1: because θ is a vector-like rotation (orientation has direction)
• Massless: because VDW does not affect the U(1) phase degree of freedom
• Speed c: the orientation waves travel at the same speed as displacement waves
• Transverse: orientation oscillations are perpendicular to propagation direction
These are precisely the properties of the photon. Light is a transverse orientation wave.
3
The coupling constant μ0
μ0 = 4π × 10−7 H/m (magnetic permeability of free space)
ε0 = 1/(μ0c²) = 8.854 × 10−12 F/m (electric permittivity)
In the lattice picture, ε0 and μ0 are the orientation-sector analogues of η and k:
• ε0 measures how easily the medium’s orientation polarises (analogue of inertial response η)
• μ0 measures how strongly orientation gradients resist distortion (analogue of stiffness k)
They satisfy c = 1/√(μ0ε0), exactly paralleling c = a√(k/η).
μ0 = 4π × 10−7 H/m (magnetic permeability of free space)
ε0 = 1/(μ0c²) = 8.854 × 10−12 F/m (electric permittivity)
In the lattice picture, ε0 and μ0 are the orientation-sector analogues of η and k:
• ε0 measures how easily the medium’s orientation polarises (analogue of inertial response η)
• μ0 measures how strongly orientation gradients resist distortion (analogue of stiffness k)
They satisfy c = 1/√(μ0ε0), exactly paralleling c = a√(k/η).
Result: Maxwell’s equations from lattice orientation waves
The orientation coupling ½κ(∇θ)² in the continuum limit becomes −(1/4μ0)FμνFμν.
Maxwell’s four equations follow automatically from the Euler–Lagrange equations of this term.
Electromagnetism is the dynamics of orientation waves in the elastic medium.
The orientation coupling ½κ(∇θ)² in the continuum limit becomes −(1/4μ0)FμνFμν.
Maxwell’s four equations follow automatically from the Euler–Lagrange equations of this term.
Electromagnetism is the dynamics of orientation waves in the elastic medium.
§5 — The Matter Sector (Dirac)
Standing waves in the lattice are particles. The binary yin–yang state of each node provides spinor structure. The lattice wave equation, combined with the yin–yang coupling, produces the Dirac equation.
Matter Lagrangian
ℒmatter = ψ̅(iγμDμ − m)ψ
Deriving the Dirac Equation from the Lattice
1
Yin–yang binary state → spinor
Each node has orientation: yin (+), yang (−), or neutral (0).
The binary subset {+, −} forms a 2-component object: the Pauli spinor.
Algebraically, the yin–yang operators generate the Clifford algebra Cl(3,0):
{σi, σj} = 2δij
where σi are the three Pauli matrices corresponding to the three spatial directions.
Each node has orientation: yin (+), yang (−), or neutral (0).
The binary subset {+, −} forms a 2-component object: the Pauli spinor.
Algebraically, the yin–yang operators generate the Clifford algebra Cl(3,0):
{σi, σj} = 2δij
where σi are the three Pauli matrices corresponding to the three spatial directions.
2
Add time → Dirac algebra
Time is the zeroth lattice direction (causal propagation). Promoting the algebra:
Cl(3,0) → Cl(3,1)
This gives the four Dirac gamma matrices satisfying:
{γμ, γν} = 2gμν
The 2-component Pauli spinor becomes a 4-component Dirac spinor ψ.
This is not a postulate — it is the unique algebraic consequence of binary states on a 3+1D lattice.
Time is the zeroth lattice direction (causal propagation). Promoting the algebra:
Cl(3,0) → Cl(3,1)
This gives the four Dirac gamma matrices satisfying:
{γμ, γν} = 2gμν
The 2-component Pauli spinor becomes a 4-component Dirac spinor ψ.
This is not a postulate — it is the unique algebraic consequence of binary states on a 3+1D lattice.
3
Lattice wave equation → Dirac equation
The wave equation on the lattice is:
(η/a³) ∂²ψ/∂t² = (k/a) ∇²ψ
For a standing wave with rest frequency ω0 = mc²/ℏ, factorising the second-order equation into two first-order equations (the standard Dirac trick) gives:
(iγμ∂μ − mc/ℏ)Ψ = 0
This is the Dirac equation. The “trick” of factorisation is forced by the yin–yang binary structure — it is the only consistent linearisation.
The wave equation on the lattice is:
(η/a³) ∂²ψ/∂t² = (k/a) ∇²ψ
For a standing wave with rest frequency ω0 = mc²/ℏ, factorising the second-order equation into two first-order equations (the standard Dirac trick) gives:
(iγμ∂μ − mc/ℏ)Ψ = 0
This is the Dirac equation. The “trick” of factorisation is forced by the yin–yang binary structure — it is the only consistent linearisation.
4
The Born rule as wave energy density
P(x) = |ψ(x)|² = wave energy density at position x
In standard quantum mechanics, the Born rule P = |ψ|² is a postulate.
In the lattice picture, it is a theorem: the energy density of a classical wave is proportional to the square of its amplitude. Where the wave is large, more energy is concentrated, and that is where “detection” (energy exchange) is most likely.
Quantum probability is classical wave energy. No measurement postulate needed.
P(x) = |ψ(x)|² = wave energy density at position x
In standard quantum mechanics, the Born rule P = |ψ|² is a postulate.
In the lattice picture, it is a theorem: the energy density of a classical wave is proportional to the square of its amplitude. Where the wave is large, more energy is concentrated, and that is where “detection” (energy exchange) is most likely.
Quantum probability is classical wave energy. No measurement postulate needed.
Result: The Dirac Lagrangian from lattice wave mechanics
ℒmatter = ψ̅(iγμDμ − m)ψ
Dμ = ∂μ + igAμ is the gauge-covariant derivative (coupling to orientation field).
The mass m comes from the standing-wave resonance condition.
The spinor structure comes from the yin–yang binary state.
The Born rule comes from wave energy density.
Three “quantum postulates” — spinors, the Dirac equation, and Born’s rule — are all mechanical consequences.
ℒmatter = ψ̅(iγμDμ − m)ψ
Dμ = ∂μ + igAμ is the gauge-covariant derivative (coupling to orientation field).
The mass m comes from the standing-wave resonance condition.
The spinor structure comes from the yin–yang binary state.
The Born rule comes from wave energy density.
Three “quantum postulates” — spinors, the Dirac equation, and Born’s rule — are all mechanical consequences.
§6 — The Higgs / Double-Well Sector
The yin–yang orientation at each node lives in a double-well potential. This is the Higgs mechanism, derived from the lattice rather than postulated.
Double-Well Potential
V(θ) = −½ k2 θ² + ¼ k4 θ4
1
This IS the Higgs potential
The Standard Model Higgs potential is:
V(φ) = −μ² |φ|² + λ |φ|4
Compare with the lattice double-well:
V(θ) = −½ k2 θ² + ¼ k4 θ4
Identification: μ² ↔ k2/2, λ ↔ k4/4, φ ↔ θ
Same functional form. Same physics. Different language.
The Standard Model Higgs potential is:
V(φ) = −μ² |φ|² + λ |φ|4
Compare with the lattice double-well:
V(θ) = −½ k2 θ² + ¼ k4 θ4
Identification: μ² ↔ k2/2, λ ↔ k4/4, φ ↔ θ
Same functional form. Same physics. Different language.
2
Symmetry breaking = choosing a well
The double well has minima at θ0 = ±√(k2/k4).
The lattice “chooses” one well: this is electroweak symmetry breaking.
The vacuum expectation value (VEV) is:
v = θ0 ⇔ v = √(μ²/λ) in SM notation
Before breaking: SU(2) × U(1)Y symmetry (both wells equivalent)
After breaking: U(1)EM symmetry (one well selected, 3 Goldstone bosons absorbed)
The double well has minima at θ0 = ±√(k2/k4).
The lattice “chooses” one well: this is electroweak symmetry breaking.
The vacuum expectation value (VEV) is:
v = θ0 ⇔ v = √(μ²/λ) in SM notation
Before breaking: SU(2) × U(1)Y symmetry (both wells equivalent)
After breaking: U(1)EM symmetry (one well selected, 3 Goldstone bosons absorbed)
3
The Higgs field = collective yin–yang amplitude
The Higgs boson is the excitation about the chosen minimum:
θ(x) = v + h(x)
where h(x) is the physical Higgs field (discovered 2012, mH = 125.1 GeV).
In the lattice picture, the Higgs is a collective oscillation of the yin–yang angle about its chosen orientation.
The Higgs boson is the excitation about the chosen minimum:
θ(x) = v + h(x)
where h(x) is the physical Higgs field (discovered 2012, mH = 125.1 GeV).
In the lattice picture, the Higgs is a collective oscillation of the yin–yang angle about its chosen orientation.
The Full Lattice Potential
4
Cosine potential from Brillouin cutoff
At the microscopic scale, the full lattice potential is not a polynomial but a cosine, uniquely determined by the Brillouin zone cutoff:
V(x) = (ka²/π²)[1 − cos(πx/a)]
The Brillouin cutoff kills all harmonics n ≥ 2, leaving only the fundamental — a unique form.
The double-well −½k2θ² + ¼k4θ4 is its Taylor expansion to quartic order.
At the microscopic scale, the full lattice potential is not a polynomial but a cosine, uniquely determined by the Brillouin zone cutoff:
V(x) = (ka²/π²)[1 − cos(πx/a)]
The Brillouin cutoff kills all harmonics n ≥ 2, leaving only the fundamental — a unique form.
The double-well −½k2θ² + ¼k4θ4 is its Taylor expansion to quartic order.
5
Barrier height: Vmax
Vmax = 2ka²/π² (at x = a, the peak of the cosine)
Numerically:
Vmax = 2 × 4.77 × 1078 × (1.616 × 10−35)² / π²
Vmax = 2 × 4.77 × 1078 × 2.611 × 10−70 / 9.8696
Vmax = 2 × 1.246 × 109 / 9.8696
Vmax = 2.524 × 108 J ≈ 0.13 EPlanck per node
Check: EPlanck = √(ℏc5/G) = 1.956 × 109 J
Vmax/EPlanck = 2.524 × 108 / 1.956 × 109 = 0.129 = 4/π³ ✓
Vmax = 2ka²/π² (at x = a, the peak of the cosine)
Numerically:
Vmax = 2 × 4.77 × 1078 × (1.616 × 10−35)² / π²
Vmax = 2 × 4.77 × 1078 × 2.611 × 10−70 / 9.8696
Vmax = 2 × 1.246 × 109 / 9.8696
Vmax = 2.524 × 108 J ≈ 0.13 EPlanck per node
Check: EPlanck = √(ℏc5/G) = 1.956 × 109 J
Vmax/EPlanck = 2.524 × 108 / 1.956 × 109 = 0.129 = 4/π³ ✓
Result: The Higgs mechanism is the yin–yang double well
V(θ) = −½k2θ² + ¼k4θ4 = SM Higgs potential
V(x) = (ka²/π²)[1 − cos(πx/a)] = exact lattice form
Vmax = 4/π³ EPlanck ≈ 0.13 EP per node
The Higgs field is not a new fundamental ingredient — it is the collective orientation of an elastic medium with binary internal structure.
V(θ) = −½k2θ² + ¼k4θ4 = SM Higgs potential
V(x) = (ka²/π²)[1 − cos(πx/a)] = exact lattice form
Vmax = 4/π³ EPlanck ≈ 0.13 EP per node
The Higgs field is not a new fundamental ingredient — it is the collective orientation of an elastic medium with binary internal structure.
§7 — The Gauge Group: SU(3) × SU(2) × U(1)
Each lattice node has 5 internal degrees of freedom in addition to its 3 spatial displacement components. These 5 + 3 = 8 DOF determine the gauge group of the Standard Model.
1
Count the degrees of freedom
Phase (1 DOF): oscillation phase of the node → U(1) → electromagnetism → 1 gauge boson (photon)
Yin–yang (2 DOF): binary orientation (up/down doublet) → SU(2) → weak force → 3 gauge bosons (W+, W−, Z)
Spatial displacement (3 DOF): x, y, z displacement directions → SU(3) → strong force → 8 gluons
Phase (1 DOF): oscillation phase of the node → U(1) → electromagnetism → 1 gauge boson (photon)
Yin–yang (2 DOF): binary orientation (up/down doublet) → SU(2) → weak force → 3 gauge bosons (W+, W−, Z)
Spatial displacement (3 DOF): x, y, z displacement directions → SU(3) → strong force → 8 gluons
2
Verify the group dimensions
U(1): 1 generator → 1 boson (photon γ) ✓
SU(2): 2² − 1 = 3 generators → 3 bosons (W+, W−, Z0) ✓
SU(3): 3² − 1 = 8 generators → 8 bosons (gluons g1…g8) ✓
Total gauge bosons: 1 + 3 + 8 = 12 ✓
U(1): 1 generator → 1 boson (photon γ) ✓
SU(2): 2² − 1 = 3 generators → 3 bosons (W+, W−, Z0) ✓
SU(3): 3² − 1 = 8 generators → 8 bosons (gluons g1…g8) ✓
Total gauge bosons: 1 + 3 + 8 = 12 ✓
The Weinberg Angle at GUT Scale
3
sin²θW(GUT) from DOF counting
At the unification scale, all couplings are equal. The hypercharge fraction is the ratio of spatial DOF to total DOF:
sin²θW(GUT) = 3 / 8 = 0.375
This is the same prediction as SU(5) grand unification — derived here from lattice geometry, not from group embedding.
At the unification scale, all couplings are equal. The hypercharge fraction is the ratio of spatial DOF to total DOF:
sin²θW(GUT) = 3 / 8 = 0.375
This is the same prediction as SU(5) grand unification — derived here from lattice geometry, not from group embedding.
Three Generations from Three Dimensions
4
Why exactly 3 generations of fermions?
The lattice has 3 spatial dimensions. Each dimension provides one independent standing-wave axis for the orientation field. The three generations are:
• Generation 1 (e, νe, u, d): standing wave along x
• Generation 2 (μ, νμ, c, s): standing wave along y
• Generation 3 (τ, ντ, t, b): standing wave along z
Higher generations would require a 4th spatial dimension, which the lattice does not have.
Ngen = 3 is derived, not postulated.
The lattice has 3 spatial dimensions. Each dimension provides one independent standing-wave axis for the orientation field. The three generations are:
• Generation 1 (e, νe, u, d): standing wave along x
• Generation 2 (μ, νμ, c, s): standing wave along y
• Generation 3 (τ, ντ, t, b): standing wave along z
Higher generations would require a 4th spatial dimension, which the lattice does not have.
Ngen = 3 is derived, not postulated.
Running to MZ (three derived terms, 0.019%)
5
sin²θW at MZ (MS-bar)
Three terms, all from d=3 geometry:
sin²θW = 15/64 − d·α/2 + α·ln(6π&sup5;)/(2d+1)
= 0.23438 − 0.01095 + 0.00780
= 0.23123
Observed (MS-bar at MZ): 0.23122. Error: 0.019%
• 15/64: tree level from d-cube vertex counting (7/8 vertices connect weakly)
• −d·α/2: one-loop VP (d=3 axes × α/2 per axis)
• +α·ln(F)/(2d+1): threshold running (7 exchange paths, same as g−2 and n−p mass difference)
Three terms, all from d=3 geometry:
sin²θW = 15/64 − d·α/2 + α·ln(6π&sup5;)/(2d+1)
= 0.23438 − 0.01095 + 0.00780
= 0.23123
Observed (MS-bar at MZ): 0.23122. Error: 0.019%
• 15/64: tree level from d-cube vertex counting (7/8 vertices connect weakly)
• −d·α/2: one-loop VP (d=3 axes × α/2 per axis)
• +α·ln(F)/(2d+1): threshold running (7 exchange paths, same as g−2 and n−p mass difference)
Result: The Standard Model gauge group from lattice DOF
SU(3) × SU(2) × U(1) = {3 spatial} × {2 transverse} × {1 longitudinal}
sin²θW(MZ) = 15/64 − 3α/2 + α·ln(6π&sup5;)/7 = 0.23123 (0.019%)
Ngen = 3 (spatial dimensions = torus axes)
Nc = 3 (colour charges = displacement directions)
Every piece of the Standard Model gauge structure is a counting exercise on the lattice.
SU(3) × SU(2) × U(1) = {3 spatial} × {2 transverse} × {1 longitudinal}
sin²θW(MZ) = 15/64 − 3α/2 + α·ln(6π&sup5;)/7 = 0.23123 (0.019%)
Ngen = 3 (spatial dimensions = torus axes)
Nc = 3 (colour charges = displacement directions)
Every piece of the Standard Model gauge structure is a counting exercise on the lattice.
§8 — The Unified Action (Complete)
Collecting all sectors, the complete GWT action is:
The GWT Action
S = ∫d4x √(−g) [ (ka/32)R − (1/4μ0)FμνFμν + ψ̅(iγμDμ − m)ψ − V(θ) ]
Every coefficient maps to the three lattice constants {k, a, η}. Let us verify each one.
Coefficient Map
1
Einstein–Hilbert coefficient
ka/32 = c4/(16πG)
Numerical check (performed in §3):
ka/32 = 4.77 × 1078 × 1.616 × 10−35 / 32 = 2.409 × 1042
c4/(16πG) = 8.077 × 1033 / 3.353 × 10−9 = 2.409 × 1042 ✓
ka/32 = c4/(16πG)
Numerical check (performed in §3):
ka/32 = 4.77 × 1078 × 1.616 × 10−35 / 32 = 2.409 × 1042
c4/(16πG) = 8.077 × 1033 / 3.353 × 10−9 = 2.409 × 1042 ✓
2
Maxwell coefficient
1/(4μ0) = c²/(4/ε0) (using μ0 = 1/(ε0c²))
μ0 comes from the lattice orientation coupling κ:
μ0 = 4π × 10−7 H/m ↔ κ/(a³ c²)
ε0 = 1/(μ0c²) = 8.854 × 10−12 F/m
1/(4μ0) = 1/(4 × 4π × 10−7) = 7.958 × 104 H−1m−1
1/(4μ0) = c²/(4/ε0) (using μ0 = 1/(ε0c²))
μ0 comes from the lattice orientation coupling κ:
μ0 = 4π × 10−7 H/m ↔ κ/(a³ c²)
ε0 = 1/(μ0c²) = 8.854 × 10−12 F/m
1/(4μ0) = 1/(4 × 4π × 10−7) = 7.958 × 104 H−1m−1
3
Fermion masses: all derived from {α, mPlanck, 6π5}
The proton mass sets the scale:
mp = 4ΛQCD (virial theorem in d=3, §5 mass gap derivation)
The electron mass derives from the proton via the 6π5 formula:
Bare: mp/me = 6π5 = 6 × 306.020 = 1836.12
With VP correction: 6π5(1 + α²/2√2) = 1836.153
Measured: 1836.153 (error: < 0.001 ppm)
[VP correction from ΣQ² = 1 quark charge identity, d=3 only]
The muon mass (Koide-extended):
mμ/me = 206.77
Measured: 206.768 (error: 0.005%)
Neutrino masses:
mν = me³/(Nc × mp²) (see-saw from lattice, zero new inputs)
The proton mass sets the scale:
mp = 4ΛQCD (virial theorem in d=3, §5 mass gap derivation)
The electron mass derives from the proton via the 6π5 formula:
Bare: mp/me = 6π5 = 6 × 306.020 = 1836.12
With VP correction: 6π5(1 + α²/2√2) = 1836.153
Measured: 1836.153 (error: < 0.001 ppm)
[VP correction from ΣQ² = 1 quark charge identity, d=3 only]
The muon mass (Koide-extended):
mμ/me = 206.77
Measured: 206.768 (error: 0.005%)
Neutrino masses:
mν = me³/(Nc × mp²) (see-saw from lattice, zero new inputs)
4
Higgs VEV from lattice constants
v = (5/2) mPlanck α8
Numerical computation:
mPlanck = √(ℏc/G) = 2.176 × 10−8 kg = 1.221 × 1019 GeV/c²
α = 1/137.042
α8 = (1/137.042)8 = (7.297 × 10−3)8
α2 = 5.325 × 10−5
α4 = 2.836 × 10−9
α8 = 8.041 × 10−18
v = (5/2) × 1.221 × 1019 × 8.041 × 10−18
v = 2.5 × 9.814 × 101
v = 245.4 GeV
v = (5/2) mPlanck α8
Numerical computation:
mPlanck = √(ℏc/G) = 2.176 × 10−8 kg = 1.221 × 1019 GeV/c²
α = 1/137.042
α8 = (1/137.042)8 = (7.297 × 10−3)8
α2 = 5.325 × 10−5
α4 = 2.836 × 10−9
α8 = 8.041 × 10−18
v = (5/2) × 1.221 × 1019 × 8.041 × 10−18
v = 2.5 × 9.814 × 101
v = 245.4 GeV
GWT Higgs VEV
245.4 GeV
Measured VEV
246.2 GeV
0.3%
5
Higgs quartic coupling λ
Primary (breather spectrum + scalar VP):
MH = m(8, 24) × πα/(d−1) = 125.28 GeV
λ = (MH/v)²/2 = (125.28/246.22)²/2
λ = 0.1295 (0.4% error)
Cross-check (tree-level): λ = 1/2d = 1/8 = 0.125 (3.1% error, before VP dressing)
Primary (breather spectrum + scalar VP):
MH = m(8, 24) × πα/(d−1) = 125.28 GeV
λ = (MH/v)²/2 = (125.28/246.22)²/2
λ = 0.1295 (0.4% error)
Cross-check (tree-level): λ = 1/2d = 1/8 = 0.125 (3.1% error, before VP dressing)
GWT λ
0.1295
Measured λ
0.129 ± 0.006
0.4%
Side-by-Side: GWT vs Standard Model
The GWT action IS the Standard Model action
Standard Model:
SSM = ∫d4x √(−g) [ c4/(16πG) R − ¼FμνFμν + ψ̅(iγμDμ − m)ψ + |Dμφ|² − V(φ) ]
GWT:
SGWT = ∫d4x √(−g) [ (ka/32)R − (1/4μ0)FμνFμν + ψ̅(iγμDμ − m)ψ − V(θ) ]
Same structure. Same equations. Same predictions.
The only difference: GWT has zero free parameters. Every coefficient is computed from {k, a, η}.
Standard Model:
SSM = ∫d4x √(−g) [ c4/(16πG) R − ¼FμνFμν + ψ̅(iγμDμ − m)ψ + |Dμφ|² − V(φ) ]
GWT:
SGWT = ∫d4x √(−g) [ (ka/32)R − (1/4μ0)FμνFμν + ψ̅(iγμDμ − m)ψ − V(θ) ]
Same structure. Same equations. Same predictions.
The only difference: GWT has zero free parameters. Every coefficient is computed from {k, a, η}.
§9 — Lattice Corrections (Testable Predictions)
The continuum limit is an approximation. The exact lattice dispersion relation introduces corrections that differ from standard QFT and are in principle testable.
The Lattice Dispersion Relation
Exact Lattice Dispersion
ω² = (4k/η) sin²(qa/2)
1
Taylor expansion for qa « 1
sin²(qa/2) = (qa/2)² − (qa/2)4/3 + …
sin²(qa/2) = q²a²/4 [1 − (qa)²/12 + (qa)4/360 − …]
ω² = (4k/η)(q²a²/4)[1 − (qa)²/12 + …]
ω² = (ka²/η)q²[1 − (qa)²/12 + …]
ω² = c²q²[1 − (qa)²/12 + …]
where we used c² = ka²/η = a²(k/η).
sin²(qa/2) = (qa/2)² − (qa/2)4/3 + …
sin²(qa/2) = q²a²/4 [1 − (qa)²/12 + (qa)4/360 − …]
ω² = (4k/η)(q²a²/4)[1 − (qa)²/12 + …]
ω² = (ka²/η)q²[1 − (qa)²/12 + …]
ω² = c²q²[1 − (qa)²/12 + …]
where we used c² = ka²/η = a²(k/η).
2
Leading correction to the speed of light
The group velocity is:
vg = dω/dq = c × cos(qa/2)
For qa « 1:
vg ≈ c [1 − (qa)²/8 + …]
The fractional speed reduction is:
Δv/c = −(qa)²/8 = −(Ea/ℏc)²/8
The group velocity is:
vg = dω/dq = c × cos(qa/2)
For qa « 1:
vg ≈ c [1 − (qa)²/8 + …]
The fractional speed reduction is:
Δv/c = −(qa)²/8 = −(Ea/ℏc)²/8
Correction Magnitudes at Various Scales
3
LHC scale: E = 13 TeV = 2.08 × 10−6 J
qa = Ea/(ℏc) = 2.08 × 10−6 × 1.616 × 10−35 / (1.055 × 10−34 × 2.998 × 108)
qa = 3.36 × 10−41 / 3.163 × 10−26
qa = 1.06 × 10−15
(qa)²/12 = (1.06 × 10−15)² / 12
Δω/ω ≈ 10−31
Completely unobservable at the LHC. The continuum limit is perfect at TeV energies.
qa = Ea/(ℏc) = 2.08 × 10−6 × 1.616 × 10−35 / (1.055 × 10−34 × 2.998 × 108)
qa = 3.36 × 10−41 / 3.163 × 10−26
qa = 1.06 × 10−15
(qa)²/12 = (1.06 × 10−15)² / 12
Δω/ω ≈ 10−31
Completely unobservable at the LHC. The continuum limit is perfect at TeV energies.
4
Cosmic ray scale: E = 1020 eV = 1.6 × 101 J
qa = Ea/(ℏc) = 16 × 1.616 × 10−35 / 3.163 × 10−26
qa = 2.584 × 10−34 / 3.163 × 10−26
qa = 8.17 × 10−9
(qa)²/12 = (8.17 × 10−9)² / 12
Δω/ω ≈ 5.6 × 10−18
Marginally detectable with future cosmic ray observatories. The GZK cutoff may encode lattice effects.
qa = Ea/(ℏc) = 16 × 1.616 × 10−35 / 3.163 × 10−26
qa = 2.584 × 10−34 / 3.163 × 10−26
qa = 8.17 × 10−9
(qa)²/12 = (8.17 × 10−9)² / 12
Δω/ω ≈ 5.6 × 10−18
Marginally detectable with future cosmic ray observatories. The GZK cutoff may encode lattice effects.
5
Planck scale: E = EPlanck = 1.956 × 109 J
qa = π (at the Brillouin zone boundary)
(qa)²/12 = π²/12 = 0.822
ωexact = 2√(k/η) = ωmax
vg = c × cos(π/2) = 0
At the Planck scale, the group velocity drops to zero.
The continuum approximation breaks down by about 4%: the exact ω is ~4% less than the continuum prediction c|q|.
qa = π (at the Brillouin zone boundary)
(qa)²/12 = π²/12 = 0.822
ωexact = 2√(k/η) = ωmax
vg = c × cos(π/2) = 0
At the Planck scale, the group velocity drops to zero.
The continuum approximation breaks down by about 4%: the exact ω is ~4% less than the continuum prediction c|q|.
| Scale | Energy | qa | Δω/ω | Observable? |
|---|---|---|---|---|
| Atomic | ~eV | ~10−28 | ~10−57 | No |
| LHC | 13 TeV | ~10−15 | ~10−31 | No |
| Cosmic ray | 1020 eV | ~10−8 | ~10−17 | Marginal |
| Planck | 1.22 × 1019 GeV | π | ~4% | In principle (BH physics) |
Unique GWT prediction:
Standard QFT predicts ω² = c²q² exactly (linear dispersion to all energies).
GWT predicts ω² = (4k/η)sin²(qa/2) (nonlinear at high energy).
The leading correction −(qa)²/12 is a falsifiable prediction: if ultra-high-energy cosmic rays show no dispersion deviation at the 10−17 level, the lattice spacing must be smaller than 10−35 m. Current data are consistent with a = lP.
Standard QFT predicts ω² = c²q² exactly (linear dispersion to all energies).
GWT predicts ω² = (4k/η)sin²(qa/2) (nonlinear at high energy).
The leading correction −(qa)²/12 is a falsifiable prediction: if ultra-high-energy cosmic rays show no dispersion deviation at the 10−17 level, the lattice spacing must be smaller than 10−35 m. Current data are consistent with a = lP.
§10 — The Equations of Motion
The Euler–Lagrange equations applied to the GWT action yield all of fundamental physics. Varying each field recovers the standard equations of motion.
Varying gμν → Einstein Field Equations
1
Vary the metric
δS/δgμν = 0 ⇒ Rμν − ½gμνR = (8πG/c4) Tμν
The coefficient works out because:
8πG/c4 = 8π × 2c4/(πka) / c4 = 16/(ka)
= 16 / (7.708 × 1043) = 2.076 × 10−43 m/N
Check: 8πG/c4 = 8π × 6.674 × 10−11 / (8.077 × 1033)
= 1.676 × 10−9 / 8.077 × 1033 = 2.076 × 10−43 m/N ✓
Einstein’s field equations emerge exactly.
δS/δgμν = 0 ⇒ Rμν − ½gμνR = (8πG/c4) Tμν
The coefficient works out because:
8πG/c4 = 8π × 2c4/(πka) / c4 = 16/(ka)
= 16 / (7.708 × 1043) = 2.076 × 10−43 m/N
Check: 8πG/c4 = 8π × 6.674 × 10−11 / (8.077 × 1033)
= 1.676 × 10−9 / 8.077 × 1033 = 2.076 × 10−43 m/N ✓
Einstein’s field equations emerge exactly.
Varying Aμ → Maxwell Equations
2
Vary the gauge potential
δS/δAν = 0 ⇒ ∂μFμν = μ0 Jν
In component form, this gives Maxwell’s four equations:
∇ · E = ρ/ε0 (Gauss)
∇ · B = 0 (no magnetic monopoles)
∇ × E = −∂B/∂t (Faraday)
∇ × B = μ0J + μ0ε0 ∂E/∂t (Ampère–Maxwell)
All four equations are orientation-wave dynamics on the lattice.
δS/δAν = 0 ⇒ ∂μFμν = μ0 Jν
In component form, this gives Maxwell’s four equations:
∇ · E = ρ/ε0 (Gauss)
∇ · B = 0 (no magnetic monopoles)
∇ × E = −∂B/∂t (Faraday)
∇ × B = μ0J + μ0ε0 ∂E/∂t (Ampère–Maxwell)
All four equations are orientation-wave dynamics on the lattice.
Varying ψ̅ → Dirac Equation
3
Vary the conjugate spinor
δS/δψ̅ = 0 ⇒ (iγμDμ − m)ψ = 0
This is the Dirac equation. Each component of ψ satisfies the Klein–Gordon equation:
(∂μ∂μ + m²c²/ℏ²)ψ = 0
which is the lattice wave equation for a standing wave with rest mass m.
δS/δψ̅ = 0 ⇒ (iγμDμ − m)ψ = 0
This is the Dirac equation. Each component of ψ satisfies the Klein–Gordon equation:
(∂μ∂μ + m²c²/ℏ²)ψ = 0
which is the lattice wave equation for a standing wave with rest mass m.
Non-Relativistic Limit → Schrödinger Equation
4
Step 1: Start from the lattice wave equation
The massive wave equation on the lattice:
(η/a³) ∂²ψ/∂t² = (k/a) ∇²ψ
In continuum notation (using c² = ka²/η):
∂²ψ/∂t² = c² ∇²ψ
The massive wave equation on the lattice:
(η/a³) ∂²ψ/∂t² = (k/a) ∇²ψ
In continuum notation (using c² = ka²/η):
∂²ψ/∂t² = c² ∇²ψ
5
Step 2: Standing wave rest frequency
A massive particle is a standing wave with rest frequency:
ω0 = mc²/ℏ
For an electron:
ω0 = (9.109 × 10−31) × (2.998 × 108)² / (1.055 × 10−34)
ω0 = 8.187 × 10−14 / 1.055 × 10−34
ω0 = 7.763 × 1020 s−1
This is the electron’s de Broglie frequency — the rate at which the standing wave oscillates.
A massive particle is a standing wave with rest frequency:
ω0 = mc²/ℏ
For an electron:
ω0 = (9.109 × 10−31) × (2.998 × 108)² / (1.055 × 10−34)
ω0 = 8.187 × 10−14 / 1.055 × 10−34
ω0 = 7.763 × 1020 s−1
This is the electron’s de Broglie frequency — the rate at which the standing wave oscillates.
6
Step 3: Slowly varying envelope
Write the total wave function as a fast oscillation times a slow envelope:
ψ(x,t) = φ(x,t) × e−iω0t
where φ varies slowly compared to ω0.
Substituting into the wave equation and dropping ∂²φ/∂t² « ω0 ∂φ/∂t:
−2iω0 ∂φ/∂t = c² ∇²φ
Multiply both sides by −ℏ/(2ω0) = −ℏ²/(2mc²):
iℏ ∂φ/∂t = −(ℏ²/2m) ∇²φ
Write the total wave function as a fast oscillation times a slow envelope:
ψ(x,t) = φ(x,t) × e−iω0t
where φ varies slowly compared to ω0.
Substituting into the wave equation and dropping ∂²φ/∂t² « ω0 ∂φ/∂t:
−2iω0 ∂φ/∂t = c² ∇²φ
Multiply both sides by −ℏ/(2ω0) = −ℏ²/(2mc²):
iℏ ∂φ/∂t = −(ℏ²/2m) ∇²φ
7
Step 4: Add the potential
An external potential V(x) shifts the local rest frequency: ω0 → ω0 + V/ℏ
This adds a Vφ term to the right side:
iℏ ∂φ/∂t = −(ℏ²/2m) ∇²φ + Vφ
This is the time-dependent Schrödinger equation.
It is not a postulate — it is the non-relativistic limit of a classical wave on the lattice.
An external potential V(x) shifts the local rest frequency: ω0 → ω0 + V/ℏ
This adds a Vφ term to the right side:
iℏ ∂φ/∂t = −(ℏ²/2m) ∇²φ + Vφ
This is the time-dependent Schrödinger equation.
It is not a postulate — it is the non-relativistic limit of a classical wave on the lattice.
All four equations of motion recovered:
• Vary gμν: Rμν − ½gμνR = (8πG/c4)Tμν (Einstein)
• Vary Aμ: ∂μFμν = μ0Jν (Maxwell)
• Vary ψ̅: (iγμDμ − m)ψ = 0 (Dirac)
• NR limit: iℏ∂φ/∂t = −(ℏ²/2m)∇²φ + Vφ (Schrödinger)
Every equation of fundamental physics is a special case of wave mechanics on the elastic lattice.
• Vary gμν: Rμν − ½gμνR = (8πG/c4)Tμν (Einstein)
• Vary Aμ: ∂μFμν = μ0Jν (Maxwell)
• Vary ψ̅: (iγμDμ − m)ψ = 0 (Dirac)
• NR limit: iℏ∂φ/∂t = −(ℏ²/2m)∇²φ + Vφ (Schrödinger)
Every equation of fundamental physics is a special case of wave mechanics on the elastic lattice.
§11 — Summary: The Standard Model from Springs
The Standard Model has 25 free parameters that must be measured experimentally. In GWT, every one of these is computed from three lattice constants {k, a, η} — which themselves reduce to {2/π, 1, 2/π} in Planck units, i.e., zero free parameters.
Complete Parameter Map
| SM Parameter | GWT Derivation | GWT Value | Measured | Accuracy |
|---|---|---|---|---|
| G (gravity) | 2c4/(πka) | 6.674 × 10−11 | 6.674 × 10−11 | exact |
| αEM | 6−1/3·exp(−48·8/(9π²)) | 1/137.042 | 1/137.036 | 0.004% |
| αs(MZ) | d²/(2dπ²)·(1+αs²·8/3) | 0.1179 | 0.1179 | 0.08% |
| sin²θW(MZ) | 15/64−3α/2+αln(F)/7 | 0.23123 | 0.23122 | 0.019% |
| me | mp/(6π5) | 0.5110 MeV | 0.5110 MeV | 0.002% |
| mμ/me | Koide extended | 206.77 | 206.768 | 0.005% |
| mp | 4ΛQCD | 938.3 MeV | 938.3 MeV | exact |
| rp | j0 cavity | 0.841 fm | 0.841 fm | exact |
| Higgs VEV (v) | (5/2)mPlα8 | 245.4 GeV | 246.2 GeV | 0.3% |
| λ (Higgs quartic) | m(8,24)+VP → (MH/v)²/2 | 0.1295 | 0.129 | 0.4% |
| fπ | pion decay from αs | 93.8 MeV | 92.1 MeV | 1.9% |
| Δm²31 (neutrino) | me³/(Ncmp²) | 2.523 × 10−3 eV² | 2.534 × 10−3 eV² | 0.4% |
| H0 | lattice growth | 66.4 km/s/Mpc | 67.4 km/s/Mpc | 1.5% |
| ΩΛ | (d−1)/d = 2/3 | 0.667 | 0.685 | 2.7% |
| SU(3) × SU(2) × U(1) | DOF: 3 + 2 + 1 | exact match | exact match | exact |
| Ngen = 3 | d = 3 spatial | 3 | 3 | exact |
The Central Result
The Standard Model of particle physics — with its 25 free parameters, its unexplained gauge group, its mysterious Higgs mechanism, and its inexplicable number of generations — is the continuum limit of a 3D elastic lattice with binary internal structure.
Every term in the Standard Model Lagrangian maps to a mechanical term in the lattice Lagrangian:
- Einstein–Hilbert R ← lattice compression energy (ka/32 = c4/(16πG))
- Maxwell FμνFμν ← orientation gradient energy
- Dirac ψ̅(iγμDμ − m)ψ ← standing wave + yin–yang coupling
- Higgs V(φ) ← double-well orientation potential
- SU(3) × SU(2) × U(1) ← 3 spatial + 2 yin–yang + 1 phase DOF
The lattice has zero free parameters (k = η = 2/π, a = 1 in Planck units). The Standard Model has 25. The discrepancy is resolved: the 25 parameters were never free — they are computable outputs of geometry.
Physics is not complicated. Physics is springs.