The Stiffness Is the Curvature Die Steifigkeit ist die Krümmung
The Skyrme Term as the Curvature of the Field's Target Sphere, and the Coefficient as the Door to the Next Mode
The Skyrme Term as the Curvature of the Field's Target Sphere, and the Coefficient as the Door to the Next Mode
Abstract
The particle paper of this series ended at one fitted constant: Skyrme's stiffness e = 5.45, the term that catches Derrick's collapse and sets the size of every knot, postulated there as it was postulated in 1961. This note removes the postulate and prices the constant. First, an identity of pure geometry, walked symbol by symbol: the Skyrme term is the square of the pulled-back curvature two-form of the target sphere. For a sphere of constant curvature, the curvature two-form is vielbein wedge vielbein, R^ab = e^a∧e^b; its pullback by the winding field is the commutator of the currents; its square is Skyrme's term, exactly. One boundary must be drawn at the outset, because it is the paper's most abusable line and an earlier draft crossed it: this is the curvature of the field's target sphere S³ — the internal manifold the winding field U maps into — computed by pullback, exactly as verified in frame_quaternion_verification.py. It is a genuine second geometric invariant of the winding, so that the two invariants of a wound field are its pulled-back target-metric (tension) and its pulled-back target-curvature (stiffness). But it is not automatically the foundation theory's spacetime frame curvature R = dω + ω∧ω (the curvature of the gravitational connection): those are two different manifolds' curvatures, and identifying them is a further claim this paper does not make. What is established is narrower and clean: Skyrme's postulated quartic term is the squared pulled-back target curvature — a standard differential-geometry identity — so within the effective field it is not an alien addition but the natural second invariant. The debt thereby shrinks from a term to a number — the ratio of the curvature-energy price to the tension — and for the one fiber whose cellar physics can enter (the pion field of QCD), that number has a known and verifiable address: the first excited mode of the tower, the ρ meson. Hidden-local-symmetry theory generates the Skyrme term by integrating out the ρ, predicting e ≈ g_ρππ; this note verifies the prediction against data: g = 5.976 from the measured ρ width, g = 5.952 from the KSRF mass relation, against the fitted e = 5.45 — agreement at 9 percent, with the stiffness length unit 1/(eF_π) = 0.281 fm landing within 10 percent of the ρ's Compton radius 0.255 fm. The pattern this establishes — the stiffness of a rung is the door to its next mode; the amputated interior returns as a tower whose ground floor sets the coefficient — is then set beside — not extended to — the weak fiber, and only as a flagged, non-viable curiosity: the lowest computed weak tone, (3/2)m_W = 120.6 GeV, lies 3.7 percent from the Higgs, but that tone is a Dirac (spinorial) eigenmode while the Higgs is a scalar, and a spin-½ mode cannot be a spin-0 particle. That categorical mismatch blocks the identification, so the pattern has exactly one demonstrated instance — the pion fiber and its ρ (a genuine scalar-isovector state) — not two; the weak comparison is a number without a mechanism until a scalar mode of the wound weak S³ is exhibited. Priority throughout belongs to the literature — Skyrme; Bando, Kugo and Yamawaki; the KSRF relation — and is stated loudly; the contribution claimed is the identification of the pair, the translation to
frame_quaternion_verification.py — frame_quaternion_verification
A. HAMILTON'S TABLE (the base algebra P1 is built on)
i^2 = -1 OK
j^2 = -1 OK
k^2 = -1 OK
ij = k OK
jk = i OK
ki = j OK
ji = -k OK
B. |U| = 1 EVERYWHERE (frame is unit-norm: metric = theta.thetabar)
max | |U|^2 - 1 | over sample points = 2.22e-16 OK
C. CURRENTS L_mu = U^dag d_mu U ARE ALGEBRA-VALUED (real part ~ 0)
max |Re tr L_mu / 2| = 1.68e-10 OK
D. MAURER-CARTAN FLATNESS d_mu L_nu - d_nu L_mu + [L_mu,L_nu] = 0
(the composite connection is flat; curvature lives on the target)
max |curvature of composite connection| = 1.85e-09 OK (flat)
E. CORE IDENTITY (stiffness note): Skyrme term = |pullback of R^{ab}|^2
verifying |[L_mu,L_nu]|^2 proportional to sum(L^a_mu L^b_nu - L^a_nu L^b_mu)^2
|[L,L]|^2 total = 3.836709e+02
pullback-curvature total = 9.591773e+01
ratio = 4.000000 (constant => the two are the SAME object up to the
fixed representation factor; the Skyrme term is curvature^2)
=> tension energy = |theta|^2, stiffness energy = |R|^2 : the pair (theta,R).
# -*- coding: utf-8 -*-
"""Verification for the Frame Quaternion Field Equations note (It Is All One).
Everything the series stands on lives in this one note: the frame theta (the
vielbein), the curvature R = d(omega) + omega^omega, and the template
"curvature of the winding = source". This script confirms, by direct numerical
computation on a real hedgehog field U(x) -- not by algebra on paper -- the
identities that the particle paper and the stiffness note quietly rely on:
(A) Hamilton's table: i^2=j^2=k^2=-1, ij=k, jk=i, ki=j (the base algebra)
(B) U is a unit quaternion everywhere: |U| = 1 (metric = theta.thetabar)
(C) the currents L_mu = U^dagger d_mu U are pure-imaginary (algebra-valued frame)
(D) Maurer-Cartan flatness: d_mu L_nu - d_nu L_mu + [L_mu, L_nu] = 0
(the COMPOSITE connection has zero curvature -- the subtlety the
stiffness note flags; the nonzero curvature is the target sphere's)
(E) THE CORE IDENTITY of the stiffness paper, in components:
|[L_mu, L_nu]|^2 = |pullback of R^{ab}=e^a^e^b|^2
i.e. the Skyrme term IS the squared pulled-back curvature two-form.
Pure standard library: quaternions as 2x2 complex matrices (Pauli rep).
"""
import cmath, math
# --- quaternions as 2x2 complex matrices: 1, i, j, k ----------------------
def mul(A, B):
return [[A[0][0]*B[0][0]+A[0][1]*B[1][0], A[0][0]*B[0][1]+A[0][1]*B[1][1]],
[A[1][0]*B[0][0]+A[1][1]*B[1][0], A[1][0]*B[0][1]+A[1][1]*B[1][1]]]
def add(A, B): return [[A[r][c]+B[r][c] for c in range(2)] for r in range(2)]
def sub(A, B): return [[A[r][c]-B[r][c] for c in range(2)] for r in range(2)]
def scal(s, A): return [[s*A[r][c] for c in range(2)] for r in range(2)]
def dag(A): # conjugate transpose
return [[A[0][0].conjugate(), A[1][0].conjugate()],
[A[0][1].conjugate(), A[1][1].conjugate()]]
def fnorm2(A): return sum(abs(A[r][c])**2 for r in range(2) for c in range(2))
def comm(A, B): return sub(mul(A, B), mul(B, A))
I2 = [[1+0j,0],[0,1+0j]]
# Pauli-based quaternion rep: i=-i*sx, j=-i*sy, k=-i*sz (so i*j=k etc.)
qi = [[0, -1j],[-1j, 0]]
qj = [[0, -1],[1, 0]]
qk = [[-1j, 0],[0, 1j]]
def quat(w, x, y, z):
return add(add(scal(w, I2), scal(x, qi)), add(scal(y, qj), scal(z, qk)))
# --- (A) Hamilton's table -------------------------------------------------
def close(A, B, tol=1e-12): return fnorm2(sub(A, B)) < tol
print("A. HAMILTON'S TABLE (the base algebra P1 is built on)")
checks = [("i^2 = -1", close(mul(qi,qi), scal(-1,I2))),
("j^2 = -1", close(mul(qj,qj), scal(-1,I2))),
("k^2 = -1", close(mul(qk,qk), scal(-1,I2))),
("ij = k ", close(mul(qi,qj), qk)),
("jk = i ", close(mul(qj,qk), qi)),
("ki = j ", close(mul(qk,qi), qj)),
("ji = -k ", close(mul(qj,qi), scal(-1,qk)))]
for name, ok in checks:
print(f" {name} {'OK' if ok else 'FAIL'}")
# --- the hedgehog field U(x) = cos f + sin f (xhat . (i,j,k)) -------------
def profile(r): # any monotone f: pi at 0 -> 0 at infinity
return math.pi * math.exp(-r)
def U_of(x, y, z):
r = math.sqrt(x*x + y*y + z*z)
f = profile(r)
if r < 1e-9:
return quat(math.cos(f), 0, 0, 0)
s = math.sin(f) / r
return quat(math.cos(f), s*x, s*y, s*z)
# --- (B) unit quaternion everywhere --------------------------------------
print("\nB. |U| = 1 EVERYWHERE (frame is unit-norm: metric = theta.thetabar)")
maxdev = 0.0
for (x,y,z) in [(0.3,0.1,-0.2),(1.0,0.0,0.0),(0.4,0.7,0.5),(-0.6,0.2,0.9)]:
U = U_of(x,y,z)
n = fnorm2(U)/2 # trace-normalised: |U|^2 for unit quaternion = 1
maxdev = max(maxdev, abs(n-1))
print(f" max | |U|^2 - 1 | over sample points = {maxdev:.2e} "
f"{'OK' if maxdev<1e-9 else 'FAIL'}")
# --- currents by central finite difference at a test point ----------------
P = (0.35, 0.15, -0.25); hstep = 1e-5
def Ushift(i, s):
p = list(P); p[i] += s; return U_of(*p)
def dU(i):
return scal(1/(2*hstep), sub(Ushift(i, hstep), Ushift(i, -hstep)))
Uc = U_of(*P); Ud = dag(Uc)
L = [mul(Ud, dU(i)) for i in range(3)] # L_x, L_y, L_z
# --- (C) currents are pure-imaginary (traceless-antihermitian) -------------
print("\nC. CURRENTS L_mu = U^dag d_mu U ARE ALGEBRA-VALUED (real part ~ 0)")
realpart = max(abs((L[i][0][0]+L[i][1][1]).real)/2 for i in range(3))
print(f" max |Re tr L_mu / 2| = {realpart:.2e} {'OK' if realpart<1e-6 else 'FAIL'}")
# --- (D) Maurer-Cartan flatness ------------------------------------------
# dL_nu/dx_mu via finite difference of the current field
def L_at(p, i):
U0 = U_of(*p); return mul(dag(U0), scal(1/(2*hstep),
sub(U_of(*[p[k]+ (hstep if k==i else 0) for k in range(3)]),
U_of(*[p[k]- (hstep if k==i else 0) for k in range(3)]))))
def dL(mu, nu):
p_plus = [P[k]+(hstep if k==mu else 0) for k in range(3)]
p_minus = [P[k]-(hstep if k==mu else 0) for k in range(3)]
return scal(1/(2*hstep), sub(L_at(p_plus, nu), L_at(p_minus, nu)))
print("\nD. MAURER-CARTAN FLATNESS d_mu L_nu - d_nu L_mu + [L_mu,L_nu] = 0")
print(" (the composite connection is flat; curvature lives on the target)")
mc_max = 0.0
for (mu,nu) in [(0,1),(0,2),(1,2)]:
F = add(sub(dL(mu,nu), dL(nu,mu)), comm(L[mu], L[nu]))
mc_max = max(mc_max, math.sqrt(fnorm2(F)))
print(f" max |curvature of composite connection| = {mc_max:.2e} "
f"{'OK (flat)' if mc_max<1e-3 else 'FAIL'}")
# --- (E) the core identity: Skyrme term = squared pulled-back curvature ----
# component form: L_mu = L^a_mu T_a. Extract L^a_mu = -1/2 tr(T_a^dag L_mu)...
# simpler: |[L_mu,L_nu]|^2 vs sum_ab (L^a_mu L^b_nu - L^a_nu L^b_mu)^2
def comps(Lmu): # (L^1,L^2,L^3) on (qi,qj,qk): L^a = -1/2 tr(q_a * L) (since q_a^dag=-q_a)
return [ -0.5*(mul(qa, Lmu)[0][0]+mul(qa, Lmu)[1][1]) for qa in (qi,qj,qk) ]
print("\nE. CORE IDENTITY (stiffness note): Skyrme term = |pullback of R^{ab}|^2")
print(" verifying |[L_mu,L_nu]|^2 proportional to sum(L^a_mu L^b_nu - L^a_nu L^b_mu)^2")
lhs_tot = rhs_tot = 0.0
for (mu,nu) in [(0,1),(0,2),(1,2)]:
lhs = fnorm2(comm(L[mu], L[nu]))
cm, cn = comps(L[mu]), comps(L[nu])
rhs = sum(abs(cm[a]*cn[b]-cm[b]*cn[a])**2 for a in range(3) for b in range(3))
lhs_tot += lhs; rhs_tot += rhs
ratio = lhs_tot/rhs_tot
print(f" |[L,L]|^2 total = {lhs_tot:.6e}")
print(f" pullback-curvature total = {rhs_tot:.6e}")
print(f" ratio = {ratio:.6f} (constant => the two are the SAME object up to the")
print(f" fixed representation factor; the Skyrme term is curvature^2)")
print("\n => tension energy = |theta|^2, stiffness energy = |R|^2 : the pair (theta,R).")
1The Question Left Standing
The particle paper built a baryon from the space quaternion and closed on one honest embarrassment: the quartic term that stops Derrick's collapse was postulated, and its dimensionless coefficient e = 5.45 fitted. Every scale of the knot — its size, its mass, the whole calibration — hangs on that one number. Two readings were possible. Either the term is an alien addition, bolted on because Derrick demands something, in which case the framework has merely relocated its ignorance; or the term is a natural invariant the formalism already owns, in which case only the coefficient is owed. This note establishes the second reading by identity, then prices the coefficient by measurement.
2The Two Objects of the Frame Note
Recall the two objects the field-equations note placed at the base of the series. The frame quaternion θ: at each point, the local clock and rulers — the vielbein, whose norm is the metric, ds² = θ·θ̄. And the curvature R = dω + ω∧ω: the rotation a frame silently accumulates when carried around a small loop. The note's template runs 'curvature of the winding equals source.' Add to this Postulate P2's anatomy of a world: every world of
3The Identity, Walked Slowly
Let U(x) be the winding field of the particle paper — at every point of space, a unit quaternion, a point on the three-sphere S³. Every symbol in what follows is built from U.
31 The currents are pulled-back vielbeins
Define the currents L_μ = U†∂_μU. They are valued in the algebra (pure quaternion-imaginary), with components L^a_μ on the units: L_μ = L^a_μ T_a, where T_a are the quaternion units in their 2×2 representation and a runs over 1, 2, 3. Geometrically the L^a_μ are the vielbeins of the target sphere, pulled back to space by the map U: they measure how far, per unit step in direction μ, the field moves along the target's a-th axis. The gradient energy of the winding is their square:
E₂ = (
— the pulled-back metric of the sphere, weighted by the tension
32 The sphere's curvature is vielbein wedge vielbein
Now P2, taken literally. For a sphere of constant curvature (unit radius), the curvature two-form is not an independent field; it is built from the vielbeins themselves:
R^{ab} = e^a ∧ e^b
This single line is the geometric content of 'a sphere has one curvature everywhere': the rotation accumulated around a small loop is proportional to the loop's area, with a universal constant of proportionality — curvature equals area form, componentwise e^a∧e^b. It is the defining property of the constant-curvature worlds that P2 assigns to every rung of
33 Pull back, and square
Pull the target's curvature two-form back to space along U — substitute the pulled-back vielbeins:
(U*R^{ab})_{μν} = L^a_μ L^b_ν − L^a_ν L^b_μ
And compare with the commutator of the currents. Since [T_a, T_b] = ε_{abc}T_c (Hamilton's table),
[L_μ, L_ν] = ε_{abc} L^a_μ L^b_ν T_c ,
whose square, traced, is componentwise identical (up to a fixed numerical factor) to the square of the pulled-back curvature:
Tr( [L_μ, L_ν][L^μ, L^ν] ) ∝ Σ_{ab} ( L^a_μL^b_ν − L^a_νL^b_μ )² = |U*R|²
The left side is Skyrme's term. The right side is the curvature energy of the wound sphere. They are the same object:
E = (tension) · |pullback θ|² + (stiffness) · |pullback R|²
The pair (tension, stiffness) is the pair (target metric, target curvature) of the field's S³ — the pulled-back |θ_target|² and |R_target|². Skyrme, in 1961, without the language, wrote down the metric energy and the curvature energy of a wound constant-curvature target sphere. The debt shrinks accordingly: not 'why is there a quartic term?' — the target geometry supplies exactly this one, as its second invariant — but 'what is the ratio of its coefficient to the tension?' One number. (What this does not do: it does not show the stiffness is the foundation theory's spacetime curvature R = dω + ω∧ω. The suggestive echo — a wound sphere carries a metric energy and a curvature energy in both the target and the spacetime settings — is a formal parallel, not an identity; asserting the identity would require deriving the pion field's target S³ from the spacetime frame, which this paper does not do.)
34 The Maurer–Cartan subtlety, stated so the note does not overclaim
One piece of craftsmanship must be recorded. The composite connection built from the full currents is pure gauge, and its total curvature vanishes identically — the Maurer–Cartan identity dL + L∧L = 0. The nonvanishing object in Section 3.3 is the target's curvature pulled back — equivalently, in the symmetric-space (coset) decomposition of the sphere, the curvature of the coset part, for which R = −[e, e] is exactly the constant-curvature statement used above. The distinction is standard differential geometry of symmetric spaces; it changes no result, but a referee will ask, and the answer is written here first.
4What Remains: One Number
With the form fixed by geometry, dimensional analysis isolates the residue. The tension carries
5The Number's Address: the One Cellar Physics Can Enter
For the pion field — the chiral field of QCD, of which the Skyrme model is the effective facade — the cellar is known, and the coefficient has been derived rather than fitted. The mechanism (hidden local symmetry; Bando, Kugo and Yamawaki) is precisely the storeys picture of the companion architecture note: the Skyrme term is not fundamental; it is generated when the tower above the pion is integrated out, and the tower's ground floor — the ρ meson at 775 MeV, the first excited mode of the same flavor quantum numbers — dominates the sum. The prediction: e ≈ g_ρππ, the ρ's coupling, with the KSRF relation m_ρ² = 2g²f_π² tying the coupling to the masses. This note verifies both against data: Two independent roads — the decay width and the mass relation — give the same g ≈ 6.0; the fitted stiffness sits 9 percent below it, inside the Skyrme model's honest 10–20 percent band; and the stiffness length is, to 10 percent, the ρ's own Compton cell. The sentence this buys, and it is the sentence of the note: The stiffness of a rung is the door to its next mode. The Rattenschwanz that the seal amputates from the facade does not vanish without trace — it returns as a tower, and the tower's ground floor sets the stiffness coefficient of the storey below it. The cellar greets the facade with the address of its first floor. Figure 1. Left: the pion fiber's stiffness number, four ways — two measurements of the first excited mode's coupling (green) against the fitted Skyrme coefficient (blue), agreeing inside the model's honest ±20% band. This is the one demonstrated instance. Right: the weak-fiber comparison — the lowest computed tone against the measured Higgs, 3.7 percent apart — shown only to be flagged NON-VIABLE: the tone is a Dirac (spin-½) eigenmode, the Higgs a scalar (spin-0); a categorical mismatch, not a bookkeeping detail.
6The Pattern, Translated Up the Ladder
If the pattern holds — stiffness coefficient of fiber n = coupling to its first excited mode — it should hold on every rung. Translate it to the weak fiber. Its harmonic tower was computed in the discrete-orbits note: the lowest Dirac tone of the wound S³ of radius
(3/2)·m_W = 120.6 GeV,
and the measured Higgs sits at 125.25 GeV — 3.7 percent away, the very curiosity that note filed under 'derive it or bury it.' It is tempting to read this as the Higgs as the ρ of the weak winding, and there is a fit of function to admire: stiffening the W's winding to its 2.5-attometre range is the Higgs's job in the standard model. But this paper must state plainly that the comparison is not yet viable, and not merely 'unresolved'. The computed weak-fiber tone is a Dirac (spinorial) eigenmode; the Higgs is a scalar. A spin-½ mode cannot be a spin-0 particle — this is a categorical mismatch, not a bookkeeping factor to be sorted out later. Until a genuinely scalar mode of the wound weak S³ is exhibited near that energy, the 3.7-percent coincidence is a number without a mechanism, and it does not constitute a second instance of the pattern. So, honestly: the stiffness-as-first-mode pattern has one demonstrated instance — the pion fiber, at 9 percent, with the ρ a genuine scalar-isovector state of the right quantum numbers. The weak 'instance' is recorded only as a flagged numerical curiosity, blocked by the spin mismatch above, and Eddington patrols every corridor where pure numbers meet measured masses.
7Falsifiers and the Program
The pattern is exposed on three sides, and says so. (i) QCD side: if improved Skyrme fits (pion mass, quantum corrections) converge on a stiffness incompatible with g_ρ beyond the model's stated band, the anchor instance fails and the translation loses its ground. (ii) Weak side: the pattern requires the spin bookkeeping of the (3/2)m_W tone to be resolvable in the fiber's favor; a demonstration that no scalar mode of the wound S³ sits near that energy kills the second instance. (iii) Lepton side — the open case: the electron's Q-ball needs its own tower; the pattern predicts its stiffness is set by the first excited lepton mode — and the muon, at 206.8 electron masses, is the obvious candidate tenant, though no calculation of a lepton-fiber tower yet exists in this series. Constructing it is the named next step, and it would convert the mass-spectrum null result of the condensation arc — masses off the
8What Is Claimed, and What Is Not
9The Sentence
The stiffness is the curvature of the winding's target sphere, squared — a clean geometric identity within the effective field, though not (this paper is careful to say) the foundation's spacetime curvature; and its coefficient is the door-price of the first mode above it — measured in the pion fiber at nine percent, with the ρ a genuine scalar-isovector state — so that at least the strong rung's rented pillar is revealed to be the next rung of the same ladder, while the weak-fiber echo remains a flagged coincidence, blocked by a spin mismatch, waiting for a scalar mode to earn or lose it.
References
T. H. R. Skyrme, Proc. Roy. Soc. A 260, 127 (1961); K. Kawarabayashi and M. Suzuki, Phys. Rev. Lett. 16, 255 (1966); Riazuddin and Fayyazuddin, Phys. Rev. 147, 1071 (1966) — the KSRF relation; M. Bando, T. Kugo and K. Yamawaki, Phys. Rep. 164, 217 (1988) — hidden local symmetry and the ρ-generated Skyrme term; G. S. Adkins, C. R. Nauenberg and E. Witten, Nucl. Phys. B 228, 552 (1983); G. H. Derrick, J. Math. Phys. 5, 1252 (1964); Particle Data Group values of m_ρ, Γ_ρ, m_π, f_π, m_W, m_H; and the papers and notes of this series (the Frame Quaternion note — θ, R, and the template; the Postulates — P2, the anatomy of a world; From Space Quaternion to Particle — Derrick's wall and the fitted e; Why the Orbits Are Discrete — the weak fiber's towers and the (3/2)m_W curiosity; The Three Storeys — the seal and the returning Rattenschwanz). Verification script: stiffness_check.py. (Citations from memory; the literature-verification pass —
10Verification
The companion scripts, with their recorded output. Each script's docstring states what it establishes and what it does not; the Source tab shows the file itself, unedited.
stiffness_check.py — stiffness_check
pion momentum in rho decay p = 361.6 MeV g from measured rho width = 5.976 g from KSRF m_rho/(sqrt2 f_pi) = 5.952 (the two roads agree to 0.4 percent) fitted Skyrme stiffness e = 5.450 (ANW 1983, massless pion) ratio e / g_KSRF = 0.916 (8 percent below) stiffness length 1/(e F_pi) = 0.281 fm (ANW calibration F_pi = 129 MeV) rho Compton radius hbar/m_rho = 0.255 fm (ratio 1.10) (3/2) m_W = 120.57 GeV vs m_H = 125.25 GeV (3.7 percent apart) -- quarantined
# -*- coding: utf-8 -*-
"""Verification for 'The Stiffness Is the Curvature' (July 2026).
Checks, with nothing but the standard library:
1. g_rho-pi-pi extracted from the measured rho width
2. g from the KSRF relation m_rho^2 = 2 g^2 f_pi^2
3. comparison with the fitted Skyrme stiffness e (ANW 1983)
4. the stiffness length unit 1/(e F_pi) vs the rho's Compton radius
5. the weak-fiber instance: (3/2) m_W vs m_H
"""
import math
# --- inputs (PDG values, MeV unless noted) --------------------------------
m_rho = 775.26 # rho(770) mass
Gamma = 149.1 # rho full width
m_pi = 139.57 # charged pion mass
f_pi = 92.1 # pion decay constant (F_pi = 2 f_pi = 184.2 in Skyrme units)
e_ANW = 5.45 # fitted Skyrme parameter, Adkins-Nauenberg-Witten 1983
hbarc = 197.327 # MeV fm
m_W = 80.377e3 # W mass, MeV
m_H = 125.25e3 # Higgs mass, MeV
# --- 1. g from the width --------------------------------------------------
# Gamma(rho -> pi pi) = g^2 p^3 / (6 pi m_rho^2), p = pion momentum in rho frame
p = math.sqrt(m_rho**2 / 4 - m_pi**2)
g_width = math.sqrt(6 * math.pi * m_rho**2 * Gamma / p**3)
print(f"pion momentum in rho decay p = {p:8.1f} MeV")
print(f"g from measured rho width = {g_width:8.3f}")
# --- 2. g from KSRF -------------------------------------------------------
g_ksrf = m_rho / (math.sqrt(2) * f_pi)
print(f"g from KSRF m_rho/(sqrt2 f_pi) = {g_ksrf:6.3f}")
print(f" (the two roads agree to {100*abs(g_width-g_ksrf)/g_ksrf:.1f} percent)")
# --- 3. against the fitted stiffness -------------------------------------
print(f"fitted Skyrme stiffness e = {e_ANW:8.3f} (ANW 1983, massless pion)")
print(f"ratio e / g_KSRF = {e_ANW/g_ksrf:8.3f} ({100*(1-e_ANW/g_ksrf):.0f} percent below)")
# --- 4. the stiffness length ---------------------------------------------
# Use the ANW fit's own calibration F_pi = 129 MeV (the particle paper's units),
# since e = 5.45 was fitted jointly with it.
F_pi_ANW = 129.0
L_stiff = hbarc / (e_ANW * F_pi_ANW) # fm
L_rho = hbarc / m_rho # fm
print(f"stiffness length 1/(e F_pi) = {L_stiff:8.3f} fm (ANW calibration F_pi = 129 MeV)")
print(f"rho Compton radius hbar/m_rho = {L_rho:8.3f} fm (ratio {L_stiff/L_rho:.2f})")
# --- 5. the weak-fiber instance (quarantined) ----------------------------
tone = 1.5 * m_W
print(f"(3/2) m_W = {tone/1000:7.2f} GeV vs m_H = {m_H/1000:7.2f} GeV"
f" ({100*abs(m_H-tone)/m_H:.1f} percent apart) -- quarantined")