Funkcja dzeta Riemanna

Z testwiki
Przejdź do nawigacji Przejdź do wyszukiwania
Wykres funkcji ζ w dziedzinie liczb rzeczywistych
Wykres funkcji ζ w dziedzinie liczb zespolonych uzyskany techniką kolorowania dziedziny.

Funkcja zeta Riemanna (funkcja dzeta Riemanna, funkcja ζ) – zespolona funkcja specjalna zdefiniowana w postaci szeregu

ζ(s)=n=11ns

dla dowolnej liczby zespolonej s o części rzeczywistej (s)>1 oraz jako przedłużenie analityczne powyższego szeregu dla pozostałych liczb zespolonych[1].

Funkcję ζ po raz pierwszy zdefiniował Leonhard Euler w XVIII w., jednak rozważał jej wartości jedynie dla zmiennych rzeczywistych. Dopiero Bernhard Riemann w artykule z listopada 1859 r. O liczbie liczb pierwszych mniejszych od zadanej wielkości (niem. Über die Anzahl der Primzahlen unter einer gegebenen Grösse) rozszerzył definicję Eulera na wszystkie liczby zespolone, udowodnił meromorficzność funkcji, przedstawił i udowodnił równanie funkcyjne opisujące funkcję na całej płaszczyźnie zespolonej i wykazał zależność między rozmieszczeniem jej miejsc zerowych a liczbą liczb pierwszych. Artykuł ten zawierał również sformułowanie hipotezy Riemanna, określanej jako najważniejszy problem otwarty matematyki[2].

Funkcja dzeta znajduje bardzo wiele zastosowań w analitycznej teorii liczb.

Postacie funkcji

Dla zmiennej o części rzeczywistej większej niż 1

Funkcja ζ jest pierwotnie definiowana za pomocą szeregu. W swojej pracy Riemann udowodnił także postać

ζ(s)=1Γ(s)0xs1ex1dx

dla (s)>1, wykorzystując przy tym odwrotność funkcji gamma.

W perspektywie teorii liczb, zdecydowanie największą rolę odgrywa iloczyn Eulera funkcji zeta, będący postaci

ζ(s)=p(11ps)1,

gdzie p oznacza iloczyn po wszystkich liczbach pierwszych[3]Szablon:Odn.

W pasie krytycznym

Często wykorzystywaną reprezentacją funkcji ζ na pasie 0<(s)<1 jest

ζ(s)=ss1s1{u}us+1du,

gdzie {x} oznacza część ułamkową. Postać tę można odczytać ze wzoru sumacyjnego Eulera[4]Szablon:Odn.

Na całej płaszczyźnie zespolonej

ζ(s)=1121sn=012n+1k=0n(1)k(nk)(k+1)s.

Równanie funkcyjne

Równaniem udowodnionym przez Riemanna, które opisuje zachowanie funkcji ζ

ζ(s)=2sπs1sin(πs2)Γ(1s)ζ(1s)

dla dowolnej liczby zespolonej s, gdzie Γ to funkcja gamma. Równanie to wiąże wartości funkcji dla zmiennych zespolonych s i 1s, symetrycznych wobec siebie względem prostej krytycznej (s)=12. Z równania tego można również odczytać, że trywialnymi miejscami zerowymi funkcji ζ2k=2,4,6, (ponieważ wtedy wartości 2sπs1, funkcji Γ(1s) i ζ(1s) są skończone, a sin(2kπ2)=sin(kπ)=0). Jednocześnie, jeśli s=2k (jest dodatnią liczbą parzystą), to ζ(s)0, ponieważ w tych miejscach występują bieguny funkcji Γ(1s). Ponadto, jeśli s0 jest nietrywialnym miejscem zerowym ζ (0<(s0)<1), to jest nim również 1s0. Jeśli s0 nie leży na prostej krytycznej, to są to dwie różne wartości, dlatego nietrywialne miejsca zerowa poza prostą krytyczną muszą występować w parach.

Przedstawimy poniżej trzy dowody prawdziwości równania funkcyjnego wg Titchmarsha (spośród aż siedmiu przedstawionych)Szablon:Odn.

Dowód 1. W pierwszym dowodzie wyprowadzamy, a następnie korzystamy z postaci funkcji ζ wykorzystującej całkę z funkcji {u}us1.

Wzór sumacyjny Eulera mówi, że dla dowolnej funkcji f o ciągłej pochodnej zachodzi

y<nxf(n)=yxf(u)du+yx{u}f(u)du{x}f(x)+{y}f(y).

Biorąc f(n)=ns, otrzymamy

1nx1ns=1xduuss1x{u}us+1du{x}xs+1.

Aby z lewej strony równania otrzymać funkcję Riemanna, chcemy aby x. Widzimy, że

1xduus=1x1ss1=1s11xs1(s1)=1s1.

Stąd

ζ(s)=1+1s1s1{u}us+1du=ss1s1{u}us+1du

dla (s)>1.

Teraz skorzystajmy z szeregu Fouriera zbieżnego do części ułamkowej. Mamy

{x}=12+n=1sin(2nπx)nπ

dla x niecałkowitych. Podstawiając pod całkę, otrzymamy

ζ(s)=sπn=11n0sin(2nπu)us+1du=sπn=1(2nπ)sn0sinuus+1du=sπ(2π)s(Γ(s))sin(πs2)ζ(1s).

Upraszczając i korzystając z Γ(s)=sΓ(1s), otrzymamy równanie.

Dowód 2. Dowód ten przeprowadzany jest ze szczególnym uwzględnieniem narzędzi analizy zespolonej.

Zacznijmy od udowodnienia szczególnej postaci funkcji ζ (przedstawionej wcześniej w artykule). Całkując przez podstawienie, pokazujemy, że

0xs1enxdx=1ns0us1eudu=Γ(s)ns,

więc dla (s)>1 mamy

Γ(s)ζ(s)=n=1Γ(s)ns=n=10xs1enxdx=0xs1n=1enxdx=0xs1ex1dx

gdzie przy trzeciej równości skorzystaliśmy z twierdzenia o całkowaniu wyraz po wyrazie, ponieważ występujący szereg funkcyjny jest jednostajnie zbieżny, a ostatnia równość zachodzi, ponieważ występujący tam szereg jest zwykłym szeregiem geometrycznym.

Rozważmy całkę

I(s)=Czs1ez1dz,
kontur Hankela

gdzie C oznacza kontur Hankela (krzywą składającą się z prostej ϵ>0 od + do ρ (0<ρ<2π), fragmentu dodatnio określonego (przeciwnie do ruchu wskazówek zegara) okręgu |z|=ρ okrążając 0 i prostej ϵ od ρ do +). Przyjęliśmy tutaj

zs1=e(s1)logs,

gdzie logarytm jest rzeczywisty na początku konturu. Na zadanym okręgu mamy

|zs1|=e((s)1)log|z|(s)arg(z)|z|(s)1e2π|(s)|

i

|ez1|>A|z|

dla pewnej stałej A>0, więc

1A|z2s|>|zs1ez1|.

Stąd, jeśli (s)>1, to przy ρ0 wartość powyższej całki na części okręgu |z|=ρ dąży do 0. Dlatego, zakładając dalej ρ0, mamy

I(s)=0xs1ex1dx+0(xe2πi)s1ex1dx=(e2πis1)Γ(s)ζ(s)=2πieπisΓ(1s)ζ(s).

Dlatego

ζ(s)=eπisΓ(1s)2πiI(s)

dla (s)>1. Jednakże całka I(s) zbiega jednostajnie na każdym skończonym obszarze płaszczyzny zespolonej. W ten sposób przedłużymy analitycznie funkcję Riemanna.

Niech Cn będzie konturem równym ϵ>0 od + do (2n+1)π, następnie dodatnio określonym fragmentem kwadratu (±1±i)(2n+1)π, a potem ϵ od (2n+1)π do +. Całkowana funkcja ma pomiędzy konturami C a Cn bieguny w punktach ±2πi,,n2πi. Residua w punktach m2πi i m2πi to w sumie

(2mπe12πi)s1+(2mπe32πi)s1=(2mπ)s1eπi(s1)2cos(π(s1)2)=2(2mπ)s1eπissin(πs2).

Dlatego z twierdzenia o residuach, mamy

I(s)=Cnzs1ez1dz+4πieπissin(πs2)m=1n(2mπ)s1.

Biorąc n i wiedząc, że funkcja (ez1)1 jest ograniczona oraz |zs1|=O(z(z)1) wnioskujemy, że całka po konturze Cn dąży do 0. Dlatego

I(s)=4πieπissin(πs2)m=1(2mπ)s1=4πieπissin(πs2)(2π)s1ζ(1s).

Upraszczając, otrzymamy równanie funkcyjne.

Dowód 3.

Jeśli (s)>0, to, całkując przez podstawienie,

0x12s1en2πxdx=Γ(s2)nsπs2.

Stąd, dla (s)>1 zachodzi

ζ(s)Γ(s2)πs2=n=1Γ(s2)nsπs2=n=10x12s1en2πxdx=0x12s1n=1en2πxdx,

gdzie, podobnie jak w poprzednim dowodzie, w ostatniej równości skorzystaliśmy z twierdzenia o całkowaniu wyraz po wyrazie, ponieważ występujący szereg funkcyjny jest jednostajnie zbieżny (stąd silniejsze założenie o części rzeczywistej s). Oznaczmy

ψ(x)=n=1en2πx.

Wówczas

ζ(s)=πs2Γ(s2)0x12s1ψ(x)dx.

Korzystając ze wzoru sumacyjnego Poissona, otrzymamy

n=en2πx=1xn=en2πx,

a stąd

2ψ(x)+1=1x(2ψ(1x)+1).

Dlatego

πs2Γ(s2)ζ(s)=01xs21ψ(x)dx+1xs21ψ(x)dx.

Powyższe wyrażenie jest równe

01xs21(1xψ(1x)+12x12)dx+1xs21ψ(x)dx,

co z kolei równa się

1s11s+01xs232ψ(1x)dx+1xs21ψ(x)dx.

Zatem

πs2Γ(s2)ζ(s)=1s(s1)+1(xs212+xs2+1)ψ(x)dx,

gdzie całka jest zbieżna dla wszystkich s zespolonych, więc wyrażenie po lewej można przedłużyć analitycznie. Ponadto, prawa strona nie zmienia wartości po zastąpieniu s przez 1s. Dlatego

πs2Γ(s2)ζ(s)=π1s2Γ(1s2)ζ(1s).

Po uproszczeniu otrzymamy równanie funkcyjne.

Wzory związane z funkcją zeta

Szablon:Osobny artykuł Na półpłaszczyźnie (s)>1 funkcja Riemanna jest wyrażona przez iloczyn Eulera ζ(s)=p(11p)1,

gdzie p oznacza iloczyn po wszystkich liczbach pierwszych[3]Szablon:Odn.

Ponadto, dla (s)>1 prawdziwe są tożsamości

1ζ(s)=p(11ps)=n=1μ(n)ns

oraz

ζ(s)ζ(2s)=p(1+1ps)=n=1|μ(n)|ns

gdzie μ to funkcja Möbiusa,

ζ(2s)ζ(s)=p(1+1ps)1=n=1λ(n)ns,

gdzie λ to funkcja Liouville’a,

Wszystkie powyższe iloczyny można otrzymać przez zwykłe przekształcenia algebraiczne na produkcie Eulera funkcji ζ, a szeregi po prawej – przez wymnażanie czynnikówSzablon:Odn.

Prawdziwe są również wzory

ζ(s)2=n=1τ(n)ns,

gdzie τ oznacza liczbę dodatnich dzielników, a także ogólniej

ζ(s)k=n=1τk(n)ns,

dla k=2,3,, gdzie τk(n) oznacza liczbę sposobów na przedstawienie liczby całkowitej dodatniej n jako iloczyn k czynników całkowitych dodatnich.

ζ(s)2ζ(2s)=n=12ω(n)ns,

gdzie ω to funkcja pierwsza omega, czyli liczba dzielników pierwszych. Ponadto

ζ(s)3ζ(2s)=n=1τ(n)2ns,
ζ(s)4ζ(2s)=n=1{τ(n)}2ns,

przy czym wszystkie powyższe wzory są prawdziwe na obszarze zbieżności szeregów, czyli (s)>1Szablon:Odn. Dodatkowo

ζ(s1)ζ(s)=n=1φ(n)ns,

gdzie φ oznacza tocjent Eulera i

121s12sζ(s1)=n=1a(n)ns,

gdzie a(n) jest największym dzielnikiem nieparzystym liczby n. Te dwa wzory są prawdziwe dla (s)>2Szablon:Odn.

Biorąc logarytm zespolony iloczynu Eulera, mamy

logζ(s)=logp(11ps)1=plog(11ps)=pk=11kpks,

przy czym ostatnia równość zachodzi, ponieważ k=1(kxk)1 jest szeregiem potęgowym funkcji log(1x1). Różniczkując wyraz po wyrazie, otrzymujemy

ζ(s)ζ(s)=pk=1logppks=n=1Λ(n)ns,

gdzie Λ oznacza funkcję von MangoldtaSzablon:Odn.

Związek z liczbami Bernoulliego:

ζ(2n)=(1)n+1B2n(2π)2n2(2n)!,

dla każdej liczby parzystej dodatniej 2n, gdzie Bk to k-ta liczba Bernoulliego. Ponadto dla liczb całkowitych ujemnych n:

ζ(n)=Bn+1n+1.

Zatem funkcja ζ przyjmuje wartość 0 dla każdej ujemnej liczby parzystej.

Logarytm funkcji zeta można wyrazić jakoSzablon:Odn:

logζ(z)=z2π(x)x(xz1)dx,

gdzie π(x) to funkcja licząca liczby pierwszeSzablon:Odn.

Znaczenie w teorii liczb

Pierwsza strona artykułu Riemanna O liczbie liczb pierwszych mniejszych od zadanej wielkości z 1859 r.

Funkcja ζ jest podstawowym obiektem badań analitycznej teorii liczb. Jest to funkcja meromorficzna, która swoim zachowaniem opisuje zjawiska dyskretne, takie jak rozmieszczenie liczb pierwszych. Mówiąc dokładniej, nietrywialne miejsca zerowe funkcji zeta występują we wzorach opisujących chociażby funkcję liczącą liczby pierwsze czy funkcję Czebyszewa.

Znając wzór

ζ(s)ζ(s)=n=1Λ(n)ns

i oznaczając drugą funkcję Czebyszewa jako ψ(x)=nxΛ(n) oraz ψ0(x)=ψ(x)12Λ(x) dla x całkowitych i ψ0(x)=ψ(x) dla wszystkich pozostałych, możemy skorzystać ze wzoru Perrona by uzyskać

ψ0(x)=12πilimT1iT1+iTζ(s)ζ(s)xssds.

Stąd otrzymujemy wzór explicite

ψ0(x)=xlog(2π)ρxρρρtxρtρt,

gdzie ρ i ρt oznaczają odpowiednio sumy po wszystkich nietrywialnych i trywialnych miejscach zerowych funkcji ζ. Widzimy, że

ρtxρtρt=k=112kx2k=log(11x2).

Zatem[5]

ψ0(x)=xlog(2π)log(11x2)ρxρρ.

Niech π będzie funkcją liczącą liczby pierwsze. Równoważnie, jeśli oznaczymy przez π0(x)=π(x)12 dla x pierwszych oraz π0(x)=π(x) dla pozostałych x, sumując po częściach, otrzymamy[6]

π0(x)=n=1μ(n)nli(x1n)ρn=1μ(n)nli(xρn),

gdzie li oznacza logarytm całkowy.

Z powyższych możemy wnioskować, że udowodnienie, że ζ nie ma żadnych zer na prostej (s)=1 jest równoważne z twierdzeniem o liczbach pierwszych, a hipoteza Riemanna jest równoważna z błędem w szacowaniu rzędu π(x)=Li(x)+O(x12+ϵ).

Miejsca zerowe

Moduły funkcji Z Riemanna-Siegela (linia przerywana) oraz funkcji zeta Riemanna (linia ciągła) na prostej krytycznej

Równanie funkcyjne mówi, że funkcja Riemanna przyjmuje wartość równą 0 dla wszystkich ujemnych liczb parzystych, czyli 2,4,6, Są to tzw. zera trywialne. Są trywialne w takim sensie, że łatwo jest udowodnić ich występowanie, ponieważ są one miejscami, w których sinus przyjmuje wartość 0. O wiele większe znaczenie mają zera nietrywialne, czyli wszystkie miejsca zerowe niebędące zerami trywialnymi.

Wiadomo, że wszystkie nietrywialne miejsca zerowe muszą leżeć na pasie krytycznym, zdefiniowanym jako {s:0<(s)<1}. Hipoteza Riemanna mówi, że wszystkie nietrywialne miejsca zerowe leżą na prostej krytycznej (s)=12.

Niech N(T) oznacza liczbę miejsc zerowych funkcji Riemanna postaci s=12+it przy 0<t<T. W 1921 Hardy i Littlewood udowodnili, że nieskończenie wiele nietrywialnych zer leży na prostej krytycznej[7], a dokładniej mówiąc wykazali, że dla każdego ϵ>0 istnieje stała T0 taka, że

N(T)>T34ϵ

dla wszystkich T>T0.

Współcześnie wiadomo, że

N(T)=T2πlogT2πT2π+O(logT)

dla T4. Ponadto, jeśli hipoteza Riemanna jest prawdziwa, toSzablon:Odn

N(T)=T2πlogT2πT2π+O(logTloglogT)

W 1989 Conrey udowodnił, że przynajmniej 25 z nich musi leżeć na tej prostej[8].

Niektóre wartości

ζ(1)=1120,0833333
ζ(2)=1+122+132+=π261,6449341Szablon:Odn
ζ(4)=1+124+134+=π4901,0823232Szablon:Odn
ζ(6)=1+126+136+=π69451,0173431Szablon:Odn
ζ(8)=1+128+138+=π894501,0040774
ζ(10)=1+1210+1310+=π10935551,0009946

Ogólnie, dla p, mamy:

ζ(2p)=(1)p+1B2p(2π)2p2(2p)!,Szablon:Odn

gdzie B2p to liczba Bernoulliego z indeksem 2p.

Zobacz też

Przypisy

Szablon:Przypisy

Bibliografia

Linki zewnętrzne

Szablon:Funkcje specjalne

Szablon:Kontrola autorytatywna