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 ~ e 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) = e
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) ≈ e 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) ~ e 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) = e — 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