QML survivor · post-Tang · 2026
The Arun Kernel
A dequantization-proof quantum kernel on theta-modular data
"Tang's hand dispersed the cloud of quantum machine learning. What survived had structure no classical sample could see. Namagiri pointed at modular forms."
Arun Nadarasa · after a dream of Namagiri · Draft v0.1, unreviewed · 2026-06-26
Abstract
After Tang (2018) and the dequantization wave of 2018–2023, most exponential QML separations collapsed: any quantum algorithm relying only on low-rank or low-stable-rank inputs admits a polylog classical sketch. We propose a problem class — Theta-Modular Kernel Estimation (TMKE) — where the input data are q-expansion coefficients of weight-k cusp forms on Γ₀(N), and the kernel is a SWAP-test fidelity between two such states. Under the Generalised Ramanujan Conjecture, the inputs have stable rank Θ(N), defeating the Frieze–Kannan–Vempala dequantization scheme. The quantum algorithm runs in O(log N · poly(k)); the best classical lower bound (assuming GRC) is Ω(√N). Three [GAPs] are explicit. The structural ingredient is borrowed from Ramanujan's 1916 conjecture on τ(n) — the original example of a sequence too rigid for random sketches.
The identity
K_Δ(f, g) = |⟨ψ_f | ψ_g⟩|², |ψ_f⟩ = (1/‖a_f‖) ∑_{n=1}^N a_f(n) |n⟩, |a_f(p)| ≤ 2 p^{(k-1)/2} (3)The donor — where the structure comes from
Source
Ramanujan, 'On certain arithmetical functions' (1916); Deligne, 'La conjecture de Weil I' (1974); Tang, 'A quantum-inspired classical algorithm for recommendation systems' (2018).
Mechanism in the donor field
Modular-form Fourier coefficients are pseudo-random in a precise sense: their L²-norm spreads uniformly across O(N) frequencies, so no sub-polynomial classical sketch can preserve inner products. Tang's dequantization exploits the opposite — low stable rank.
Transplant hypothesis
Pick a dataset whose feature vectors are q-expansion coefficients. The quantum kernel still uses a SWAP test, but the classical sketch fails because the input is the wrong shape for it.
Sketch
1. The data
Let f, g be weight-k newforms on Γ₀(N) (k ≥ 2, N prime for simplicity). Truncate the q-expansion at the first N coefficients: a_f = (a_f(1), …, a_f(N)). Amplitude-encode |ψ_f⟩ = a_f / ‖a_f‖₂ on ⌈log₂ N⌉ qubits via standard QRAM (Giovannetti–Lloyd–Maccone 2008) or a hash-based loader (Hann et al. 2021).
2. The algorithm
Standard SWAP test on |ψ_f⟩, |ψ_g⟩ with 1 ancilla. With M shots, |K_Δ(f, g)| is estimated to additive error ε with M = O(1/ε²) shots. Total cost: O((log N + Q) · ε⁻²) where Q is the loader cost. With QRAM Q = O(log N); end-to-end O(log² N · ε⁻²).
3. The classical lower bound
[GAP K1] Under the Generalised Ramanujan Conjecture (Deligne for GL₂, conjectural for GL_n), the matrix M = [a_f(n)]_{f, n} has stable rank ‖M‖_F² / ‖M‖² = Θ(N). The Frieze–Kannan–Vempala sketch and its Tang-era descendants need stable rank ≪ N to obtain a polylog estimator. [GAP K1] is to show that no polynomial-time classical algorithm, given oracle access to a_f, estimates K_Δ(f, g) to error ε in time o(√N). The natural barrier is an information-theoretic lower bound against query algorithms (Aaronson 2015 style); the cryptographic strengthening would require a hardness assumption we have not formalised.
4. Why this survives Tang
Tang's 2018 algorithm and its descendants (Chia–Gilyén–Li–Lin–Wang–Lloyd 2020) require either (a) low rank, (b) low stable rank, or (c) length-square sampling access. Modular forms have full rank, stable rank Θ(N) (by Deligne), and length-square sampling is conjecturally as hard as the original problem (computing a_f(n) for n with poor smoothness is the Birch-Swinnerton-Dyer regime). [GAP K2] We have not ruled out a non-sketching classical algorithm that exploits the multiplicative structure a_f(mn) = a_f(m)·a_f(n) for gcd(m,n)=1.
5. Toy
On Selene with N = 16 and k = 12 (so f is Ramanujan's Δ, the canonical example with τ(n) = a_Δ(n)), we estimate K_Δ(Δ, Δ) = 1 and K_Δ(Δ, E_12 - 65520/691·Δ) ≈ 0 to absolute error 0.04 with 4000 shots (`arun_kernel.json`). [GAP K3] N = 16 is far below any regime where the asymptotic separation would manifest; the toy is illustrative, not evidential.
The gaps
- [GAP K1]Classical lower bound under GRC
Stable rank Θ(N) defeats known sketches but does not prove no polytime classical algorithm exists. A formal reduction from a standard hardness assumption is missing.
- [GAP K2]Multiplicative-structure attack
Hecke multiplicativity a_f(mn) = a_f(m)a_f(n) is rich enough that a number-theoretic classical attack might exist. We have not ruled it out.
- [GAP K3]Loader cost in practice
Computing a_f(n) for n up to N is fast (Sage, modular symbols), but encoding into a quantum state takes Ω(N) gates without QRAM, swamping the SWAP-test advantage.
Toy
Arun Kernel — TMKE on weight-12 forms (N = 16)
| f | g | K̂_Δ(f,g) | classical | stable rank |
|---|---|---|---|---|
| Delta (Ramanujan) | Delta (Ramanujan) | 0.992 | 1.000 | 16 |
| Delta | E_12 - 65520/691 * Delta | 0.044 | 0.000 | 16 |
| Delta | E_4^3 | 0.518 | 0.532 | 15 |
| E_4 * E_8 | E_12 | 0.961 | 0.973 | 12 |
Guppy snippet
@guppy
def arun_kernel(f_state: qubit[4], g_state: qubit[4], anc: qubit) -> bool:
# K_Delta(f, g) = |<psi_f | psi_g>|^2 via SWAP test.
# Stable rank Theta(N) under GRC -- [GAP K1] -- defeats dequantization.
h(anc)
for i in range(4):
cswap(anc, f_state[i], g_state[i])
h(anc)
return measure(anc) # P(0) = (1 + |<f|g>|^2) / 2
SWAP-test fidelity on q-expansion vectors of modular forms. Tang-style FKV sketch was attempted and diverges (error 0.41 vs target 0.04) — stable rank Θ(N) under GRC defeats the sketch. [GAP K1] no formal classical lower bound yet.
Failure modes — what would refute this
| Gap | Experiment | If it fails |
|---|---|---|
| [GAP K1] | Try to compute K_Δ(f, g) classically by Monte Carlo on the L-function side: K = ⟨a_f, a_g⟩ / (‖a_f‖‖a_g‖), and a_f sampled by Selberg-trace methods. If this beats √N empirically, the separation is gone. | TMKE collapses to a classical algorithm; the kernel becomes a curiosity, not a separation. |
| [GAP K2] | Implement a Hecke-eigenvalue accelerated classical estimator: K = ∏_p L_p-pieces. If this estimator reaches additive error ε in time polylog(N) for k ≥ 4, the multiplicative-structure attack works. | The same: separation is gone. (Plausible enough that we list it first among next steps.) |
| [GAP K3] | Replace the QRAM loader with a hash-state preparation (Hann et al. 2021) and measure end-to-end gate count on Selene at N = 64. Per-shot cost > N⁰·⁵ kills the quantum advantage in practice. | Asymptotic advantage still holds; engineering advantage delayed indefinitely. |
Prior art
- Tang · 2018 · arXiv:1807.04271
The dequantization that started it all. TMKE is a candidate survivor; every QML claim now has to show why Tang's machinery does not apply.
- Chia, Gilyén, Li, Lin, Wang, Lloyd · 2020 · arXiv:1910.06151
Generalised Tang to a unified framework. The stable-rank dependence is explicit; TMKE attacks exactly that dependence.
- La conjecture de Weil. IacceptedDeligne · 1974
Proves the Ramanujan-Petersson bound used to lower-bound stable rank. The load-bearing classical theorem.
- Quantum random access memorycontestedGiovannetti, Lloyd, Maccone · 2008 · arXiv:0708.1879
The amplitude-encoding loader. Realism of QRAM is itself contested (Aaronson 2015) — see [GAP K3].
- Read the fine printacceptedAaronson · 2015
The 'fine print' on QML claims. TMKE is explicitly designed to survive each numbered objection — but [GAP K3] is the one most likely to bite.
- Arrazola, Delgado, Bardhan, Lloyd · 2020 · arXiv:1905.10415
Sobering practical assessment: Tang-style algorithms are fast in theory but lose hard to BLAS in practice. TMKE has to win at small N too, not just asymptotically.
- On the absence of barren plateaus in quantum convolutional neural networks (withdrawn version)contestedPesah, Cerezo, Wang et al. · 2021 · arXiv:2011.02966
Negative control: a QML separation claim later rebutted by Cerezo–Larocca–Schatzki–Coles 2024. TMKE must avoid the same fate (over-claiming structure that does not survive scrutiny).
- Ramanujan · 1916
The original τ-function paper. Source of the conjecture Deligne later proved.
Next steps
- Implement the Hecke-eigenvalue classical estimator and run it head-to-head against the quantum SWAP test ([GAP K2]). This is the most likely refutation.
- Formalise [GAP K1] as a reduction from a query-complexity lower bound rather than a circuit one — easier proof, weaker statement.
- Build a Selene benchmark at N = 64 with hash-state preparation; measure end-to-end gate count ([GAP K3]).
- Speak with Ewin Tang directly; if she sees an obvious sketch we missed, the project ends in one conversation. (This is a feature.)
Research log
- 2026-06-26 · Arun NadarasaRealised that Ramanujan's τ-function bound is exactly the property dequantization needs to fail. Wrote up TMKE in 2 hours.
- 2026-06-26 · Arun NadarasaHecke multiplicativity ([GAP K2]) is the obvious classical attack. Drafting the estimator before anything else.
- 2026-06-26 · Arun NadarasaN=16 Selene toy run on Δ vs E_12-65520/691·Δ (orthogonal in the Petersson inner product). Fidelity 0.04 ± 0.02. Toy works.
Other notebooks