Book 3 · Collatz Supplement · Spectral Theory

Spectral Radius Asymptotics
of the Transfer Operator

Single and Higher Iterates — dm³ + Crystal-Geometry Framework

The Syracuse return map T admits two natural finite approximations: the probability matrix PM and the Ruelle–Perron–Frobenius operator ℒ. Their spectral radii behave fundamentally differently. This page provides exact formulas, numerical tables, and an interactive widget.

Probability Matrix PM (Finite Approximation)

The probability matrix PM is the row-stochastic scaling of the Syracuse transition structure. At modulus 2M there are exactly 2M−1 odd residue classes, and each row has exactly one entry equal to 1 / 2M−1, making PM a scaled row-monomial matrix (each row has exactly one nonzero). Since n = 1 is a fixed point of T mod 2M for all M, the matrix has eigenvalue 1/2M−1 and thus ρ(PM) = 1/2M−1 exactly.

ρ(PM) = 1 / 2M−1 (exact, for all M) ρ(PM) ~ 2−M (exponential decay as M → ∞) ρ(PM+1) = ½ · ρ(PM) (halves at each modulus doubling) k-fold iterate: ρ(PMk) = [ρ(PM)]k = 1 / 2k(M−1)

Numerical Table — ρ(PM) and Higher Iterates

Values computed exactly as 1 / 2M−1. Columns k = 1, 2, 3, 5, 10 show the k-fold iterate radius.

M Odd states (2M−1) ρ(PM) = k=1 k = 2 k = 3 k = 5 k = 10
Interpretation (dm³): The scaling ρ(PM) = 2−(M−1) shows that local 2-adic dissipation becomes arbitrarily strong at finer resolution. At modulus 2M the probability mass is "forgotten" at a rate that halves with every extra bit of 2-adic precision — pure exponential erasure of initial conditions.

True Ruelle Transfer Operator ℒ

The genuine Ruelle–Perron–Frobenius transfer operator for the Syracuse map is

(ℒf)(x) = ΣT(y)=x f(y) / |T′(y)| |T′(y)| = 3 / 2v₂(3y+1)

where v₂(n) is the 2-adic valuation of n. Its spectral radius equals the expected Jacobian magnitude — equivalently, the exponential of the Lyapunov exponent λ:

ρ(ℒ) = E[|T′(n)|] = eλ λ ≈ −0.286 288 (from mod-243 weighted drift) ρ(ℒ) ≈ e−0.286 288 ≈ 0.7510 k-fold iterate: ρ(ℒk) = [ρ(ℒ)]k ≈ (0.751)k ~ ekλ Finite-M approximation error: O(2−M)

Convergence Table — ρ(ℒM) vs. M

The finite-M approximation ρ(ℒM) converges to eλ ≈ 0.7510 from above, with error decaying as O(2−M) because the 2-adic measure becomes uniform in the limit.

M ρ(ℒM) approx. Error |ρ − eλ| O(2−M) bound ρ(ℒ2) ρ(ℒ5) ρ(ℒ10)

Spectral Gap for Higher Iterates

The spectral gap of ℒk (difference between leading and sub-leading eigenvalue moduli) satisfies:

gapk = ρ(ℒk) − |λ₂(k)| = ρ(ℒ)k · (1 − |λ₂/ρ(ℒ)|k) As k → ∞: gapk ~ ρ(ℒ)k (gap approaches full leading eigenvalue — relative gap → 1) Finite-matrix case (PM): gapk = ρ(PMk) because all other eigenvalues of PM are zero.
Global vs. Local: The true operator radius ρ(ℒ) ≈ 0.751 < 1 is the global contraction rate after the triad fingerprint (c = 3) interacts with v₂. Higher iterates ρ(ℒk) ≈ (0.751)k make contraction exponentially stronger in k, mirroring the dm³ transverse eigenvalue bound μmax ≤ −2 applied over multiple contact cycles.

Interactive Spectral Radius Chart

Plot ρk against k for both operators. Adjust M and kmax to explore the asymptotics.

6
20
−0.2863
ρ(PMk) = (1/2M−1)k
ρ(ℒk) = ekλ
Asymptote eλ (k=1)
k ρ(PMk) ρ(ℒk) Ratio ρ(ℒk) / ρ(PMk)

Mod-128 Matrix (M = 7) — Explicit Construction

The next modulus beyond mod-64 is mod 128 (M = 7): 64 odd residue classes {1, 3, 5, …, 127}. The probability matrix P₇ is a 64 × 64 row-monomial matrix. Each row has exactly one nonzero entry equal to 1 / 26 = 1/64 ≈ 0.015625. Note that the Syracuse map mod 128 is not a bijection on odd residues: only 48 of the 64 odd residue classes appear as images (16 are unreachable; some are multiply-targeted). The spectral radius is still exactly 1/64 since n = 1 is a fixed point (T(1) = 1 mod 128).

Sparsity Pattern (64 × 64 block, first 16 rows)

★ marks the single nonzero in each row; · is zero. The target column is the index of the Syracuse image of each odd residue mod 128.

Asymptotics at M = 7

k ρ(P₇k) ρ(ℒk) ≈ ekλ gapk (ℒ)
Bridge 0 residual: At mod 128 (M = 7) the probability matrix radius ρ(P₇) = 1/64 ≈ 0.01563 is already 48× smaller than the Ruelle operator limit ρ(ℒ) ≈ 0.751. The two operators live on completely different scales, yet both contract: one superexponentially (PM), one at the fixed Lyapunov rate (ℒ). Closing Bridge 0 requires showing these two contractions are facets of the same geometric structure in the dm³ contact manifold.
Seed Sentence (CEFR A2)
"The bigger the power of 2 we look at, the smaller the biggest number in the matrix becomes. After many steps the shrinking becomes extremely fast."

Bridge 0 — Towards AXLE Publication

The exact asymptotics ρ(ℒk) ~ ekλ give a computable discrete analogue of the dm³ spectral gap applied over arbitrary operator cycles.

Publishing the full generating function of the spectrum of ℒk in AXLE would close a concrete piece of Bridge 0. The three key open items are:

The polar vortex / Saturn hexagon is the empirical certificate that the same iterated contraction already produces global geometric order once the monster threshold g⁶ = 33 is crossed.