What Happens When You Fire an Electron at a Proton?
Three Energy Regimes, Three Answers, One Algebra
Three Energy Regimes, Three Answers, One Algebra
April 2026
We take one physical system — an electron approaching a proton — and ask the simplest possible question: what comes out? The answer depends entirely on how hard the electron hits. At 13.6 eV, nothing flies out at all: the electron is captured, and you get a hydrogen atom. At 100 MeV, the electron bounces off and the proton recoils, but both remain intact. At 4 GeV, the proton gets excited and spits out a pion. At 20 GeV, the proton shatters, and you discover that it was built from quarks all along. Every number in this document is computed from first principles. No magic, no secret constants — just masses,
1The Setup: One Electron, One Proton, One Question
Imagine you have a proton sitting on a table. You pick up an electron and throw it at the proton. What happens? The answer depends on one thing: how fast you throw. And “how fast” is really “how much energy.” In physics, we measure this energy in electron-volts (eV). One electron-volt is the energy an electron gains when it falls through a one-volt battery — about 1.6 × 10⁻¹⁹ joules — the natural unit for atomic and particle physics. The beautiful thing about this experiment is that as you increase the energy, you peel back layers of reality like an onion. At low energy, you see the proton as a point. At medium energy, you see it as a ball. At high energy, you see inside the ball and discover it is made of smaller things. The mathematics changes at each layer — from quaternions to octonions — and this document shows you exactly where and why.
11 The two numbers that govern everything
Two numbers control all that follows.
Proton radius: rₚ = 0.8414 fm (1 femtometre = 10⁻¹⁵ m, a million times smaller than an atom)
Electron de Broglie wavelength: λ = h/p = 2πħ
Think of the electron as a flashlight: its wavelength is the size of the smallest detail it can illuminate. The fundamental rule of wave physics is that you cannot resolve anything smaller than your wavelength.
(1)λ = 2πħ
When λ ≫ rₚ: the electron sees a featureless dot. When λ ≈ rₚ: it starts to see shape. When λ ≪ rₚ: it sees inside. The electron does not change. The proton does not change. What changes is the resolution. The resolution criterion λ ≫ rₚ (= 0.84 fm): proton looks like a point → Regime 1 and 2 λ ≈ rₚ: proton’s surface visible → Regime 3 λ ≪ rₚ: interior visible → Regime 4 The transition energy (where λ = rₚ) is p = 2πħ
12 Fundamental constants used throughout
ħ
mᵉ = 0.511 MeV/
Mₚ = 938.272 MeV/
Mπ⁰ = 134.977 MeV/
Mπ⁺ = 139.570 MeV/
MΔ = 1232 MeV/
2Regime 1: The Gentle Touch (13.6 eV)
Nothing flies out — you get a hydrogen atom
At room temperature, an electron has about 0.025 eV of kinetic energy. Even at 13.6 eV, it is barely crawling by particle physics standards. Let us compute its wavelength from equation (1). At 13.6 eV, the electron momentum (non-relativistic, since 13.6 eV ≪ mᵉ
21 The hydrogen energy levels
The energy levels of hydrogen are one of physics’ most beautiful results. In quaternion language, the electron’s state is described by:
Q = E + Lₓ·ι1 + Lᵧ·ι2 + Sᵣ·ι3
where E is the energy, Lₓ, Lᵧ are angular momentum components, Sᵣ is spin, and ι1, ι2, ι3 are the three quaternion imaginary units. The energy eigenvalues come from the requirement that the quaternion norm be preserved under rotation: Eₙ = −mᵉ
3Regime 2: The Hard Bounce (100 MeV)
Elastic scattering — the proton stays intact
Now throw the electron much harder: 100 MeV, about 7.4 million times the energy of the hydrogen case. The electron is now ultra-relativistic. We compute the Lorentz factor: γ = E/mᵉ
λ = 2πħ
The electron’s flashlight is 15 times wider than the proton — still much too coarse to see inside. The proton stays intact. The electron bounces off the proton’s Coulomb field.
31 Elastic scattering kinematics
In elastic scattering, both particles are the same before and after. Only the momenta change direction. The scattered electron energy E′ at angle θ is given exactly by:
(5)E' = E / ( 1 + (2E/Mₚ) sin²(θ/2) ) [elastic, proton at rest]
Let us compute E′ and the kinetic energy given to the proton Tₚ = E − E′ for several angles: Even at backscatter (180°), the proton receives only 17.7 MeV and moves at 0.19c. The 1836:1 mass ratio means the electron is like a marble bouncing off a bowling ball.
32 The Mott cross section
Not all angles are equally likely. The differential cross section for relativistic electron scattering from a point charge is the Mott formula:
(2)dσ/dΩ = (
The sin⁴(θ/2) factor in the denominator makes forward scattering overwhelmingly more probable. Let us compute dσ/dΩ in units of fm²/sr at 100 MeV:
33 Centre-of-mass energy and the pion threshold
Can we break the proton at 100 MeV? We need to compute the centre-of-mass energy √s — the total energy available for particle creation. For a beam electron of energy E hitting a stationary proton: √s = √( (E + Mₚ)² − pᵉ² ) = √( Mₚ² + 2MₚE ) [target at rest] (3) At E = 100 MeV: √s = √(938.3² + 2×938.3×100) = √(896,660) = 947.0 MeV. To produce even the lightest new particle (a neutral pion, π⁰), we need:
(4)√sₜℍ = Mₚ + Mπ = 938.3 + 135.0 = 1073.3 MeV
We have 947.0 MeV and need 1073.3 MeV. We are 126.3 MeV short. The proton is safe. The minimum beam energy to produce a single pion is found by setting √s = 1073.3 MeV:
Eₜℎʳ = (√sₜℎʳ² − Mₚ²) / (2Mₚ) = (1073.3² − 938.3²)/(2×938.3) = (1,152,013 − 880,447)/1876.6 = 144.7 MeV
Below 145 MeV beam energy, elastic scattering is the only available process. This makes the 100 MeV regime clean: two particles in, two particles out, nothing created. Regime 2 summary: 100 MeV Wavelength: λ = 12.4 fm (λ/rₚ = 14.7) Lorentz factor: γ = 195.7 (β = 0.999987) What comes out: e⁻ (deflected) + p (recoils) Energy to proton: 0.16 to 17.7 MeV depending on angle Cross section: Mott formula eq.(2) — forward strongly favoured Pion threshold: needs 145 MeV. At 100 MeV, proton is safe.
4Regime 3: The First Cracks (4 GeV)
The proton gets excited — and breaks
At 4 GeV the electron has 40 times more energy than in Regime 2. Its Lorentz factor γ = 4000/0.511 = 7,828. Its de Broglie wavelength:
λ = 2πħ
For the first time λ < rₚ. The electron can now see structure inside the proton. The centre-of-mass energy from equation (3):
√s = √(938.3² + 2×938.3×4000) = √(8,347,633) = 2890 MeV = 2.89 GeV
41 Reaction channels open at 4 GeV
Every channel whose threshold √sₜℎʳ ≤ 2.89 GeV is available:
42 The Δ(1232) resonance: the proton’s first excited state
The dominant process is excitation of the Δ⁺(1232) resonance. Just as hydrogen can absorb a photon and jump to a higher level, the proton absorbs energy and becomes the Δ. The mass difference:
(11)ΔE = MΔ − Mₚ = 1232 − 938.3 = 293.7 MeV
The Δ has spin 3/2 (compared to 1/2 for the proton), isospin 3/2, and decays in τ = 5.6 × 10⁻²⁴ s. At the speed of light, it travels only cτ = 1.69 fm — barely twice the proton radius — before decaying. It is born and dies within the proton’s volume. Decay modes: Δ⁺ → p + π⁰ (33%): Proton survives. Neutral pion flies out. π⁰ then decays: π⁰ → γγ (lifetime 8.4×10⁻¹⁷ s). Two photons hit the calorimeter. Δ⁺ → n + π⁺ (67%): Proton converts to neutron. Charged pion escapes. π⁺ is long-lived (τ = 26.0 ns) and traverses the detector.
43 A specific event: 4 GeV, θ = 10°
Let us trace one event in full. A 4 GeV electron hits a stationary proton and scatters at 10°. First compute the elastic scattered energy from equation (4):
E'ᵉˡ = 4000 / (1 + (2×4000/938.3)×sin²(5°)) = 4000 / (1 + 8.521×0.00760) = 4000/1.0648 = 3757 MeV
The momentum transfer squared from equations (5) and (6):
Q² = 4 × 4000 × 3757 × sin²(5°) = 60,112,000 × 0.00760 = 0.457 GeV²
If instead the proton gets excited into the Δ(1232), the virtual photon must carry exactly MΔ − Mₚ = 293.7 MeV of invariant mass. The electron comes out with E′ = 3757 − 293.7 = 3463 MeV, losing an extra 294 MeV compared to elastic scattering. This 294 MeV “missing energy” is the experimental signature of the resonance. In the Δ’s rest frame, the decay pion gets:
p_π* = √[ (MΔ² - (Mₚ+Mπ)²)(MΔ² - (Mₚ-Mπ)²) ] / (2MΔ) = 229 MeV/
T_π = √(p_π² + Mπ²) - Mπ = √(52441 + 18225) - 135 = 266 - 135 = 131 MeV The Δ is boosted forward in the lab with βΔ = 0.571, γΔ = 1.22. Boosting the pion forward:
p_π(lab) ≈ γΔ(p_π* + βΔE_π*) = 1.22 × (229 + 0.571×266) = 1.22 × 381 = 465 MeV/
Regime 3 summary: 4 GeV Wavelength: λ = 0.310 fm (λ/rₚ = 0.368) Centre-of-mass energy: √s = 2.89 GeV All pion channels open. Δ(1232) dominant. Elastic E′ at 10°: 3757 MeV | Inelastic (via Δ) E′: 3463 MeV Signature: 294 MeV missing energy + forward pion at ~465 MeV/
5Regime 4: Shattering the Proton (20 GeV)
Deep inelastic scattering — the electron sees quarks
In 1968, at SLAC, physicists fired 20 GeV electrons at stationary protons. The de Broglie wavelength:
λ = 2πħ
The electron now has 1/0.074 ≈ 14 “pixels” across the proton. It does not see a ball. It sees individual point-like objects inside. This was the experimental discovery of quarks.
51 The kinematic variables of deep inelastic scattering
In deep inelastic scattering (DIS), the electron exchanges a virtual photon with the proton. Four key variables:
(6)Q² = 4EE′ sin²(θ/2)
(7)ν = E − E′
(8)x = Q² / (2Mₚν)
(9)W² = Mₚ² + 2Mₚν − Q²
Q² is the momentum transfer squared: higher Q² means finer resolution (smaller wavelength). ν is the energy transferred to the proton. x (the Bjorken variable) is the fraction of the proton’s momentum carried by the struck quark. W is the invariant mass of the hadronic debris.
52 A single event in full detail: θ = 10°, E′ = 12 GeV
Before: electron E = 20 GeV rightward. Proton at rest E = 0.938 GeV. Total energy 20.938 GeV. Total forward momentum 20.000 GeV/
Q² = 4 × 20 × 12 × sin²(5°) = 960 × 0.00760 = 7.30 GeV²
ν = 20 − 12 = 8.0 GeV
x = 7.30 / (2 × 0.938 × 8.0) = 7.30 / 15.01 = 0.487
W = √(0.938² + 2×0.938×8.0 − 7.30) = √(0.880 + 15.008 − 7.30) = √8.59 = 2.93 GeV
Interpretation: the electron struck a quark carrying 48.7% of the proton’s momentum. At this x value, it is most likely an up quark (the proton contains two up quarks and one down quark; up quarks are favoured at high x). The hadronic debris has invariant mass W = 2.93 GeV and total energy Eℎʳˡ = ν + Mₚ = 8.938 GeV, flying forward at angle θℎʳˡ ≈ arctan(p⊥/p∥). Average particle multiplicity at W = 2.93 GeV is approximately 5–6 hadrons. Typical particle list from this event: 1 proton or neutron (2–6 GeV) The two “spectator” quarks that were not struck, carrying most of the forward momentum. 2–3 π⁺ pions (1–3 GeV each) Quark-antiquark pairs formed when the struck quark was ripped out. 1–2 π⁻ pions (0.5–2 GeV each) Charge conservation requires roughly equal π⁺ and π⁻ production. 1–2 π⁰ pions Decay instantly to two photons; hit the electromagnetic calorimeter. Occasionally 1 kaon (K⁺ or K⁰) If a strange quark pair was created from the vacuum. Conservation check at this event: Regime 4 summary: 20 GeV Wavelength: λ = 0.062 fm (λ/rₚ = 0.074) Resolution: 14 “pixels” across the proton Bjorken x = Q²/(2Mν): fraction of proton momentum carried by struck quark At x = 0.487: most likely hitting a valence up quark Hadronic debris: W = 2.93 GeV, 5-6 particles, all forward What comes out: e⁻ (deflected) + hadronic jet
ep_collision.py — ep_collision
===========================================================================
ELECTRON-PROTON COLLISION: REAL NUMBERS
Electron beam hits stationary proton (fixed target)
===========================================================================
===========================================================================
SCENARIO A: 100 MeV electron on stationary proton
===========================================================================
--- Incoming electron ---
Kinetic energy: T = 100.0 MeV
Total energy: E = 100.51 MeV
Momentum: p = 100.51 MeV/c
Lorentz gamma: 196.7
Speed: beta = 0.999987 c
de Broglie lam: 12.34 fm
lam / r_proton: 14.7 -> cannot resolve proton interior
--- Centre-of-mass frame ---
sqrt(s) = 1033.91 MeV = 1.0339 GeV
Available energy above rest masses: 95.13 MeV
--- Pion production threshold ---
Need sqrt(s) >= 1073.8 MeV to produce even one pi0
That requires E_beam = 145.3 MeV (T = 144.8 MeV)
We have sqrt(s) = 1033.9 MeV -> BELOW threshold
--- What happens (elastic scattering) ---
The proton STAYS INTACT. Only two particles come out:
1. Scattered electron (deflected)
2. Recoiling proton
--- Scattered electron (lab frame) ---
theta_lab Ee' (MeV) Te' (MeV) dE (MeV) E_recoil Tp (MeV)
--------------------------------------------------------------------
10 deg 100.35 99.84 0.16 938.44 0.16
20 deg 99.87 99.35 0.65 938.92 0.65
30 deg 99.09 98.58 1.42 939.69 1.42
45 deg 97.45 96.94 3.06 941.33 3.06
60 deg 95.40 94.89 5.11 943.38 5.11
90 deg 90.79 90.27 9.73 948.00 9.73
120 deg 86.60 86.09 13.91 952.19 13.91
150 deg 83.77 83.26 16.74 955.02 16.74
180 deg 82.78 82.27 17.73 956.01 17.73
--- Recoiling proton ---
At theta_e = 90 deg:
Electron: E' = 90.79 MeV, theta = 90 deg
Proton: T = 9.73 MeV, phi = 42.1 deg (forward)
Proton speed: beta = 0.1429 c
-> Proton barely moves. It is 1836x heavier than the electron.
--- Cross section (Mott, point proton) ---
theta = 10 deg: dsig/dOmega = 0.8811 fm^2/sr = 8.8108 mbarn/sr
theta = 30 deg: dsig/dOmega = 0.0105 fm^2/sr = 0.1052 mbarn/sr
theta = 90 deg: dsig/dOmega = 0.0001 fm^2/sr = 0.0009 mbarn/sr
===========================================================================
SCENARIO B: 4 GeV electron on stationary proton
===========================================================================
--- Incoming electron ---
Kinetic energy: T = 4.0 GeV
Momentum: p = 4.0005 GeV/c (ultra-relativistic: E ~ pc)
Lorentz gamma: 7829
de Broglie lam: 0.3099 fm
lam / r_proton: 0.3683 -> RESOLVES proton interior\!
--- Centre-of-mass ---
sqrt(s) = 2.896 GeV
Available hadronic mass: W = 2896 MeV
--- Production thresholds vs our sqrt(s) = 2896 MeV ---
[OPEN ] e + p -> e + p (elastic)
threshold: sqrt(s) = 939 MeV = 0.939 GeV
[OPEN ] e + p -> e + p + pi0
threshold: sqrt(s) = 1074 MeV = 1.074 GeV
[OPEN ] e + p -> e + n + pi+
threshold: sqrt(s) = 1080 MeV = 1.080 GeV
[OPEN ] e + p -> e + Delta(1232) -> e + p + pi
threshold: sqrt(s) = 1233 MeV = 1.233 GeV
[OPEN ] e + p -> e + p + pi+ + pi-
threshold: sqrt(s) = 1218 MeV = 1.218 GeV
[OPEN ] e + p -> e + p + pi0 + pi0
threshold: sqrt(s) = 1209 MeV = 1.209 GeV
[OPEN ] e + p -> e + p + rho0(->pipi)
threshold: sqrt(s) = 1714 MeV = 1.714 GeV
[OPEN ] e + p -> e + p + 3pi
threshold: sqrt(s) = 1357 MeV = 1.357 GeV
[OPEN ] e + p -> e + p + p + pbar
threshold: sqrt(s) = 2815 MeV = 2.815 GeV
--- Scattered electron at selected angles ---
For ELASTIC scattering (proton stays intact):
theta_lab Ee' (MeV) Q^2 (GeV^2) lam_probe (fm)
--------------------------------------------------
6 deg 3909.2 0.171 2.995
10 deg 3757.1 0.457 1.835
20 deg 3182.3 1.536 1.001
30 deg 2546.1 2.729 0.750
45 deg 1778.9 4.169 0.607
--- The Delta(1232) resonance ---
The first thing the electron 'breaks off' the proton.
p + gamma* -> Delta++, Delta+, Delta0, or Delta-
Delta(1232) has spin 3/2, isospin 3/2
Decays in ~5.6 x 10^-24 s (barely exists\!)
Delta+ -> p + pi0 (BR ~33%%)
-> n + pi+ (BR ~67%%)
Lifetime: tau = hbar/Gamma, where Gamma ~ 117 MeV
tau = hbar/Gamma = 5.63e-24 s
Travels: c*tau = 1.69 fm before decaying
-> It decays INSIDE the proton\! You never see the Delta directly.
You see its decay products: a proton/neutron + pion(s).
--- What the detector sees at 4 GeV ---
FORWARD (small angle):
- Scattered electron (most of the beam energy)
- Hadronic debris: proton or neutron + pions, mostly forward
LARGE ANGLE (theta > 30 deg):
- Hard-scattered electron (lost a LOT of energy)
- Multiple pions spraying forward
Typical event at theta = 10 deg, elastic:
Electron out: E' = 3757 MeV at theta = 10 deg
Proton out: T = 243 MeV at phi = 65.3 deg
Proton speed: beta = 0.6079 c
Typical inelastic event (Delta production) at theta = 10 deg:
Electron out: E' = 3438 MeV at theta = 10 deg (lost 562 MeV\!)
Q^2 = 0.418 GeV^2
The 562 MeV goes into creating the Delta(1232)
Delta -> proton (938 MeV) + pi0 (135 MeV)
or Delta -> neutron (940 MeV) + pi+ (140 MeV)
In Delta rest frame: pion gets T = 131 MeV, proton gets T = 28 MeV
Pion momentum (Delta frame): p = 229 MeV/c
Delta boost: gamma = 1.22, beta = 0.5710
Forward pion in lab: p ~ 464 MeV/c
-> Pion flies FORWARD (same direction as beam)
===========================================================================
SCENARIO C: 20 GeV electron on stationary proton (SLAC 1968)
===========================================================================
--- Incoming electron ---
Energy: E = 20.001 GeV
Momentum: p = 20.001 GeV/c
de Broglie lam: 0.0620 fm
lam / r_proton: 0.0737
-> Resolves to 1/14 of proton radius\!
This electron sees INSIDE the proton.
--- Centre-of-mass ---
sqrt(s) = 6.20 GeV
Enough energy for: ~38 pions
--- Deep inelastic scattering kinematics ---
The electron hits a SINGLE QUARK inside the proton.
[output truncated at 160 lines — run the script for the rest]
#\!/usr/bin/env python3
"""
Real-world electron-proton collision: full kinematics.
What goes in, what comes out, at what angle, with what energy.
"""
import numpy as np
# CONSTANTS
m_e = 0.51099895 # electron mass, MeV/c2
m_p = 938.272088 # proton mass, MeV/c2
m_pi0 = 134.977 # neutral pion mass, MeV/c2
m_pip = 139.570 # charged pion mass, MeV/c2
m_n = 939.565378 # neutron mass, MeV/c2
m_Delta = 1232.0 # Delta(1232) resonance mass, MeV/c2
m_rho = 775.26 # rho meson mass, MeV/c2
alpha = 1/137.036 # fine structure constant
hbar_c = 197.3269804 # hbar*c in MeV*fm
r_p = 0.8414 # proton charge radius, fm
hbar_s = 6.582119e-22 # hbar in MeV*s
def de_broglie(p_MeV):
return 2 * np.pi * hbar_c / p_MeV
def lorentz_gamma(T, m):
return (T + m) / m
def beta_func(T, m):
g = lorentz_gamma(T, m)
return np.sqrt(1 - 1/g**2)
def cm_energy_sq(E_lab, m_beam, m_target):
return m_beam**2 + m_target**2 + 2 * m_target * E_lab
print("=" * 75)
print("ELECTRON-PROTON COLLISION: REAL NUMBERS")
print("Electron beam hits stationary proton (fixed target)")
print("=" * 75)
# ============================================================
# SCENARIO A: 100 MeV ELECTRON - ELASTIC
# ============================================================
print("\n" + "=" * 75)
print("SCENARIO A: 100 MeV electron on stationary proton")
print("=" * 75)
T_e = 100.0
E_e = T_e + m_e
p_e = np.sqrt(E_e**2 - m_e**2)
lam = de_broglie(p_e)
gamma_e = lorentz_gamma(T_e, m_e)
beta_e = beta_func(T_e, m_e)
print("\n--- Incoming electron ---")
print(" Kinetic energy: T = %.1f MeV" % T_e)
print(" Total energy: E = %.2f MeV" % E_e)
print(" Momentum: p = %.2f MeV/c" % p_e)
print(" Lorentz gamma: %.1f" % gamma_e)
print(" Speed: beta = %.6f c" % beta_e)
print(" de Broglie lam: %.2f fm" % lam)
resolves = "cannot" if lam/r_p > 1 else "CAN"
print(" lam / r_proton: %.1f -> %s resolve proton interior" % (lam/r_p, resolves))
s = cm_energy_sq(E_e, m_e, m_p)
sqrt_s = np.sqrt(s)
print("\n--- Centre-of-mass frame ---")
print(" sqrt(s) = %.2f MeV = %.4f GeV" % (sqrt_s, sqrt_s/1000))
print(" Available energy above rest masses: %.2f MeV" % (sqrt_s - m_e - m_p))
threshold_sqrt_s = m_e + m_p + m_pi0
E_threshold_lab = (threshold_sqrt_s**2 - m_e**2 - m_p**2) / (2 * m_p)
T_threshold = E_threshold_lab - m_e
print("\n--- Pion production threshold ---")
print(" Need sqrt(s) >= %.1f MeV to produce even one pi0" % (m_e + m_p + m_pi0))
print(" That requires E_beam = %.1f MeV (T = %.1f MeV)" % (E_threshold_lab, T_threshold))
below = "BELOW" if sqrt_s < threshold_sqrt_s else "ABOVE"
print(" We have sqrt(s) = %.1f MeV -> %s threshold" % (sqrt_s, below))
print("\n--- What happens (elastic scattering) ---")
print(" The proton STAYS INTACT. Only two particles come out:")
print(" 1. Scattered electron (deflected)")
print(" 2. Recoiling proton")
print("\n--- Scattered electron (lab frame) ---")
print(" %8s %12s %12s %10s %12s %10s" % ("theta_lab", "Ee' (MeV)", "Te' (MeV)", "dE (MeV)", "E_recoil", "Tp (MeV)"))
print(" " + "-" * 68)
angles = [10, 20, 30, 45, 60, 90, 120, 150, 180]
for theta_deg in angles:
theta = np.radians(theta_deg)
sin2 = np.sin(theta/2)**2
E_prime = E_e / (1 + 2*E_e*sin2/m_p)
T_prime = E_prime - m_e
delta_E = E_e - E_prime
T_p = delta_E
E_p_total = m_p + T_p
print(" %6d deg %12.2f %12.2f %10.2f %12.2f %10.2f" % (theta_deg, E_prime, T_prime, delta_E, E_p_total, T_p))
print("\n--- Recoiling proton ---")
print(" At theta_e = 90 deg:")
theta_e = np.pi/2
E_prime_90 = E_e / (1 + 2*E_e*np.sin(theta_e/2)**2/m_p)
T_p_90 = E_e - E_prime_90
p_p = np.sqrt(T_p_90**2 + 2*T_p_90*m_p)
p_e_prime = np.sqrt(E_prime_90**2 - m_e**2)
phi = np.arctan2(p_e_prime * np.sin(theta_e), p_e - p_e_prime * np.cos(theta_e))
print(" Electron: E' = %.2f MeV, theta = 90 deg" % E_prime_90)
print(" Proton: T = %.2f MeV, phi = %.1f deg (forward)" % (T_p_90, np.degrees(phi)))
print(" Proton speed: beta = %.4f c" % beta_func(T_p_90, m_p))
print(" -> Proton barely moves. It is %.0fx heavier than the electron." % (m_p/m_e))
print("\n--- Cross section (Mott, point proton) ---")
for theta_deg in [10, 30, 90]:
theta = np.radians(theta_deg)
sin2 = np.sin(theta/2)**2
sin4 = sin2**2
cos2 = np.cos(theta/2)**2
recoil = 1 / (1 + 2*E_e*sin2/m_p)
dsdo = (alpha * hbar_c / (2*E_e))**2 * cos2 / sin4 * recoil
barn = dsdo / 100
mbarn = barn * 1000
print(" theta = %3d deg: dsig/dOmega = %.4f fm^2/sr = %.4f mbarn/sr" % (theta_deg, dsdo, mbarn))
# ============================================================
# SCENARIO B: 4 GeV ELECTRON - RESONANCE REGION
# ============================================================
print("\n\n" + "=" * 75)
print("SCENARIO B: 4 GeV electron on stationary proton")
print("=" * 75)
T_e = 4000.0
E_e = T_e + m_e
p_e = np.sqrt(E_e**2 - m_e**2)
lam = de_broglie(p_e)
gamma_e = E_e / m_e
print("\n--- Incoming electron ---")
print(" Kinetic energy: T = %.1f GeV" % (T_e/1000))
print(" Momentum: p = %.4f GeV/c (ultra-relativistic: E ~ pc)" % (p_e/1000))
print(" Lorentz gamma: %.0f" % gamma_e)
print(" de Broglie lam: %.4f fm" % lam)
print(" lam / r_proton: %.4f -> RESOLVES proton interior\!" % (lam/r_p))
s = cm_energy_sq(E_e, m_e, m_p)
sqrt_s = np.sqrt(s)
W_available = sqrt_s - m_e
print("\n--- Centre-of-mass ---")
print(" sqrt(s) = %.3f GeV" % (sqrt_s/1000))
print(" Available hadronic mass: W = %.0f MeV" % W_available)
print("\n--- Production thresholds vs our sqrt(s) = %.0f MeV ---" % sqrt_s)
channels = [
("e + p -> e + p (elastic)", m_e + m_p),
("e + p -> e + p + pi0", m_e + m_p + m_pi0),
("e + p -> e + n + pi+", m_e + m_n + m_pip),
("e + p -> e + Delta(1232) -> e + p + pi", m_e + m_Delta),
("e + p -> e + p + pi+ + pi-", m_e + m_p + 2*m_pip),
("e + p -> e + p + pi0 + pi0", m_e + m_p + 2*m_pi0),
("e + p -> e + p + rho0(->pipi)", m_e + m_p + m_rho),
("e + p -> e + p + 3pi", m_e + m_p + 3*m_pip),
("e + p -> e + p + p + pbar", m_e + 3*m_p),
]
for label, threshold in channels:
status = "OPEN " if sqrt_s >= threshold else "closed"
print(" [%s] %s" % (status, label))
print(" threshold: sqrt(s) = %.0f MeV = %.3f GeV" % (threshold, threshold/1000))
print("\n--- Scattered electron at selected angles ---")
print(" For ELASTIC scattering (proton stays intact):")
print(" %8s %12s %12s %14s" % ("theta_lab", "Ee' (MeV)", "Q^2 (GeV^2)", "lam_probe (fm)"))
print(" " + "-" * 50)
for theta_deg in [6, 10, 20, 30, 45]:
theta = np.radians(theta_deg)
sin2 = np.sin(theta/2)**2
E_prime = E_e / (1 + 2*E_e*sin2/m_p)
Q2 = 4 * E_e * E_prime * sin2
Q2_GeV = Q2 / 1e6
lam_probe = 2 * np.pi * hbar_c / np.sqrt(Q2)
print(" %6d deg %12.1f %12.3f %14.3f" % (theta_deg, E_prime, Q2_GeV, lam_probe))
print("\n--- The Delta(1232) resonance ---")
print(" The first thing the electron 'breaks off' the proton.")
print(" p + gamma* -> Delta++, Delta+, Delta0, or Delta-")
print(" Delta(1232) has spin 3/2, isospin 3/2")
print(" Decays in ~5.6 x 10^-24 s (barely exists\!)")
print(" Delta+ -> p + pi0 (BR ~33%%)")
print(" -> n + pi+ (BR ~67%%)")
print(" Lifetime: tau = hbar/Gamma, where Gamma ~ 117 MeV")
tau_Delta = hbar_s / 117.0
print(" tau = hbar/Gamma = %.2e s" % tau_Delta)
print(" Travels: c*tau = %.2f fm before decaying" % (tau_Delta * 3e23))
print(" -> It decays INSIDE the proton\! You never see the Delta directly.")
print(" You see its decay products: a proton/neutron + pion(s).")
print("\n--- What the detector sees at 4 GeV ---")
print(" FORWARD (small angle):")
print(" - Scattered electron (most of the beam energy)")
print(" - Hadronic debris: proton or neutron + pions, mostly forward")
print(" LARGE ANGLE (theta > 30 deg):")
print(" - Hard-scattered electron (lost a LOT of energy)")
print(" - Multiple pions spraying forward")
print("")
print(" Typical event at theta = 10 deg, elastic:")
theta = np.radians(10)
sin2 = np.sin(theta/2)**2
E_prime = E_e / (1 + 2*E_e*sin2/m_p)
T_p = E_e - E_prime
p_p = np.sqrt(T_p**2 + 2*T_p*m_p)
p_ep = np.sqrt(E_prime**2 - m_e**2)
phi_p = np.degrees(np.arctan2(p_ep*np.sin(theta), p_e - p_ep*np.cos(theta)))
print(" Electron out: E' = %.0f MeV at theta = 10 deg" % E_prime)
print(" Proton out: T = %.0f MeV at phi = %.1f deg" % (T_p, phi_p))
print(" Proton speed: beta = %.4f c" % beta_func(T_p, m_p))
print("\n Typical inelastic event (Delta production) at theta = 10 deg:")
theta = np.radians(10)
sin2 = np.sin(theta/2)**2
E_prime_D = (2*m_p*E_e - (m_Delta**2 - m_p**2)) / (2*m_p + 4*E_e*sin2)
nu = E_e - E_prime_D
Q2_D = 4 * E_e * E_prime_D * sin2
print(" Electron out: E' = %.0f MeV at theta = 10 deg (lost %.0f MeV\!)" % (E_prime_D, nu))
print(" Q^2 = %.3f GeV^2" % (Q2_D/1e6))
print(" The %.0f MeV goes into creating the Delta(1232)" % nu)
print(" Delta -> proton (%.0f MeV) + pi0 (%.0f MeV)" % (m_p, m_pi0))
print(" or Delta -> neutron (%.0f MeV) + pi+ (%.0f MeV)" % (m_n, m_pip))
E_pi_cm = (m_Delta**2 - m_p**2 + m_pi0**2) / (2*m_Delta)
p_pi_cm = np.sqrt(E_pi_cm**2 - m_pi0**2)
T_pi_cm = E_pi_cm - m_pi0
E_p_cm = (m_Delta**2 + m_p**2 - m_pi0**2) / (2*m_Delta)
T_p_cm = E_p_cm - m_p
print(" In Delta rest frame: pion gets T = %.0f MeV, proton gets T = %.0f MeV" % (T_pi_cm, T_p_cm))
print(" Pion momentum (Delta frame): p = %.0f MeV/c" % p_pi_cm)
gamma_D = (m_p + nu) / m_Delta
beta_D = np.sqrt(1 - 1/gamma_D**2)
print(" Delta boost: gamma = %.2f, beta = %.4f" % (gamma_D, beta_D))
p_pi_lab_fwd = gamma_D * (p_pi_cm + beta_D * E_pi_cm)
print(" Forward pion in lab: p ~ %.0f MeV/c" % p_pi_lab_fwd)
print(" -> Pion flies FORWARD (same direction as beam)")
# ============================================================
# SCENARIO C: 20 GeV ELECTRON - DEEP INELASTIC
# ============================================================
print("\n\n" + "=" * 75)
print("SCENARIO C: 20 GeV electron on stationary proton (SLAC 1968)")
print("=" * 75)
T_e = 20000.0
E_e = T_e + m_e
p_e = np.sqrt(E_e**2 - m_e**2)
lam = de_broglie(p_e)
print("\n--- Incoming electron ---")
print(" Energy: E = %.3f GeV" % (E_e/1000))
print(" Momentum: p = %.3f GeV/c" % (p_e/1000))
print(" de Broglie lam: %.4f fm" % lam)
print(" lam / r_proton: %.4f" % (lam/r_p))
print(" -> Resolves to 1/%.0f of proton radius\!" % (r_p/lam))
print(" This electron sees INSIDE the proton.")
s = cm_energy_sq(E_e, m_e, m_p)
sqrt_s = np.sqrt(s)
print("\n--- Centre-of-mass ---")
print(" sqrt(s) = %.2f GeV" % (sqrt_s/1000))
print(" Enough energy for: ~%.0f pions" % ((sqrt_s - m_p - m_e)/m_pip))
print("\n--- Deep inelastic scattering kinematics ---")
print(" The electron hits a SINGLE QUARK inside the proton.")
print(" Key variables:")
print(" Q^2 = 4EE'sin^2(theta/2) -- resolution (momentum transfer squared)")
print(" nu = E - E' -- energy transferred to proton")
print(" x = Q^2/(2M*nu) -- Bjorken x (quark momentum fraction)")
print(" W^2 = M^2 + 2M*nu - Q^2 -- invariant mass of hadronic debris")
print(" y = nu/E -- inelasticity")
print("\n %5s %8s %9s %7s %6s %7s %5s %30s" % ("theta", "E'(GeV)", "Q2(GeV2)", "nu(GeV)", "x", "W(GeV)", "y", "What happens"))
print(" " + "-" * 85)
cases = [
(6, 18.0), (6, 15.0), (6, 10.0), (6, 5.0),
(10, 17.0), (10, 12.0), (10, 5.0),
(20, 15.0), (20, 10.0), (20, 5.0),
(34, 10.0), (34, 5.0),
]
for theta_deg, Ep_GeV in cases:
theta = np.radians(theta_deg)
sin2 = np.sin(theta/2)**2
Ep = Ep_GeV * 1000
E = E_e
Q2 = 4 * E * Ep * sin2 / 1e6
nu_val = (E - Ep) / 1000
x_val = Q2 / (2 * m_p/1000 * nu_val) if nu_val > 0 else 0
W2 = (m_p/1000)**2 + 2*(m_p/1000)*nu_val - Q2
W = np.sqrt(max(W2, 0))
y_val = nu_val / (E/1000)
if x_val > 0.95 and W < 1.0:
what = "elastic (proton intact)"
elif W < 1.5:
what = "quasi-elastic / Delta(1232)"
elif x_val > 0.5:
what = "valence quark hit"
elif x_val > 0.2:
what = "quark hit, multi-pion out"
elif x_val > 0.05:
what = "sea quark / soft, many pi"
else:
what = "very soft, huge hadron shower"
if x_val > 1.01 or x_val < 0:
continue
print(" %4d deg %7.1f %9.2f %7.1f %6.3f %7.2f %5.2f %30s" % (theta_deg, Ep_GeV, Q2, nu_val, x_val, W, y_val, what))
print("\n--- What the detector sees at 20 GeV ---")
print("\n Example event: theta = 6 deg, E' = 15 GeV")
theta = np.radians(6)
sin2 = np.sin(theta/2)**2
Ep = 15000.0
Q2 = 4 * E_e * Ep * sin2 / 1e6
nu_val = (E_e - Ep) / 1000
x_val = Q2 / (2 * m_p/1000 * nu_val)
W2 = (m_p/1000)**2 + 2*(m_p/1000)*nu_val - Q2
W = np.sqrt(W2)
lam_Q = 2*np.pi*hbar_c/1000/np.sqrt(Q2)
print(" Q^2 = %.2f GeV^2 -> probing at lam = %.3f fm" % (Q2, lam_Q))
print(" nu = %.1f GeV (energy dumped into proton)" % nu_val)
print(" x = %.3f (hit a quark carrying %.1f%% of proton momentum)" % (x_val, x_val*100))
print(" W = %.2f GeV (mass of hadronic debris)" % W)
print("")
print(" PARTICLES OUT:")
print(" 1. SCATTERED ELECTRON")
print(" E' = 15 GeV, theta = 6 deg, still ultra-relativistic")
print(" -> Detected in electromagnetic calorimeter")
print("")
print(" 2. HADRONIC JET (the broken proton)")
n_pions_approx = int((W*1000 - m_p) / m_pip)
print(" Total mass W = %.2f GeV" % W)
print(" Contains ~%d pions + remnant nucleon" % n_pions_approx)
print(" Typical particles:")
print(" - 1 proton or neutron (the 'spectator' quarks)")
print(" - %d-%d pions (pi+, pi-, pi0)" % (max(1,n_pions_approx-1), n_pions_approx+2))
print(" - Occasionally: K+, K-, K0 (if strange quarks produced)")
print(" - Rarely: eta, rho, omega mesons")
print(" All flying FORWARD (same direction as beam)")
boost_g = (m_p/1000 + nu_val)/W
print(" Boost: gamma ~ %.1f (the debris is highly boosted)" % boost_g)
print("\n--- Bjorken scaling: the smoking gun for quarks ---")
print(" If the proton were a smooth blob, the cross section would")
print(" depend on BOTH Q^2 and nu independently.")
print(" If the proton contains point-like quarks, the cross section")
print(" depends only on x = Q^2/(2M*nu). This is Bjorken scaling.")
print("")
print(" Same x ~ 0.33 at different Q^2:")
print(" %5s %8s %8s %6s %6s %8s" % ("theta", "E'(GeV)", "Q^2", "nu", "x", "Same x?"))
print(" " + "-" * 45)
for theta_deg, Ep_GeV in [(6, 17.3), (10, 14.5), (20, 8.6), (34, 4.7)]:
theta = np.radians(theta_deg)
sin2 = np.sin(theta/2)**2
Ep = Ep_GeV * 1000
Q2 = 4 * E_e * Ep * sin2 / 1e6
nu_val = (E_e - Ep) / 1000
x_val = Q2 / (2 * m_p/1000 * nu_val)
same = "~ yes" if abs(x_val-0.33)<0.05 else "no"
print(" %4d deg %7.1f %8.2f %6.1f %6.3f %8s" % (theta_deg, Ep_GeV, Q2, nu_val, x_val, same))
print(" -> If F2(x,Q^2) ~ F2(x) regardless of Q^2, quarks are POINT-LIKE.")
print(" This is exactly what SLAC found in 1968.")
print("\n--- Where is the proton's momentum? ---")
print(" Sum over all quarks: integral x*[f_u(x) + f_d(x) + ...] dx ~ 0.54")
print(" But the proton's total momentum is integral x*f(x) dx = 1.00")
print(" Missing 46%% -> carried by GLUONS (colourless, chargeless)")
print(" The electron can't scatter off gluons directly.")
print(" But the missing momentum PROVES they exist.")
# ============================================================
# SUMMARY TABLE
# ============================================================
print("\n\n" + "=" * 75)
print("SUMMARY: WHAT FLIES OUT AT EACH ENERGY")
print("=" * 75)
lam20 = 2*np.pi*hbar_c/20000
pixels = r_p/lam20
print("""
Energy What goes in What comes out
------------- ------------------- ----------------------------------------
13.6 eV e- (thermal) NOTHING flies out -- bound state\!
-> hydrogen atom, photon emission
100 MeV e- (v ~ c) e- (deflected) + p (barely moves)
Proton stays INTACT.
Like bouncing a marble off a bowling
ball -- you learn it's round, not
what's inside.
4 GeV e- (gamma ~ 8000) e- (lost energy) + p + pi0
or e- + n + pi+
Proton gets EXCITED (Delta resonance)
then decays. First cracks appear.
20 GeV e- (gamma ~ 40000) e- (hard scatter) + HADRON JET:
(SLAC) proton/neutron + 3-8 pions
occasionally kaons, etas
Proton SHATTERED. Electron hit a
single quark. Quarks rip apart ->
new quark-antiquark pairs form ->
pions (= quark + antiquark bound
states). You NEVER see a free quark.
""")
print(" KEY INSIGHT: At 20 GeV, the electron wavelength is %.4f fm." % lam20)
print(" The proton radius is %.4f fm." % r_p)
print(" The electron sees %.0f 'pixels' across the proton." % pixels)
print(" Each pixel contains: one quark (a point, as far as we can tell).")
print("")
print(" The quarks carry only 54%% of the proton's momentum.")
print(" The other 46%% is carried by gluons -- the non-associative algebra itself.")
print(" You can't scatter an electron off the algebra. But the missing momentum")
print(" proves it's there.")
# ============================================================
# ONE SPECIFIC EVENT
# ============================================================
print("\n" + "=" * 75)
print("ONE SPECIFIC EVENT: 20 GeV electron, theta = 10 deg, E' = 12 GeV")
print("=" * 75)
E_beam = 20000.0
theta_deg = 10
theta = np.radians(theta_deg)
sin2 = np.sin(theta/2)**2
E_prime = 12000.0
Q2 = 4 * E_beam * E_prime * sin2 / 1e6
nu_MeV = E_beam - E_prime
nu_GeV = nu_MeV / 1000
x_val = Q2 / (2 * m_p/1000 * nu_GeV)
W2 = (m_p/1000)**2 + 2*(m_p/1000)*nu_GeV - Q2
W = np.sqrt(W2)
print("\n BEFORE:")
print(" Electron: E = 20.000 GeV, p = 20.000 GeV/c, direction = ->")
print(" Proton: E = 0.938 GeV (at rest), p = 0")
print(" Total: E = %.3f GeV, p = %.3f GeV/c" % ((E_beam + m_p)/1000, E_beam/1000))
quark_type = "u (up)" if x_val > 0.25 else "d (down) or sea"
print("\n INTERACTION:")
print(" Virtual photon exchanged: Q^2 = %.2f GeV^2" % Q2)
print(" Energy transferred: nu = %.1f GeV" % nu_GeV)
print(" Bjorken x = %.3f -> hit a quark with %.1f%% of proton momentum" % (x_val, x_val*100))
print(" This is a %s quark (most likely)" % quark_type)
print("\n AFTER (what the detector records):")
print(" +-- SCATTERED ELECTRON")
print(" | Energy: E' = %.1f GeV" % (E_prime/1000))
print(" | Angle: theta = %d deg from beam axis" % theta_deg)
print(" | Momentum: p = %.1f GeV/c" % (E_prime/1000))
p_e_x = E_prime * np.cos(theta) / 1000
p_e_y = E_prime * np.sin(theta) / 1000
print(" | p_x = %.2f GeV/c, p_y = %.2f GeV/c" % (p_e_x, p_e_y))
print(" |")
p_had_x = (E_beam - E_prime*np.cos(theta)) / 1000
p_had_y = -E_prime*np.sin(theta) / 1000
E_had = (m_p + nu_MeV) / 1000
p_had = np.sqrt(p_had_x**2 + p_had_y**2)
theta_had = np.degrees(np.arctan2(abs(p_had_y), p_had_x))
n_pions = int((W*1000 - m_p) / m_pip)
n_avg = 2 + 1.5 * np.log(W**2) if W > 1 else 1
print(" +-- HADRONIC JET (the shattered proton)")
print(" Total mass: W = %.2f GeV" % W)
print(" Total energy: E = %.2f GeV" % E_had)
print(" Total momentum: p = %.2f GeV/c" % p_had)
print(" Direction: theta ~ %.1f deg from beam (mostly forward)" % theta_had)
print(" Avg multiplicity: ~%.0f hadrons" % n_avg)
print("")
print(" Typical particle list:")
print(" - 1 proton (E ~ 3-6 GeV, forward)")
print(" - 2-3 pi+ (E ~ 1-3 GeV each)")
print(" - 1-2 pi- (E ~ 0.5-2 GeV each)")
print(" - 1-2 pi0 (each -> 2 gamma, E_gamma ~ 0.3-1 GeV)")
print(" - Total visible: ~6-8 particles")
print("\n CONSERVATION CHECK:")
print(" Energy: %.3f + %.3f = %.3f GeV in" % (E_beam/1000, m_p/1000, (E_beam+m_p)/1000))
print(" %.1f + %.2f = %.3f GeV out [check]" % (E_prime/1000, E_had, E_prime/1000 + E_had))
print(" Momentum (x): %.3f in -> %.2f + %.2f = %.3f out [check]" % (E_beam/1000, p_e_x, p_had_x, p_e_x+p_had_x))
print(" Momentum (y): 0 in -> %.2f + %.2f = %.4f out [check]" % (p_e_y, p_had_y, p_e_y+p_had_y))
print(" Charge: -1 + 1 = 0 in -> -1 + (total = 0) out [check]")
print("\n THE PUNCHLINE:")
print(" The electron bounced off a POINT-LIKE object inside the proton.")
print(" That object carried x = %.3f = %.1f%% of the proton's momentum." % (x_val, x_val*100))
print(" The proton was destroyed in the process -- but NO free quarks emerged.")
print(" Instead, the energy created new quark-antiquark pairs from vacuum,")
print(" which instantly bound into pions and other hadrons.")
print(" Confinement is absolute: you put in %.0f GeV trying to free a quark," % nu_GeV)
print(" and all you get is more bound states.")
6The Missing Momentum: Gluons
If you add up the momentum carried by all the quarks the electron can scatter from, you get a surprise. The momentum sum rule says:
(14)∫₀¹ x [ u(x) + ū(x) + d(x) + d̅(x) + g(x) ] dx = 1
But when you measure the quark contributions alone: ∫₀¹ x [ u(x) + d(x) ] dx ≈ 0.54 → neutral constituents carry ∼46% (gluons ≈ 41%, strange sea ≈ 5%) (15)
The proton’s momentum budget:
7Confinement: Why You Cannot Free a Quark
The most striking fact about the 20 GeV experiment is what does not come out: a free quark. The electron smashes a single quark with 8 GeV of energy transfer. The quark flies away. But it never emerges alone. The quark-quark potential has two terms. At short range, a weak attractive well (from gluon exchange). At long range, a linearly rising string tension:
(12)V(r) = −
The string tension κ ≈ 0.18 GeV² ≈ 0.9 GeV/fm means the force between separated quarks is approximately constant at about 15 tonnes — independent of distance. As the struck quark flies away, the energy stored in the stretched colour field grows at 0.9 GeV per femtometre. When the field energy exceeds the rest mass energy of a quark-antiquark pair (about 300 MeV for up/down quarks), the field snaps and a new pair is created from the vacuum. The new quark joins the departing quark to form a meson (pion). The process repeats. You put in 8 GeV trying to free a quark and all you get is more bound states. The distance at which the first string breaks:
rₛₜʳʳʳ = 2mᵠ / κ ≈ 300 MeV / (0.18 GeV²) = 1.67 GeV⁻¹ = 0.33 fm
After 0.33 fm of separation, a new quark-antiquark pair is created. How many pions does the snapping string produce? Hadron multiplicity grows logarithmically with W; empirically, W = 2.93 GeV yields 5–6 hadrons on average.
Confinement from the octonion algebra Observable states must be associative: (AB)C = A(BC) Single quark: carries colour e₅, e₆, or e₇ — non-associative — unobservable Meson (q q̅): colour × anticolour — associative — observable Baryon (qqq): R×
8Bjorken Scaling: Proof That Quarks Are Points
If the proton were a continuous blob of charge, its structure function F₂(x, Q²) would change as you changed Q²: zooming in would reveal finer and finer texture. But SLAC found something different: F₂(x, Q²) → F₂(x) [Bjorken scaling: independent of Q²] (13) F₂ depends only on the ratio x = Q²/(2Mν), not on Q² alone. Whether you probe at Q² = 1 GeV² or Q² = 10 GeV², you see the same pattern. This is Bjorken scaling, and it means the objects inside the proton are point-like — they have no internal structure of their own. The scaling data from the original SLAC experiment (approximate): The ratio of maximum to minimum F₂ at any fixed x is less than 10% across a factor-of-five range in Q². This is scaling to within experimental precision. The objects inside the proton have no size. They are points. They are quarks.
9Summary: One System, Four Answers
The transition energy
The boundary between quaternion and octonion physics occurs where λ = rₚ:
pₜʳʳʳₛ = 2πħ
Below this momentum, the proton’s second quaternion is invisible. The proton looks like a point with a charge. Above this momentum, the second quaternion is revealed: confinement energy in e₄, colour charges in e₅, e₆, e₇, gluon couplings in the non-associative products. SLAC did not discover “quarks.” SLAC discovered the proton’s second quaternion.
Equation Index
(1) λ = 2πħ
(2) dσ/dΩ = (
(3) √s = √( (E + Mₚ)² − pᵉ² ) = √( Mₚ² + 2MₚE ) [target at rest]
(4) √sₜℍ = Mₚ + Mπ = 938.3 + 135.0 = 1073.3 MeV
(5) E' = E / ( 1 + (2E/Mₚ) sin²(θ/2) ) [elastic, proton at rest]
(6) Q² = 4EE′ sin²(θ/2)
(7) ν = E − E′
(8) x = Q² / (2Mₚν)
(9) W² = Mₚ² + 2Mₚν − Q²
(10) Tₚ⁻⁺ˣ = 2Mₚpᵉ² / ( Mₚ² + 2MₚE + mᵉ² ) ≈ 2Mₚpᵉ² / Mₚ² for E ≫ Mₚ
(11) ΔE = MΔ − Mₚ = 1232 − 938.3 = 293.7 MeV
(12) V(r) = −
(13) F₂(x, Q²) → F₂(x) [Bjorken scaling: independent of Q²]
(14) ∫₀¹ x [ u(x) + ū(x) + d(x) + d̅(x) + g(x) ] dx = 1
(15) ∫₀¹ x [ u(x) + d(x) ] dx ≈ 0.54 → neutral constituents carry ∼46% (gluons ≈ 41%, strange sea ≈ 5%)