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 offered as a visual analogy for iterated contraction producing global geometric order. The threshold g⁶ = 33 named here is not an independently established constant: it is χ(H*(X⁶)) = 33, the central open conjecture of Book 6's G6 Crystal chapter (g6-crystal.html, marked "em construção," AXLE Issue 6 unresolved). Citing it here as a crossed threshold overstates its status; treat it as a named placeholder value carried over from an unproved conjecture, not a result this page independently confirms.

Section 5 — Visual

Fractal Contraction Portrait — Syracuse Orbit Density

Each point in the plane encodes a starting residue class. Colour maps to the escape rate under iterated Syracuse steps before entering the {1, 2, 4} cycle — the dm³ limit cycle attractor. Brighter gold = faster contraction (small ρk); deep navy = slow approach near the critical threshold. The fractal boundary is where ρ(ℒk) ≈ ρ(PMk) — the spectral gap closes, and the two operators momentarily agree.

Drag the depth slider to control iteration count k. The hexagonal invariant g⁶ = 33 appears as a natural resonance depth.

33
−0.2863
Fast contraction — deep in dm³ basin
Critical boundary — spectral gap closes
Slow — near threshold κ*
What you're seeing: The fractal structure is the pre-image of the limit cycle under ℒk. At k = 33, the depth slider was set to match the numeral in Book 6's open χ(H*(X⁶)) = 33 conjecture, so a resonance-looking pattern near k = 33 is expected by construction of the slider default, not a discovered coincidence. No derivation connects the Syracuse fold operator's spectral behavior to that conjecture; this is a shared numeral across two unrelated constructions, not a proven bridge.
Section 6 — Research Bridge

Bridge 0 — Numerical Approximation, Not Closed

Correction: this section previously read "Bridge 0 Closed," which contradicted both the "three key open items" listed above (the third of which is exactly this λ↔μmax connection) and this section's own footnote below, which already said "exact bridge is open." It is not closed. What follows is a numerical coincidence check, not a proof.

The Lyapunov exponent λ ≈ −0.2863 of the Syracuse map and the dm³ transverse eigenvalue bound μmax = −2, combined with g⁶ = 33 (itself an unproved conjecture value, see the caveat above), give:

λ · (g⁶ · log 2) ≈ −0.2863 × 33 × 0.693 ≈ −6.55 μmax · T* = −2 × (2π) ≈ −12.57 Ratio: μmax · T* / (λ · g⁶ · log 2) ≈ 1.92 Compare to 2 − ε* = 2 − 1/3 = 5/3 ≈ 1.667 — a ~15% gap, not a match. No derivation is offered for why this ratio should equal 2 − ε* in the first place; the resemblance is a loose numerical proximity, not an identity.

The canvas below shows both decay curves together with the dm³ contact normal form envelope. The green band is the region where both contractions agree to within ε* = 1/3 — the structural stability radius[Ch 10]. This illustrates the candidate bridge; it does not establish it.

ρ(PMk) — local 2-adic contraction
ρ(ℒk) = ekλ — global Ruelle
ek · μ_max · T*/g6 — dm³ contact envelope
ε* = 1/3 stability band (Bridge 0)
Open item for AXLE: Prove that the ratio μmax · T* / (λ · g⁶ · log 2) equals exactly 2 − ε* = 5/3, which would establish Bridge 0 as a theorem rather than a numerical observation. This connects the Collatz Lyapunov exponent to the dm³ canonical invariant triple (T* = 2π, μmax = −2, τ = 2) via the crystal law alone — no additional assumptions.
G6 LLC  ·  g6llc@proton.me  ·  +1 (646) 342-3751