Richard Hoekstra Papers · Atlas · Data

The drift rate of the Hecke automaton at conductor 163163

Setup and prior work

We use the Hecke automaton 𝒜163\mathcal A_{163} defined and studied in the companion paper . Briefly: EE is the elliptic curve of conductor 163163 and rank 11, with Mordell–Weil generator P=(0,0)P=(0,0), and f=anqnf=\sum a_n q^n its weight-22 newform. The six small Hecke operators have eigenvalues a3=0,a7=2,a11=6,a19=6,a43=7,a67=2.a_3=0,\quad a_7=2,\quad a_{11}=-6,\quad a_{19}=-6,\quad a_{43}=7,\quad a_{67}=-2. Their gate types on the formal-group quotient at the prime 33 are papv3(ap)u3(ap)mod3type30annihilator7202swap11611shift (+1)19611shift (+1)43701identity67201identity\begin{array}{c|cccc} p & a_p & v_3(a_p) & \mathrm{u}_3(a_p)\bmod 3 & \text{type} \\\hline 3 & 0 & - & - & \text{annihilator} \\ 7 & 2 & 0 & 2 & \text{swap} \\ 11 & -6 & 1 & 1 & \text{shift }(+1) \\ 19 & -6 & 1 & 1 & \text{shift }(+1) \\ 43 & 7 & 0 & 1 & \text{identity} \\ 67 & -2 & 0 & 1 & \text{identity} \end{array} The transition rule on a unit state (r,k)𝔽3××0(r,k)\in\mathbb{F}_3^\times\times\mathbb{Z}_{\geq 0} under a non-annihilator gate TpT_p is (r,k)(u3(ap)rmod3,k+v3(ap)).(r,k)\;\longmapsto\;\bigl(\mathrm{u}_3(a_p)\cdot r\bmod 3,\;k+v_3(a_p)\bigr). The shell coordinate kk is therefore a deterministic counter on each trajectory, additively driven by the values v3(ap)v_3(a_p) of the applied gates .

The drift theorem

Throughout this section, let X1,X2,X_1,X_2,\ldots be i.i.d. uniform random variables on the five-element set Σ*={T7,T11,T19,T43,T67}\Sigma^*=\{T_7,T_{11},T_{19},T_{43},T_{67}\} of non-annihilator gates of 𝒜163\mathcal A_{163}. Given an initial state (r0,k0)𝔽3××0(r_0,k_0)\in\mathbb{F}_3^\times\times\mathbb{Z}_{\geq 0}, define the random trajectory {(rn,kn)}n0\{(r_n,k_n)\}_{n\geq 0} by iterating the transition rule with XiX_i. Set ξi=v3(ap(Xi)){0,1},Kn=k0+i=1nξi.\xi_i \;=\; v_3\!\bigl(a_{p(X_i)}\bigr) \;\in\;\{0,1\}, \qquad K_n \;=\; k_0 + \sum_{i=1}^n \xi_i.

[lem:single] Each ξi\xi_i is a Bernoulli random variable with parameter p=2/5p=2/5.

Reading off the gate table, exactly two of the five gates in Σ*\Sigma^* have v3(ap)=1v_3(a_p)=1 (namely T11T_{11} and T19T_{19}), and the remaining three have v3(ap)=0v_3(a_p)=0. Since XiX_i is uniform on Σ*\Sigma^*, ξi=1\xi_i=1 with probability 2/52/5 and ξi=0\xi_i=0 with probability 3/53/5.

[thm:drift] For every initial state (r0,k0)(r_0,k_0) and every n1n\geq 1, 𝔼[Knk0]=2n5,Var[Knk0]=6n25.\mathbb{E}[K_n-k_0] \;=\; \frac{2n}{5}, \qquad \mathrm{Var}[K_n-k_0] \;=\; \frac{6n}{25}. The rate μ=2/5\mu=2/5 is independent of the state and of the residue coordinate.

By Lemma [lem:single], ξ1,,ξn\xi_1,\ldots,\xi_n are i.i.d. Bin(1,2/5)\mathrm{Bin}(1,2/5). Independence of state follows from the fact that the increment ξi\xi_i depends only on the gate XiX_i, not on (ri1,ki1)(r_{i-1},k_{i-1}). Linearity of expectation gives 𝔼[Knk0]=n2/5\mathbb{E}[K_n-k_0]=n\cdot 2/5, and independence gives Var[Knk0]=n(2/5)(3/5)=6n/25\mathrm{Var}[K_n-k_0]=n\cdot(2/5)(3/5)=6n/25.

[thm:law] For every n1n\geq 1, the random variable Knk0K_n-k_0 is binomial: Knk0Bin(n,2/5).K_n-k_0 \;\sim\; \mathrm{Bin}(n,2/5). In particular, the moment generating function is 𝔼[etKn]=etk0(35+25et)n\mathbb{E}[\,e^{tK_n}\,]=e^{tk_0}\bigl(\tfrac{3}{5}+\tfrac{2}{5}e^t\bigr)^{\!n}.

Knk0K_n-k_0 is a sum of nn i.i.d. Bernoulli(2/5)(2/5) variables by Lemma [lem:single].

[thm:hoeff] For every ϵ>0\epsilon>0 and every n1n\geq 1, $$\Pr\!\Bigl[\,\bigl|K_n-k_0-\tfrac{2n}{5}\bigr|\geq \epsilon n\,\Bigr] \;\leq\; 2\exp\!\bigl(-2n\epsilon^2\bigr).$$ Equivalently, Kn/n2/5K_n/n\to 2/5 a.s. exponentially fast in nn.

Hoeffding’s inequality applied to the bounded i.i.d. sum Knk0=ξiK_n-k_0=\sum\xi_i with ξi[0,1]\xi_i\in[0,1].

[cor:clt] $\displaystyle \frac{K_n-k_0-\frac{2n}{5}}{\sqrt{6n/25}} \;\xrightarrow{d}\; \mathcal N(0,1).$

Where the constant 2/52/5 comes from

[prop:formula] Let SS be any finite set of gate primes for an elliptic newform ff, and let S*SS^*\subseteq S be the subset on which ap0a_p\neq 0. Then under uniform i.i.d. sampling from S*S^*, the mean shift increment of the associated Hecke automaton is μ(S)=1|S*|pS*v3(ap).\mu(S) \;=\; \frac{1}{|S^*|}\sum_{p\in S^*} v_3(a_p).

Linearity of expectation, exactly as in the proof of Theorem [thm:drift].

For E=163a1E=163a1 with S={3,7,11,19,43,67}S=\{3,7,11,19,43,67\}, we have |S*|=5|S^*|=5 and pS*v3(ap)=0+1+1+0+0=2\sum_{p\in S^*}v_3(a_p)=0+1+1+0+0=2, giving μ(S)=2/5\mu(S)=2/5.

The constant μ(S)\mu(S) is independent of the residue coordinate, of the initial state, of which permutation gate appears (swap vs. identity), and of the Atkin–Lehner sign of the level. It depends only on the multiset {v3(ap)}\{v_3(a_p)\} for the chosen gate set.

Connection to the trace inequality

In  the trace Gram matrix G(d,d)=Tr(TdTdS2(Γ0(163)))G(d,d')=\mathop{\mathrm{Tr}}(T_dT_{d'}\mid S_2(\Gamma_0(163))) on the seven-element Heegner set S={3,7,11,19,43,67,163}S=\{3,7,11,19,43,67,163\} was shown to satisfy |G(d,d)|2min(G(d,d),G(d,d)).|G(d,d')| \;\leq\; 2\min\bigl(G(d,d),G(d',d')\bigr). We record one quantitative consequence which the present paper makes visible.

[prop:trace-asymp] Let S*={7,11,19,43,67}SS^*=\{7,11,19,43,67\}\subset S be the non-annihilator subset. The arithmetic mean of Tr(Tp2S2(Γ0(163)))/p\mathop{\mathrm{Tr}}(T_p^2\mid S_2(\Gamma_0(163)))/p over pS*p\in S^* is bounded above and below by absolute constants independent of which Heegner prime pp one picks; numerically, with the data of , minpS*G(p,p)p=640679.55,maxpS*G(p,p)p=2821914.84.\min_{p\in S^*}\frac{G(p,p)}{p}\;=\;\frac{640}{67}\approx 9.55, \qquad \max_{p\in S^*}\frac{G(p,p)}{p}\;=\;\frac{282}{19}\approx 14.84. The diagonals are roughly linear in pp across the gate set (within a factor of 1.55\approx 1.55), while the off-diagonals are bounded by twice the smaller diagonal — a nontrivial constraint compatible with, but not implied by, the linear growth on the diagonal.

Using G(7,7)=92G(7,7)=92, G(11,11)=136G(11,11)=136, G(19,19)=282G(19,19)=282, G(43,43)=583G(43,43)=583, G(67,67)=640G(67,67)=640 from , the ratios G(p,p)/pG(p,p)/p are 92/713.1492/7\approx 13.14, 136/1112.36136/11\approx 12.36, 282/1914.84282/19\approx 14.84, 583/4313.56583/43\approx 13.56, 640/679.55640/67\approx 9.55. Direct evaluation.

We do not claim a closed form for the diagonal growth, nor a conceptual link between the drift rate μ=2/5\mu=2/5 and the trace constants. The two phenomena live on different objects: the drift is a property of the formal-group automaton, while the trace inequality is a property of S2(Γ0(163))S_2(\Gamma_0(163)). Whether the gate-level identity pv3(ap)=2\sum_p v_3(a_p)=2 controls anything about the modular trace remains open.

Family-level statement

[prop:family] For any imaginary quadratic field of class number one with prime discriminant N-N, N{3,7,11,19,43,67,163}N\in\{3,7,11,19,43,67,163\}, and any rational elliptic newform fNf^N at level NN (when one exists), the same construction yields a Hecke automaton 𝒜N\mathcal A_N with a well-defined drift rate μ(N)=1|SN*|pSN*v3(ap(fN))\mu(N) \;=\; \frac{1}{|S^*_N|}\sum_{p\in S^*_N} v_3(a_p(f^N)) under uniform sampling from any chosen non-annihilator gate set SN*S^*_N. The map Nμ(N)N\mapsto\mu(N) is a finite, computable invariant of the level once a gate set has been chosen.

Proposition [prop:formula] applied at each level.

We do not evaluate μ(N)\mu(N) for N163N\neq 163 here. Doing so requires a choice of gate set SN*S^*_N at each level and the corresponding Sage computation of ap(fN)a_p(f^N). At levels N{3,7}N\in\{3,7\} the space S2(Γ0(N))S_2(\Gamma_0(N)) is zero so no fNf^N exists, and at N{11,19}N\in\{11,19\} the space is one-dimensional and the construction reduces trivially. The interesting cases are N{43,67,163}N\in\{43,67,163\}. We state Proposition [prop:family] only as an existence result; numerical evaluation across the family is deferred.

What is not proved here

[prop:weighted] Under sampling proportional to ap2a_p^2 on S*={7,11,19,43,67}S^*=\{7,11,19,43,67\}, with weights {4,36,36,49,4}\{4,36,36,49,4\} summing to 129129, the mean shift increment is μw=40+361+361+490+40129=72129=2443.\mu_w \;=\; \frac{4\cdot 0+36\cdot 1+36\cdot 1+49\cdot 0+4\cdot 0}{129} \;=\; \frac{72}{129} \;=\; \frac{24}{43}. Equivalently, μw=24/43\mu_w=24/43 where 4343 is the cofactor in pap2=129=343\sum_p a_p^2=129=3\cdot 43 from .

Weighted average of v3(ap)v_3(a_p) with weights ap2a_p^2.

The two natural drift rates are μunif=25,μap2=2443.\mu_{\mathrm{unif}} \;=\; \frac{2}{5}, \qquad \mu_{a_p^2} \;=\; \frac{24}{43}. Their ratio is (24/43)/(2/5)=60/43(24/43)/(2/5)=60/43. We record this without interpretation.

Reproducibility

The five-line verification of all numerical claims in this paper:

ap = {3:0, 7:2, 11:-6, 19:-6, 43:7, 67:-2}
nonzero = {p:a for p,a in ap.items() if a != 0}
v3 = lambda n: 0 if n%3 else 1+v3(n//3)
mu = sum(v3(abs(a)) for a in nonzero.values()) / len(nonzero)
assert mu == 2/5

9

R. Hoekstra, The Hecke automaton at conductor 163163: monotonicity, /2\mathbb{Z}/2\mathbb{Z} symmetry, and confinement on the formal-group quotient, companion paper.

R. Hoekstra, A trace inequality for Hecke operators at class-number-one Heegner primes of level 163163.