p6.py — src/27-fields

1# framed-rational status: MIXED -- exact framed-rational core; continuum constructs are confined to clearly-labelled [approx] / degenerate-idealisation comparisons (Tier 2/3/5 or measured-data), never to an exact claim.2"""3P6: asymptotic freedom and Lambda_QCD.45(a) SIGN of asymptotic freedom from the gluon self-coupling (P5):6    b0 = (11/3) C_A - (2/3) n_f T_F  =  11 - (2/3) n_f   for SU(3) (C_A=3, T_F=1/2).7    The +11 is anti-screening from the gluon self-coupling (non-abelian, P5);8    the -2/3 n_f screens.  b0>0 (asymptotic freedom) for the derived matter content;9    an abelian U(1) has b0 = -(4/3) n_f < 0 (screening, no AF).10(b) Lambda_QCD by dimensional transmutation: Lambda = mu exp(-2pi/(b0 alpha_s(mu))),11    exponentially below the substrate (Planck) scale -> the ~20-order QCD/Planck gap.12(c) m_proton ~ Lambda_QCD, so gravity's weakness for ordinary matter13    (m_p/m_P)^2 = exp(-4pi/(b0 alpha_s)) is the SQUARE of the dimensional-transmutation14    factor: O2's hierarchy for protons IS QCD dimensional transmutation, squared.15"""16import math17from fractions import Fraction as Fr1819print("="*68)20print("P6  asymptotic freedom and the QCD scale")21print("="*68)2223# ---- (a) the sign: b0 from the self-coupling ----24def b0_SU3(nf): return Fr(11,1) - Fr(2,3)*nf      # = (11/3)*3 - (2/3) nf25def b0_U1(nf):  return -Fr(4,3)*nf                 # abelian QED-like (screening)26print("\n(a) beta-function coefficient  b0 = (11/3)C_A - (2/3) n_f  (SU(3): 11 - (2/3)n_f)")27print("    the +11 is anti-screening from the GLUON SELF-COUPLING (P5, C_A=3)")28for nf in (0, 2, 6, 16, 17):29    b = b0_SU3(nf)30    print(f"    n_f={nf:2d}:  b0(SU3) = {b} = {float(b):+.2f}  ->  {'asymptotically FREE' if b>0 else 'not free'}"31          f"     [abelian U(1): b0 = {b0_U1(nf)} (screening)]")32print("    one generation n_f=2 (u,d): b0=29/3>0 ; three generations n_f=6: b0=7>0.")3334# ---- (b) Lambda_QCD by dimensional transmutation ----35M_P = 1.220890e19          # Planck mass, GeV36M_Z = 91.187637print("\n(b) Lambda_QCD = mu * exp(-2 pi / (b0 alpha_s(mu)))   (dimensional transmutation)")38# using the measured anchor alpha_s(M_Z)=0.1179, b0(n_f=5)=23/339b0_5 = float(b0_SU3(5))40als_MZ = 0.117941Lam = M_Z*math.exp(-2*math.pi/(b0_5*als_MZ))42print(f"    anchor alpha_s(M_Z=91 GeV)=0.118, b0(n_f=5)={b0_5:.2f}:"43      f"  Lambda_QCD ~ {Lam*1000:.0f} MeV  (observed ~200-300 MeV)")44# the hierarchy from the substrate scale with an O(1) bare coupling45print("    from the substrate scale with an O(1) bare coupling:")46for als_MP, b0v in [(1/40, 7.0), (1/50, 7.0), (1/45, 7.67)]:47    lnratio = 2*math.pi/(b0v*als_MP)48    ratio = math.exp(-lnratio)49    print(f"      alpha_s(M_P)={als_MP:.4f}, b0={b0v}:  Lambda/M_P = exp(-{lnratio:.1f}) = {ratio:.1e}"50          f"   -> Lambda ~ {ratio*M_P*1000:.0f} MeV")5152# ---- (c) the gravity connection: (m_p/m_P)^2 = (dim transmutation)^2 ----53m_p = 0.938354print("\n(c) the proton mass is the QCD scale, exponentially below the Planck mass:")55print(f"    m_p/m_P = {m_p/M_P:.2e} = exp(-{math.log(M_P/m_p):.1f})   (dimensional transmutation)")56print(f"    gravitational coupling  (m_p/m_P)^2 = {(m_p/M_P)**2:.2e} = exp(-{2*math.log(M_P/m_p):.1f})")57print("    => O2's '10^-38 gravitational coupling of protons' IS the SQUARE of the QCD")58print("       dimensional-transmutation factor.  Gravity is weak for ordinary matter")59print("       because the proton mass is exp(-2pi/(b0 alpha_s)) below the Planck scale.")60print("="*68)61print("AF sign derived from the self-coupling; the exponential QCD/Planck hierarchy is")62print("dimensional transmutation; the precise b0 reduces to the matter content (sect 7),")63print("the precise Lambda to the bare coupling (P2/P8) + running.  sigma ~ Lambda^2 (QCD-1).")64print("="*68)65