Appendix

Five-Term Functional and Robustness Scans

The free-energy functional that produces stable condensations, and the robustness scans that show the 3+e configuration is the minimum-energy outcome across parameter space.

Appendix

The Full Five-Term Functional and Robustness Scans

1. The Full Five-Term Functional

The Big Flare-Up Theory (BFUT) condensation energy functional has five terms. Each term encodes a distinct physical mechanism. Together they determine which configuration of substrate units is energetically preferred.

The full functional evaluated for a configuration of n units with k co-rotating:

E = coop + imb + geom + c3ph

1.1 Term by Term

Term 1 - Cooperation (coop)

coop = -J x pairs_sum(s)

pairs_sum(s) = sum of s_i x s_j over all distinct pairs i < j

Physical meaning: Co-rotating units attract each other by the Bernoulli mechanism. High substrate velocity at the shared interface between two co-rotating regions creates low pressure, drawing them together. The cooperation energy grows with the number of co-rotating pairs.

k co-rotating Pairs Binding energy
1 0 0 (no pairs, unstable)
2 1 -J = -1.0 (marginal)
3 3 -3J = -3.0 (qualitative jump - first stable nucleus)
4 6 -6J = -6.0

The jump from k=2 (one pair, -J) to k=3 (three pairs, -3J) is qualitative not gradual. This is why the 3-core is the first stable nucleus. Below k=3 the core cannot survive substrate fluctuations.

Term 2 - Imbalance Penalty (imb)

imb = λ x (sum(s))^2

Physical meaning: A net circulation asymmetry costs energy. If all units circulate in the same direction, sum(s) = n and the penalty is large. The balanced 2+2 configuration has sum(s) = 0 and zero penalty. The 3+1 configuration has sum(s) = 3-1 = 2, giving a moderate penalty λ x 4 = 2.4.

Term 3 - Geometric Cost (geom)

geom = (k-3)^2 + α x (n-k)

Physical meaning: Two independent geometric costs. First, (k-3)^2 penalises deviation of the primary group size from 3 - the three-sphere close-packing geometry. Second, α x (n-k) penalises each counter-circulating unit for the geometric asymmetry it introduces. When the expelled unit has mass fraction μ, this term scales as α x μ.

Term 4 - Circulation Phase Reward (c3ph)

c3ph = D_s x cos(3 x φ)

Physical meaning: The fifth term explicitly encodes the topology of the three-sphere packing into the energy functional. The phase φ measures the circulation configuration:

Config φ cos(3φ) c3ph = D_s x cos(3φ) Effect
3+1 π/3 -1 -D_s = -1.5 Rewarded
2+2 π/2 0 0 Neutral
4+0 0 +1 +D_s = +1.5 Penalised

D_s must be positive. If D_s were negative, the 4+0 configuration would be the energy minimum and no stable charged matter would form. The fifth term gives E(4+0) = 6.10, with a 1D robustness of 97.56%, and explicitly encodes the three-sphere topology into the functional.

1.2 Parameters

Symbol Value Name Physical role
J 1.0 Cooperation strength Bernoulli binding per co-rotating pair
λ 0.6 Imbalance penalty Cost of net circulation asymmetry
α 0.5 Geometric asymmetry Cost per counter-circulating unit (scales with μ for expelled unit)
D_s 1.5 Phase reward cos(3φ) circulation topology reward. Must be positive.

2. The Per-Unit Energy Scan (Code 1)

Code 1 answers: for n co-rotating substrate units, which n minimises energy per unit E(n)/n? All units are at the primary phase φ = π/3, so cos(3φ) = -1 and c3ph = -D_s for all n.

n E(n) E(n)/n Note
1 3.1000 3.1000 No pairs. Unstable.
2 0.9000 0.4500 One pair. Marginal.
3 0.9000 0.3000 MINIMUM E/unit. The 3-core attractor.
4 3.1000 0.7750
5 7.5000 1.5000 Rising steeply
6-12 ... ... Continues to rise

Figure A. Total condensation energy E(n) for n=1 to 12. n=3 highlighted.

Figure B. Energy per unit E(n)/n for n=1 to 12. n=3 is the unambiguous minimum.

3. The Four-Unit Partition (Post 3-Core Formation)

The moment n=3 forms, the three-sphere packing geometry simultaneously creates the interstitial region. The four-unit bound system forms at E=1.400 model units. The three configurations and their full-functional energies:

Config E (model units) cos(3φ) Status
3+1 1.4000 -1 (rewarded) MINIMUM. SELECTED.
2+2 4.0000 0 (neutral) Symmetric. No net charge.
4+0 6.1000 +1 (penalised) All co-rotating. Penalised by D_s.

Figure C. Four-unit partition energies. 3+1 is the clear minimum.

Figure D. Energy through the three stages of proton formation.

4. Three-Sphere Packing Geometry

Three substrate condensations of radius r_q in close-packed contact. The three centres form an equilateral triangle of side 2r_q. The outer radius of the assembly equals r_p, the measured proton charge radius. This is the only measured input.

r_outer = r_q x (1 + 2/sqrt(3)) = 2.1547 x r_q = r_p

r_q = r_p / (1 + 2/sqrt(3)) = 0.8414 / 2.1547 = 0.3905 fm

V_gap / V_q = (2*sqrt(3) - π) / (4*π/3) = 0.0770

Both results are universal geometric constants. No free parameters.

Figure E. Three-sphere packing geometry. Green arrows: co-rotating quarks. Red: interstitial unit (counter-rotates). Yellow dashed: outer radius = r_p.

5. Why the Interstitial Unit Counter-Rotates

The counter-rotation is mechanically imparted, not assumed. When the interstitial substrate exits through the gap between any two quarks, it encounters two co-rotating surfaces - one on each side. Both quarks rotate in the same direction. Each imparts a tangential force in the opposite direction to the passing substrate. Together they impart a net counter-clockwise torque.

This is the gear analogy: a gear placed between two co-rotating gears of the same handedness always rotates in the opposite direction. The result is the same regardless of which gap the substrate exits from, because all three quarks rotate in the same direction.

Counter-rotation in BFUT is the definition of opposite charge. The negative charge of the expelled unit is therefore not assigned or assumed. It is mechanically imparted during expulsion by the same co-rotation that defines the quarks as positively charged.

6. Proton Energy with Interstitial Unit

When the detached unit has mass fraction μ relative to a core unit, the P16 functional is evaluated with s = [+1, +1, +1, -μ]. The geometric asymmetry term scales with μ because the geometric displacement is proportional to the actual mass of the detached unit:

E(μ) = -J x pairs_sum([1,1,1,-μ])

+ λ x (3-μ)^2

+ (3-3)^2 + α x μ

+ D_s x (-1)

At μ=1.0 (standard 3+1): E = 1.4000 model units (baseline confirmed)

Minimum: E = 0.8958 model units at μ = 0.083

Driving force for expulsion: -0.504 model units

μ E(μ) Reduction from 1.400 Note
1.000 1.4000 0.0000 Standard 3+1 baseline
0.500 1.0000 -0.4000
0.230 0.9087 -0.4913 P19 upper bound
0.083 0.8958 -0.5042 MINIMUM
0.077 0.8959 -0.5041 P19 lower bound

The result is robust across the full physical range μ = 0.077 to 0.23. Proton stability does not depend on fine-tuning the interstitial mass fraction.

Figure F. E(μ) vs μ across the full range 0 to 1. Minimum at μ=0.083. Green band: physical range.

Figure G. Proton energy across the physical range μ=0.077 to 0.230. Result is robust.

7. The Connecting Identity

P19 Section 21.6 establishes an algebraic identity connecting the interstitial volume fraction, the electron mass, and the compression energy. E_unit cancels, so the ratio is purely geometric:

E_gap = E_unit x V_gap/V_q = 298.661 x 0.0770 = 22.997 MeV

m_e_vss = E_unit / (6 x π^4) = 298.661 / 584.45 = 0.511009 MeV

E_gap / m_e_vss = 6 x π^4 x V_gap/V_q = 584.45 x 0.076996 = 44.9991, approximately 45.0

Physical meaning: The expelled substrate dissipates approximately 44/45 of the compression energy into the surrounding substrate during proton stabilisation. The remaining approximately 1/45 is retained as the stable counter-rotating condensate whose mass is m_e_vss. The chain from V_gap to E_gap to m_e_vss is one identity with E_unit as the common factor that cancels.

8. Modular Organisation: The Universal Structural Unit

For n = 3m co-rotating core units with m > 1, m independent three-core structures are always energetically preferred over a single large co-rotating condensate. The energy advantage grows with n. At n=24 the gap is 501.900 model units.

n m E single E modular Gap Winner
3 1 0.900 0.900 0.000 EQUAL
6 2 14.100 1.800 +12.300 MODULAR
9 3 47.100 2.700 +44.400 MODULAR
12 4 99.900 3.600 +96.300 MODULAR
24 8 509.100 7.200 +501.900 MODULAR

The comparison counts co-rotating core units only. The three-core unit is the preferred structural unit for matter at all scales, and its electron is generated subsequently by interstitial ejection.

Figure H. Total energy: single condensate vs modular units for n = 3m, from n = 3 to 24. The two coincide at n = 3.

Figure I. Energy advantage of modular organisation. ΔE = (n − 3)(1.1n − 2.5). Positive = modular wins. Gap grows with n.

9. Robustness: Parameter Space Analysis

The 3+1 selection is not a fragile result at a single parameter point. Scanning λ and α across [0.2, 1.2] with J=1.0 fixed:

Scan 3+1 wins Parameters varied
1D 97.56% λ in [0.2, 1.2]
2D 95.95% λ x α in [0.2, 1.2]^2
3D 90.43% λ x α x D_s in [0.2,1.2]^2 x [0.5,2.5]

Figure J. Parameter space map. Green: 3+1 is minimum energy. Red: other configuration wins. P16 working point marked.

10. All Key Results at a Glance

Quantity Value Source
r_p 0.8414 fm CODATA 2018 measurement
r_q 0.3905 fm Three-sphere geometry
V_gap / V_q 0.0770 Universal geometric constant
E_unit = m_p/π 298.661 MeV Proton mass formula
m_e_vss = E_unit/(6*π^4) 0.511009 MeV Electron mass (measured: 0.510999)
E_gap/m_e_vss 45.00 Connecting identity - E_unit cancels
E(3+1) baseline 1.4000 model units P16 functional, standard 3+1
E minimum (μ=0.083) 0.8958 model units P16 functional, interstitial expelled
Driving force -0.504 model units Energetic basis for expulsion
2D robustness 95.95% Five-term functional scan
Confinement F 0.574 GeV/fm vs QCD 0.900 GeV/fm (36%)