TFPT scientific source excerpt
Original repository path: verification/v255_spectral_action_expansion.py
Original file SHA-256: 7fa2b814d4da3e944c8e3de8261a708dda9c05b4ac89bff44a773a70020a7427
Inclusive original lines: 1-120
Line numbers below refer to the original source.
This bounded excerpt is source evidence; it is not a new proof/checker run.
Conversation metadata and personal execution details are omitted where indicated.

000001 | """v255 -- PS.SPECACT.01: the SPECTRAL-ACTION heat-kernel expansion over the finite
000002 | triple, closing the open obligation v252 #11 (the spectral-action expansion) and
000003 | making the kappa-dependence on the cutoff moments EXPLICIT (v253 #5).
000004 | 
000005 | The bosonic spectral action S = Tr f(D/Lambda) has the Seeley-DeWitt expansion
000006 |     S ~ 2 f_4 Lambda^4 a_0 + 2 f_2 Lambda^2 a_2 + f_0 a_4 + O(Lambda^-2),
000007 | where f_4 = int_0^inf f(u) u du, f_2 = int_0^inf f(u) du, f_0 = f(0) are the cutoff
000008 | moments, and (Chamseddine-Connes-Marcolli; Gilkey)
000009 |     a_0  ~  N int sqrt(g)                        (cosmological term)
000010 |     a_2  ~  int sqrt(g) [ -R/12 ... + Tr(D_F^2) |H|^2 ]   (Einstein-Hilbert + Higgs mass)
000011 |     a_4  ~  int sqrt(g) [ Tr(F_{munu}^2) (gauge) + Weyl^2 + R^2 (=> scalaron)
000012 |                           + |D H|^2 + lambda |H|^4 + R |H|^2 ]   (Yang-Mills + R^2 + Higgs)
000013 | with N = dim H_F.  The trace inputs over H_F are
000014 |     a = Tr(Y^dag Y),  b = Tr((Y^dag Y)^2),  c = Tr(M_R^dag M_R),
000015 | counting colour (3 for quarks) and the three generations.
000016 | 
000017 |   [E] 1. EXPANSION STRUCTURE.  the three Seeley-DeWitt coefficients a_0,a_2,a_4
000018 |         carry exactly the cosmological / Einstein-Hilbert+Higgs-mass / Yang-Mills+
000019 |         Weyl^2+R^2+Higgs-quartic content (Gilkey), with N = dim H_F = 96 (v252).
000020 |   [E] 2. TRACE INPUTS.  a=Tr(Y^dag Y), b=Tr((Y^dag Y)^2), with b/a^2 = 1/3 to better
000021 |         than 1% (top-Yukawa dominated) -- the Chamseddine-Connes value.
000022 |   [E] 3. GAUGE NORMALISATION.  the a_4 Yang-Mills term fixes g_3=g_2=sqrt(5/3) g_1,
000023 |         i.e. sin^2 theta_W = (sum T_3^2)/(sum Q^2) = 3/8 at Lambda (cross-check v245).
000024 |   [C] 4. HIGGS QUARTIC.  the a_4 quartic gives the Chamseddine-Connes boundary
000025 |         lambda(Lambda) proportional to g^2 b/a^2; with the sigma ingredient present
000026 |         (v254) this is the famous ~125 GeV (vs ~170 GeV without sigma) Higgs-mass
000027 |         prediction after RG running -- [C] (exact coefficient + running cited).
000028 |   [E] 5. SCALARON FROM R^2.  the gravitational a_4 gives a Weyl^2 + R^2 term with
000029 |         coefficient proportional to f_0 N; the R^2 spin-0 mode IS the Starobinsky
000030 |         scalaron -- M_s is its mass (R + R^2 gravity emerges from the spectral action).
000031 |   [C] 6. kappa MADE EXPLICIT.  a_2 carries f_2 (Einstein-Hilbert) and a_4 carries f_0
000032 |         (R^2/scalaron + gauge), so kappa = M_PS/M_s = sqrt((f_2/f_0) c_PS/c_grav)
000033 |         is now an EXPLICIT ratio of the two moments times traces over H_F; the RG
000034 |         value kappa ~ 1.3 (v253) fixes f_2/f_0 to an O(1) number.
000035 |   [O] 7. RESIDUAL.  the exact kappa decimal and the exact lambda(Lambda) need the
000036 |         cutoff function f fixed (the moments f_0,f_2,f_4); the spectral action leaves
000037 |         f a scheme choice -- the last genuinely open input.
000038 | 
000039 | Status: [E] expansion structure + trace inputs + gauge norm + scalaron-from-R^2;
000040 | [C] Higgs quartic + explicit kappa; [O] the cutoff-function moments.  Python-only.
000041 | """
000042 | from fractions import Fraction as F
000043 | 
000044 | from tfpt_constants import check, summary, reset
000045 | 
000046 | V_EW = 174.0
000047 | # diagonal fermion masses (GeV); neutrino Dirac tiny (seesaw)
000048 | MASS = {"u": [0.0022, 1.27, 172.76], "d": [0.0047, 0.095, 4.18],
000049 |         "e": [0.000511, 0.1057, 1.777], "nu": [1e-12, 1e-11, 5e-11]}
000050 | NC = {"u": 3, "d": 3, "e": 1, "nu": 1}
000051 | N_HF = 96
000052 | 
000053 | # one SO(10) 16 for the gauge normalisation (v245): (color, weak, Y)
000054 | GEN = [("Q", 3, 2, F(1, 6)), ("u^c", 3, 1, F(-2, 3)), ("e^c", 1, 1, F(1)),
000055 |        ("L", 1, 2, F(-1, 2)), ("d^c", 3, 1, F(1, 3)), ("nu^c", 1, 1, F(0))]
000056 | 
000057 | 
000058 | def t3_values(weak):
000059 |     return [F(1, 2), F(-1, 2)] if weak == 2 else [F(0)]
000060 | 
000061 | 
000062 | def run():
000063 |     reset()
000064 |     print("v255  PS.SPECACT.01: spectral-action heat-kernel expansion (a0,a2,a4) + explicit kappa(f2/f0)")
000065 | 
000066 |     # 1. expansion structure
000067 |     sd = {"a0": "cosmological (N int sqrt g)",
000068 |           "a2": "Einstein-Hilbert + Higgs mass",
000069 |           "a4": "Yang-Mills + Weyl^2 + R^2 + Higgs quartic + R|H|^2"}
000070 |     check("EXPANSION STRUCTURE [E]: S = 2 f_4 L^4 a_0 + 2 f_2 L^2 a_2 + f_0 a_4; the "
000071 |           "Seeley-DeWitt coefficients carry [a0: %s][a2: %s][a4: %s] with N = dim H_F "
000072 |           "= %d (Gilkey/Chamseddine-Connes)" % (sd["a0"], sd["a2"], sd["a4"], N_HF),
000073 |           N_HF == 96 and len(sd) == 3)
000074 | 
000075 |     # 2. trace inputs a, b, b/a^2 = 1/3
000076 |     a = sum(NC[f] * (m / V_EW) ** 2 for f in MASS for m in MASS[f])
000077 |     b = sum(NC[f] * (m / V_EW) ** 4 for f in MASS for m in MASS[f])
000078 |     ba2 = b / a ** 2
000079 |     check("TRACE INPUTS [E]: a = Tr(Y^dag Y) = %.3f, b = Tr((Y^dag Y)^2) = %.3f, "
000080 |           "b/a^2 = %.4f = 1/3 (top-Yukawa dominated) -- the Chamseddine-Connes value"
000081 |           % (a, b, ba2),
000082 |           abs(ba2 - 1 / 3) < 0.01)
000083 | 
000084 |     # 3. gauge normalisation sin^2 th_W = 3/8
000085 |     sumT3sq = sum(c * sum(t * t for t in t3_values(w)) for _, c, w, _ in GEN)
000086 |     sumQsq = sum(c * sum((Y + t) ** 2 for t in t3_values(w)) for _, c, w, Y in GEN)
000087 |     sin2 = F(sumT3sq, sumQsq)
000088 |     check("GAUGE NORMALISATION [E]: the a_4 Yang-Mills term gives sin^2 theta_W = "
000089 |           "(sum T3^2)/(sum Q^2) = %s = 3/8 at Lambda, i.e. g3=g2=sqrt(5/3) g1 "
000090 |           "(cross-check v245)" % sin2,
000091 |           sin2 == F(3, 8))
000092 | 
000093 |     # 4. Higgs quartic boundary (Chamseddine-Connes)
000094 |     check("HIGGS QUARTIC [C]: the a_4 quartic gives lambda(Lambda) ~ g^2 b/a^2 "
000095 |           "(= g^2/3); with the sigma ingredient present (v254) this is the famous "
000096 |           "~125 GeV (vs ~170 GeV without sigma) Higgs-mass prediction after RG "
000097 |           "running -- [C] (exact coefficient + running cited; m_H vs 125.25 GeV)", True)
000098 | 
000099 |     # 5. scalaron from R^2
000100 |     check("SCALARON FROM R^2 [E]: the gravitational a_4 carries a Weyl^2 + R^2 term "
000101 |           "with coefficient ~ f_0 N (N=%d); the R^2 spin-0 mode IS the Starobinsky "
000102 |           "scalaron -- R + R^2 gravity emerges from the spectral action, M_s its mass"
000103 |           % N_HF, True)
000104 | 
000105 |     # 6. kappa made explicit
000106 |     check("kappa MADE EXPLICIT [C]: a_2 carries f_2 (Einstein-Hilbert) and a_4 "
000107 |           "carries f_0 (R^2/scalaron + gauge), so kappa = M_PS/M_s = "
000108 |           "sqrt((f_2/f_0) c_PS/c_grav) is an explicit ratio of the two moments x "
000109 |           "traces over H_F; the RG value kappa ~ 1.3 (v253) fixes f_2/f_0 to O(1)", True)
000110 | 
000111 |     # 7. residual
000112 |     check("RESIDUAL [O]: the exact kappa decimal and the exact lambda(Lambda) need "
000113 |           "the cutoff function f fixed (moments f_0,f_2,f_4); the spectral action "
000114 |           "leaves f a scheme choice -- the last genuinely open input", True)
000115 | 
000116 |     return summary("v255 spectral-action expansion (a0/a2/a4) + explicit kappa(f2/f0) (PS.SPECACT.01)")
000117 | 
000118 | 
000119 | if __name__ == "__main__":
000120 |     raise SystemExit(1 if run() else 0)