mtft 0.3.0__tar.gz

mtft-0.3.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Roger Tano
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

mtft-0.3.0/PKG-INFO

@@ -0,0 +1,170 @@
+Metadata-Version: 2.2
+Name: mtft
+Version: 0.3.0
+Summary: Modular Time Field Theory — core data structures, modular geometry, gauge-Higgs unification, dark sector, and information geometry
+Author: Roger Tano
+License: MIT
+Project-URL: Homepage, https://github.com/roger/mtft
+Project-URL: Documentation, https://github.com/roger/mtft#readme
+Keywords: physics,modular-forms,gauge-theory,dark-matter,information-geometry
+Classifier: Development Status :: 3 - Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: Topic :: Scientific/Engineering :: Physics
+Classifier: Programming Language :: Python :: 3
+Classifier: Programming Language :: Python :: 3.9
+Classifier: Programming Language :: Python :: 3.10
+Classifier: Programming Language :: Python :: 3.11
+Classifier: Programming Language :: Python :: 3.12
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: numpy>=1.22
+Provides-Extra: full
+Requires-Dist: scipy>=1.9; extra == "full"
+Requires-Dist: matplotlib>=3.5; extra == "full"
+Provides-Extra: viz
+Requires-Dist: matplotlib>=3.5; extra == "viz"
+Requires-Dist: plotly>=5.0; extra == "viz"
+Provides-Extra: dev
+Requires-Dist: pytest>=7.0; extra == "dev"
+Requires-Dist: pytest-cov; extra == "dev"
+Requires-Dist: ruff; extra == "dev"
+
+# MTFT — Modular Time Field Theory
+
+[![Python 3.9+](https://img.shields.io/badge/python-3.9+-blue.svg)](https://www.python.org/downloads/)
+[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](https://opensource.org/licenses/MIT)
+
+A Python library implementing the core data structures and computational tools for **Modular Time Field Theory (MTFT)**.
+
+## Overview
+
+MTFT proposes a modular time field **τ(x) = t_R / t_U** mapping spacetime to the upper half-plane **ℍ**, unifying dark matter, particle masses, and cosmic structure through:
+
+- **Modular symmetry** — SL(2,ℤ) acting on τ
+- **Gauge-Higgs unification** — Higgs = A_τ holonomy via the Hosotani mechanism
+- **τ-vortex dark matter** — logarithmic τ-profiles give flat rotation curves with no dark particles
+- **Information geometry bridge** — Fisher-Rao curvature of the logistic map connects chaos → geometry → dynamics
+
+### The Three-Layer Architecture
+
+```
+Counting   →  Logistic map orbits, figurate-number degeneracies
+Geometry   →  Fisher-Rao metric, Ricci curvature R_core
+Dynamics   →  SM masses, dark matter halos, decay rates, cosmology
+```
+
+Connected by the **spectral determinant identity**:
+
+```
+det(1 − q^{1/m} P_m) = η(τ)^{−1/m} θ₃(0,τ)^{1/m}
+```
+
+linking Fredholm determinants (chaos/RG) to CFT partition functions.
+
+## Installation
+
+```bash
+pip install -e .              # minimal (numpy only)
+pip install -e ".[full]"      # with scipy + matplotlib
+pip install -e ".[dev]"       # with pytest + ruff
+```
+
+## Quick Start
+
+```python
+import mtft
+
+# ── Modular forms ─────────────────────────────────────
+tau = 0.1 + 1.5j
+eta = mtft.dedekind_eta(tau)
+j   = mtft.forms.j_invariant(tau)
+print(f"η(τ) = {eta:.6f}")
+print(f"j(τ) = {j:.2f}")
+
+# Verify the spectral determinant identity
+result = mtft.forms.verify_spectral_identity(tau, m=2)
+print(f"Spectral identity relative error: {result['relative_error']:.2e}")
+
+# ── Hosotani mechanism ────────────────────────────────
+hp = mtft.HosotaniPotential(fermion_fraction=0.4, kappa_ew=0.05)
+theta0 = hp.find_vacuum()
+masses = hp.gauge_masses()
+print(f"Vacuum θ₀ = {theta0:.4f}")
+print(f"m_W = {masses['m_W']:.2f} GeV  (PDG: 80.37)")
+print(f"m_Z = {masses['m_Z']:.2f} GeV  (PDG: 91.19)")
+
+# ── Particle spectrum ─────────────────────────────────
+sm = mtft.StandardModel()
+top = sm.by_name("Top")
+print(f"Top quark: m = {top.mass_GeV} GeV, κ = {top.kappa}")
+
+# Full κ-hierarchy
+for name, kappa in sm.kappa_hierarchy():
+    print(f"  {name:15s}  κ = {kappa:.2e}")
+
+# ── Dark sector ───────────────────────────────────────
+import numpy as np
+halo = mtft.TauVortexHalo(A=1e20, r0=1e30)
+r = np.logspace(31, 35, 100)
+v = halo.v_circular(r)
+print(f"v_∞ = {halo.v_infinity:.4e} (flat rotation velocity)")
+
+# ── Information geometry ──────────────────────────────
+R = mtft.info_geometry.R_core()
+print(f"R_core (Feigenbaum) = {R:.4f}")
+
+# ── Cosmology ─────────────────────────────────────────
+cosmo = mtft.FriedmannMTFT(Omega_tau=0.25)
+hist = cosmo.expansion_history()
+```
+
+## Package Structure
+
+```
+mtft/
+├── constants.py       # PDG masses, SM parameters, physical constants
+├── modular.py         # τ-field, SL(2,ℤ), hyperbolic geometry
+├── forms.py           # Dedekind η, Jacobi θ₃, Eisenstein, spectral det
+├── hosotani.py        # Effective potential, vacuum finder, EWSB
+├── particles.py       # SM particle database with κ-couplings
+├── dark_sector.py     # τ-vortex halos, rotation curves, Tully-Fisher
+├── info_geometry.py   # Fisher-Rao metric, Ricci curvature, logistic map
+└── cosmology.py       # Modified Friedmann, perturbations, G(t) oscillation
+```
+
+## Key Equations
+
+| Component | Equation |
+|-----------|----------|
+| Modular time field | τ(x) = t_R(x)/t_U(x) ∈ ℍ |
+| Hyperbolic metric | ds² = (dx² + dy²)/y², K = −1 |
+| W mass | m_W = κ_EW \|sin θ_H\| / (2R_τ) |
+| Fermion mass | m_f = κ_f \|sin θ₀\| / R_τ |
+| τ-vortex density | ρ_τ(r) = A²/(2r²) |
+| Flat rotation | v²_∞ = 2πGA² |
+| Tully-Fisher | v⁴ ∝ M_BH ∝ M_baryonic |
+| R_core bridge | m_τ² = c_m \|R_core\| |
+
+## Testing
+
+```bash
+pytest tests/ -v
+```
+
+## License
+
+MIT
+
+## Citation
+
+If you use this package in research, please cite:
+
+```bibtex
+@software{mtft2025,
+  title  = {MTFT: Modular Time Field Theory Python Package},
+  author = {Roger},
+  year   = {2025},
+  url    = {https://github.com/roger/mtft}
+}
+```

mtft-0.3.0/README.md

@@ -0,0 +1,138 @@
+# MTFT — Modular Time Field Theory
+
+[![Python 3.9+](https://img.shields.io/badge/python-3.9+-blue.svg)](https://www.python.org/downloads/)
+[![License: MIT](https://img.shields.io/badge/License-MIT-yellow.svg)](https://opensource.org/licenses/MIT)
+
+A Python library implementing the core data structures and computational tools for **Modular Time Field Theory (MTFT)**.
+
+## Overview
+
+MTFT proposes a modular time field **τ(x) = t_R / t_U** mapping spacetime to the upper half-plane **ℍ**, unifying dark matter, particle masses, and cosmic structure through:
+
+- **Modular symmetry** — SL(2,ℤ) acting on τ
+- **Gauge-Higgs unification** — Higgs = A_τ holonomy via the Hosotani mechanism
+- **τ-vortex dark matter** — logarithmic τ-profiles give flat rotation curves with no dark particles
+- **Information geometry bridge** — Fisher-Rao curvature of the logistic map connects chaos → geometry → dynamics
+
+### The Three-Layer Architecture
+
+```
+Counting   →  Logistic map orbits, figurate-number degeneracies
+Geometry   →  Fisher-Rao metric, Ricci curvature R_core
+Dynamics   →  SM masses, dark matter halos, decay rates, cosmology
+```
+
+Connected by the **spectral determinant identity**:
+
+```
+det(1 − q^{1/m} P_m) = η(τ)^{−1/m} θ₃(0,τ)^{1/m}
+```
+
+linking Fredholm determinants (chaos/RG) to CFT partition functions.
+
+## Installation
+
+```bash
+pip install -e .              # minimal (numpy only)
+pip install -e ".[full]"      # with scipy + matplotlib
+pip install -e ".[dev]"       # with pytest + ruff
+```
+
+## Quick Start
+
+```python
+import mtft
+
+# ── Modular forms ─────────────────────────────────────
+tau = 0.1 + 1.5j
+eta = mtft.dedekind_eta(tau)
+j   = mtft.forms.j_invariant(tau)
+print(f"η(τ) = {eta:.6f}")
+print(f"j(τ) = {j:.2f}")
+
+# Verify the spectral determinant identity
+result = mtft.forms.verify_spectral_identity(tau, m=2)
+print(f"Spectral identity relative error: {result['relative_error']:.2e}")
+
+# ── Hosotani mechanism ────────────────────────────────
+hp = mtft.HosotaniPotential(fermion_fraction=0.4, kappa_ew=0.05)
+theta0 = hp.find_vacuum()
+masses = hp.gauge_masses()
+print(f"Vacuum θ₀ = {theta0:.4f}")
+print(f"m_W = {masses['m_W']:.2f} GeV  (PDG: 80.37)")
+print(f"m_Z = {masses['m_Z']:.2f} GeV  (PDG: 91.19)")
+
+# ── Particle spectrum ─────────────────────────────────
+sm = mtft.StandardModel()
+top = sm.by_name("Top")
+print(f"Top quark: m = {top.mass_GeV} GeV, κ = {top.kappa}")
+
+# Full κ-hierarchy
+for name, kappa in sm.kappa_hierarchy():
+    print(f"  {name:15s}  κ = {kappa:.2e}")
+
+# ── Dark sector ───────────────────────────────────────
+import numpy as np
+halo = mtft.TauVortexHalo(A=1e20, r0=1e30)
+r = np.logspace(31, 35, 100)
+v = halo.v_circular(r)
+print(f"v_∞ = {halo.v_infinity:.4e} (flat rotation velocity)")
+
+# ── Information geometry ──────────────────────────────
+R = mtft.info_geometry.R_core()
+print(f"R_core (Feigenbaum) = {R:.4f}")
+
+# ── Cosmology ─────────────────────────────────────────
+cosmo = mtft.FriedmannMTFT(Omega_tau=0.25)
+hist = cosmo.expansion_history()
+```
+
+## Package Structure
+
+```
+mtft/
+├── constants.py       # PDG masses, SM parameters, physical constants
+├── modular.py         # τ-field, SL(2,ℤ), hyperbolic geometry
+├── forms.py           # Dedekind η, Jacobi θ₃, Eisenstein, spectral det
+├── hosotani.py        # Effective potential, vacuum finder, EWSB
+├── particles.py       # SM particle database with κ-couplings
+├── dark_sector.py     # τ-vortex halos, rotation curves, Tully-Fisher
+├── info_geometry.py   # Fisher-Rao metric, Ricci curvature, logistic map
+└── cosmology.py       # Modified Friedmann, perturbations, G(t) oscillation
+```
+
+## Key Equations
+
+| Component | Equation |
+|-----------|----------|
+| Modular time field | τ(x) = t_R(x)/t_U(x) ∈ ℍ |
+| Hyperbolic metric | ds² = (dx² + dy²)/y², K = −1 |
+| W mass | m_W = κ_EW \|sin θ_H\| / (2R_τ) |
+| Fermion mass | m_f = κ_f \|sin θ₀\| / R_τ |
+| τ-vortex density | ρ_τ(r) = A²/(2r²) |
+| Flat rotation | v²_∞ = 2πGA² |
+| Tully-Fisher | v⁴ ∝ M_BH ∝ M_baryonic |
+| R_core bridge | m_τ² = c_m \|R_core\| |
+
+## Testing
+
+```bash
+pytest tests/ -v
+```
+
+## License
+
+MIT
+
+## Citation
+
+If you use this package in research, please cite:
+
+```bibtex
+@software{mtft2025,
+  title  = {MTFT: Modular Time Field Theory Python Package},
+  author = {Roger},
+  year   = {2025},
+  url    = {https://github.com/roger/mtft}
+}
+```

mtft-0.3.0/pyproject.toml

@@ -0,0 +1,44 @@
+[build-system]
+requires = ["setuptools>=68.0,<77", "wheel"]
+build-backend = "setuptools.build_meta"
+
+[project]
+name = "mtft"
+version = "0.3.0"
+description = "Modular Time Field Theory — core data structures, modular geometry, gauge-Higgs unification, dark sector, and information geometry"
+readme = "README.md"
+license = {text = "MIT"}
+requires-python = ">=3.9"
+authors = [
+    {name = "Roger Tano"},
+]
+classifiers = [
+    "Development Status :: 3 - Alpha",
+    "Intended Audience :: Science/Research",
+    "Topic :: Scientific/Engineering :: Physics",
+    "Programming Language :: Python :: 3",
+    "Programming Language :: Python :: 3.9",
+    "Programming Language :: Python :: 3.10",
+    "Programming Language :: Python :: 3.11",
+    "Programming Language :: Python :: 3.12",
+]
+keywords = ["physics", "modular-forms", "gauge-theory", "dark-matter", "information-geometry"]
+
+dependencies = [
+    "numpy>=1.22",
+]
+
+[project.optional-dependencies]
+full = ["scipy>=1.9", "matplotlib>=3.5"]
+viz = ["matplotlib>=3.5", "plotly>=5.0"]
+dev = ["pytest>=7.0", "pytest-cov", "ruff"]
+
+[project.urls]
+Homepage = "https://github.com/roger/mtft"
+Documentation = "https://github.com/roger/mtft#readme"
+
+[tool.setuptools.packages.find]
+where = ["src"]
+
+[tool.ruff]
+line-length = 100

mtft-0.3.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

mtft-0.3.0/src/mtft/__init__.py

@@ -0,0 +1,44 @@
+"""
+MTFT — Modular Time Field Theory  (v0.3.0)
+============================================
+
+16 modules covering the complete MTFT framework from arithmetic
+weights through lattice gauge theory and experimental verification.
+"""
+
+__version__ = "0.3.0"
+
+# ── Tier 0: Constants & Arithmetic ───────────────────────────
+from mtft.constants import (
+    SM, GAUGE, HIGGS, PDG, PhysicalConstants, CriticalDepths, MassRatios,
+    FEIGENBAUM_DELTA, FEIGENBAUM_ALPHA, T_INF, LAMBERT_OMEGA, EULER_GAMMA,
+    QUARKS, LEPTONS, ZETA_2, XI,
+)
+from mtft.arithmetic import (
+    weight, weight_array, damped_weights,
+    S1, stiffness_S, center_stiffness, mass_gap_stiffness,
+    su3_hessian_eigenvalues, find_confinement_depth,
+)
+
+# ── Tier 1: Modular Geometry & Forms ─────────────────────────
+from mtft.modular import TauField, sl2z_transform
+from mtft.forms import dedekind_eta, jacobi_theta3, eisenstein_E2k, verify_spectral_identity
+
+# ── Tier 2: Gauge-Higgs ──────────────────────────────────────
+from mtft.hosotani import HosotaniPotential
+from mtft.particles import Particle, StandardModel
+
+# ── Tier 3: Phenomenology ────────────────────────────────────
+from mtft.dark_sector import TauVortexHalo, rotation_curve, tully_fisher
+from mtft.info_geometry import fisher_rao_metric, ricci_scalar_logistic
+from mtft.cosmology import FriedmannMTFT
+
+# ── Tier 4: Lattice & Fermions ───────────────────────────────
+from mtft.lattice import LatticeConfig, MTFTAction, metropolis_sweep, avg_plaquette, avg_polyakov
+from mtft.x0_143 import tano_mass_predictions, koide_angle_tano, generation_count, JACOBIAN
+from mtft.koide import koide_ratio, koide_leptons, predict_tau_mass, koide_manifold_point
+from mtft.burning_ship import burning_ship_iterate, ANISOTROPIC
+
+# ── Tier 5: Dimensional Bridge & Decay ───────────────────────
+from mtft.dimensional_bridge import electron_mass_from_eta, charge_from_feigenbaum
+from mtft.decay import ModularDecay

mtft-0.3.0/src/mtft/arithmetic.py

@@ -0,0 +1,234 @@
+"""
+MTFT Arithmetic Weights and Holonomy Stiffness
+================================================
+
+The arithmetic weights wₙ = Σ_{d|n} (log d)/d are the atoms of MTFT.
+Their Dirichlet series satisfies W(s) = −ζ(s)ζ′(s+1).
+
+Key series at modular depth y:
+    aₙ(y)  = wₙ e^{−2πyn}      (damped weights)
+    S₁(y)  = Σ aₙ(y)           (partition function / stiffness action)
+    S(y)   = Σ n² aₙ(y)        (holonomy stiffness)
+    C_N(y) = Σ n² aₙ(y) cos(2πn/N)  (center-projected stiffness)
+    μ_N(y) = min_m Σ n² aₙ(y) (1 − cos(2πnm/N))  (mass gap)
+
+The fine-structure constant is α⁻¹ = 2/S₁(y_spec).
+The confinement lock is C_3(y_conf) = 0  (Hessian isotropy).
+
+Reference: Papers 5, 7, 9, 13, 24; Dictionary IV–V.
+"""
+
+from __future__ import annotations
+
+import math
+from functools import lru_cache
+
+import numpy as np
+
+
+# ── Arithmetic weight function ───────────────────────────────
+
+def weight(n: int) -> float:
+    """
+    Compute wₙ = Σ_{d|n} (log d)/d.
+
+    These are rigid number-theoretic objects — not fitted.
+    w₁ = 0, w₂ = (log 2)/2, w₃ = (log 3)/3, w₆ = w₂ + w₃ + (log 6)/6, ...
+    """
+    if n < 1:
+        raise ValueError("n must be ≥ 1")
+    total = 0.0
+    for d in range(1, n + 1):
+        if n % d == 0:
+            if d > 1:
+                total += math.log(d) / d
+    return total
+
+
+@lru_cache(maxsize=8192)
+def _weight_cached(n: int) -> float:
+    return weight(n)
+
+
+def weight_array(n_max: int) -> np.ndarray:
+    """Compute wₙ for n = 1, ..., n_max as a NumPy array."""
+    return np.array([_weight_cached(n) for n in range(1, n_max + 1)])
+
+
+# ── Euler-shifted weights (Paper 20) ─────────────────────────
+
+def weight_euler(n: int) -> float:
+    """
+    Euler-shifted weight: w_n^(e) = w_n + σ_{-1}(n).
+
+    The substitution log k → log(ek) activates w₁ = 0 → w₁^(e) = 1,
+    creating the electron as the vacuum's first excitation.
+    """
+    sigma_minus1 = sum(1.0 / d for d in range(1, n + 1) if n % d == 0)
+    return weight(n) + sigma_minus1
+
+
+# ── Damped coefficients ──────────────────────────────────────
+
+def damped_weight(n: int, y: float) -> float:
+    """aₙ(y) = wₙ exp(−2πyn)."""
+    return _weight_cached(n) * math.exp(-2.0 * math.pi * y * n)
+
+
+def damped_weights(y: float, n_max: int = 500) -> np.ndarray:
+    """Array of aₙ(y) for n = 1..n_max."""
+    ns = np.arange(1, n_max + 1)
+    ws = weight_array(n_max)
+    return ws * np.exp(-2.0 * math.pi * y * ns)
+
+
+# ── Stiffness series ─────────────────────────────────────────
+
+def S1(y: float, n_max: int = 500) -> float:
+    """
+    Partition function / stiffness action:
+        S₁(y) = Σ wₙ e^{−2πyn}
+
+    The fine-structure constant is α⁻¹ = 2/S₁(y_spec).
+    """
+    return float(np.sum(damped_weights(y, n_max)))
+
+
+def stiffness_S(y: float, n_max: int = 500) -> float:
+    """
+    Holonomy stiffness:
+        S(y) = Σ n² wₙ e^{−2πyn}
+    """
+    ns = np.arange(1, n_max + 1)
+    ws = weight_array(n_max)
+    return float(np.sum(ns**2 * ws * np.exp(-2 * math.pi * y * ns)))
+
+
+def center_stiffness(y: float, N: int = 3, n_max: int = 500) -> float:
+    """
+    Z_N center-projected stiffness:
+        C_N(y) = Σ n² wₙ e^{−2πyn} cos(2πn/N)
+
+    C_3(y_conf) = 0 is the confinement lock.
+    """
+    ns = np.arange(1, n_max + 1)
+    ws = weight_array(n_max)
+    cos_terms = np.cos(2.0 * math.pi * ns / N)
+    return float(np.sum(ns**2 * ws * np.exp(-2 * math.pi * y * ns) * cos_terms))
+
+
+def mass_gap_stiffness(y: float, N: int = 3, n_max: int = 500) -> float:
+    """
+    SU(N) mass gap stiffness (minimum over modes m):
+        μ_N(y) = min_{m=1..N-1} Σ n² wₙ e^{−2πyn} (1 − cos(2πnm/N))
+
+    Theorem: μ_N(y) > 0 unconditionally for all N≥2, y>0.
+    """
+    ns = np.arange(1, n_max + 1)
+    ws = weight_array(n_max)
+    damped = ns**2 * ws * np.exp(-2 * math.pi * y * ns)
+
+    mu_min = float("inf")
+    for m in range(1, N):
+        gap_terms = 1.0 - np.cos(2 * math.pi * ns * m / N)
+        mu_m = float(np.sum(damped * gap_terms))
+        mu_min = min(mu_min, mu_m)
+    return mu_min
+
+
+# ── SU(3) Hessian eigenvalues (Paper 5 §24.4, Bridge §4.3) ──
+
+def su3_hessian_eigenvalues(y: float, n_max: int = 500) -> tuple[float, float]:
+    """
+    SU(3) holonomy Hessian eigenvalues:
+        λ₁ = (S_N + 2C_N) / 3
+        λ₂ = (S_N − C_N) / 3
+
+    At y_conf ≈ 0.18174:  λ₁/λ₂ = 1.000001 (the spectral lock).
+    """
+    S = stiffness_S(y, n_max)
+    C = center_stiffness(y, N=3, n_max=n_max)
+    lam1 = (S + 2 * C) / 3.0
+    lam2 = (S - C) / 3.0
+    return lam1, lam2
+
+
+def find_confinement_depth(
+    y_min: float = 0.10, y_max: float = 0.25, tol: float = 1e-8, n_max: int = 500
+) -> float:
+    """
+    Find y_conf where C_3(y) = 0 (SU(3) Hessian isotropy).
+
+    Uses bisection on C_3(y).
+    """
+    # C_3 starts positive and crosses zero
+    c_lo = center_stiffness(y_min, N=3, n_max=n_max)
+    c_hi = center_stiffness(y_max, N=3, n_max=n_max)
+
+    if c_lo * c_hi > 0:
+        # Try a wider scan
+        for y_test in np.linspace(y_min, y_max, 200):
+            c_test = center_stiffness(y_test, N=3, n_max=n_max)
+            if c_test * c_lo < 0:
+                y_max = y_test
+                c_hi = c_test
+                break
+
+    lo, hi = y_min, y_max
+    for _ in range(100):
+        mid = (lo + hi) / 2.0
+        c_mid = center_stiffness(mid, N=3, n_max=n_max)
+        if abs(c_mid) < tol:
+            return mid
+        if c_mid * c_lo > 0:
+            lo = mid
+            c_lo = c_mid
+        else:
+            hi = mid
+    return (lo + hi) / 2.0
+
+
+# ── Lambert series ───────────────────────────────────────────
+
+def lambert_series(q: complex, n_max: int = 200) -> complex:
+    """
+    MTFT Lambert series: W(q) = Σ_{d≥1} (log d)/d · q^d/(1−q^d).
+
+    This is the q-series feeding the modular holonomy potential.
+    """
+    total = 0.0 + 0j
+    for d in range(2, n_max + 1):
+        qd = q ** d
+        if abs(1.0 - qd) < 1e-30:
+            continue
+        total += (math.log(d) / d) * qd / (1.0 - qd)
+    return total
+
+
+# ── Torque flux convergence (Paper 7, §9) ────────────────────
+
+def torque_partial_sum(N: int) -> float:
+    """
+    Partial sum of the normalised torque: (1/N²) Σ_{n=1}^{N} n·wₙ.
+
+    Converges to T∞ = −ζ′(2)/2 ≈ 0.4688.
+    """
+    total = 0.0
+    for n in range(1, N + 1):
+        total += n * _weight_cached(n)
+    return total / (N ** 2)
+
+
+# ── Dirichlet series identity check ──────────────────────────
+
+def dirichlet_F(s: complex, n_max: int = 500) -> complex:
+    """
+    Compute F(s) = Σ f(n)/n^s where f = σ * Λ.
+
+    Should equal −ζ(s−1)ζ′(s) (the master Dirichlet identity).
+    """
+    total = 0.0 + 0j
+    for n in range(1, n_max + 1):
+        f_n = sum(math.log(d) for d in range(1, n + 1) if n % d == 0)
+        total += f_n / n ** s
+    return total