MathLabs

几何学

黎曼几何

研究在每点都配备长度和角度度量方式的空间上的微分几何,是广义相对论的数学语言。

直观直觉:没有直尺的几何学

设想一只生活在丘陵地形上的蚂蚁,或一枚在弯曲时空中航行的火箭:它们都没有直尺,也没有平坦的方格纸,但无论站在哪里,都仍能测量短距离和小角度。黎曼几何仅凭这种平滑地附着在空间每一点上的局部、无穷小的量尺,构建出关于长度、角度、笔直性与曲率的整套理论。它是广义相对论——时空本身是弯曲的——背后的数学语言,也是从机械臂到机器学习模型等一切自然配置空间是弯曲而非平坦的系统背后的语言。

展示负曲率的三维鞍形曲面图。
鞍形曲面 z=x2−y2z = x^2 - y^2 是负曲率空间最简单的局部模型:在这里以平行方向出发的测地线会越来越快地相互远离,这与正曲率球面上发生的情况恰好相反。

大学黎曼度量:在每一点测量长度和角度

定义: 黎曼度量

光滑流形 MM 上的黎曼度量gg在每一点光滑地给出切空间上的一个内积,从而可以计算切向量的长度以及它们之间的夹角。在局部坐标 x1,…,xnx^1,\dots,x^n 下,它由一个对称正定的函数矩阵gijg_{ij}编码,无穷小距离的平方写作ds2=gij dxidxjds^2 = g_{ij}\,dx^i dx^j(对重复指标使用爱因斯坦求和约定)。

ds2=gij dxidxjds^2 = g_{ij}\,dx^i dx^j

曲线γ\gamma从 aa 到 bb 的长度即为L(γ)=∫abgγ(t)(γ′(t),γ′(t)) dtL(\gamma) = \int_a^b \sqrt{g_{\gamma(t)}(\gamma'(t),\gamma'(t))}\,dt:这与把无穷小的直线段加起来的想法完全一致,只是每一段是用局部的量尺gg来测量,而不是用固定的全局量尺。这一个公式,在gg为常数时给出普通的欧几里得长度,而根据代入的度量不同,还能给出球面上的大圆距离、庞加莱圆盘中的双曲距离,或相对论中沿轨迹的固有时。

L(γ)=∫abgγ(t)(γ′(t),γ′(t)) dtL(\gamma) = \int_a^b \sqrt{g_{\gamma(t)}(\gamma'(t),\gamma'(t))}\,dt

大学Levi-Civita联络与测地线

定义: Levi-Civita联络

要在弯曲空间中沿曲线对向量场求导,需要一个比较邻近点处切向量的规则,称为联络∇\nabla。Levi-Civita联络是唯一无挠且与度量相容(平行移动保持长度与角度)的联络。在坐标中,它由克里斯托费尔符号Γijk=12gkl(∂igjl+∂jgil−∂lgij)\Gamma^k_{ij} = \frac{1}{2}g^{kl}\left(\partial_i g_{jl} + \partial_j g_{il} - \partial_l g_{ij}\right)给出,该式纯粹由度量及其一阶导数构成,其中 gklg^{kl} 是 gklg_{kl} 的逆矩阵。

Γijk=12gkl(∂igjl+∂jgil−∂lgij)\Gamma^k_{ij} = \frac{1}{2}g^{kl}\left(\partial_i g_{jl} + \partial_j g_{il} - \partial_l g_{ij}\right)

测地线是能够平行移动自身速度的曲线,也就是说,它在给定度量下 尽可能笔直 地前进;它是直线的自然推广,并在局部上使两端点之间的长度最短。用坐标写出,这一条件就是测地线方程d2xkdt2+Γijkdxidtdxjdt=0\dfrac{d^2 x^k}{dt^2} + \Gamma^k_{ij}\dfrac{dx^i}{dt}\dfrac{dx^j}{dt} = 0,这是一个二阶常微分方程组,其解由初始点和初始速度唯一确定,恰如牛顿第二定律,只是克里斯托费尔符号在此充当一个依赖速度的力。

d2xkdt2+Γijkdxidtdxjdt=0\frac{d^2 x^k}{dt^2} + \Gamma^k_{ij}\frac{dx^i}{dt}\frac{dx^j}{dt} = 0

进阶曲率:黎曼曲率、里奇曲率与数量曲率

定义: 黎曼曲率张量

黎曼曲率张量R(X,Y)Z=∇X∇YZ−∇Y∇XZ−∇[X,Y]ZR(X,Y)Z = \nabla_X\nabla_Y Z - \nabla_Y\nabla_X Z - \nabla_{[X,Y]}Z衡量二阶协变导数不可交换的程度,即一个向量沿由 XX 和 YY 张成的小闭合回路平行移动后会发生多大变化。将它再与两个向量缩并,便得到截面曲率K(X,Y)=⟨R(X,Y)Y,X⟩∣X∣2∣Y∣2−⟨X,Y⟩2K(X,Y) = \dfrac{\langle R(X,Y)Y,X\rangle}{|X|^2|Y|^2 - \langle X,Y\rangle^2},即与由 X,YX,Y 张成的平面相切的曲面的二维切片的高斯曲率——它是把普通微分几何中曲面曲率直接推广到任意维度的结果。

R(X,Y)Z=∇X∇YZ−∇Y∇XZ−∇[X,Y]ZR(X,Y)Z = \nabla_X\nabla_Y Z - \nabla_Y\nabla_X Z - \nabla_{[X,Y]}Z

定义: 里奇曲率与数量曲率

里奇曲率Ric(Y,Z)=tr(X↦R(X,Y)Z)\mathrm{Ric}(Y,Z) = \mathrm{tr}\big(X \mapsto R(X,Y)Z\big)通过对黎曼张量的一对指标求迹,将包含某个固定方向的所有二维平面的截面曲率取平均;它和度量本身一样,是一个对称双线性型。再用度量求一次迹,便得到数量曲率R=trg(Ric)R = \mathrm{tr}_g(\mathrm{Ric}),它是每一点上把各个方向的平均曲率汇总成的一个数。里奇曲率正是广义相对论爱因斯坦场方程中直接出现的量,其符号决定了小球体积相对于平坦空间是变大还是变小。

Ric(Y,Z)=tr(X↦R(X,Y)Z),R=trg(Ric)\mathrm{Ric}(Y,Z) = \mathrm{tr}\big(X \mapsto R(X,Y)Z\big), \qquad R = \mathrm{tr}_g(\mathrm{Ric})
nn 维黎曼流形的三种曲率
曲率对象类型衡量的内容
截面曲率 K(X,Y)K(X,Y)每个二维平面对应一个数与该平面相切的二维切片的高斯曲率
里奇 Ric(Y,Z)\mathrm{Ric}(Y,Z)一个对称双线性型沿某方向的所有二维平面上截面曲率的平均值
数量 RR每一点上的一个数所有方向上的平均曲率,例如常曲率 kk 时 R=n(n−1)kR=n(n-1)k

进阶关键定理:Hopf–Rinow 与 Myers

设 (M,g)(M,g) 为连通的黎曼流形。以下条件等价:(i) MM 测地完备,即每条测地线都能延拓为定义在整个 R\mathbb{R} 上的解;(ii) MM 关于由 gg 诱导的距离是完备度量空间;(iii) MM 的每个闭且有界的子集都是紧的。此外,当这些条件成立时,MM 中任意两点都可由一条最小测地线连接。

为什么成立?

完备性是局部(只能保证测地线在短时间内存在的常微分方程)与整体(一个性质良好的空间,其中距离的表现就像普通空间中一个闭且有界的区域,并且你凭直觉画出的最短路径总是确实存在)之间的桥梁。

证明

(i) ⇒\Rightarrow 最小测地线的存在性:固定 p,q∈Mp,q \in M,令 r=d(p,q)r = d(p,q)。小的度量球面 Sε(p)S_\varepsilon(p) 是紧的,因此连续函数 x↦d(x,q)x \mapsto d(x,q) 在其上某点 x0=γ(ε)x_0 = \gamma(\varepsilon) 处取得最小值,其中 γ\gamma 是从 pp 经过 x0x_0 的单位速度测地线(由测地完备性,它对所有时间都有定义)。可以证明 d(γ(ε),q)=r−εd(\gamma(\varepsilon),q) = r-\varepsilon,再将其推广到集合 A={t∈[0,r]:d(γ(t),q)=r−t}A = \{t \in [0,r] : d(\gamma(t),q) = r-t\}:由连续性,AA 非空且闭;对满足 t0<rt_0<r 的任意 t0∈At_0 \in A,在 γ(t0)\gamma(t_0) 处重复同样的最小球面论证,可证明 AA 在 [0,r][0,r] 中也是开的,故 A=[0,r]A=[0,r],从而 γ(r)=q\gamma(r)=q。

(i) ⇒\Rightarrow (ii):对 MM 中的柯西列 (xn)(x_n),上述论证表明任意两个足够接近的点都可由一条长度等于其间距离的最小测地线相连;固定点周围的闭度量球是指数映射下 TpMT_pM 中闭球的连续像,因而是紧的,所以柯西列最终落在某个紧集内并收敛。

(ii) ⇒\Rightarrow (iii):度量完备性连同局部紧性(每个黎曼流形都是局部紧的)蕴含闭球是完全有界且完备的,因而根据度量空间中紧性的标准刻画是紧的;一般的闭且有界集合是某个闭球的闭子集,因而作为紧集的闭子集也是紧的。

(iii) ⇒\Rightarrow (i):假设单位速度测地线 γ\gamma 只在极大区间 [0,T)[0,T)(其中 T<∞T<\infty)上有定义。当 t→T−t \to T^- 时,γ(t)\gamma(t) 始终位于以 γ(0)\gamma(0) 为中心、半径为 TT 的闭球内,由 (iii) 知该闭球是紧的;于是满足 tn→Tt_n \to T 的序列 γ(tn)\gamma(t_n) 有收敛子列,而测地线常微分方程(具有光滑系数的二阶方程组)的标准存在性理论表明,解可以光滑地延拓超过 TT,这与 TT 的极大性矛盾。因此每条测地线都能延拓到整个 R\mathbb{R}。

定理: Myers定理

设 (Mn,g)(M^n,g) 是完备黎曼流形,满足 Ric≥(n−1)k g\mathrm{Ric} \geq (n-1)k\,g(其中 k>0k>0 为常数,即每个方向的里奇曲率都至少为 (n−1)k(n-1)k)。那么 MM 是紧的,diam(M)≤πk\mathrm{diam}(M) \leq \dfrac{\pi}{\sqrt{k}},且其万有覆盖也是紧的,因此 π1(M)\pi_1(M) 是有限的。

为什么成立?

正的里奇曲率意味着测地线平均而言会相互汇聚,就像在球面上一样;如果这种汇聚在每个方向上都足够强,一条测地线就不可能永远保持距离最短——它最终会遇到一个共轭点,附近的测地线在那里追上它,这就限制了任意两点之间可能的最大距离。

证明

根据Hopf–Rinow定理,由于 MM 完备,任意两点都可由一条最小测地线相连。为反证,假设某条单位速度最小测地线 γ:[0,L]→M\gamma:[0,L]\to M 的长度满足 L>π/kL > \pi/\sqrt{k}。

由于 γ\gamma 使长度最小,对沿 γ\gamma 且在两端点处为零的任意分段光滑向量场 VV,第二变分(指标型)满足 I(V,V)≥0I(V,V) \geq 0,其中 I(V,V)=∫0L(∣V′∣2−⟨R(V,γ′)γ′,V⟩)dtI(V,V) = \int_0^L \left(|V'|^2 - \langle R(V,\gamma')\gamma', V\rangle\right)dt。沿 γ\gamma 选取一个处处与 γ′\gamma' 正交的平行标准正交标架 E1,…,En−1E_1,\dots,E_{n-1},并用向量场 Vi(t)=sin⁡(πt/L) Ei(t)V_i(t) = \sin(\pi t/L)\,E_i(t)(i=1,…,n−1i=1,\dots,n-1)进行检验。

由于 EiE_i 是平行的,Vi′=πLcos⁡(πt/L)EiV_i' = \frac{\pi}{L}\cos(\pi t/L)E_i,故 I(Vi,Vi)=∫0L[(πL)2cos⁡2(πt/L)−sin⁡2(πt/L)⟨R(Ei,γ′)γ′,Ei⟩]dtI(V_i,V_i) = \int_0^L\left[\left(\frac{\pi}{L}\right)^2\cos^2(\pi t/L) - \sin^2(\pi t/L)\langle R(E_i,\gamma')\gamma',E_i\rangle\right]dt。对 ii 求和并利用 ∑i⟨R(Ei,γ′)γ′,Ei⟩=Ric(γ′,γ′)≥(n−1)k\sum_i \langle R(E_i,\gamma')\gamma',E_i\rangle = \mathrm{Ric}(\gamma',\gamma') \geq (n-1)k,得到 ∑iI(Vi,Vi)≤(n−1)∫0L[(πL)2cos⁡2(πt/L)−ksin⁡2(πt/L)]dt\sum_i I(V_i,V_i) \leq (n-1)\int_0^L\left[\left(\frac{\pi}{L}\right)^2\cos^2(\pi t/L) - k\sin^2(\pi t/L)\right]dt。

由于 ∫0Lcos⁡2(πt/L) dt=∫0Lsin⁡2(πt/L) dt=L/2\int_0^L\cos^2(\pi t/L)\,dt = \int_0^L\sin^2(\pi t/L)\,dt = L/2,右边等于 (n−1)⋅L2[(πL)2−k](n-1)\cdot\frac{L}{2}\left[\left(\frac{\pi}{L}\right)^2 - k\right],只要 L>π/kL > \pi/\sqrt{k} 该式就严格为负。于是某个 I(Vi,Vi)<0I(V_i,V_i) < 0,与 γ\gamma 使长度最小矛盾。因此每条最小测地线都满足 L≤π/kL \leq \pi/\sqrt{k},从而 diam(M)≤πk\mathrm{diam}(M) \leq \dfrac{\pi}{\sqrt{k}};由Hopf–Rinow定理,直径有限的完备流形是紧的,把同样的论证用于万有覆盖(它继承相同的里奇曲率下界)可知它也是紧的,从而 π1(M)\pi_1(M) 必为有限。

进阶实际应用与典型例题

黎曼几何是广义相对论的语言:爱因斯坦场方程把时空的里奇曲率(经洛伦兹符号调整后)与其中的物质和能量联系起来,而自由下落的物体沿着这个弯曲度量的测地线运动。机器人学与运动规划把关节角度或形状空间当作黎曼流形,以避开障碍物、使能量最小化的测地线来规划运动。信息几何在统计模型空间上放置黎曼度量(Fisher信息度量),由此产生的测地线给出机器学习中自然梯度下降之类算法所用的最短统计路径。计算机视觉中的形状分析把物体轮廓当作黎曼流形上的点,用测地距离来比较不同的形状。

例题: 庞加莱上半平面的常曲率

庞加莱上半平面模型的度量为 ds2=λ(x,y)2(dx2+dy2)ds^2 = \lambda(x,y)^2(dx^2+dy^2),其中共形因子 λ(x,y)=1/y\lambda(x,y) = 1/y(y>0y>0)。利用平面上共形度量高斯曲率的标准公式 K=−1λ2Δ(ln⁡λ)K = -\dfrac{1}{\lambda^2}\Delta(\ln\lambda),计算 KK。

解答

写出共形因子的对数:由于 λ(x,y)=1/y\lambda(x,y) = 1/y,有 ln⁡λ=ln⁡(1/y)=−ln⁡y\ln\lambda = \ln(1/y) = -\ln y,这是一个仅依赖于 yy 的函数。

对 yy 求两次导数:∂y(−ln⁡y)=−1/y\partial_y(-\ln y) = -1/y,再求 ∂y2(−ln⁡y)=1/y2\partial_y^2(-\ln y) = 1/y^2。由于 ln⁡λ\ln\lambda 不依赖 xx,拉普拉斯算子为 Δ(ln⁡λ)=∂x2(ln⁡λ)+∂y2(ln⁡λ)=0+1/y2=1/y2\Delta(\ln\lambda) = \partial_x^2(\ln\lambda) + \partial_y^2(\ln\lambda) = 0 + 1/y^2 = 1/y^2。

代入曲率公式:K=−1λ2Δ(ln⁡λ)K = -\dfrac{1}{\lambda^2}\Delta(\ln\lambda),取 λ2=1/y2\lambda^2 = 1/y^2,得到 K=−y2⋅1y2=−1K = -y^2 \cdot \dfrac{1}{y^2} = -1。

yy 的依赖性恰好完全抵消,因此在上半平面上处处 K=−1K=-1:双曲平面具有恒定的负曲率,不同于球面(恒定正曲率)或平面(曲率为零)。

例题: 球面的里奇曲率与数量曲率,及Myers定理检验

利用截面曲率为常数 kk 的空间具有黎曼张量 Rijkl=k(gikgjl−gilgjk)R_{ijkl} = k(g_{ik}g_{jl}-g_{il}g_{jk}) 这一事实,计算单位球面 S2S^2 的里奇曲率与数量曲率,然后检验Myers定理预测的直径上界。

解答

半径为 11 的圆球面 S2S^2 具有常截面曲率 k=1k=1,因此其黎曼张量取常曲率形式 Rijkl=k(gikgjl−gilgjk)R_{ijkl} = k(g_{ik}g_{jl}-g_{il}g_{jk})(n=2n=2)。

求一次迹得到里奇张量 Ric=(n−1)k g\mathrm{Ric} = (n-1)k\,g:当 n=2n=2、k=1k=1 时为 Ric=(2−1)⋅1⋅g=g\mathrm{Ric} = (2-1)\cdot 1\cdot g = g,故每个切方向的里奇曲率恰为 11。

再求一次迹得到数量曲率 R=n(n−1)kR = n(n-1)k:当 n=2n=2、k=1k=1 时为 R=2⋅1⋅1=2R = 2\cdot 1\cdot 1 = 2。

由于 k=1k=1 时 Ric≥(n−1)k g\mathrm{Ric} \geq (n-1)k\,g 成立,Myers定理给出的 diam(M)≤πk\mathrm{diam}(M) \leq \dfrac{\pi}{\sqrt{k}} 变为 diam(S2)≤π/1=π\mathrm{diam}(S^2) \leq \pi/\sqrt{1} = \pi。这恰好吻合:沿大圆测量,单位球面上两个对径点之间的测地距离恰好是 π\pi。

某完备黎曼流形满足 Ric≥(n−1)k g\mathrm{Ric} \geq (n-1)k\,g(k=4k=4)。Myers定理给出的直径上界是多少?

哪个应用数学领域在统计模型空间上放置黎曼度量(Fisher信息度量),并在自然梯度下降等优化算法中使用其测地线?

某三维空间具有常截面曲率 k=2k=2。利用 R=n(n−1)kR = n(n-1)k,其数量曲率 RR 是多少?

在黎曼流形上所有可能的联络中,唯一刻画Levi-Civita联络 ∇\nabla 的两个性质是什么?

参考文献

  1. Manfredo P. do Carmo (1992). Riemannian Geometry
  2. Grigori Perelman (2002). The entropy formula for the Ricci flow and its geometric applications · arXiv:math/0211159 [预印本,未经同行评审]
  3. Peter Topping (2006). Lectures on the Ricci Flow