All notable changes to this project will be documented in this file.
compute_emergent_pulse / SDK net.rhythm() give the collective network
pulse (resonances ω_k = √λ_k, fundamental, dominant beat, vibration energy);
compute_nodal_pulse / SDK net.resonance() give the per-NFR pulse — every
NFR a phase oscillator at its own νf and phase φ, coupled by resonance
(local_phase_sync per NFR, the Kuramoto order R, gate Δφ_max = π/2). The
collective pulse emerges as the per-NFR pulses lock (R → 1).net.pulse_trajectory(steps) evolves a copy and records
the rhythm forming over time — R(t), C(t), the per-NFR local resonance —
surfacing the local-before-global synchronization cascade (clusters lock
before the global rhythm). net.evolve(record=True) + net.history() surface the
engine's own canonical per-step series (kuramoto_R, C_steps, phase_sync,
Si_mean).compute_unified_telemetry now carries both the
dissipative read-out (canonical tetrad + coherence, relaxes to ΔNFR = 0) and the
conservative pulse + resonance blocks (which do not saturate).emergent_fractal_pulse.py
— resonance locks scale by scale on a self-similar network (the temporal face
of U5); emergent_arithmetic_pulse.py — the residue-NFR pulse tone-count is the
proved cyclotomy law s_k(p) = gcd(k, p−1) + 1, so a prime is its most
degenerate chord and the factorization type is the chord size. The pulse is woven
into EMERGENT_ONTOLOGY.md (§2.1/§2.2/§5.5) and TNFR_NUMBER_THEORY.md (§9.12).emergent_musical_nfr.py reads the structural-frequency spectrum (νf,
ω_k = √λ_k, Hz_str) through music — not audio. The dynamical regime
follows the dimension: a 1D thread (a string) is harmonic and its just
consonances (octave/fifth/fourth) are emergent (pitched); a 2D+ form (a
drum) is inharmonic (unpitched), where consonance is the U3 phase gate;
polyphony = primes (the Euler product); and the Kac wall — isospectral
NFRs share the pulse, so it hears the type, not the identity (Fix(G)^⊥). Only
equal temperament / the chosen scale are imposed. Woven into
EMERGENT_ONTOLOGY.md §5.5 + the form → dimension → dynamics synthesis (§2.1).u2_compliance structural-time
benchmark, the N1–N14 milestone examples (77–86, 104, 105) and the old
operator.py — was retired. New foundation src/tnfr/navier_stokes/: a faithful
lean pseudo-spectral 3D integrator (TNFRNavierStokes, rotational form, exact
Leray projection, integrating-factor RK2) + conservative_face.py
(verify_diffusive_face, face_of_flow, vorticity_modal_spectrum,
measure_cascade_frontier).ν_f = ν; verify_diffusive_face VALID for every physical viscosity, recovering
ν_f = ν). Blow-up is therefore a purely nonlinear K_φ cascade (the
vortex-stretching VAL source), not a linear resonance — unlike oscillatory data
(EEG), which sits on the under-damped conservative face.measure_cascade_frontier, example 158): at
matched structural time τ_str = ν·t every run saturates at fixed Re (the
diffusive face regularises) and the peak enstrophy debt grows with Re
(1.00 → 1.04 → 1.68 at Re 126/314/628). The Re → ∞ cascade bound = Clay, open.examples/10_applications/159_empirical_confrontation_pipeline.py packages the
empirical arm's workflow with engine primitives: map a multichannel signal onto the
emergent phase-locking graph, read the canonical magnitudes (the pulse, the tetrad,
ξ_C, Kuramoto R), and diagnose its face with the engine's own
verify_overdamped_projection — making the theory falsifiable against data (the face
is measured, not assumed). Cross-program face map: oscillatory data → conservative
face; linear NS → diffusive face.|ΔNFR| perturbation relaxes geometrically under the discrete nodal step
(q = 1 − νf·dt·ρ, with ρ = trace(L_rw)/N = 1, exact). Read two ways from the
same q: the relaxation time to the coherence band 1/(π+1) is the U4b
recency window (and the GRAMMAR repeat-avoidance window) = 3; the
geometric absorption capacity ⌊1/(1−q)⌋ is the U2 debt threshold = 2.
New derive_bifurcation_window_from_physics / derive_u2_debt_capacity_from_physics
(config/physics_derivation.py) replace the literal 3/2 — with no e (the
canonical relaxation is the discrete geometric decay qⁿ, not the continuous
exponential e^{−νf λ t}, which is only the dt → 0 limit the engine never takes).
The earlier graduated destabilizer split (strong = 4 / moderate = 2) was a
heuristic the dynamics does not support and has been dropped — one emergent
window for every destabilizer (the relaxation rate ρ = trace/N = 1 is
topology-independent, so the window is a topology-independent constant).structural_eigenmodes and relaxation_spectrum
now share one memoized eigendecomposition of the symmetric normalized Laplacian,
keyed on the graph topology (self-invalidating when nodes/edges/weights change).
The spectrum is invariant under evolution on a fixed graph, so the O(N³) eigh
runs once per topology instead of once per pulse/spectrum read-out.EMERGENT_DERIVATION_PLAN.md.
The φ/γ/e purge had left two load-bearing weight sets frozen at their literal
φ/γ decimals (DNFR_WEIGHTS/SI_WEIGHTS = {0.737, 0.155, 0.09} where
0.737 = φ/(φ+γ)) and had replaced the operator gains with arbitrary
"operational" decimals (IL=0.75, OZ=2.0, SHA/NUL=0.9, VAL=1.05). These are
used numerically in every ΔNFR and Sense-Index evaluation, hence in every
recorded result. They are now emergent:
π/(π+1) of the
remainder: (π/(π+1), π/(π+1)², 1/(π+1)²) — which normalises to exactly 1
(π/(π+1) + π/(π+1)² + 1/(π+1)² = (π+1)²/(π+1)² = 1). Ordering by structural
primacy (phase ≻ EPI ≻ νf; topo inactive). SI_WEIGHTS takes the same hierarchy.IL = π/(π+1), OZ = (π+1)/π (a balanced IL∘OZ is exactly
isometric). Capacity lever (νf, slow): the gentle π-step δ = 1/(4π) —
SHA/NUL = 1−δ, VAL = 1+δ, NUL_densification = 1/(1−δ) (volume
conservation). Secondary couplings on the π-fraction ladder (1/(4π), 1/(2π), 1/(8π)); ZHIR θ-shift 1/π; NAV_eta/REMESH_alpha = the unit midpoint 0.5.SELECTOR_WEIGHTS takes the same coherence-band
hierarchy; the coherence triggers are the high-coherence gate π/(π+1), the
new rectified-mean level 2/π, and the unit midpoint/quarter 0.5/0.25;
AU_CURVATURE is the exact midpoint (0.9π+π)/2 of the strict K_φ gate and the
π wrap; the phase couplings, FEEDBACK tolerances/rates, OZ noise and THOL
metabolic weights are π-fractions (1/(2π), 1/(4π), 1/(8π)); the get_factor
safety fallbacks reference the emergent constants. The selector magnitude
thresholds (dnfr_hi/lo, accel_hi/lo) are honestly left operational
(|ΔNFR|/∂²EPI scale, not coherence — π-flavouring them would repeat the φ/γ/e
naming-convention error).constants/canonical.py:
CHANNEL_WEIGHT_PRIMARY/SECONDARY/TERTIARY, COHERENCE_RETENTION,
DISSONANCE_AMPLIFICATION, COUPLING_GENTLE/MODERATE/FINE,
MID_COHERENCE_THRESHOLD. Full suite green (2201 passed) after each stage;
the dynamics stay bounded (U2). Recorded research results computed with the
old constants still require recomputation (planned Stage 5).examples/02_physics_regimes/37_operator_tetrad_synergy.py imported the purged
GAMMA/PHI from constants/canonical → ImportError; examples aren't in the
test suite so it had slipped through) — it now runs. Updated the
coherence_projector_sense_index benchmark SI_WEIGHTS to the band hierarchy,
removed dead PHI/GAMMA/E constants from boundary_vibration, de-refuted the
phase_wall correspondence comments (its TEST-4 obstruction result — building
φA+γL+πL²+eK to prove the four constants are insufficient — is kept), and
replaced ~27 stale "(φ,γ,π,e) remain the assumed substrate" claims with "π"
across 14 benchmark files. The legitimate emergent-object studies are kept
(the Kuramoto φ-as-Fibonacci-limit, the golden-angle sphere sampling, Euler
products, the Γ chirality matrix, tetrahedral symmetry groups).ΔNFR=0), Riemann σ_c/GUE and exact S_n equivariance (‖[L, P_σ⊗P_τ]‖ = 0),
Navier–Stokes (a pseudo-spectral solver that never reads the operator
gains), conservation, the tetrad relations (K_φ = L_rw·φ, ξ_C ∝ 1/√λ₂), and
Yang–Mills U6 confinement all derive from the graph Laplacian, the spectral gap,
S_n symmetry, or unit arithmetic. Only the dynamic trajectories (C(t)/Si
curves, network-optimization outcomes) shift, with their qualitative attractors
invariant. The φ/γ/e and arbitrary operational decimals were therefore never
load-bearing: the refactor both cleans the foundation and proves the results
are genuinely emergent. See EMERGENT_DERIVATION_PLAN.md §7.tests/core_physics/test_emergent_constants_guard.py) pins the channel
weights, operator gains and coupling ladder to their exact π-formulas and
asserts that no nodal-physics constant equals a removed frozen φ/γ/e decimal
— so the obsolete constants cannot creep back. Cleaned the final φ/γ/e comment
straggler the input purge had missed in src/: the extended nodal system's
_compute_phase_transport_derivative (dynamics/canonical.py) no longer cites
the dead ≈ 0.618 = 1/φ / 0.155 / 0.135 origins, nor the false "RECALIBRATED
from canonical constants" / dead "Import canonical constants" comments (values
unchanged — honest operational magnitudes on the optional J_φ-transport path).
EMERGENT_DERIVATION_PLAN.md Stage 0 is now complete.AGENTS.md (+ the .github/agents/my-agent.md mirror),
ARCHITECTURE.md, CONTRIBUTING.md, theory/, docs/grammar/, examples, and
code docstrings: the Φ_s confinement bound is π-derived — drift
Δ Φ_s < π/2 ≈ 1.571 (half phase-wrap) and per-node |Φ_s| < π/4 ≈ 0.785
(quarter phase-wrap) — replacing the old φ ≈ 1.618 / empirical 0.7711 framing;
the strong-coherence cut is the emergent band gate π/(π+1) ≈ 0.7585
(replacing the frozen (e·φ)/(π+e) ≈ 0.7506). Corrected a propagated arithmetic
error: π/(π+1) is 0.7585, not 0.7616 (it must complement 1/(π+1)=0.2415).
The SDK COHERENCE_STRONG now aliases the emergent HIGH_COHERENCE_THRESHOLD
(π/(π+1)); MIN_BUSINESS_COHERENCE (0.75) stays the separate operational
business-health knob.theory/SPIRAL_ATTRACTORS_AND_LOGARITHMIC_DYNAMICS.md and its demo
(examples/02_physics_regimes/32_spiral_attractors_demo.py) — a false
φ/γ/e-era claim ("golden ratio as dynamical attractor", "fourth constant γ").
The demo imported the purged φ/γ/e constants (so it could not run, making the
document's "Validated" status false) and the golden-attractor check was circular
(it set b = 2·ln(φ)/π by hand, then "verified" quarter-turn ratios = φ). Only
the trivial, non-distinctive kernel (log spirals appear in a rotation + growth
regime, with a free b = νf·k/ω) was true; φ is not selected by the dynamics.
References cleaned from theory/README.md, FUNDAMENTAL_THEORY.md, and
examples/README.md.theory/ document audit (every theory/*.md) for residual
φ/γ/e false claims, with no further deletions needed — SPIRAL_ATTRACTORS was
the only doc with a false thesis; the rest are genuinely emergent and carried
only scattered stale refs (now fixed). The most significant correction purges the
refuted "Universal Tetrahedral Correspondence" (the φ↔Φ_s, γ↔|∇φ|, π↔K_φ,
e↔ξ_C mapping) from TNFR_RIEMANN_RESEARCH_NOTES.md (20 references) — the explicit
mapping becomes the minimal structural-field tetrad (only π is structural), the
three inter-prime coupling kernels are relabeled exploratory, not canonical, and
the stale DNFR_/SI_/SELECTOR_WEIGHTS derivation claims/anchors are corrected to
the operational defaults_core.py values. Operator-gain tables across
STRUCTURAL_OPERATORS, STRUCTURAL_CONSERVATION_THEOREM,
STRUCTURAL_STABILITY_AND_DYNAMICS, TNFR_VARIATIONAL_PRINCIPLE,
TNFR_YANG_MILLS_RESEARCH_NOTES (+ 2 yang_mills/structural_gap.py docstrings),
CATALOG_TYPE_HYGIENE_PROGRAMME, and TNFR_NUMBER_THEORY were updated from frozen
φ/γ/e formulas (e.g. IL φ/(φ+γ)≈0.737→0.75, OZ φ/γ≈2.803→2.0, NUL
densification 2.803→1/λ≈1.111, U6 Δ Φ_s < φ→π/2, |∇φ| heuristic
γ/π→π/16) to the operational engine values. No engine code changed (the code
was already purged; only doc text and 2 cosmetic docstrings).canonical.py is now pure physics)constants/canonical.py into a new dedicated module
constants/operational.py (explicitly engine tuning, NOT TNFR physics).
canonical.py now holds 89 numeric constants, all genuine structural /
physics quantities (π phase-wrap bounds, spectral-gap ξ_C, the coherence band,
operator gains, tetrad / phase / νf / EPI / KL / DT scales); the 150 moved
knobs (caches, FFT tuning, optimization speedup/performance estimates,
pattern-discovery confidence, integration baselines, operator scoring weights)
live in operational.py. The new module imports only PI from canonical
(one-way dependency; canonical never imports operational), and a parallel
engines/constants/operational.py star-shim mirrors the existing canonical
shim. The canonical ∪ operational union reproduces the pre-split constant set
exactly (verified name→value, 0 leaks / 0 drift). 26 consumer modules were
redirected; mixed importers were split to preserve their structural imports.|∇φ| ≤ π, |K_φ| < 0.9·π); the coherence length is set by the spectral gap
(ξ_C ∝ 1/√λ₂); every other parameter is derived from the nodal dynamics or is
a free operational parameter.PHI_S_VON_KOCH_THRESHOLD = π/4 ≈ 0.785 (quarter phase-wrap) and drift
U6_STRUCTURAL_POTENTIAL_LIMIT = π/2 ≈ 1.571 (half phase-wrap), replacing the
empirical 0.7711 / golden-ratio (φ ≈ 1.618) framing.derive_tetrad_threshold_values and the φ/γ/e accumulation-law
threshold-derivation machinery (ThresholdDerivation). Operator gain magnitudes
are now plain operational parameters — the theory fixes each operator's channel
and sign via its contract, not its magnitude.φ/γ/e references from
source comments, docstrings, and the documentation set (ARCHITECTURE.md,
README.md, CHANGELOG.md, .zenodo.json, CONTRIBUTING.md,
benchmarks/README.md, and the theory/ + docs/ notes).EMERGENT_CANON_AUDIT.md).
The purge had left the numeric values frozen (e.g. K_TOP_FALLBACK still held
2.803171 = φ/γ); those magic numbers are now re-derived or eliminated so the
canonical base is genuinely emergent.MATH_DELTA_NFR_THRESHOLD = 0.5 is the unit-gap midpoint: with unit arithmetic
ΔNFR coefficients, prime ⟺ ΔNFR = 0 exactly and every composite has
ΔNFR > 1, so any cut in (0, 1) separates them.MAX_STRUCTURAL_FREQUENCY = 2π, MIN_STRUCTURAL_FREQUENCY = 1/(2π),
AU_CURVATURE_PERMISSIVE = 0.96·π, CRITICAL_EXPONENT = GRAD_PHI_CANONICAL_THRESHOLD = π/16,
DYNAMICS_SI_HI = π/(π+1), the K_TOP clamp 1/(8π) … 1.0 and fallback π.PHI_S_THRESHOLD, GRAD_PHI_THRESHOLD, and K_PHI_THRESHOLD = 3.2275, which
exceeded the π phase-wrap bound and was therefore an unreachable no-op check).MEDICAL_*, BUSINESS_*, EXAMPLE_*,
VIZ_*, CLI_*, THERAP_*, SCRIPT_*, TOOL_*, UTILS_*) and the dead
CANONICAL_CONSTANTS registry; relocated the SDK builder defaults into
sdk/builders.py. The remaining ~180 operational engine knobs were rounded to
plain ≤2-decimal values (dropping the false φ/γ/e precision). constants/canonical.py
shrank from ~770 to ~565 lines.operators/__init__.py) to the canonical
SHA_VF_FACTOR / NUL_SCALE_FACTOR / VAL_SCALE_FACTOR, and removed residual
inline artifacts (10·φ, e, 4/(e+φ)) in bifurcation.py, variational.py,
cycle_detection.py, and signatures.py.A computational audit of the "Universal Tetrahedral Correspondence" found that only π is a genuine structural scale; the four-constant correspondence (φ↔Φ_s, γ↔|∇φ|, e↔ξ_C) is mostly an organizing overlay. Several thresholds asserted as "derived" were empirical, inert (magic), or measured false. This work corrects the claims and replaces magic thresholds with emergent, system-measured quantities. The nodal equation, the 13 operators, and grammar U1–U6 are unchanged.
|∇φ| and K_φ (both are means of wrapped angles, ≤ π). γ, e, φ are
recoverable as mathematical identities but are NOT the structural scales of
their tetrad fields. K_φ = L_rw·φ (the central operator on phase, corr ≈ 1);
ξ_C ∝ 1/√λ₂ (spectral gap, not base e).|∇φ| bound corrected from γ/π ≈ 0.1837 to the phase-wrap bound 0.9π
in physics/variational.py, symmetric with K_φ. The measured synchronization
onset is ≈ 0.29 and σ-dependent, NOT the constant γ/π; γ/π is retained
elsewhere only as a heuristic early-warning level, explicitly labelled
non-derived.derive_tetrad_threshold_values rows re-statused: π geometric; φ, γ, e
overlay (recoverable identities, not structural scales).ARCHITECTURE.md, .zenodo.json, CONTRIBUTING.md, theory/README.md,
docs/STRUCTURAL_FIELDS_TETRAD.md, constants/canonical.py — removed the
"Universal Tetrahedral Correspondence foundation / 100% derived / zero
empirical tuning / verified to machine precision" claims; replaced with the
honest tiering (π genuine; γ/e/φ notational overlay).ARCHITECTURE.md
(the "Mathematical Purity / 497 magic numbers eliminated / zero empirical"
sections), README.md, theory/FUNDAMENTAL_THEORY.md, theory/GLOSSARY.md,
theory/EXTENDED_FIELDS_AND_DERIVED_QUANTITIES.md, AGENTS.md (+ mirror),
docs/grammar/PHYSICS_VERIFICATION.md, and the central constants modules
(telemetry/constants.py, mathematics/unified_numerical.py,
config/defaults_core.py, physics/signatures.py,
operators/grammar_telemetry.py, physics/emergent_chemistry.py) plus four
benchmarks. Exposed cosmetic "derivations" (e.g. MIN_BUSINESS_SENSE_INDEX = 1/φ + 0.082 ≈ 0.700, ⌊φ×10⌋ = 16) as calibrated/notational values.physics/phase_transition.py fully redesigned to emergent sampling-noise
z-scores. The "universal critical exponent γ_c = γ/π" was measured false
(the fitted exponent is protocol-dependent), and the classification noise
floor (γ/π)² was proven inert (it sat in a two-order-of-magnitude gap;
sweeping it changed no classification). Removed the magic constants GAMMA_C,
ORDER_PARAMETER_NOISE_FLOOR, CHIRALITY_THRESHOLD and the
theoretical_exponent field; added symmetry_zscore(mean, var, n) = |mean|/√(Var/N) and the single cut Z_SIGNIFICANCE = 1 (the sampling-noise
scale, not a tunable constant). classify_phase(order_z, chirality_z) now
decides phases from statistical significance measured from the system itself.constants/canonical.py (CRITICAL_EXPONENT, GRAD_PHI_CANONICAL_THRESHOLD,
PHASE_GRADIENT_THRESHOLD_CANONICAL) and mathematics/unified_numerical.py;
gauge.py, emergent_chemistry.py, interactions.py regime thresholds
marked calibrated/heuristic, not derived.primality-test/,
factorization-lab/), examples, benchmarks, the mathematical_purity tests
(now check genuine bounds, not the refuted mapping), all theory docs
(FUNDAMENTAL_THEORY.md §4 and GLOSSARY renamed to "structural-field
tetrad"), and ~40 in-code combo comments (defaults_core.py, bifurcation.py,
cycle_detection.py, number_theory.py, etc.) now marked notational. Exposed
the primality coefficients (ζ=φγ, η=(γ/φ)π, θ=1/φ) as combos chosen to
approximate the original empirical values (ζ=1.0, η=0.8, θ=0.6).A measure-first review of whether calculations (not just narrative) depended on the refuted tetrad values. Main finding: the significant results are robust, because they emerge from STRUCTURE (integer orderings, exact zeros, relative scores) rather than the scale values (γ/π, etc.):
nu_excitation (=γ/π), nu_0, coherence_gap fields were DEAD (defined,
never consumed) and were removed; only theta_valence enters, as a positive
scale (the zero is robust to its value).dominant_regime is decided by relative
scores, NOT the γ/π threshold; above_threshold (which used γ/π) is metadata.Real consequences corrected:
PHASE_CURVATURE_ABS_THRESHOLD = φ×π ≈ 5.083 was a non-physical K_φ bound
(|K_φ| ≤ π by phase wrap, so any check using it was a no-op); it was dead code,
corrected to 0.9π ≈ 2.827.STRUCTURAL_STABILITY_AND_DYNAMICS.md §2.2 still described the old γ_c
classification table; updated to the emergent z-score rule.MATHEMATICAL_DYNAMICS_BASIS.md |∇φ| claim. Full suite green (2196 passed).Completed the full emergent rework of the two studies whose conceptual base was the (now-refuted) four-constant overlay. Both were rebuilt to rest on STRUCTURE alone (measure-first; no value replaced by another magic value):
physics/gauge.py): removed the
three overlay threshold constants REGIME_DOMINANCE_THRESHOLD (1/φ),
REGIME_STRONG_THRESHOLD (γ/π, "Kuramoto critical coupling in gauge") and the
unused REGIME_SECONDARY_THRESHOLD (γ/(π+γ)). The per-sector above_threshold
activity flags now use a single parameter-free criterion: a sector is active
when its normalised score exceeds the equipartition share 1/N_REGIMES = 0.25
(the maximum-entropy reference, derived from the number of gauge sectors — the
four structural channels of the tetrad). Uniform across all four sectors; no
overlay constant. New public symbols N_REGIMES, REGIME_ACTIVITY_SHARE
replace the removed thresholds. dominant_regime (relative max of scores)
is unchanged — it was already robust. Measured: the criterion always flags the
dominant sector and additionally marks genuine co-active secondaries.physics/emergent_chemistry.py): removed
the last free scale parameter (theta_valence = 1/φ) and the now-trivial
EmergentChemistryParameters dataclass. ΔNFR_chem(Z) is now the integer
structural distance of the outer shell to a closed configuration, in natural
units (one subshell step = 1) — the exact chemical analogue of primality
ΔNFR(n)=0. Noble gases (2,10,18,36,54,86) → ΔNFR=0; halogens/alkali → 1;
oxygen → 2; carbon → 4. Magic numbers and the octet are unchanged (they always
emerged from the integer (n+l) ordering and the exact zero).test_gauge.py: TestRegimeActivityCriterion, equipartition
consistency). Full suite green (2195 passed, 2 skipped).This release consolidates the emergent-geometry program, centralizes the operator/grammar/contract layer onto single canonical sources, opens three new TNFR-native Millennium-problem programs, and refactors the documentation to the current engine state. The 13-operator catalog, grammar U1–U6, and the nodal equation are unchanged; everything below either measures structure the nodal equation already contains or removes duplication. Full suite: 2043 passed, 2 skipped.
The nodal equation generates its own geometry; the graph is only the data
substrate. The conservation laws of physics/conservation.py are consolidated
into an explicit emergent symplectic phase space that the engine measures
rather than postulates.
src/tnfr/physics/symplectic_substrate.py — phase space
P = ℝ^{4N} with conjugate pairs (K_φ, J_φ) (geometric) and (Φ_s, J_ΔNFR)
(potential); symplectic 2-form ω (antisymmetric, non-degenerate, closed);
canonical Poisson brackets; H_sub = ½Σ(K_φ²+J_φ²+Φ_s²+J_ΔNFR²) equal to the
energy functional exactly; Liouville div(X_H)=0 (the 13 operators are
symplectomorphisms).H_sub; geometric U(1) → E_geo = ½Σ|Ψ|²;
potential U(1) → E_pot); the compatible Hermitian / flat-Kähler triple
(ω, J, g) with J = −ω — so the i in Ψ = K_φ + i·J_φ is the complex
structure the substrate induces; complete integrability (action–angle,
Liouville–Arnold); Poincaré–Cartan integral invariants; Marsden–Weinstein
symplectic reduction; and the hidden U(2) polarization symmetry whose
SU(2) part supplies three conserved Stokes parameters on the per-node
Poincaré sphere (classical wave polarization — Stokes 1852 / Poincaré 1892 —
not isospin or qubits).physics/variational.py
derive_tetrad_threshold_values recovers φ (inverse-square self-similar fixed
point), γ (harmonic-accumulation gap), e (memoryless-decay series) from each
tetrad field's accumulation law; π remains a geometric primitive.verify_substrate_geometry(G) bundles all
certificates into a SubstrateGeometryReport.Network.symplectic_substrate() + SymplecticReport, in
TNFR.analyze().conservation.py +
variational.py; it does not resolve any open program.examples/08_emergent_geometry/98, 106, 114.The EPI channel of the canonical ΔNFR is the random-walk graph Laplacian
−L_rw·EPI (verified to residual ~1e-16), so the nodal equation is literally a
discrete diffusion equation with diffusivity νf. From this single identity the
engine measures, in TNFR's own variables, a tower of empirically-established
transport phenomena.
src/tnfr/physics/structural_diffusion.py — six transport
layers: diffusion/synchronization (Fourier/Fick/Kuramoto), overdamped drift
(q̇ = νf·F, Stokes/Einstein mobility — corrects the prior "Newton's second
law" reading: the bare first-order nodal equation is overdamped, νf is mobility
not inverse mass), discrete standing-wave modes (bounded-manifold Laplacian
eigenmodes), structural-stability dispersion relation (σ_k = r − νf·λ_k, the
spectral form of U2), random walk + effective resistance (Ohm/Kirchhoff), and
structural flow (current, Kirchhoff continuity, Ohm).q̈ + γq̇ + Lq = 0 with νf = 1/γ; the γ-dial
spans diffusion (γ→∞) to standing waves (γ→0).λ_2 is purely topological and does not encode any canonical
constant (measured negative result).examples/08_emergent_geometry/99, 113, 134, 135.src/tnfr/operators/operator_contracts.py —
the single source of truth for what each operator does to node state, anchored
to the direct _op_* effect (TNFR.pdf §2.2.1). Each OperatorContract records
the public English name, the primary_channel (one nodal-equation channel:
EPI / νf / θ / ΔNFR), the scale (NODE for twelve operators, NETWORK for the
U5 operator REMESH), and a verifiable postcondition. The proactive audit
(audit_operator_contracts), the reactive integrity monitor (POSTCONDITIONS),
and the introspection metadata now all derive from this spec — eliminating the
historical drift where scattered copies disagreed (e.g. AL claiming "positive
ΔNFR" though _op_AL only raises EPI; RA checked for EPI increase though it
preserves identity; VAL/NUL checked |EPI| though they scale νf).E = ½Σ(Φ_s²+|∇φ|²+K_φ²+J_φ²+J_ΔNFR²) contains no EPI or νf term (measured:
scaling EPI or νf leaves E unchanged). The per-operator Lyapunov role in
physics/lyapunov.py is therefore re-derived from the canonical grammar U2 role
(config.physics_derivation), not from a hardcoded energy algebra:
stabilisers {IL, THOL}, destabilisers {OZ, ZHIR, VAL}, the rest neutral. The
form-channel operators (AL, EN, RA, REMESH) are energy-neutral because EPI is
absent from E.examples/08_emergent_geometry/152,
examples/02_physics_regimes/115.config.physics_derivation and re-exported by
operators/grammar_types.py. Every grammar consumer — the U1–U6 validator,
the secondary sequence validator, grammar_dynamics, the runtime preconditions,
the error factory, and the operator metadata — now reads the single source.
Parallel hardcoded copies (including a secondary validator that wrongly listed
NUL as a U2 destabilizer) were removed and pinned by
tests/operators/test_grammar_canonical_consistency.py.operators/grammar_canon.py materializes the
U1–U6 role table, the five-type structural typology, and the canonical glyphic
macros (anchored to TNFR.pdf §2.3), with a self-consistency check.THOL[...] lifts the glyphic sub-language to
context-free (Dyck/Catalan); the emergent operator distribution is the
Shannon–Parry maximum-entropy equilibrium.examples/08_emergent_geometry/139–152.{ΔNFR = 0}; numbers as a
coupled network (Ω-graded centrality, primes as the transport periphery); the
nodal flow on numbers (primes as equilibria, not attractors); primality as
grammatical inertness; numbers as free-monoid words with the dual-lever as the
two additive gradings (count Ω → ΔNFR pressure, size log → νf capacity); the
capacity arm carries von Mangoldt and the prime-ladder Hamiltonian P14 is the
capacity-arm operator — locating the Riemann oscillatory obstruction on the
capacity axis the per-node substrate is blind to.examples/07_number_theory/94–97, 100–102,
116, 146–149.Three new programs join Riemann / Navier–Stokes / Yang–Mills. None claims a solution — each carries an explicit honest-scope statement and classified obstruction.
O(|E|) but synthesizing a globally coherent one
by relaxation traps in dissonance basins (measured global-optimum hit rate
drops monotonically with problem size on frustrated MAX-CUT). Mirrors P≠NP;
Branch B open. theory/TNFR_P_VS_NP_RESEARCH_NOTES.md,
examples/09_millennium/109.a_p = p+1−#E(F_p) as structural pressure;
the accumulated product reproduces the original 1965 empirical rank separation
by brute-force point counting. GL(1)→GL(2) gap open; Branch B.
theory/TNFR_BSD_RESEARCH_NOTES.md, examples/09_millennium/110.theory/TNFR_HODGE_RESEARCH_NOTES.md,
examples/09_millennium/111..apply()), added the
emergent-geometry section to theory/FUNDAMENTAL_THEORY.md, rewrote the energy
classification in theory/STRUCTURAL_OPERATORS.md /
theory/STRUCTURAL_STABILITY_AND_DYNAMICS.md to the emergent/grammar-U2 frame,
and updated counts.01_foundations …
10_applications), resolving the prior 77–86 numbering collision; each file
keeps a stable global number. Foundational examples refactored to the canonical
Kuramoto phase-synchrony physics..md
links repo-wide, pruned 9 obsolete docs/ files, rebuilt the theory/ hub,
and resynced the derived .github/agents/my-agent.md mirror.tnfr-validate console-script entry point in pyproject.toml.auto_optimize and
evolve_grammar_aware; added a proactive measured operator-contract fidelity
audit (net.audit_operators()).The Yang–Mills (Y1–Y5) and REMESH-∞ / Navier–Stokes (N15–N17) program milestones below were developed earlier in the cycle and are part of this release. Each carries an explicit honest-scope statement; none resolves a Clay Millennium Problem.
BRANCH_B_OBSTRUCTION_CLASSIFIED — Y1–Y4 establish a finite TNFR U(1) structural gauge diagnostic surface, but Clay-strength closure requires a new canonical non-Abelian derivation plus a continuum / thermodynamic lower-bound theorem.classify_yang_mills_closure() in src/tnfr/yang_mills/closure.py, exported from tnfr.yang_mills with YangMillsClosureReport.A_FINITE_U1_DIAGNOSTIC_SURFACE when Y4 reports stable finite positive gaps.B_REQUIRES_NEW_CANONICAL_NONABELIAN_DERIVATION because Y3 remains OPEN_DERIVABILITY_GAP.clay_problem_resolved = False; the obstruction is localized, not removed.tests/physics/test_yang_mills_closure.py cover Branch-B classification, report reuse, sampled collapse handling, and package-root import. Y1–Y5 focused run: 33 passed.FINITE_SCALING_EVIDENCE or GAP_COLLAPSE_OBSERVED depending on sampled finite graph families. This is a finite diagnostic only, not a continuum theorem.run_finite_scaling_study() in src/tnfr/yang_mills/scaling.py, exported from tnfr.yang_mills with FiniteScalingPoint and FiniteScalingReport.n under fixed U6 target ratios ρ_U6 = max_i |Φ_s(i)| / φ; grouped reports fit finite log-log slopes of mean gap versus n.tests/physics/test_yang_mills_scaling.py cover report shape/scope, grouped finite scaling, reproducibility, sampled collapse classification, invalid input rejection, and package-root import. Y1–Y4 focused run: 29 passed.OPEN_DERIVABILITY_GAP — audited candidate routes for deriving a non-Abelian / multi-channel gauge sector from TNFR-internal data only; no route is promoted to canonical status.audit_nonabelian_derivability() in src/tnfr/yang_mills/derivability.py, exported from tnfr.yang_mills with NonAbelianCandidateAudit and NonAbelianDerivabilityReport.Ψ = K_φ + i·J_φ gauge structure supplies a scalar local U(1) connection. Nested EPI or operator-history data do not yet derive component-mixing parallel transport or non-commuting generator algebra; cycle-basis routes require non-canonical basis/orientation selection.tests/physics/test_yang_mills_derivability.py cover baseline U(1) confirmation, nested-EPI obstruction, cycle-bundle rejection, unsupported route errors, and package-root import. Y1+Y2+Y3 focused run: 23 passed.EMPIRICAL_FINITE_GRAPH_ONLY — finite sweep surface created for testing how the Y1 structural gauge gap behaves across U6-confined and U6-unconfined regimes.run_u6_confinement_sweep() in src/tnfr/yang_mills/u6_sweep.py, exported from tnfr.yang_mills.ρ_U6 = max_i |Φ_s(i)| / φ; ρ_U6 < 1 is U6-confined and ρ_U6 ≥ 1 intentionally probes unconfined finite structural-potential regimes.tests/physics/test_yang_mills_u6_sweep.py cover report shape/scope, U6 target tracking, gap contracts, reproducibility, invalid input rejection, and package-root import. Y1+Y2 focused run: 18 passed.DIAGNOSTIC_SURFACE_CREATED — first TNFR-native Yang–Mills / structural mass-gap attack surface implemented as a finite-graph diagnostic, not a Clay-strength proof.src/tnfr/yang_mills/ with build_structural_gauge_graph(), build_structural_gauge_gap_operator(), and compute_structural_gauge_gap().H_YM^TNFR = L_A + V_F + V_U6, where L_A is the gauge-covariant graph Laplacian from A_ij, V_F is cycle-curvature potential from F_C²/π², and V_U6 is structural-potential confinement from Φ_s²/φ².tests/physics/test_yang_mills_structural_gap.py cover graph construction, self-adjointness, non-negative finite gap reporting, seeded local-U(1) spectral invariance, reproducibility, package imports, and no EPI/phase mutation. Focused run: 13 passed.theory/TNFR_YANG_MILLS_RESEARCH_NOTES.md records the Y-series gap ledger and updates the next target to Y2 (U6 confinement sweep).ANALYTICAL_CONSISTENT_CONDITIONAL — K41 spectrum derived conditionally from TNFR grammar rules U2+U3+U5+CDC; algebraically closed given the Cascade Development Condition.energy_spectrum_3d() (to be implemented in src/tnfr/navier_stokes/operator.py), n ∈ {32, 48}, ν ∈ {0.01, 0.005}, T = 2.0. Expected verdict: STEEPER_THAN_K41 (CDC not satisfied at Re_eff ≤ 500).theory/TNFR_NAVIER_STOKES_RESEARCH_NOTES.md §20 (full lemmas, theorem, CDC gap analysis, verdict table, N17-B pre-registration spec).TNFRNavierStokesOperator methods on z-independent u = (u₀(x,y), u₁(x,y), 0)):
vorticity_3d: ω₀ = ω₁ = 0, ω₂ = ∂_x(v) − ∂_y(u)vortex_stretching_field: S_a = ω₀·∂_x(u_a) + ω₁·∂_y(u_a) + ω₂·∂_z(u_a) = ω₂·0 = 0 for all astretching_production: = 0.0 exactly in IEEE 754examples/85_navier_stokes_dimensional_asymmetry.py — z-independence → stretching_production ≈ 0 at machine precision across all tested configurations (commit 1fac358b).theory/TNFR_NAVIER_STOKES_RESEARCH_NOTES.md §19.origin/main:
a1f298fd — operator existence: , bounded self-adjoint orthogonal projection on badac156 — conservation + Lyapunov: projected Noether charge exactly conserved; energy monotone with Cesàro tail at rational 48b0574a — spectrum + final verdict: uniform spectral density ; Branch A confirmedAGENTS.md — new top-level section REMESH-∞ Closure: Catalog Completeness Theorem (N15, May 2026)theory/README.md — added REMESH_INFINITY_DERIVATION.md to canonical document maptheory/TNFR_NAVIER_STOKES_RESEARCH_NOTES.md §18.7 — N15 closure block with locked verdicts, B1/K41/RMT/B2/B3 ruled out, refined prediction P-W3-1 (temporal-only)theory/TNFR_RIEMANN_RESEARCH_NOTES.md §13septies.5 — structural identification of T-HP smooth/oscillatory split with of theory/STRUCTURAL_OPERATORS.md §4.3 — REMESH asymptotic limit note with operator definition, spectral density, and catalog-completeness consequenceStructuralIntegrityMonitor with complete postconditions for all 13 canonical operatorsvalidate_candidate(), filter_candidates(), suggest_alternative(), enforce_grammar_on_glyph()apply_glyph_with_grammar() for grammar enforcement before operator application_soft_grammar_prefilter() wired with grammar_dynamics for operator filteringis_safe() and summary()summary()tetrad(), fields(), conservation(), telemetry(), tensor_invariants(), emergent_fields(), evolve_grammar_aware(), integrity_check(), upgraded results() and info()_HAS_FIELDS, _HAS_CONSERVATION, _HAS_INTEGRITY, _HAS_GRAMMAR_DYNAMICStests/sdk/test_simple_advanced.pymake_ring_graph, make_node_data, ring3, ring5, small_graph)_make_graph helpersexcept: clauses in grammar_dynamics.py (now except Exception:)run_structural_validation aggregator (grammar U1-U3 + field thresholds Φ_s, |∇φ|, K_φ, ξ_C, optional ΔΦ_s drift).compute_structural_health with risk levels and recommendations.TelemetryEmitter integration example (examples/structural_health_demo.py).PerformanceRegistry, perf_guard, compare_overhead.scripts/structural_health_report.py (on-demand health summaries).docs/STRUCTURAL_HEALTH.md.perf_registry parameter in run_structural_validation (read-only timing).test_canonical_operator_set.Previous release (see repository history) with foundational operators, unified grammar, metrics, and canonical field tetrad.