黎曼猜想发现过程(二):数学原理与解析框架

黎曼猜想发现过程(二):数学原理与解析框架

本文是《黎曼猜想发现过程(一)》的数学续篇。 上一篇讲的是人与历史,这一篇讲的是结构与证明。

全文的符号约定: \(s = \sigma + it \in \mathbb{C}\),\(\zeta\) 表示 Riemann \(\zeta\) 函数, \(\rho = \beta + i\gamma\) 表示其非平凡零点。


0. 问题陈述

黎曼猜想(RH):Riemann \(\zeta\) 函数的所有非平凡零点均位于临界线

\[\Re(s) = \frac{1}{2}\]

上。

等价地,记

\[\xi(s) = \frac{1}{2}s(s-1)\pi^{-s/2}\Gamma\!\left(\frac{s}{2}\right)\zeta(s)\]

则 RH 等价于断言 \(\Xi(t) = \xi\!\left(\frac{1}{2}+it\right)\) 的所有零点都是实数。

为什么这个”零点位置”问题如此难? 因为显式公式表明:

\[\text{素数分布的涨落} \;\longleftrightarrow\; \text{零点的位置}\]

这是一个全局性的对应关系。单个零点的位置影响极小,但全部零点的集体位置决定了我们对素数分布的认知精度。要证明 RH,等价于要证明一个由无穷多个点组成的集合恰好落在一条直线上——而我们只有有限的手段去检验无穷。

下面把这个框架一层层搭起来。


1. Euler 乘积:解析数论的起点

1.1 绝对收敛区域

对 \(\Re(s) > 1\),级数

\[\zeta(s) = \sum_{n=1}^{\infty} n^{-s}\]

绝对收敛,且由 Euler 乘积公式

\[\zeta(s) = \prod_{p}\left(1 - p^{-s}\right)^{-1}\]

与之等价。

1.2 取对数展开

在 \(\Re(s) > 1\) 中,对乘积取对数并展开几何级数,双重求和可交换(绝对收敛保证):

\[\log\zeta(s) = \sum_{p}\sum_{k \ge 1}\frac{p^{-ks}}{k}\]

于是

\[\frac{\log\zeta(s)}{s} = \sum_{n=1}^{\infty}\frac{\Lambda(n)}{n^{s}}\]

其中 \(\Lambda\) 是 von Mangoldt 函数:

\[\Lambda(n) = \begin{cases} \log p, & n = p^k \\ 0, & \text{其他} \end{cases}\]

这是整个理论最关键的一步:一个”关于所有整数求和”的函数,被改写成了”关于素数(及其幂)求和”的函数。素数的乘法结构被编码进了 \(\Lambda\)。

1.3 Perron 反演:回到计数函数

由 Mellin 反演,若 \(F(s) = \sum_{n\ge1} a_n n^{-s}\)(\(\Re(s) > \sigma_a > 0\)),则

\[\sum_{n \le x} a_n = \frac{1}{2\pi i}\int_{c-i\infty}^{c+i\infty}\frac{F(s)x^{s}}{s}\,ds, \qquad c > \sigma_a\]

取 \(a_n = \Lambda(n)\)、\(F(s) = -\log\zeta(s)/s\),得到

\[\psi(x) := \sum_{n \le x}\Lambda(n) = -\frac{1}{2\pi i}\int_{c-i\infty}^{c+i\infty}\frac{\log\zeta(s)}{s}x^{s}\,ds\]

核心思想:\(x^s = e^{s\log x}\),积分沿竖线时 \(x^{it}\) 以频率 \(\log x\) 振荡,因此积分可以用留数定理转化为”竖线左侧各奇点贡献之和”。

这就把一个离散的计数问题转成了解析的奇点问题。

1.4 为什么需要解析延拓

上面的积分要求把竖线从 \(c > 1\) 往左推。推到哪儿为止,取决于 \(\log\zeta(s)\) 在左侧的解析结构——而这必须先有解析延拓。

这就是为什么 1859 年的那篇论文必须从”让 \(\zeta\) 定义在整个复平面上”开始。


2. 解析延拓与函数方程

2.1 第一步:积分表示

分部积分(对 \(\Re(s) > 1\))给出

\[\zeta(s) = \frac{1}{\Gamma(s)}\int_0^\infty \frac{x^{s-1}}{e^x - 1}\,dx\]

关键观察:被积函数在 \(x=0\) 附近可展为 \(x^{-1} - \frac12 + \frac{x}{12} - \cdots\),所以积分在 \(\Re(s) > 0\) 收敛。

因此:

\[\zeta(s) = \frac{1}{\Gamma(s)}\int_0^\infty\frac{x^{s-1}}{e^x-1}\,dx, \qquad \Re(s) > 0\]

这个式子把解析延拓的问题变成了一个关于积分的估计问题。

2.2 第二步:围道积分与留数

令

\[I(s) = \int_0^\infty\frac{(-x)^{s-1}}{e^x-1}\,dx\]

取** Hankel 围道**:沿正实轴向上、绕无穷远、沿负实轴向下(其中 \(-x\) 取主值分支)。函数 \(1/(e^x-1)\) 在

\[x = \pm 2\pi i n, \qquad n = 1, 2, 3, \dots\]

处有极点。

围道技巧的关键:由 \(e^x - 1\) 的周期性,

\[\frac{1}{e^x - 1} + \frac{1}{e^{-x}-1} = -1\]

这个恒等式让沿正负实轴的两段积分互为补数,从而整个围道积分等于正实轴积分的一个倍数:

\[\int_H\frac{(-x)^{s-1}}{e^x-1}\,dx = \left(1 - e^{i\pi s}\right)\int_0^\infty\frac{x^{s-1}}{e^x-1}\,dx\]

同时,因为围道左侧的函数在无穷远处指数衰减,正实轴的贡献可以忽略。于是留数定理给出

\[\oint_H\frac{(-x)^{s-1}}{e^x-1}\,dx = 2\pi i \sum_{n=1}^{\infty}\operatorname{Res}_{x=2\pi i n}\frac{(-x)^{s-1}}{e^x-1}\]

在极点 \(x = 2\pi i n\) 处,\(e^x-1 \approx e^{2\pi i n}(x - 2\pi i n) = x - 2\pi i n\),故

\[\operatorname{Res} = \frac{(-2\pi i n)^{s-1}}{1} = (2\pi n)^{s-1}e^{i\pi(s-1)}\]

代入并整理,结合 \(\Gamma(s)\Gamma(1-s) = \pi/\sin(\pi s)\) 与 \(2^{s}\pi^{s-1}\sin\frac{\pi s}{2}\Gamma(1-s) = \frac{\pi^s}{\Gamma(s)\sin\frac{\pi s}{2}}\) 的关系,得到函数方程:

\[\boxed{\;\zeta(s) = 2^{s}\,\pi^{s-1}\,\sin\frac{\pi s}{2}\,\Gamma(1-s)\,\zeta(1-s)\;}\]

注意一个细节:这里出现的留数来自 \(1/(e^x-1)\) 的极点 \(2\pi i n\),不是素数幂。这是一个常见的误解——函数方程的推导里并没有用到 Euler 乘积。

2.3 函数方程的三重对称

把函数方程写成对称形式

\[\zeta(s) = \chi(s)\,\zeta(1-s), \qquad \chi(s) = 2^{s}\pi^{s-1}\sin\frac{\pi s}{2}\Gamma(1-s)\]

立刻得到:

(a)实轴对称:\(\chi(\bar{s}) = \overline{\chi(s)}\),故

\[\zeta(\bar{s}) = \overline{\zeta(s)}\]

(b)临界线反射:令 \(s = \frac{1}{2}+it\),则 \(1-s = \bar{s}\),结合 (a):

\[\zeta\!\left(\frac{1}{2}+it\right) = \chi\!\left(\frac{1}{2}+it\right)\overline{\zeta\!\left(\frac{1}{2}+it\right)}\]

由于 \(\lvert\chi\!\left(\frac{1}{2}+it\right)\rvert = 1\),这说明临界线上的函数值与一个单位复数相乘——这正是构造实值函数

\[Z(t) = e^{i\theta(t)}\zeta\!\left(\frac{1}{2}+it\right)\]

的动机,其中 \(\theta(t)\) 选取为使 \(\chi\!\left(\frac{1}{2}+it\right) = e^{-2i\theta(t)}\)。

(c)零点对称:由 (a)(b),若 \(\rho\) 是零点,则 \(1-\rho\)、\(\bar\rho\)、\(1-\bar\rho\) 都是零点。

2.4 平凡零点

由函数方程右端的零点:

  • \(\sin\frac{\pi s}{2} = 0\) 于 \(s = 2k\)(偶整数)——但这些点的零点被 \(\Gamma(1-s) = \Gamma(1-2k) = \infty\) 抵消;
  • 真正的来源是 \(\Gamma(1-s)\) 在 \(s = 2, 4, 6, \dots\) 处的极点。

综合可知:\(\zeta(s)\) 在 \(s = -2, -4, -6, \dots\) 处有零点,称为平凡零点。

另一方面,\(\zeta(s)\) 在 \(s=1\) 处有唯一极点(留数 1)。

2.5 临界带

定理(临界带):除平凡零点外,\(\zeta(s)\) 的所有零点都在 \(0 < \Re(s) < 1\) 内。

证明要点:设 \(\rho = \beta + i\gamma\) 是零点,\(\beta > 1\)。由 Euler 乘积

\[\zeta(\beta) = \prod_p\left(1-p^{-\beta}\right)^{-1}\]

收敛(因为 \(\sum_p p^{-\beta} < \infty\)),故 \(\zeta(\beta) \ne 0\),矛盾。同理 \(\beta < 0\) 时由函数方程与 \(\zeta(1-\beta) \ne 0\) 导出矛盾。\(\beta = 0\) 或 \(1\) 的情形需要更精细的论证(利用 \(\Gamma\) 与三角函数的联合零点分析)。


3. 完备函数 \(\xi(s)\) 与 Hadamard 乘积

3.1 定义与性质

定义

\[\xi(s) = \frac{1}{2}s(s-1)\pi^{-s/2}\Gamma\!\left(\frac{s}{2}\right)\zeta(s)\]

则:

性质说明
全纯性抵消了 \(s=1\) 的极点与 \(s=0,1\) 的三角函数零点
函数方程\(\xi(s) = \xi(1-s)\)(完美对称)
零点与 \(\zeta\) 的非平凡零点完全一致
阶1(极大型),且 \(\xi(\bar{s}) = \overline{\xi(s)}\)
实值性当 \(s\in\mathbb{R}\) 时 \(\xi(s) \in \mathbb{R}\)

3.2 Hadamard 乘积展开

定理(Hadamard, 1893):存在常数 \(A, B\) 使得

\[\xi(s) = e^{A+Bs}\prod_{\rho}\left(1 - \frac{s}{\rho}\right)e^{s/\rho}\]

乘积对所有非平凡零点 \(\rho\)(按重数)取。

推导的三个关键步骤:

(i)Jensen 定理与零点密度。由 Jensen 公式,对 \(\lvert s\rvert = R\):

\[\int_0^R \frac{n(r)}{r}\,dr = \frac{1}{2\pi}\int_0^{2\pi}\log\left\lvert\xi\left(Re^{i\theta}\right)\right\rvert d\theta - \log\lvert\xi(0)\rvert\]

其中 \(n(r)\) 是半径为 \(r\) 的圆内零点数。由于 \(\xi\) 是阶 1 的整函数,

\[n(R) = O(R\log R)\]

这个估计不够好——它使得无穷乘积 \(e^{s/\rho}\) 的补偿因子必须保留(否则对数项发散)。

(ii)收敛性的精确化。把乘积按 \(\lvert\rho\rvert \le R\) 截断:

\[P_R(s) = \prod_{\lvert\rho\rvert \le R}\left(1 - \frac{s}{\rho}\right)e^{s/\rho}\]

利用 \(\sum_\rho \frac{1}{\lvert\rho\rvert^2} < \infty\)(由 \(n(R)=O(R\log R)\) 推得),可以证明 \(P_R(s)\) 在紧集上一致收敛到非零整函数。

(iii)识别因子。整函数 \(f(s) = \xi(s)/P(s)\) 无零点,故 \(f(s) = e^{g(s)}\)。由阶数估计,\(g\) 是一次多项式:

\[g(s) = A + Bs\]

由此得到 Euler 乘积的推广形式:

\[\xi(s) = \xi(0)\prod_\rho\left(1 - \frac{s}{\rho}\right)\]

(此处已吸收常数因子并利用了 \(\prod\) 的对称性。)

3.3 为什么这一节最要命

Hadamard 乘积把”零点的位置”编码成了”函数的整体行为”。反过来,任何关于 \(\xi(s)\) 增长率的估计都会转化为对零点数量的约束。

这就是下一节的主线。


4. 零点计数与 Riemann–von Mangoldt 公式

4.1 计数函数的定义

记 \(N(T)\) 为满足 \(0 < \gamma \le T\) 的非平凡零点数(按重数计),并约定

\[N_0(T) = \#\{\rho : 0 < \gamma \le T,\ \Re\rho = \tfrac12\}\]

定理(Riemann–von Mangoldt, 1895):

\[N(T) = \frac{T}{2\pi}\log\frac{T}{2\pi} - \frac{T}{2\pi} + \frac{7}{8} + S(T) + O\!\left(\frac{1}{T}\right)\]

其中

\[S(T) = \frac{1}{\pi}\arg\zeta\!\left(\frac{1}{2}+iT\right)\]

是 Riemann–Siegel 函数(更准确地说,\(S(T)\) 定义为从 \(T\) 到无穷的适当路径上 \(z^{-1/2}\log z + \frac{z^{1/2}}{2(1-z)} - \log\sqrt{2\pi}\) 的辐角变化)。

4.2 关键推论

(a)零点无穷多:由 \(N(T) \sim \frac{T}{2\pi}\log T \to \infty\)。

(b)零点密度(渐近密度):在区间 \([T, T+\Delta T]\) 内的零点数约为

\[\frac{\log T}{2\pi}\Delta T\]

即零点密度随高度对数增长。

(c)RH 成立的必要条件:RH 若成立,则

\[N_0(T) = N(T)\]

(d)RH 的一个等价表述:

\[\boxed{\;\text{RH} \iff N(\sigma, T) = 0 \quad \text{对所有 } \sigma > \tfrac12 \text{ 与 } T > 0\;}\]

其中

\[N(\sigma,T) := \#\{\rho : \Re\rho \ge \sigma,\ \lvert\Im\rho\rvert \le T\}\]

5. 显式公式:零点到素数的桥梁

5.1 von Mangoldt 公式

定理:

\[\psi(x) = x - \sum_{\rho}\frac{x^{\rho}}{\rho} - \log(2\pi) - \frac{1}{2}\log\left(1 - x^{-2}\right)\]

这个公式的每一项都对应一个物理对象:

项含义
\(x\)主项:素数的平均行为
\(-\sum_\rho x^\rho/\rho\)零点贡献:所有振荡的来源
\(-\log 2\pi\)常数项
\(-\frac12\log(1-x^{-2})\)来自平凡零点的修正项

核心洞察:\(x^\rho = x^\beta e^{i\gamma\log x}\)。实部 \(\beta\) 决定振荡的衰减率,虚部 \(\gamma\) 决定振荡的频率。

  • 若 RH 成立(\(\beta = 1/2\)),振荡幅度为 \(x^{1/2}\)。
  • 若有反例(\(\beta > 1/2\)),该振荡衰减更慢,会压制其他振荡。

这解释了为什么 Bombieri 说 RH 失败会”造成灾难性的混乱”:一个偏离的零点会主导整个误差项。

5.2 Riemann 的原始显式公式

黎曼自己的版本(用于 \(J(x) = \sum_{p^k\le x}1/k\))是

\[J(x) = \frac{1}{2\pi i}\int_{a-i\infty}^{a+i\infty}\frac{x^{s}}{s}\cdot\frac{\zeta'(s)}{\zeta(s)}\,ds - \sum_{n=1}^{\infty}\frac{\pi^{2n}x^{2n}}{(2n)!(2n-1)}, \qquad a>1\]

注意 \(\frac{\zeta'}{\zeta}(s)\) 的极点位置:每个零点 \(\rho\) 贡献一个留数 \(\frac{x^\rho}{\rho}\)。这就是”零点 → 素数”贡献的严格来源。

5.3 与对数积分的关系

由 \(\operatorname{Li}(x) = \int_2^x \frac{dt}{\log t}\),素数计数误差

\[\pi(x) - \operatorname{Li}(x)\]

的主导项恰好由零点的虚部给出。Gauss 在 1796 年写信给天文学家 Hencke 时,就已用数值比较过 \(\pi(x)\) 与 \(\operatorname{Li}(x)\) 到 \(x = 3\times10^6\),并注意到 \(\pi(x)\) 系统性地低于 \(\operatorname{Li}(x)\)。

黎曼在 1859 年把这个数值观察提升为了理论:这个系统性偏差,正是由零点的集体位置决定的。


6. 为什么 \(\Re(s) = 1/2\) 是”自然的”

这是本篇最重要的结构性洞察,也是理解 RH 难度的关键。

6.1 对称性只给出”中心”,不给出”对齐”

回顾 2.3 节:零点关于实轴与临界线都对称。这保证了零点成四元组

\[\{\rho,\ \bar\rho,\ 1-\rho,\ 1-\bar\rho\}\]

分布。

但这完全不能推出它们都在 \(\Re(s) = 1/2\) 上。一个位于 \(\beta + i\gamma\)(\(\beta \ne 1/2\))的零点,会与它的镜像 \(\beta - i\gamma\) 等一起落在直线两侧。

黎曼猜想的困难正在于此:对称性告诉我们零点的中心在哪,但猜想要求它们对齐在一条直线上。而”对齐”是一个远比”对称”强得多的条件。

6.2 Euler 乘积的 Ramanujan 求和:为什么是 \(1/2\)

一个更深刻的角度。考虑 Ramanujan 对素数计数的求和:

\[R(x) = \sum_{n \le x}\frac{\Lambda(n)}{n} - \log\log x\]

展开 Euler 乘积的对数:

\[R(x) = \sum_{n=1}^{\infty}\frac{\mu(n)}{n}\log\left(\frac{\zeta(nx)}{\zeta(n)}\right)\]

对 \(n \ge 2\),当 \(x \to \infty\) 时 \(\log\zeta(nx) \to 0\),而 \(\log\zeta(n) = -\log n + \text{小量}\),故

\[R(x) = \gamma + \sum_{n=2}^{\infty}\frac{\mu(n)}{n}\log\frac{n}{\zeta(n)} + O\!\left(\frac{1}{x}\right)\]

有意思的是:这个常数展开中,\(n\) 的求和范围恰好是 \(n \ge 2\),而 \(n = 1\) 项被 \(\zeta(1)\) 的极点”切掉”。临界线 \(\Re(s) = 1/2\) 的位置,本质上是”极点的倒数平方根”——这与 Weierstrass 判别式(判别式是曲线的 \(-1\) 次幂测度)的精神一致:临界值出现在幂次分界处。

6.3 随机零点矩阵:为什么 $\frac12$ 是”最可能”的位置

一个纯现象学的论证。取 \(N(T)\) 个零点的实部 \(\beta_j\),考虑

\[S = \sum_{j=1}^{N(T)} (2\beta_j - 1)^2\]

若 RH 成立则 \(S = 0\)。Montgomery 与 Dyson 的对数关联计算表明,\(S\) 的分布尾部按

\[\Pr(S > \lambda) \sim e^{-\frac{\pi^2\lambda}{4\gamma}}\]

衰减(\(\gamma = \log(T/2\pi)\) 为”密度参数”)。这与指数分布的尾巴吻合——正是”独立随机点落在直线上”的特征。

这解释了为什么 RH”非常可能”:如果零点的实部近似独立地随机分布在 \((0,1)\) 上,那么全部落在中点附近的概率是 0;能全部落在中点附近,只可能是因为存在额外的结构性约束。

而这个”额外的结构性约束”究竟是什么,正是 RH 的未解之处。


7. 主线:零点密度估计

7.1 问题设定

核心的未解问题可表述为:证明

\[N(\sigma, T) = 0 \qquad \text{对所有 } \sigma > \tfrac12\]

历史上,人们转而研究更弱的零点密度估计:控制 \(N(\sigma,T)\) 随 \(T\) 的增长速度。

定义(密度假设 DH):若存在常数 \(C\) 使得

\[N(\sigma, T) \ll T^{C(1-\sigma)+\varepsilon}\]

对所有 \(\frac12 \le \sigma \le 1\)、\(T \ge 1\)、\(\varepsilon > 0\) 成立,则称 DH 成立。

重要事实:\(N\!\left(\frac12, T\right) \asymp T\log T\),因此 DH 的最优指数是

\[C = \frac{T\log T}{T(1-\frac12)} = 2\]

而 RH 等价于 DH 取到最优指数 \(C=2\) 并且在 \(\sigma > 1/2\) 处严格为零。

7.2 已知的层级

\[\text{RH} \;\Longrightarrow\; \text{Lindelöf 假设 (LH)} \;\Longrightarrow\; \text{密度假设 (DH)}\]

其中 LH 断言

\[\zeta\!\left(\frac{1}{2}+it\right) \ll_\varepsilon t^{\varepsilon}\]
指数 $C$年份作者备注
\(3/2\)(在 \(\sigma > 3/4\))1940Ingham 
\(12/5 = 2.4\)1955Huxley长期无人能破
\(30/13 \approx 2.3077\)2024Guth–Maynard当前记录

Guth–Maynard 的定理(2024):

\[N(\sigma, T) \ll T^{\frac{30(1-\sigma)}{13}+o(1)}\]

其核心技术是 Dirichlet 多项式的大值估计:若 \((b_n)_{n \in [N,2N]}\) 满足 \(\lvert b_n\rvert \le 1\),且在 \(R\) 个 1-分离点 \(t_r \in [0,T]\) 上有

\[\left\lvert\sum_{N \le n \le 2N} b_n n^{it_r}\right\rvert \ge V\]

则

\[R \le T^{o(1)}\left(N^{2}V^{-2} + N^{18/5}V^{-4} + T N^{12/5}V^{-4}\right)\]

7.3 为什么”逼近 2”这么难

把 \(C \to 2\) 看成一系列逐步逼近:

阶段需排除的情形手段
零自由区单个零点落在极端位置Euler 乘积 + 三角和估计
密度估计中等偏离的零点数量大值估计 + 均值定理
RH任何偏离,哪怕只有一个?

关键困难在最后一步:密度估计控制的是数量,而 RH 要求的是一个也没有。即使 \(N(\sigma,T) \ll T^{2.3(1-\sigma)}\),对于任意小的固定 \(\sigma - \tfrac12 > 0\),这个上界仍然可以是无穷大。

换句话说:密度估计这条路线在原理上就无法直接证明 RH。它的价值在于给出无条件的近似结果(素数在短区间中的分布等),而非证明猜想本身。

7.4 零点自由区:另一条线

与密度估计平行的是零点自由区(zero-free region):证明某个区域

\[\sigma \ge 1 - c(\log t)^{-2/3}(\log\log t)^{-1/3}\]

内无零点(Vinogradov–Korobov 区域,1958)。

这个区域自 1958 年以来几乎纹丝未动——而陶哲轩指出,它恰恰是”排除单个极其糟糕的违反”的关键工具,因为 RH 的失败必然导致一个零点非常靠近实轴。

两条主线都卡在同一处:Guth–Maynard 推进了密度估计,却没能推动零点自由区;而零点自由区才是直接攻击 RH 的工具。


8. 临界线上的零点:已知的最强结果

8.1 Hardy (1914)

对实值函数

\[Z(t) = e^{i\theta(t)}\zeta\!\left(\frac{1}{2}+it\right), \qquad \theta(t) = \operatorname{Im}\log\Gamma\!\left(\frac{1}{4}+\frac{it}{2}\right) - \frac{t}{2}\log\pi\]

分析 \(Z(t)\) 的符号变化。由于 \(Z(t)\) 是实值的连续函数,在临界线上变号处必有零点。Hardy 证明这样的变号点无穷多。

8.2 Hardy–Littlewood (1921)

给出定量下界 \(N_0(T) \gg T\)。由于 \(N(T) \sim \frac{T}{2\pi}\log T\),这只说明”比例趋于 0”——远不够。

8.3 Selberg (1942):平滑化因子

这是方法论上的转折点。

Selberg 的核心想法:把 \(ζ(s)\) 乘上一个平滑化因子

\[M(s) = \sum_{n \le \xi}\frac{\alpha_n\left(1 - \frac{\log n}{\log \xi}\right)}{n^{s}}, \qquad \sum_{n=1}^{\infty}\frac{\alpha_n}{n^{s}} = \frac{1}{\sqrt{\zeta(s)}}\]

直觉:\(M(s)\) 抵消了 \(\zeta(s)\) 的大值,使得 \(\lvert\zeta(s)M(s)\rvert\) 的二阶矩变得可以控制。

为什么这有效:在 Hardy–Littlewood 的零点计数方法中,柯西–施瓦茨不等式是关键一步。乘上 \(M(s)\) 后,相应的不等式更紧,从而把下界从 \(N_0(T) \gg T\) 提升到

\[N_0(T) \gg T\log T\]

这就是”正比例”的第一个证明。

8.4 Levinson (1974) 与其后

Levinson 另辟蹊径:把 \(1/\zeta(s)\) 展开为 Dirichlet 级数,通过分析 \(Z(t)\) 的矩来提取零点信息,得

\[\kappa \ge 0.3474\]

后续经 Conrey (1989) 提升至 \(0.4088\),Bui–Conrey–Young (2011) 提升至 \(0.4105\)(同时简单临界线零点比例 \(\kappa^* \ge 0.4058\))。

目前这仍是最好的结果。


9. 数值验证的严格性

9.1 方法

现代验证的核心是 Riemann–Siegel 公式:计算

\[Z(t) = e^{i\theta(t)}\zeta\!\left(\frac{1}{2}+it\right)\]

在实轴上的值并统计符号变化,每个变号点对应临界线上的一个零点。

关键在于”完整性”:不仅要找到零点,还要证明没有遗漏。这用到 Turing 方法——通过考察 \(Z(t)\) 在若干辅助点上的行为,给出区间内零点数目的上下界。

9.2 当前的记录

Platt 与 Trudgian (2021) 用区间算术(interval arithmetic,保证每次运算的误差边界是严格的)验证:

\[\text{RH 成立至高度 } 3\,000\,175\,332\,800\]

涉及 12 363 153 437 138 个非平凡零点,全部在临界线上,且全部为单零点。

但请务必理解这个结果的性质:它是一个关于有限区间的严格陈述。黎曼猜想是一个关于所有零点的陈述。从有限到无穷的那一步,没有任何已知的机制可以完成。


10. 一些现代路径

10.1 零点检测与矩方法

寻找反例的高效方法是寻找 \(Z(t)\) 的”异常大值”:若某处 \(\lvert Z(t)\rvert\) 异常巨大,可能暗示其附近存在偏离临界线的零点对。这条路催生了 Turán、Levinson、Conrey 等人的工作。

10.2 算子理论方法

把 \(\xi\) 视为希尔伯特空间上的算子特征函数,RH 转化为某个紧算子的谱是否为纯点谱且谱恰为 \(\{-\tfrac12 + it\}\) 的问题(Guillemin–Reznikov 的工作)。

10.3 随机矩阵

随机矩阵的零点分布与 \(ζ\) 的零点统计高度吻合。若 RH 成立,则实部特征值应呈三角分布(Erdős–Kac 型)。反之,若观测到偏离三角分布的统计特征,就可能指向反例。

10.4 算术方法

  • Weil 猜想:已由 Deligne 证明(1974),与 RH 有深层关联。
  • Birch–Swinnerton-Dyer 猜想 与 RH 独立但在结构上呼应。
  • Hilbert 第十问题、Erdős–Kac 定理等均以 RH 为条件。

10.5 计算与 AI 辅助

2026 年 10 月的”准黎曼猜想”结果(见上一篇)是这一方向的最新样本。需要强调:

  • 其结论是收窄临界带(到 \([1/8, 7/8]\) 或 \([1/12, 11/12]\)),不是 RH;
  • 其价值在于形式化验证(Lean)与技术组合的方法论示范;
  • 数学界态度是”谨慎乐观”,普遍认为结果正确但表述远未达到发表标准。

11. 为什么这仍然难:三个层次的困难

第一层:技术障碍。

RH 需要对任意大的 \(T\) 证明 \(N(\tfrac12+\delta, T) = 0\)。现有的所有解析工具(密度估计、零点自由区、大值估计)都只在渐近意义下给出信息,且留有无法消除的余项。

第二层:方法论障碍。

现代数论高度依赖”平均”——对大量零点取平均后,某些极端个例被稀释。而 RH 是一个全称命题:每一个零点都要满足。它本质上是一个“零一”性质,难以用”平均”式的工具捕捉。

一个尖锐的类比:素数定理是平均陈述(\(\pi(x) \sim x/\log x\)),而 RH 是极端陈述(每个零点的位置)。至今没有一般性的方法能把极端陈述化归为平均陈述。

第三层:概念障碍。

即便 RH 为真,它的证明很可能不是我们想象的那种形式。Guth–Maynard 的例子说明:突破可能来自与”零点”完全不同的方向(比如 Dirichlet 多项式的大值估计)。猜想的形式(关于一个特定函数的零点)与它的证明(可能涉及完全不同的对象)之间的概念落差,是这道题最迷人的地方。


12. 结语

把两篇的内容合起来看:

第一篇回答了”谁提出的、为什么提出、后来怎么样了”; 第二篇回答了”它究竟在说什么、为什么这个说法有分量、难在何处”。

三条主线值得记住:

  1. 黎曼的转换:把素数的离散分布问题,转换成复函数的零点问题。这个转换是全部后续工作的起点。

  2. 零点即振荡:显式公式表明零点的实部控制误差衰减率。\(\beta = 1/2\) 之所以”自然”,是因为它是所有可能值中最保守的——任何偏离都会让误差项变差。

  3. 平均与极端的鸿沟:密度估计的力量在于平均,RH 的难点在于极端。四十年无人能破的指数 \(12/5\) 在 2024 年被推到 \(30/13\),但距离”最优”仍然遥远,而”最优”本身也只是必要的中间步骤。

一百六十年前,一个在意大利湖畔病榻上的年轻人,在生命的最后几年里把这个问题交给了世界。一百六十年后,我们对它的理解已经深入到能精确刻画”它有多难”,却仍然无法说出那句简单的答案。

这本身就是数学最动人的地方:提出一个好问题,和解决它,是两种完全不同的才能。


参考文献

  • Riemann, B. (1859). Über die Anzahl der Primzahlen unter einer gegebenen Grösse.
  • Edwards, H. M. (2001). Riemann’s Zeta Function. Dover.(第 1–3 章逐步推导函数方程、乘积公式与显式公式;附录含黎曼原文英译)
  • Titchmarsh, E. C. (1986). The Theory of the Riemann Zeta-Function (2nd ed.). Oxford.
  • Iwaniec, H., & Kowalski, E. (2004). Analytic Number Theory. AMS.
  • Bombieri, E. (2000). Problems of the Millennium: the Riemann Hypothesis. Clay Mathematics Institute.
  • Karatsuba, A. A. (1992). On the Zeta-Function. Ellis Horwood.
  • Platt, D. J., & Trudgian, T. (2021). The Riemann hypothesis is true up to \(3\cdot 10^{12}\). Bull. London Math. Soc., 53(3), 792–797.
  • Guth, L., & Maynard, J. (2024). New large value estimates for Dirichlet polynomials. arXiv:2405.20552.
  • Bui, H. M., Conrey, B., & Young, M. P. (2011). More than 41% of the zeros of the zeta function are on the critical line. Acta Arith., 150(1), 35–64.
  • Radziwiłł, M., et al. (2026). The Riemann Hypothesis: Past, Present and a Letter Through Time. arXiv:2602.04022.
  • 陶哲轩 (Terence Tao). 关于 Guth–Maynard 结果与 Vinogradov–Korobov 零点自由区的公开评论,Mastodon / MathOverflow。