关于黎曼猜想

前几日关于黎曼猜想（Riemann Hypothesis）被证明的消息在网络上炒得沸沸扬扬，菲尔兹奖、阿贝尔奖双料得主Michael Francis Atiyah爵士声称他证明了这个困扰了数学界160年的问题，且其目的是为了研究精细结构常数。然而，在他公布了他的Simple Proof后，却不被学界所看好，甚至遭到了一些尖锐的批评。毕竟，这篇只有五页纸的文章，刨去介绍等内容，并没有展示太多的细节；其证明过程，又依赖于一个存疑的Todd function。如果如证明中所说，黎曼猜想真的指向了精细结构常数，这对于由第一性原理出发构建物理学体系将有极大的启发性，其结果无疑是令人激动的。但是经过计算，却发现文中所声称的极限存在问题，整个理论的正确性难以得到证实。
有鉴于此，博主研究该证明过程的计划搁置了（毕竟没有什么价值，大概可以用来提升读Paper的能力#捂脸），还是写篇文章介绍下黎曼猜想本身吧。

猜想的提出

时间回到1859年，著名的德国数学家波恩哈德·黎曼（Bernhard Riemann）在他的论文《Über die Anzahl der Primzahlen unter einer gegebenen Größe》中提及了这个著名的猜想。黎曼本人在数学界颇有建树，但他却并没有试图对论文中的猜想给出证明。
这个猜想是关于黎曼$\zeta$函数的，该函数的定义为：
$$\zeta(s)=\frac{1}{1^s}+\frac{1}{2^s}+\frac{1}{3^s}+\frac{1}{4^s}+\cdots$$
注意，$\zeta$函数必须经过解析延拓，才能在全复平面上有定义，稍后将给出解析延拓的过程。
而黎曼猜想就是：$\zeta(s)$的所有非平凡零点的实部均为$\frac{1}{2}$。

$\zeta$函数的解析延拓

前言

在中学数学中，我们所接触到的函数都是实变函数，即从实数集到实数集的映射，写为$f:\mathbb{R}\rightarrow\mathbb{R}$。在接触了虚数单位i和复数后，就可以开始定义复变函数，即复数集到复数集的映射，写为$f:\mathbb{C}\rightarrow\mathbb{C}$。实变函数的定义域被限制在了实轴上，而为了表示出一个与之对应的复变函数，就需要扩充其定义域到复平面上。这样，可以定义出复数域上的各种函数：指数、对数、三角函数等。由于复数的性质，还需要根据枝点规定割线、单值分支等，来确定出函数值。而原本定义域只在复平面上某一区域的复变函数，经过严格的数学推导（例如级数展开等），也可以扩充其定义域。这时，解析延拓的办法应运而生。对于$\zeta(s)$函数而言，非常明显，当$s$取在实轴上（即作为实变函数时），若$s\le 1$，其函数值发散；而$s>1$时收敛。特别的，$s=1$时为调和级数，而$s=2$时由高斯最早计算出了其函数值$\frac{\pi^2}{6}$。
在复数域上，可以证明，$\zeta(s)$在$Re(s)>1$上解析。接下来，我们需要将其定义域扩展到全复平面；之前的级数定义将不在对$Re(s)<1$的情况使用，因为它不收敛，无法采用Cauchy和。

梅林变换

梅林（Mellin）变换是一个积分变换，它把一个正实数域$\mathbb{R}^+$上的函数变换为一个复数域$\mathbb{C}$上的函数。设$f:\mathbb{R}^+\rightarrow\mathbb{C}$是一个函数，定义：
$$\mathcal{M}f(z)=\int_0^{+\infty}f(s)s^z\frac{ds}{s}$$
称$\mathcal{M}f:\mathbb{C}\rightarrow\mathbb{C}$为$f$的梅林变换。
取$f(s)$为$f_\lambda (s)=e^{-\lambda s}$，则有：
\begin{align*}
\mathcal{M}f(z)&=\int_0^{+\infty}e^{-\lambda s}s^z\frac{ds}{s}\\
&=\lambda^{-z}\int_0^{+\infty}e^{-s}s^{z-1}ds\\
&=\lambda^{-z}\Gamma(z)
\end{align*}
这里的$\Gamma(z)$就是一般数理方程中会涉及到的$\Gamma$函数，其性质不再赘述。
记住以上结论，取$\lambda=\pi n^2$，继续构造函数：
$$f(s)=\sum_{n=1}^{\infty}f_{\pi n^2}(s)=\sum_{n=1}^{\infty}e^{-\pi n^2 s}$$
注意到梅林变换是线性的，这时，
\begin{align*}
\mathcal{M}f(z)&=\sum_{n=1}^{\infty}\mathcal{M}f_{\pi n^2}(z)\\
&=\sum_{n=1}^{\infty}(\pi n^2)^{-z}\Gamma(z)\\
&=\pi^{-z}\sum_{n=1}^{\infty}n^{-2z}\Gamma(z)\\
&=\pi^{-z}\zeta(2z)\Gamma(z)
\end{align*}
最后一个等号使用了$\zeta$函数的原始定义。考虑到：
$$\mathcal{M}f(z)=\int_0^{+\infty}(\sum_{n=1}^{\infty}e^{-\pi n^2 s})s^{z-1}ds$$
于是，利用$\mathcal{M}f(z)$的两个等价表达式，我们得到了一个重要结论：
$$\pi^{-z}\zeta(2z)\Gamma(z)=\int_0^{+\infty}(\sum_{n=1}^{\infty}e^{-\pi n^2 s})s^{z-1}ds$$
注意，这里的推导并非十分严谨，还需要证明其收敛性，才能交换积分与极限。这要求$Re(z)>\frac{1}{2}$。

泊松求和

再来看一下Fourier变换
$$\mathcal{F}f(x)=\int_{-\infty}^{+\infty}f(t)e^{ixt}dt$$
对于足够『好』的函数$f$，可以证明泊松（Poisson）求和公式
$$\sum_{n\in\mathbb{Z}}f(n)=\sum_{n\in\mathbb{Z}}\mathcal{F}f(n)$$
这部分我们就证明这个公式。
假设$f$足够『好』（要多好有多好），定义
$$F(t)=\sum_{n\in\mathbb{Z}}f(n+t)$$
则$F$是$\mathbb{R}$上周期为$1$的函数，于是可以构造$F$的Fourier级数：令
$$c_k=\int_0^1F(t)e^{-2\pi ikt}dt$$
是第$k$个Fourier系数，则
$$F(t)\sim\sum_{k\in\mathbb{Z}}c_ke^{2\pi ikt}$$
而这个系数
\begin{align*}
c_k&=\int_0^1F(t)e^{-2\pi ikt}dt\\
&=\int_0^1\sum_{n\in\mathbb{Z}}f(n+t)e^{-2\pi ik(n+t)}dt\\
&=\int_{-\infty}^{+\infty}f(t)e^{-2\pi ikt}dt\\
&=\mathcal{F}f(k)
\end{align*}
从而
$$\sum_{n\in\mathbb{Z}}f(n)=F(0)\sim\sum_{k\in\mathbb{Z}}c_k=\sum_{k\in\mathbb{Z}}\mathcal{F}f(k)$$
就得到了Poisson求和公式形式上的『证明』。
那么$f$究竟要满足什么条件呢？
如果周期函数$F\in C^2[0,1]$，就有$F(t)=\sum_{k\in\mathbb{Z}}c_ke^{2\pi ikt}$。因此，我们要求$f$二阶连续可导，并且
$$\sum_{n\in\mathbb{Z}}||f||_{[n.n+1],\infty}+||f’||_{[n.n+1],\infty}+||f’’||_{[n.n+1],\infty}<+\infty$$
使用导数一致收敛的判别法则，就得到$F\in C^2[0,1]$啦。
也就是说，如果$f\in C^2(\mathbb{R})$，并且满足以上条件，就有Poisson求和公式
$$\sum_{n\in\mathbb{Z}}f(n)=\sum_{n\in\mathbb{Z}}\mathcal{F}f(n)$$

解析延拓

设$\lambda>0$，令
\begin{eqnarray*}
\theta(\lambda)&=&\sum_{n\in\mathbb{Z}}e^{-\pi n^2\lambda}\\
\psi(\lambda)&=&\sum_{n=1}^\infty e^{-\pi n^2\lambda}
\end{eqnarray*}
很显然$\theta(\lambda)=2\psi(\lambda)+1$。
考虑$f(t)=e^{-\pi t^2\lambda}$，计算得$\mathcal{F}f(x)=\frac{1}{\sqrt{\lambda}}e^{-\pi\frac{x^2}{\lambda}}$，而且满足相应条件，因此用Poisson求和公式，
$$\sum_{n\in\mathbb{Z}}e^{-\pi n^2\lambda}=\frac{1}{\sqrt{\lambda}}\sum_{n\in\mathbb{Z}}e^{-\pi\frac{n^2}{\lambda}}$$
即$\theta(\lambda)=\frac{1}{\sqrt{\lambda}}\theta\left(\frac{1}{\lambda}\right)$。
换成$\psi$就是$\psi\left(\frac{1}{\lambda}\right)=\sqrt{\lambda}\psi(\lambda)+\frac{\sqrt{\lambda}-1}{2}$。
第二部分的最后，得到了
\begin{align*}
\pi^{-z}\zeta(2z)\Gamma(z)&=\int_0^{+\infty}\left(\sum_{n=1}^\infty e^{-\pi n^2s}\right)s^{z-1}ds\\
&=\int_0^{+\infty}\psi(s)s^{z-1}ds
\end{align*}
因此
\begin{align*}
\pi^{-z}\zeta(2z)\Gamma(z)&=\int_0^{+\infty}\psi(s)s^{z-1}ds\\
&=\int_1^{+\infty}\psi(s)s^{z-1}ds+\int_0^1\psi(s)s^{z-1}ds\\
&=\int_1^{+\infty}\psi(s)s^{z-1}ds+\int_1^{+\infty}\psi\left(\frac{1}{s}\right)s^{1-z}\frac{ds}{s^2}\\
&=\int_1^{+\infty}\psi(s)s^{z-1}ds+\int_1^{+\infty}\left(\psi(s)\sqrt{s}+\frac{\sqrt{s}-1}{2}\right)s^{-1-z}ds\\
&=\int_1^{+\infty}\psi(s)(s^{z-1}+s^{-z-\frac{1}{2}})ds+\int_1^{+\infty}\frac{s^{-z-\frac{1}{2}}-s^{-z-1}}{2}ds\\
&=\int_1^{+\infty}\psi(s)(s^{z-1}+s^{-z-\frac{1}{2}})ds+\frac{1}{2z(2z-1)}
\end{align*}
以上计算的前提是$Re(z)>\frac{1}{2}$。
令$w=2z$，则当$Re(w)>1$时，
$$\sqrt{\pi}^{-w}\zeta(w)\Gamma\left(\frac{w}{2}\right)=\int_1^{+\infty}\psi(s)(\sqrt{s}^{w-2}+\sqrt{s}^{-w-1})ds-\frac{1}{w(1-w)}$$
现在观察：等号右边是一个亚纯函数（亚纯函数就是『能表示成两个解析函数的商』的函数，或者『没有本性奇点』的函数，两个亚纯函数经过四则运算仍然是亚纯函数），而等号左边的$\sqrt{\pi}^{-w}$、$\Gamma\left(\frac{w}{2}\right)$也都是亚纯函数。这样，上式相当于给出了$\zeta(w)$在$w\in\mathbb{C}$上的定义！

所有自然数的『和』

这还没完，试试在等号右边把$w$换为$1-w$，你会发现式子根本没变。这就说明：
$$\sqrt{\pi}^{-w}\zeta(w)\Gamma\left(\frac{w}{2}\right)=\sqrt{\pi}^{w-1}\zeta(1-w)\Gamma\left(\frac{1-w}{2}\right)$$
现在让$w=-1$，就有
$$\sqrt{\pi}\zeta(-1)\Gamma\left(-\frac{1}{2}\right)=\sqrt{\pi}^{-2}\zeta(2)\Gamma(1)$$
我们知道
$$\zeta(2)=\sum_{n=1}^\infty\frac{1}{n^2}=\frac{\pi^2}{6}$$
并且，
\begin{eqnarray*}
\Gamma(1)&=&1\\
\Gamma\left(\frac{1}{2}\right)&=&\sqrt{\pi}\\
\Gamma\left(-\frac{1}{2}\right)&=&\frac{\Gamma\left(\frac{1}{2}\right)}{-\frac{1}{2}}=-2\sqrt{\pi}
\end{eqnarray*}
因此：
$$\zeta(-1)=\frac{\sqrt{\pi}^{-2}\zeta(2)\Gamma(1)}{\sqrt{\pi}\Gamma\left(-\frac{1}{2}\right)}=\frac{\frac{\pi^2}{6}}{\sqrt{\pi}^3(-2)\sqrt{\pi}}=-\frac{1}{12}$$
我们便得出了全体自然数的『和』为$-\frac{1}{12}$这一结论。

黎曼猜想与素数

欧拉乘积公式

从$\zeta$函数的原始定义出发：
$$\zeta(s)=\frac{1}{1^s}+\frac{1}{2^s}+\frac{1}{3^s}+\cdots$$
在等式两边同时乘以第二项：
$$\frac{1}{2^s}\times\zeta(s)=\frac{1}{2^s}+\frac{1}{4^s}+\frac{1}{6^s}+\cdots$$
两式相减得：
$$(1-\frac{1}{2^s})\times\zeta(s)=\frac{1}{1^s}+\frac{1}{3^s}+\frac{1}{5^s}+\cdots$$
这时，等号右边的分母上已经没有2的倍数了。重复这一过程：
$$\frac{1}{3^s}\times(1-\frac{1}{2^s})\zeta(s)=\frac{1}{3^s}+\frac{1}{9^s}+\frac{1}{15^s}+\cdots$$
以上两式相减，又有：
$$(1-\frac{1}{3^s})(1-\frac{1}{2^s})\zeta(s)=\frac{1}{1^s}+\frac{1}{5^s}+\frac{1}{7^s}+\cdots$$
等号右边也没有3的倍数了。无限重复此过程，就有：
$$\cdots(1-\frac{1}{11^s})(1-\frac{1}{7^s})(1-\frac{1}{5^s})(1-\frac{1}{3^s})\zeta(s)=1$$
即：
$$\zeta(s)=\frac{1}{1-\frac{1}{2^s}}\times\frac{1}{1-\frac{1}{3^s}}\times\frac{1}{1-\frac{1}{5^s}}\times\frac{1}{1-\frac{1}{7^s}}\times\cdots$$
可以写成：
$$\zeta(s)=\sum_n\frac{1}{n^s}=\prod_p\frac{1}{1-p^{-s}}$$
这就是欧拉乘积公式，该表达式首先出现在1737年一篇题为《Variae observationes circa series infinitas》（无穷级数的各种观察）的论文中。
该恒等式展示了素数与$\zeta$函数间的联系。

莫比乌斯函数

奥古斯特·费迪南德·莫比乌斯之后重写了欧拉乘积公式，创造了一个新的求和。首先引入莫比乌斯函数：
$$\mu(n)=
\begin{cases}
1 & n是偶数个不同的素数相乘\\
-1 & n是奇数个不同的素数相乘\\
0 & n被某个素数的平方整除
\end{cases}
$$
和式用$\mu(n)$表示如下：
$$\frac{1}{\zeta(s)}=\sum_{n\in\mathbb{N}}\frac{\mu(n)}{n^s}$$

素数计数函数

定义为
$$\pi(x)=\sum_{p\le x}1$$
表示不超过x的素数个数。
在间断点处用左右平均代替更好，重新定义：
$$\pi(x)=\frac{1}{2}\left(\sum_{p\le x}1+\sum_{p< x}1\right)$$
令
$$\Pi(x)=\sum_n\frac{1}{n}\pi(x^{\frac{1}{n}})$$
有莫比乌斯反转
$$\pi(x)=\sum_n\frac{\mu(n)}{n}\Pi(x^{\frac{1}{n}})$$

素数定理

素数定理由高斯和勒让德独立地阐述：
$$\mathop{\lim}_{x\rightarrow\infty}\frac{\pi(x)}{\frac{x}{\ln(x)}}=1$$
目前人类证明的素数定理的误差界是
$$\pi(x)=Li(x)+O\left(x\exp\left({-\frac{A(\ln x)^{0.6}}{(\ln\ln x)^{0.2}}}\right)\right)$$
从概率的角度来说，素数定理说明若你随机选择一个自然数$x$，那么该数字为素数的概率约为$\frac{1}{\ln(x)}$。这意味着前$x$个整数中相邻素数之间的平均间隙约为$\ln(x)$。

对数积分函数

定义
$$Li(x)=\int_2^x\frac{1}{\ln(t)}dt$$
称为对数积分函数。
那么
$$\Pi(x)=Li(x)-\sum_{\rho}Li(x^{\rho})-\ln 2+\int_x^\infty\frac{dt}{t(t^2-1)\ln t}$$
其中求和号的$\rho$取$\zeta(s)$的非平凡零点，求和之后通常表现为震荡。至于剩下的常数和积分，对较大的$x$可以忽略。