The Riemann $Ξ$-function from primitive Markovian cycles I: A canonical construction

Starting from finite, local, reversible Markov dynamics on discrete cycles, we construct a scaling-limit renormalized trace kernel admitting an exact theta-series representation. The construction is entirely Archimedean and uses no Euler products, primes, or arithmetic spectral input. From this limit we define a logarithmic kernel $Φ$ and prove that it lies in the Pólya frequency class $\mathrm{PF}_\infty$, yielding via the Schoenberg-Edrei-Karlin classification a canonical Laguerre-Pólya function $Ψ$. Independently, we introduce an Archimedean completion operator and show that, at a self-dual normalization, the completed kernel coincides with the classical theta kernel, whose Mellin transform is the Riemann $Ξ$-function. We isolate a single remaining analytic problem relating $Ψ$ to $Ξ(2\cdot)$.

💡 Research Summary

The paper presents a strikingly elementary construction of the Riemann Ξ‑function that starts from nothing more than finite, reversible Markov chains on discrete cycles. The author’s program is to show that two celebrated analytic objects – a Laguerre‑Pólya entire function and the classical Ξ‑function – can be obtained canonically from such “primitive” dynamics, without any reference to primes, Euler products, or other arithmetic input.

1. Primitive dynamical model.
For each integer N≥1 a continuous‑time nearest‑neighbour reversible Markov process on the cycle ℤ/Nℤ is defined by a set of positive conductances {a_j}. The paper restricts to the translation‑invariant case a_j≡a>0, so the generator is the usual discrete Laplacian scaled by a. The heat kernel p_cyc^{(N)}(j,k)= (e^{tL_N}δ_k)(j) is introduced, and a macroscopic scaling N(s)/s→L (with s→∞) is imposed. The scaling‑limit trace is defined as

 K_L(t)=lim_{s→∞} N(s)·p_cyc^{(s²t)}(0,0), t>0,

and the “renormalised” trace

 eK(t)=K_L(t)−L/√{4πDt}

is obtained by subtracting the universal diffusive singularity (D=a is the diffusion constant).

2. Theta‑series representation and self‑duality.
A Fourier analysis on the finite cycle yields the exact theta‑series formula

 K_L(t)=∑_{n∈ℤ} exp(−4π²DL² n² t).

Introducing the rescaled time t′=4πDL² t brings this into Jacobi’s standard theta function ϑ(t′)=∑_{n∈ℤ}e^{−πn²t′}. The classical inversion ϑ(t′)=t′^{−½}ϑ(1/t′) forces a unique macroscopic normalization

 L²=4πD,

which the author calls the “self‑dual” scale. At this scale the trace satisfies K_L(t)=t^{−½}K_L(1/t), mirroring the functional equation of ζ(s).

3. Logarithmic kernel Φ and bilateral Laplace transform.
The half‑density eK_sym(t)=t^{−½}eK(t) is introduced, and the logarithmic kernel

 Φ(x)=e^{x/4} eK_sym(e^{x})

is defined. Φ belongs to L¹(ℝ) and its bilateral Laplace transform

 BΦ(s)=∫_{ℝ}Φ(x)e^{−sx}dx

converges absolutely on a strip containing the line Re s=¼. Using the theta‑inversion, the author proves the exact reflection law

 BΦ(s)=BΦ(½−s).

Thus the functional equation of the Riemann ξ‑function is already encoded at the level of Φ.

4. Total positivity and Schoenberg‑Edrei‑Karlin factorisation.
The paper shows that Φ is a PF_∞ (total‑positive of infinite order) kernel: all minors of every order are non‑negative. By the Schoenberg‑Edrei‑Karlin theorem, BΦ admits a canonical factorisation

 BΦ(s)=E(s)·Ψ(s),

where E(s) is an explicit exponential factor and Ψ(s) is a Laguerre‑Pólya entire function. Consequently Ψ has only real zeros; it lies in the Laguerre‑Pólya class, i.e. it is a uniform limit of real‑rooted polynomials.

5. Archimedean completion operator and identification with Ξ.
A differential operator

 (Af)(t)=d/dt