Richard Hoekstra Papers · Atlas · Data

Confinement from the conductor: the defect spectrum of a singular Hecke order

The curve and its Hecke node

Let N=163N = 163 and let E:y2+y=x32x+1E\colon y^2 + y = x^3 - 2x + 1 be the unique elliptic curve of conductor NN over \mathbb{Q}. It has rank 11, trivial torsion, $|\text{\textcyr{Sh}}| = 1$, jj-invariant 963/163-96^3/163, and generator P=(1,0)P = (1, 0) with canonical height ĥ(P)=0.1899\hat{h}(P) = 0.1899\ldots, satisfying L(E,1)/(Ωĥ(P))=1L'(E,1)/(\Omega \cdot \hat{h}(P)) = 1 (BSD, verified numerically).

Write f=anqnf = \sum a_n q^n for the weight-22 newform attached to EE by modularity. The characteristic polynomial of T2T_2 on S2(Γ0(163))S_2(\Gamma_0(163)) factors over \mathbb{Q} as xg5(x)g7(x)x \cdot g_5(x) \cdot g_7(x) with g5,g7g_5, g_7 irreducible of degrees 55 and 77. The root x=0x = 0 corresponds to ff (since a2(E)=0a_2(E) = 0).

Over 3\mathbb{Q}_3, the factor g5g_5 acquires one linear root α5\alpha_5 with $\alpha_5 \equiv 12 \pmod{27}$, and g7g_7 acquires one linear root α7\alpha_7 with $\alpha_7 \equiv 3 \pmod{27}$. Together with αf=0\alpha_f = 0, these three roots give eigenforms fV1=ff_{V_1} = f, fV5f_{V_5}, fV7f_{V_7} spanning the singular 33-adic Hecke factor.

The singular Hecke order is the rank-33 subring 𝒪sing=313T33T1133\mathcal{O}_{\mathrm{sing}} = \mathbb{Z}_3 \cdot 1 \oplus \mathbb{Z}_3 \cdot T_3 \oplus \mathbb{Z}_3 \cdot T_{11} \;\subset\; \mathbb{Z}_3^3 where the embedding into 33=3V7×3V1×3V5\mathbb{Z}_3^3 = \mathbb{Z}_3^{V_7} \times \mathbb{Z}_3^{V_1} \times \mathbb{Z}_3^{V_5} sends each operator to its eigenvalue triple. The pair factor is the rank-22 projection $$R = \pi_{V_1,V_5}(\mathcal{O}_{\mathrm{sing}}) = \{(\sigma_1, \sigma_2) \in \mathbb{Z}_3^2 : \sigma_1 \equiv \sigma_2 \pmod{3}\} \;\cong\; \mathbb{Z}_3[\eta]/(\eta^2 - 3\eta),$$ with conductor 𝔠=3R\mathfrak{c}= 3R and normalisation R̃=3×3\tilde{R} = \mathbb{Z}_3 \times \mathbb{Z}_3. The conductor exact sequence is 0R3×3𝔽30.0 \to R \to \mathbb{Z}_3 \times \mathbb{Z}_3 \to \mathbb{F}_3 \to 0.

The congruence underlying the singularity is $a_p(f_{V_1}) \equiv a_p(f_{V_5}) \pmod{3}$ for all primes pp. The form fV1=ff_{V_1} = f is rational; the form fV5f_{V_5} is defined over the number field [x]/(g5(x))\mathbb{Q}[x]/(g_5(x)) of discriminant 6565765657.

The lattice action

Let Λ=(/L)d\Lambda = (\mathbb{Z}/L\mathbb{Z})^d be a dd-dimensional periodic lattice. Assign Hecke operators Tp1,,Tp2dT_{p_1}, \ldots, T_{p_{2d}} to the 2d2d lattice directions, choosing from {T3,T7,T11,T19,T43,T67}\{T_3, T_7, T_{11}, T_{19}, T_{43}, T_{67}\}.

The single-defect energy of sR/3kRs \in R/3^k R is E(s)=2ds2+pdirsTps2,E(s) = 2d \cdot \|s\|^2 + \sum_{p \in \mathrm{dirs}} \|T_p \cdot s\|^2, where (σ1,σ2)2=|σ1|M2+|σ2|M2\|(\sigma_1, \sigma_2)\|^2 = |\sigma_1|_M^2 + |\sigma_2|_M^2 with ||M|\cdot|_M the centered residue modulo M=3kM = 3^k.

In the normalisation, TpT_p acts as (σ1,σ2)(λ1(p)σ1,λ2(p)σ2)(\sigma_1, \sigma_2) \mapsto (\lambda_1(p)\,\sigma_1,\, \lambda_2(p)\,\sigma_2) where λ1(p)=ap(fV1)\lambda_1(p) = a_p(f_{V_1}) and λ2(p)=ap(fV5)\lambda_2(p) = a_p(f_{V_5}) are the eigenvalues on the two branches.

E(s)E(s) is independent of the assignment of operators to directions.

E(s)=2ds2+pTps2E(s) = 2d\,\|s\|^2 + \sum_{p} \|T_p \cdot s\|^2 sums over all six operators regardless of pairing.

The vacuum sector

The residue classes R/3RR/3R decompose into three sectors: 0=𝔠/3𝔠\mathcal{B}_0 = \mathfrak{c}/3\mathfrak{c} (vacuum) and 1,2\mathcal{B}_1, \mathcal{B}_2 (matter). Since λ1(p)=ap(E)\lambda_1(p) = a_p(E) \in \mathbb{Z} for all pp, elements of the conductor 𝔠=3R\mathfrak{c}= 3R have the form s=(3a,3b)s = (3a, 3b) with a,b3a, b \in \mathbb{Z}_3.

[thm:vac] For s=(3,0)𝔠s = (3, 0) \in \mathfrak{c} the defect energy is $$E_0 = 9\Bigl(2d + \sum_{p} a_p(E)^2\Bigr) = 9 \cdot 135 = 1215 = 5 \cdot 3^5.$$ This value is independent of kk for k4k \ge 4. More generally, E(3n,0)=9n135E(3^n, 0) = 9^n \cdot 135 for 1nk31 \le n \le k - 3.

E(3,0)=69+p|3λ1(p)|2E(3,0) = 6 \cdot 9 + \sum_p |3\lambda_1(p)|^2. Since λ1(p)=ap(E)\lambda_1(p) = a_p(E) are ordinary integers, |3ap|M=3|ap||3a_p|_M = 3|a_p| once M>6max|ap|M > 6\max|a_p|, which holds for k4k \ge 4. The values 3ap(E)3a_p(E) for p{3,7,11,19,43,67}p \in \{3,7,11,19,43,67\} are {0,6,18,18,21,6}\{0, 6, -18, -18, 21, -6\}, giving (3ap)2=1161\sum (3a_p)^2 = 1161 and E0=54+1161=1215E_0 = 54 + 1161 = 1215. The geometric tower follows by scaling.

The Hecke invariant is p{3,7,11,19,43,67}ap(E)2=0+4+36+36+49+4=129=343,\sum_{p \in \{3,7,11,19,43,67\}} a_p(E)^2 = 0 + 4 + 36 + 36 + 49 + 4 = 129 = 3 \cdot 43, factoring through the conductor prime 33 and the Heegner prime 4343.

Confinement

[thm:conf] For sR\𝔠s \in R \setminus \mathfrak{c}, write s=(1+3a,1+3b)s = (1 + 3a, 1 + 3b). Then the optimal energy satisfies Emin(k)=Θ(32k)E_{\min}(k) = \Theta(3^{2k}) as kk \to \infty.

The leading term of TpsT_p \cdot s is λ2(p)1\lambda_2(p) \cdot 1, which has |λ2(p)|3=O(3k)|\lambda_2(p)|_3 = O(3^k) since λ2(p)\lambda_2(p) is an algebraic integer of degree 55 over \mathbb{Q} (not rational). The free parameters (a,b)(a, b) provide 22 degrees of freedom to minimise 1212 squared terms. The system is overconstrained and the minimum is Θ(32k)\Theta(3^{2k}).

[thm:z2] The spectra of 1\mathcal{B}_1 and 2\mathcal{B}_2 are identical at every kk.

(σ1,σ2)(σ1,σ2)(\sigma_1,\sigma_2) \mapsto (-\sigma_1,-\sigma_2) sends 12\mathcal{B}_1 \leftrightarrow \mathcal{B}_2 and preserves EE. This realises the suspension Σ2=Id\Sigma^2 = \mathrm{Id} of Dsg(R)D_{\mathrm{sg}}(R).

The vacuum sector lifts to 3\mathbb{Z}_3 (the element (3,0)(3,0) is a fixed integer). The matter sector does not: the optimal element at precision kk changes at every kk. The conductor 𝔠=3R\mathfrak{c}= 3R is the sharp boundary: Hensel lifting converges on 𝔠\mathfrak{c} and fails on R\𝔠R \setminus \mathfrak{c}.

Two-defect interaction

[thm:two] Two vacuum excitations sa=sb=(3,0)s_a = s_b = (3, 0) at distance rr on Λ\Lambda:

  1. V(r)=0V(r) = 0 for r2r \ge 2. The defects are free.

  2. V(1)=72=832V(1) = -72 = -8 \cdot 3^2 for nearest neighbours, independent of the direction assignment.

For r2r \ge 2, no link has both endpoints at defect sites, so Epair=2E0E_{\mathrm{pair}} = 2E_0 and V=0V = 0. For r=1r = 1, the two shared directed links change from s2+Tps2\|s\|^2 + \|T_p s\|^2 to sTps2\|s - T_p s\|^2 and vice versa. The difference sums to V=72V = -72 over all six operators.

The vacuum algebra

In the normalisation, ring multiplication is componentwise: (σ1,σ2)(σ1,σ2)=(σ1σ1,σ2σ2)(\sigma_1, \sigma_2) \cdot (\sigma_1', \sigma_2') = (\sigma_1\sigma_1', \sigma_2\sigma_2').

The boundary class is preserved: products of vacuum elements are vacuum elements. No product of two vacuum excitations produces matter.

σ10\sigma_1 \equiv 0 and σ10\sigma_1' \equiv 0 implies $\sigma_1\sigma_1' \equiv 0 \pmod{3}$.

(3,0)(3,0)=(9,0)(3,0) \cdot (3,0) = (9, 0) with E(9,0)=9E(3,0)=10935E(9,0) = 9 \cdot E(3,0) = 10935, exceeding the matter threshold. Two light vacuum excitations produce a heavy vacuum excitation.

(3,0)(0,6)=(0,0)(3, 0) \cdot (0, 6) = (0, 0). A σ1\sigma_1-excitation and a σ2\sigma_2-excitation annihilate to the vacuum.

The boundary grammar

On the boundary R/3R𝔽3R/3R \cong \mathbb{F}_3, each Hecke operator acts as multiplication by ap(E)mod3a_p(E) \bmod 3: p3711194367apmod3020011\begin{array}{c|cccccc} p & 3 & 7 & 11 & 19 & 43 & 67 \\ \hline a_p \bmod 3 & 0 & 2 & 0 & 0 & 1 & 1 \end{array} Three operators kill (×0\times 0: T3,T11,T19T_3, T_{11}, T_{19}), one swaps (×2\times 2: T7T_7), and two are the identity (×1\times 1: T43,T67T_{43}, T_{67}).

On a plaquette with corners b00,b10,b01,b11𝔽3b_{00}, b_{10}, b_{01}, b_{11} \in \mathbb{F}_3, the boundary action cost S=links|bi(apmod3)bj|2S_{\partial} = \sum_{\mathrm{links}} |b_i - (a_p \bmod 3) \cdot b_j|^2 defines a weighted code on 𝔽34\mathbb{F}_3^4. The unique zero-cost pattern is the vacuum (0,0,0,0)(0,0,0,0). The cost function is invariant under the bulk precision kk: the boundary grammar is a property of the conductor, not the bulk.

The Hecke automaton

The Heegner point P=(1,0)E()P = (1,0) \in E(\mathbb{Q}) has formal logarithm z(P)33z(P) \in 3\mathbb{Z}_3 with $z(P)/3 \equiv 1 \pmod{3}$ (matter sector). In the formal group E1(3)3E_1(\mathbb{Q}_3) \cong \mathbb{Z}_3, the six Hecke operators act by multiplication by ap(E)a_p(E).

[thm:automaton] The action of the Hecke eigenvalues {0,2,6,6,7,2}\{0, 2, {-6}, {-6}, 7, {-2}\} on 3\mathbb{Z}_3 by multiplication defines a one-counter automaton with states (v,r)1×𝔽3×(v, r) \in \mathbb{Z}_{\ge 1} \times \mathbb{F}_3^\times and transitions: $$\begin{array}{ll} \textup{ANNIHILATE}\; (a_3 = 0): &(v, r) \to \bot \\ \textup{SWAP}\; (a_7 = 2): &(v, r) \to (v,\, 2r) \\ \textup{SHIFT}\; (a_{11} = a_{19} = -6): &(v, r) \to (v{+}1,\, r) \\ \textup{ID}\; (a_{43} \equiv a_{67} \equiv 1): &(v, r) \to (v, r) \end{array}$$ The counter vv is monotone: it can only increase. The reachable set from the Heegner state (1,1)(1, 1) is exactly 33=𝔠3\mathbb{Z}_3 = \mathfrak{c}, the conductor ideal. The matter sector 3×\mathbb{Z}_3^\times is unreachable.

The unit gates {2,7,2}\{2, 7, -2\} generate (/3k)×(\mathbb{Z}/3^k\mathbb{Z})^\times at every kk (verified for k5k \le 5), giving full permutation within each shell. The shift gate 6=23-6 = -2 \cdot 3 increments vv by 11 (the factor 2-2 is a unit, leaving rr unchanged modulo 33). The annihilator sends everything to \bot. No gate decrements vv, so the matter sector (v=0v = 0) is unreachable.

The Heegner point PE()P \in E(\mathbb{Q}) is a confined excitation of the Hecke lattice: it generates E()E(\mathbb{Q}) but lives in the matter sector of the defect spectrum, with energy diverging as Θ(32k)\Theta(3^{2k}).

Numerical verification

kk 3k3^k Evac,1E_{\mathrm{vac},1} EmatterE_{\mathrm{matter}} Δ\Delta
3 27 324 657 333
4 81 1215 3 087 1 872
5 243 1215 23 445 22 230
6 729 1215 185 877 184 662

The vacuum energy stabilises at k=4k = 4. The matter energy grows as Em0.3532kE_m \approx 0.35 \cdot 3^{2k}. The /2\mathbb{Z}/2\mathbb{Z} symmetry is exact at every kk.

nn E(3n,0)E(3^n, 0) E/9nE / 9^n
1 1 215 135
2 10 935 135
3 98 415 135

The vacuum tower is exactly geometric with constant 135=6+343135 = 6 + 3 \cdot 43.

The probabilistic automaton

With uniform gate selection (probability 1/61/6 each), the automaton becomes a Markov chain.

[thm:boltzmann] The steady-state depth distribution under uniform gate selection with continuous injection at v=1v = 1 is geometric: ρ(v)qv,q=SA+S,\rho(v) \propto q^v, \qquad q = \frac{S}{A + S}, where AA is the number of annihilator gates and SS the number of shift gates. The permutation count PP does not enter. For (A,S,P)=(1,2,3)(A, S, P) = (1, 2, 3): q=2/3q = 2/3.

Conditioned on survival (probability (S+P)/(A+S+P)(S{+}P)/(A{+}S{+}P) per step), each step shifts with probability S/(S+P)S/(S{+}P) and permutes otherwise. The depth increment is Bernoulli with parameter S/(S+P)S/(S{+}P). Combined with geometric killing (rate A/(A+S+P)A/(A{+}S{+}P)), the steady-state depth is geometric with ratio q=(survival)×(shift rate)/(1(survival)×(1-shift rate))=S/(A+S)q = \text{(survival)} \times \text{(shift rate)} / (1 - \text{(survival)} \times \text{(1-shift rate)}) = S/(A+S).

The mean depth at death is v=1+S/A\langle v \rangle = 1 + S/A. For (1,2,3)(1, 2, 3): v=3\langle v \rangle = 3, corresponding to the 33-shell ring R/33R=R/27RR/3^3 R = R/27R.

The emergent temperature is T=1/ln((A+S)/S)=1/ln(3/2)T = 1/\ln((A{+}S)/S) = 1/\ln(3/2) for the (1,2,3)(1,2,3) class. No temperature parameter is imposed: it follows from the gate statistics.

The light cone

The six Hecke gates classify the lattice directions: Timelike (permutation, same shell):x,±ySpacelike (shift, deeper shell):±zNull (annihilation):+x\begin{array}{ll} \text{Timelike (permutation, same shell):} & -x,\, \pm y \\ \text{Spacelike (shift, deeper shell):} & \pm z \\ \text{Null (annihilation):} & +x \end{array} A vacuum soliton (3,0)(3, 0) propagates freely in the three timelike directions, irreversibly deepens in the two spacelike directions, and dies in the null direction. The signature (3,2,1)(3, 2, 1) is determined by the distribution of v3(ap(E))v_3(a_p(E)) over the six Heegner primes. The chirality (+x+x null, x-x timelike) arises from the supersingularity a3(E)=0a_3(E) = 0.

Lattice verification

A single defect (3,0)(3, 0) on a 24324^3 lattice at mod729\bmod\,729 spreads with initial depth growth rate ΔvΔn|n=0=0.40=25=SS+P,\frac{\Delta\langle v \rangle}{\Delta n}\bigg|_{n=0} = 0.40 = \frac{2}{5} = \frac{S}{S+P}, matching the automaton exactly. The depth saturates at v3.5\langle v \rangle \approx 3.5 due to mod\bmod-MM wrapping.

Universality

Among elliptic curves with a3(E)=0a_3(E) = 0 and prime conductor N<1000N < 1000, five automaton classes arise:

(A,S,P)(A,S,P) qq v\langle v\rangle TT sig. curves
(2,0,4)(2,0,4) 00 11 \infty (4,0,2)(4,0,2) 17, 811
(1,0,5)(1,0,5) 00 11 \infty (5,0,1)(5,0,1) 73
(1,1,4)(1,1,4) 12\tfrac12 22 1/ln21/\!\ln 2 (4,1,1)(4,1,1) 109
(1,2,3)(1,2,3) 23\tfrac23 33 1/ln321/\!\ln\tfrac32 (3,2,1)(3,2,1) 163,179,197,269,739
(2,1,3)(2,1,3) 13\tfrac13 32\tfrac32 1/ln31/\!\ln 3 (3,1,2)(3,1,2) 307

Discussion

The lattice model is a sigma model with singular target. The target space Spec(R)\mathrm{Spec}(R) is a nodal curve: two copies of Spec(3)\mathrm{Spec}(\mathbb{Z}_3) meeting at a single point (the conductor). The field s:ΛSpec(R)s\colon \Lambda \to \mathrm{Spec}(R) assigns to each lattice site a section of the structure sheaf of the node.

The two branches are:

The conductor 𝔠=3R\mathfrak{c}= 3R is the singular point where the branches meet. The vacuum sector consists of sections supported at the node. The matter sectors consist of sections supported on a single branch. Confinement is the obstruction to extending a branch section to a global section of the nodal curve over the full 33-adic ring.

The chain of identifications is: $$e^{\pi\sqrt{163}} \approx \mathbb{Z} \;\longrightarrow\; h(-163) = 1 \;\longrightarrow\; E/\mathbb{Q} \;\longrightarrow\; f \equiv g \pmod{3} \;\longrightarrow\; R \;\longrightarrow\; \text{confinement}.$$

The five numbers

The model is determined by five arithmetic constants: $$\begin{aligned} \textstyle\sum a_p(E)^2 &= 129 = 3 \cdot 43 & &\text{(vacuum stiffness, Heegner prime)} \\ E_0 &= 1215 = 5 \cdot 3^5 & &\text{(vacuum energy)} \\ V(1) &= -72 = -8 \cdot 3^2 & &\text{(nearest-neighbour attraction)} \\ E_m/M^2 &\to 0.35\ldots & &\text{(confinement constant)} \\ \Sigma^2 &= \mathrm{Id} & &\text{($\mathbb{Z}/2\mathbb{Z}$ of $D_{\mathrm{sg}}(R)$)} \end{aligned}$$

The discrete Poincaré disc

The automaton state (v,r)(v, r) has a natural polar structure: v1v \ge 1 is the radius (conductor depth), r𝔽3×r \in \mathbb{F}_3^\times the angle (branch choice). The Boltzmann factor (2/3)v(2/3)^v defines a discrete metric ds2=(23)vdr2+dv2.ds^2 = \bigl(\tfrac{2}{3}\bigr)^v\,dr^2 + dv^2. The angular cost decreases exponentially with depth: swapping branches is cheap in the bulk, expensive near the boundary. Geodesics prefer shallow orbits, but the shift gates push the soliton deeper at rate 2/52/5 per step. The conductor is a repulsive force — the node at v=0v = 0 pushes outward.

In the 33-adic limit kk \to \infty: the soliton drifts to v=v = \infty, angular motion costs nothing, and the soliton becomes a pure radial ray. Confinement is radial escape on the Poincaré disc.

The soliton

On the lattice, the vacuum excitation (3,0)(3, 0) is immortal: it never annihilates (the gate T3T_3 with a3=0a_3 = 0 contributes zero to the additive lattice sum, unlike the multiplicative automaton where it absorbs). It oscillates with period 33 in amplitude and breathes with period 9=329 = 3^2 in spatial extent, contracting to 5{\sim}\,5 sites every 99 steps. The breathing period is the conductor squared.

The topological charge is K0(Dsg(R))=/2K_0(D_{\mathrm{sg}}(R)) = \mathbb{Z}/2\mathbb{Z}: the σ1\sigma_1-soliton (3,0)(3,0) carries charge [B0]=+1[B_0] = +1, the σ2\sigma_2-soliton (0,6)(0,6) carries [B3]=1[B_3] = -1, and ring multiplication gives annihilation: (3,0)(0,6)=(0,0)(3,0) \cdot (0,6) = (0,0).

The path integral

The amplitude to reach lattice site xx after nn steps is the sum over all gate-sequences (paths) of length nn from the origin to xx, weighted by the product of eigenvalues: G(x,n)=paths γ:0xi=1napi(E).G(x, n) = \sum_{\text{paths } \gamma: 0 \to x} \prod_{i=1}^{n} a_{p_i}(E). Each eigenvalue apa_p contributes a direction-dependent mass: m(p)=log|ap(E)|log3.m(p) = \frac{\log |a_p(E)|}{\log 3}. The six masses are:

direction +x+x x-x +y+y y-y +z+z z-z
apa_p 00 2-2 22 77 6-6 6-6
|ap||a_p| 00 22 22 77 66 66
v3(ap)v_3(a_p) \infty 00 00 00 11 11
gate A P P P S S
m(p)m(p) \infty 0.630.63 0.630.63 1.771.77 1.631.63 1.631.63

One table. Three readings. The gate type (v3v_3) determines the automaton. The mass (log|ap|\log|a_p|) determines the path integral. The eigenvalue itself (apa_p) determines the lattice amplitude. All three are the same six integers {0,2,6,6,7,2}\{0, 2, -6, -6, 7, -2\}, which are the Fourier coefficients of the elliptic curve E:y2+y=x32x+1E\colon y^2 + y = x^3 - 2x + 1 at the six Heegner primes.

Open questions

  1. Is ap(E)2=343\sum a_p(E)^2 = 3 \cdot 43 a special value of the symmetric square LL-function L(sym2E,s)L(\mathrm{sym}^2 E, s)?

  2. The Adèlic product over supersingular primes p{2,3,17}p \in \{2, 3, 17\} gives a multi-tape automaton with constraint v3v2v_3 \le v_2 (the 33-adic world nested inside the 22-adic world). Is this constraint universal or specific to E=𝟷𝟼𝟹𝚊𝟷E = \texttt{163a1}?

  3. The quantum automaton (Hamiltonian on the state space) has ground state at v=Vmax/2v = V_{\max}/2. The Heegner point multiple 43v1P4 \cdot 3^{v-1} \cdot P connects the quantum ground state to the formal group. Is this a Gross–Zagier relation?

  4. The path integral gives direction-dependent mass m(p)=log|ap|/log3m(p) = \log|a_p|/\log 3. Is the mass spectrum {0.63,0.63,1.63,1.63,1.77}\{0.63, 0.63, 1.63, 1.63, 1.77\} related to a spectral measure on the Hecke algebra?

  5. The gate lens on neural networks. Classifying the weights of a trained network (e.g. GPT-2) by magnitude into kill/shift/permute bins reveals depth-dependent structure: shallow layers kill more (AA decreasing from 0.240.24 to 0.120.12), deep layers shift more (SS increasing from 0.290.29 to 0.430.43). Is this a universal property of trained deep networks?

  6. The five automaton classes (A,S,P)(A, S, P) for curves with a3=0a_3 = 0 and prime conductor N<1000N < 1000: do they exhaust all possibilities, or are there further classes at larger NN?