subhorizon_resolution.py โ€” src/35-hadrons

The file src/35-hadrons/subhorizon_resolution.py of the validation package of the paper 35-hadr, rendered with an anchor per line; the file itself is beside this page.

35-hadrons ledgersrc/35-hadrons/subhorizon_resolution.py ยท 4.5 KBsubhorizon_resolution.py
1#!/usr/bin/env python32# Terminal resolution of E1 and E4 by the Omega-stability discriminator.3#4# FRC binary: every predicate resolves to a sub-horizon small residue (T) or Omega-hard bedrock5# (Omega). The discriminator (as in the 28-flavour delta_0 push) is Omega-stability: a quantity6# is sub-horizon iff it is the same on every admissible shell Omega = 5 (mod 12), i.e. it does not7# depend on the carrier factorisation of Omega-1 / the cross-scale running. Framed-rational; exact.89from fractions import Fraction as F1011shells = [53, 173, 389, 1373]1213print("== the discriminator: are the E1/E4 building blocks Omega-stable? ==")14print("  Every dimensionless baryon quantity is built from three ingredients:")15print("   (i)   the colour rank N = 3;")16print("   (ii)  the gauge character sums c_1(beta), c_adj(beta);")17print("   (iii) the framed flavour ratio m_s/m_l (the chi_3 current-mass ratio, 28-flavour).")18print("  None may carry Omega for the quantity to be sub-horizon.\n")1920print("  (i) colour rank fixed on every shell (triality 3 | Omega+1, for Omega = 5 mod 12):")21for Om in shells:22    assert Om % 12 == 5 and (Om + 1) % 3 == 023    N = 324    print(f"     Omega = {Om:5d}:  3 | Omega+1 = {Om+1:5d}  ->  colour rank N = {N}  (Omega-independent)")25print("     -> N = 3 on every admissible shell: the gauge group SU(3,F_q) is Omega-independent.\n")2627print("  (ii) the character sums depend only on N, not on Omega:")28def c1_lead(N):  return F(1, 2*N*N)        # fundamental plaquette: beta/(2N^2)29def cadj_lead(N):return F(1, (2*N)**2)     # adjoint plaquette:    beta^2/(2N)^230print(f"     c_1   leading = 1/(2N^2) beta   = {c1_lead(3)} beta   (= beta/18 at N=3)")31print(f"     c_adj leading = 1/(2N)^2 beta^2 = {cadj_lead(3)} beta^2 (= beta^2/36 at N=3)")32print("     these are finite-group character sums (D6b: no Haar, no Omega->infinity): Omega-independent.\n")3334print("  (iii) m_s/m_l is a framed (chi_3) current-mass ratio of 28-flavour: sub-horizon, Omega-stable.")35print("  => every dimensionless baryon quantity is Omega-stable -> SUB-HORIZON (T).")36print("     Omega enters ONLY the absolute scale sqrt(sigma)=Lambda_QCD (dimensional transmutation,")37print("     the cross-scale running / carrier factorisation): the lone Omega-hard residue (E5/D6c), with alpha (E7).")3839# ---------------------------------------------------------------------------------------------40print("\n== E4 resolution: A_light/sqrt(sigma) is sub-horizon (T) ==")41print("  A_light/sqrt(sigma) = c_mag(beta) / sqrt(-ln c_1(beta))  [a ratio of Omega-independent")42print("  character sums and the contact wavefunction, all in units of sqrt(sigma)]. Omega-stable,")43print("  like every hadron ratio m_hadron/sqrt(sigma). Terminal: T. Only the absolute A_light = ")44print("  sqrt(sigma)*c_mag carries the Omega-hard scale (E5).  [approx] A_light/sqrt(sigma) ~ 195/440 ~ 0.44.")4546# ---------------------------------------------------------------------------------------------47print("\n== E1 resolution: the residuals are sub-horizon (T), corroborated by power-counting ==")48# SU(3)-breaking expansion in the sub-horizon parameter eps_s ~ (m_s - m_l)/scale.49# first-order strange splitting (octet, per two seeds): M_Xi - M_N ; decuplet spacing.50mN = (F("938.272")+F("939.565"))/2; mXi = (F("1314.86")+F("1321.71"))/251first_oct = (mXi - mN)                                  # ~ 2 strange units52gmo_res = abs((mN+mXi)/2 - (3*F("1115.683")+(F("1189.37")+F("1192.642")+F("1197.449"))/3)/4)53mD=F("1232.0"); mO=F("1672.45")54spacing = (mO - mD)/355third_diff = abs(mD - 3*((F("1382.8")+F("1383.7")+F("1387.2"))/3) + 3*((F("1531.80")+F("1535.0"))/2) - mO)56r_gmo = float(gmo_res/first_oct)57r_third = float(third_diff/spacing)58print(f"  GMO 27-plet residual / first-order octet splitting = {float(gmo_res):.1f}/{float(first_oct):.0f} = {r_gmo:.3f}  ~ eps_s^2  (eps_s ~ {r_gmo**0.5:.2f})")59print(f"  decuplet 3rd diff   / decuplet spacing             = {float(third_diff):.1f}/{float(spacing):.0f} = {r_third:.3f}  ~ eps_s^2")60print(f"  both residuals are O(eps_s^2) of the leading splitting, eps_s ~ 0.2 a sub-horizon (framed chi_3)")61print( "  parameter: a CONVERGENT small-parameter expansion. (An Omega-hard quantity does not so organise;")62print( "  cf. the delta_0 drive-orientation Gauss sum, which equidistributes.) Terminal: T.")6364print("\nTERMINAL RESOLUTION: E1 -> T, E4 -> T (both sub-horizon, Omega-stable). The baryon sector's")65print("ONLY Omega-hard residue is the absolute scale sqrt(sigma) = Lambda_QCD (E5/D6c), with alpha (E7).")66print("No predicate rests at 'open': each is sub-horizon (T) or the one Omega-hard scale.")67