Hudba prvočísel · díl 1 ze 3
Prvočísla krotí kvantový chaos
Výprava přes 17 řádů výšky — od prvních nul Riemannovy zeta funkce až k nulám, jejichž pořadová čísla mají 22 cifer. Všude na ně narážíme na otisky prstů prvočísel, které žádná náhoda nenapodobí.
Představte si, že chaos uvnitř atomového jádra a dokonalý řád prvočísel tančí na stejném parketu. Zní to jako scénář ke sci-fi filmu, ale přesně tohle ukazuje naše nejnovější numerická studie.
Prošli jsme nuly Riemannovy zeta funkce v obrovském rozpětí: od prvních tisíců až po výšky kolem 1021, kde nuly nesou pořadová čísla s 22 ciframi. Otázka zněla jednoduše — chovají se nuly jako náhodný chaos, nebo si s sebou nesou paměť na prvočísla? Odpověď je jednoznačná: prvočísla zanechávají v rozložení nul otisky, které žádná náhoda nedokáže napodobit.
Zeta funkce ζ(s) = 1 + 1/2s + 1/3s + … je s prvočísly svázána tak těsně, že Bernhard Riemann už v roce 1859 ukázal, jak z jejích nul přesně „přehrát“ rozložení prvočísel. Netriviální nuly mají tvar ½ + i·γ — leží na takzvané kritické přímce (Riemannova hypotéza tvrdí, že všechny). Číslo γ je jejich „výška“ t. Prvních deset vidíte v záhlaví článku.
01Kvantový chaos a atomová jádra: proč fyziky zajímá čistá matematika?
Na začátku 70. let spočítal matematik Hugh Montgomery, jak často se dvě nuly zeta funkce vyskytují v dané vzdálenosti od sebe. Když vzorec ukázal fyzikovi Freemanu Dysonovi, ten v něm okamžitě poznal známou tvář: přesně tak vypadá statistika vlastních čísel velkých náhodných matic — konkrétně souboru GUE (Gaussian Unitary Ensemble).
Teorie náhodných matic se přitom proslavila popisem energetických hladin těžkých atomových jader, kde vládne kvantový chaos. Nuly zeta funkce se tedy na malých vzdálenostech chovají jako hladiny chaotického kvantového systému: navzájem se odpuzují a téměř nikdy se neocitnou těsně vedle sebe.
Zní to skvěle. Jenže má to háček.
Náhodné matice popisují lokální chaos — jak se nuly odpuzují zblízka. Jakmile se ale podíváme na dlouhé vzdálenosti a začneme počítat, kolik nul se vejde do širšího okna, teorie náhodných matic selhává. Podle ní by kolísání počtu nul mělo růst bez omezení. Ve skutečnosti se růst zastaví a nasytí.
Stojíte na mostě nad dálnicí a počítáte auta, která projedou úsekem dlouhým L. Řidiči drží odstup (to je ono „odpuzování“). Pokaždé napočítáte trochu jiný počet — a právě rozptyl tohoto počtu značíme Σ²(L). U náhodných matic roste s délkou úseku donekonečna. U nul zeta funkce se v jistou chvíli zarazí, jako by provoz řídil skrytý dirigent. Tím dirigentem jsou prvočísla.
02Berryho kouzlo: prvočísla jako oběžné dráhy
V roce 1988 přišel britský fyzik Sir Michael Berry s elegantní myšlenkou. Nuly zeta funkce si můžeme představit jako energetické hladiny (zatím neznámého) kvantového systému — a prvočísla jako jeho uzavřené periodické orbity. Každé prvočíslo p a jeho mocniny přispívají „orbitou“ s periodou ln pk.
Z této představy Berry odvodil semiklasický vzorec pro rozptyl počtu nul Σ²(L), který nemá žádné volné parametry. Nic se nedolaďuje, nic se nepřizpůsobuje datům. Vzorec předpovídá, že rozptyl bude nejprve oscilovat a pak se nasytí.
Studie ho systematicky otestovala na 18 datových blocích přes 17 řádů výšky:
- Přesnost: Berryho vzorec sedí na reálná data s odchylkou pouhých 0,26–1,21 %.
- Podpis prvočísla 2: v grafu rozptylu jsou vidět jemné vlnky. Vytváří je nejkratší prvočíselná orbita — prvočíslo 2 — s periodou ≈ H/ln 2. Ve výšce t ≈ 104 to vychází zhruba na 10,6 průměrné mezery (viz horní graf výše).
03Velký test: skutečná prvočísla vs. náhodný podvrh
Možná vás napadne námitka: „Není to celé jen tím, že prvočísel je určitý počet a jsou čím dál řidší?“ Abychom ji vyvrátili, postavili jsme Berryho vzorec proti Cramérovu modelu náhodných prvočísel — fiktivním „prvočíslům“ se stejnou hustotou jako ta pravá, ale s náhodnými polohami.
Poučení: na přesných polohách prvočísel záleží. Nuly zeta funkce si i ve výšce 1021 „pamatují“ konkrétní aritmetickou posloupnost 2, 3, 5, 7, 11, … — nestačí jim jen správný počet.
04Nasycení není konstantní: výška mění pravidla
Dlouho se tiše předpokládalo, že po vhodném přeškálování (takzvaném unfoldingu, kdy průměrnou mezeru mezi nulami nastavíme na 1) se nuly na všech výškách chovají stejně. Pro rozptyl na dlouhých vzdálenostech to neplatí.
Proč? Čím výš vystoupáte, tím víc prvočísel se vejde pod takzvaný Heisenbergův čas H = ln(t/2π) — tedy tím víc orbit má periodu ln p kratší než H. A každé další prvočíslo přidá kousek k výšce plató. Roste to ovšem neuvěřitelně pomalu, zhruba jako logaritmus logaritmu výšky: 17 řádů přidá jen asi dvě desetiny.
05Čtyřstupňové síto: jak poznat pravé spektrum zeta funkce?
Abychom poznatky podrobili zátěžovému testu, sestavili jsme čtyřstupňové síto — něco jako detektor lži pro kandidátní spektra:
- 1 · Hustota a tuhost: sedí celkový počet hladin a jejich stabilita?
- 2 · GUE mezery: chovají se sousední hladiny lokálně jako náhodné matice?
- 3 · Prvočíselný obsah: obsahuje spektrum oscilace prvočísel?
- 4 · Nasycení: zastaví se růst rozptylu na dlouhých vzdálenostech?
| Model spektra | 1 · Hustota a tuhost | 2 · GUE mezery | 3 · Prvočíselné oscilace | 4 · Nasycení |
|---|---|---|---|---|
| Čistá náhodná matice GUE žádná aritmetika | ✓ projde | ✓ projde | ✗ selže | ✗ selže |
| Skutečné nuly zeta funkce | ✓ projde | ✓ projde | ✓ projde | ✓ projde |
| Hladiny modelu xp + prvočíselné orbity | ✓ projde | ✓ projde | ✓ projde | ✓ projde |
Čistá náhodná matice projde prvními dvěma stupni, ale na prvočíslech a nasycení selže — nemá v sobě žádnou aritmetiku. Skutečné nuly i teoretické hladiny Berryho–Keatingova modelu xp doplněné o prvočíselné orbity projdou všemi čtyřmi. To je cenný návod, jak by musel vypadat hledaný kvantový „operátor“ — pokud podle Hilbertovy–Pólyovy domněnky vůbec existuje.
06Co z toho plyne (a co ne)
✓ Co studie ukázala
Berryho semiklasická teorie popisuje chování nul přes 17 řádů výšky bez jediného volného parametru — a pravá prvočísla v ní nelze nahradit náhodnými. Všechna data a skripty jsou zveřejněny s otevřenou licencí a otisky SHA-256, takže si výpočty může kdokoli zopakovat.
✗ Co z toho neplyne
Práce nic netvrdí o platnosti Riemannovy hypotézy. Ukazuje však, že kdo chce pochopit tajemství zeta funkce, musí vzít fyziku chaosu a prvočísla zároveň — jedno bez druhého nestačí.
Slovníček pojmů
- GUE
- Gaussovský unitární soubor — náhodné hermitovské matice. Jejich vlastní čísla se odpuzují stejně jako nuly zeta funkce na krátkých vzdálenostech.
- Unfolding (přeškálování)
- Přepočet nul tak, aby průměrná mezera mezi sousedními byla 1. Teprve pak lze srovnávat nuly z různých výšek.
- Σ²(L) — rozptyl počtu
- Jak moc kolísá počet nul v okně délky L kolem svého průměru L.
- Heisenbergův čas H
- Zde H = ln(t/2π). Orbity (prvočísla) s periodou ln p kratší než H určují dlouhodosahové chování nul.
- Cramérův model
- Náhodná „prvočísla“: každé číslo n je vybráno s pravděpodobností 1/ln n. Hustota sedí, polohy jsou náhodné.
- Hilbertova–Pólyova domněnka
- Myšlenka, že nuly zeta funkce jsou vlastními čísly nějakého samoadjungovaného operátoru — kvantového systému, který zatím nikdo nenašel.
Studie I · Berryho nasycení nul Riemannovy zeta funkce
Celá série na Zenodu: Studie I · Studie II · Studie III · licence CC BY 4.0 (text, data) a MIT (kód)
Grafy jsou interaktivní — najeďte myší na křivky a body.
The music of the primes · part 1 of 3
Primes tame quantum chaos
An expedition across 17 orders of magnitude in height — from the first zeros of the Riemann zeta function to zeros whose index numbers have 22 digits. Everywhere we find the fingerprints of the primes, which no randomness can imitate.
Imagine the chaos inside an atomic nucleus and the perfect order of the prime numbers dancing on the same floor. It sounds like a sci-fi script, yet this is exactly what our latest numerical study shows.
We examined the zeros of the Riemann zeta function over an enormous range: from the first few thousand up to heights around 1021, where zeros carry index numbers with 22 digits. The question was simple — do the zeros behave like random chaos, or do they carry a memory of the primes? The answer is clear: the primes leave fingerprints in the distribution of the zeros that no randomness can reproduce.
The zeta function ζ(s) = 1 + 1/2s + 1/3s + … is tied to the primes so tightly that in 1859 Bernhard Riemann showed how its zeros can exactly “play back” the distribution of primes. The non-trivial zeros have the form ½ + i·γ — they lie on the so-called critical line (the Riemann Hypothesis says all of them do). The number γ is their “height” t. You can see the first ten in the header of this article.
01Quantum chaos and atomic nuclei: why do physicists care about pure maths?
In the early 1970s the mathematician Hugh Montgomery computed how often two zeta zeros occur at a given distance from each other. When he showed the formula to the physicist Freeman Dyson, Dyson immediately recognised a familiar face: this is exactly the eigenvalue statistics of large random matrices — specifically of the GUE (Gaussian Unitary Ensemble).
Random-matrix theory made its name describing the energy levels of heavy atomic nuclei, where quantum chaos rules. So at short distances the zeta zeros behave like the levels of a chaotic quantum system: they repel each other and almost never sit right next to one another.
Sounds great. But there is a catch.
Random matrices describe local chaos — how zeros repel at close range. As soon as we look at long distances and count how many zeros fit into a wider window, random-matrix theory fails. It predicts that the fluctuation in the number of zeros keeps growing without limit. In reality the growth stops and saturates.
You stand on a bridge over a motorway and count the cars passing through a stretch of length L. Drivers keep their distance (that is the “repulsion”). Each time you get a slightly different count — and the variance of that count is what we call Σ²(L). For random matrices it grows forever with the length of the stretch. For zeta zeros it levels off at some point, as if the traffic were run by a hidden conductor. That conductor is the primes.
02Berry's magic: primes as orbits
In 1988 the British physicist Sir Michael Berry came up with an elegant idea. We can think of the zeta zeros as the energy levels of a (still unknown) quantum system — and of the primes as its closed periodic orbits. Each prime p and its powers contribute an “orbit” with period ln pk.
From this picture Berry derived a semiclassical formula for the number variance Σ²(L) that has no free parameters. Nothing is tuned, nothing is fitted to the data. The formula predicts that the variance first oscillates and then saturates.
The study tested it systematically on 18 data blocks across 17 orders of height:
- Accuracy: Berry's formula matches the real data with a deviation of just 0.26–1.21%.
- The signature of the prime 2: the variance chart shows fine ripples. They come from the shortest prime orbit — the prime 2 — with period ≈ H/ln 2. At height t ≈ 104 that is about 10.6 mean spacings (see the top chart above).
03The big test: real primes vs. a random fake
You might object: “Isn't it all just because there are a certain number of primes and they get sparser?” To rule this out we pitted Berry's formula against the Cramér model of random primes — fictitious “primes” with the same density as the real ones but random positions.
The lesson: the exact positions of the primes matter. Even at height 1021 the zeta zeros “remember” the specific arithmetic sequence 2, 3, 5, 7, 11, … — the right count alone is not enough.
04Saturation is not constant: height changes the rules
It was long quietly assumed that after suitable rescaling (so-called unfolding, which sets the mean spacing between zeros to 1) the zeros behave the same at every height. For long-range variance this is not true.
Why? The higher you climb, the more primes fit below the so-called Heisenberg time H = ln(t/2π) — that is, the more orbits have a period ln p shorter than H. Each extra prime adds a little to the height of the plateau. But it grows incredibly slowly, roughly like the logarithm of the logarithm of the height: 17 orders of magnitude add only about two tenths.
05A four-stage sieve: how to recognise the true zeta spectrum?
To stress-test the findings we built a four-stage sieve — a kind of lie detector for candidate spectra:
- 1 · Density and rigidity: do the total level count and its stability match?
- 2 · GUE spacings: do neighbouring levels behave locally like random matrices?
- 3 · Prime content: does the spectrum contain prime oscillations?
- 4 · Saturation: does the growth of the variance stop at long range?
| Spectral model | 1 · Density & rigidity | 2 · GUE spacings | 3 · Prime oscillations | 4 · Saturation |
|---|---|---|---|---|
| Pure GUE random matrix no arithmetic | ✓ passes | ✓ passes | ✗ fails | ✗ fails |
| Actual zeta zeros | ✓ passes | ✓ passes | ✓ passes | ✓ passes |
| Levels of the xp model + prime orbits | ✓ passes | ✓ passes | ✓ passes | ✓ passes |
A pure random matrix passes the first two stages but fails on primes and saturation — it contains no arithmetic. Both the actual zeros and the theoretical levels of the Berry–Keating xp model supplemented with prime orbits pass all four. That is a valuable guide to what the sought-after quantum “operator” would have to look like — if, as the Hilbert–Pólya conjecture suggests, it exists at all.
06What follows (and what does not)
✓ What the study showed
Berry's semiclassical theory describes the behaviour of the zeros across 17 orders of height without a single free parameter — and the real primes in it cannot be replaced by random ones. All data and scripts are published under an open licence with SHA-256 checksums, so anyone can reproduce the computations.
✗ What does not follow
The work makes no claim about the Riemann Hypothesis. It does show, however, that anyone wanting to understand the zeta function has to take chaos physics and the primes together — one without the other is not enough.
Glossary
- GUE
- Gaussian Unitary Ensemble — random Hermitian matrices. Their eigenvalues repel each other exactly like zeta zeros at short distances.
- Unfolding
- Rescaling the zeros so that the mean spacing between neighbours is 1. Only then can zeros from different heights be compared.
- Σ²(L) — number variance
- How much the number of zeros in a window of length L fluctuates around its mean L.
- Heisenberg time H
- Here H = ln(t/2π). Orbits (primes) with period ln p shorter than H govern the long-range behaviour of the zeros.
- Cramér model
- Random “primes”: each number n is chosen with probability 1/ln n. The density is right, the positions are random.
- Hilbert–Pólya conjecture
- The idea that the zeta zeros are the eigenvalues of some self-adjoint operator — a quantum system nobody has found yet.
Study I · Berry saturation of the Riemann zeros
The whole series on Zenodo: Study I · Study II · Study III · licence CC BY 4.0 (text, data) and MIT (code)
The charts are interactive — hover over the curves and points.
2 komentáře