系列 · 泛函分析 · 第 9 篇

泛函分析(九):无界算子 —— 当有界性失效时

闭算子、对称与自伴的区别、亏指数、Friedrichs 扩张、无界自伴算子的谱定理以及 Stone 定理。

当算子拒绝在整个空间上定义#

我第一次写下“位置算子 $X$$L^2(\mathbb{R})$ 上”这句话时没多想——直到我试图把它平方。$(Xf)(x) = xf(x)$ 看起来人畜无害,但 $f(x) = 1/(1+|x|)$ 满足 $f \in L^2(\mathbb{R})$ ,而 $Xf(x) = x/(1+|x|)$ 在无穷远处不衰减为零,根本不在 $L^2$ 中。也就是说,$X$$L^2$ 中的某些元素送出 $L^2$ 之外。它甚至不是“有界但范数大”——它是“在某些点根本没定义”。整个之前关于有界算子的理论在这里失效。

这个尴尬在物理里几乎是常态。位置算子 $X$ 、动量算子 $P = -i\hbar \, d/dx$ 、Hamilton 算子 $H = -\Delta + V$ 、Laplace 算子——量子力学中所有“真正的”可观测量都不在整个 $L^2$ 上定义。它们只在某个稠密子空间 $\mathcal{D}(T) \subset L^2$ 上有意义。这个子空间不是装饰,而是数据的一部分:换一个 $\mathcal{D}(T)$ ,谱可能完全不同(比如 $-d^2/dx^2$$[0,1]$ 上配 Dirichlet 还是 Neumann 给出不同的谱),物理意义可能完全不同(比如自由粒子还是被困粒子)。

这一篇要把无界算子的理论写明白:定义域为什么是数据的一部分、闭性如何替代有界性、对称性和自共轭性的真实差距是什么、亏指数如何检测自共轭延拓的存在、Friedrichs 延拓如何从能量形式自动构造一个自共轭算子、Kato-Rellich 微扰定理如何处理 Schrödinger 算子。最重要的是:谱定理在无界设定下依然成立。Stone 定理把自共轭算子和强连续幺正群一一对应起来——这是量子动力学的数学骨架。

历史背景很重要。冯·诺伊曼在20世纪20年代末发展了无界算子理论,正是为了给量子力学奠定坚实的基础。当时薛定谔方程和海森堡矩阵力学已经在使用,但还没有一个清晰的“自伴”概念来描述微分算子——也没有适用于这类算子的谱定理——因此数学基础还不清楚。冯·诺伊曼的工作(1929-1932年发表,最终成果是他的《量子力学的数学基础》)通过引入闭算子、亏指标和自伴延拓等概念解决了这个问题。这里发展的理论基本上就是他的框架,经过后人的整理和完善,但核心思想没有改变。

无界算子的挑战#

前两篇文章中,我谈到了谱理论是量子力学的线性代数基础。问题在于,物理学家真正关心的几乎所有算子——位置算子、动量算子、拉普拉斯算子、薛定谔哈密顿算子——都是无界的。它们没有定义在整个希尔伯特空间上。它们只在某些稠密子集上定义,这些子集依赖于输入函数的正则性和衰减性质。之前的谱分析工具无法直接应用。我们需要扩展这些工具。

这种扩展非常微妙。对于无界算子,简单地写 “$T = T^*$ ” 已经不再明确,因为两边可能有不同的定义域。对称算子(在共同定义域上满足 $\langle Tx, y\rangle = \langle x, Ty\rangle$ )和自伴算子(此外还要求 $T$ 的定义域等于 $T^*$ 的定义域)之间有明显的区别。对于有界算子,这两者是一致的;但对于无界算子,它们在细微之处有所不同。这个差距正是数学物理中的大部分困难所在。如果处理得当,谱定理、泛函演算和斯通定理都可以扩展——这样我们就能严格地进行量子力学研究。

历史背景很重要。冯·诺伊曼在20世纪20年代末发展了无界算子理论,正是为了给量子力学奠定坚实的基础。当时薛定谔方程和海森堡矩阵力学已经在使用,但还没有一个清晰的“自伴”概念来描述微分算子——也没有适用于这类算子的谱定理——因此数学基础还不清楚。冯·诺伊曼的工作(1929-1932年发表,最终成果是他的《量子力学的数学基础》)通过引入闭算子、亏指标和自伴延拓等概念解决了这个问题。这里发展的理论基本上就是他的框架,经过后人的整理和完善,但核心思想没有改变。

定义域及其为何编码物理#

Hilbert空间$H$ 上的无界算子是一个线性映射$T: D(T) \to H$ ,其中定义域$D(T)$$H$ 的一个稠密线性子空间。这个映射不必在$H$ 的所有点上都有定义;定义域是数据的一部分。两个公式相同但定义域不同的算子是不同的算子——它们可能有不同的谱和不同的物理解释。

无界算子及其稠密但真子空间的定义域

典型例子:考虑$L^2[0,1]$ 上的$T = -i\,d/dx$ 。有几种自然的选择:

  • $D_{\max} = \{f \in L^2 : f$ 绝对连续,$f' \in L^2\}$ —— 没有边界条件。这给出一个非对称算子。
  • $D_{\text{Dir}} = \{f \in D_{\max} : f(0) = f(1) = 0\}$ —— Dirichlet条件。对称但不是自共轭。
  • $D_{\text{per}} = \{f \in D_{\max} : f(0) = f(1)\}$ —— 周期边界条件。自共轭,谱为$\{2\pi n : n \in \mathbb{Z}\}$
  • $D_{C_c^\infty} = C_c^\infty(0,1)$ —— 光滑紧支集函数。对称,可闭,但远非自共轭。

这四种选择给出了四个不同的算子,尽管它们共享公式$-id/dx$ 。周期算子的特征值为$2\pi n$ ,特征函数为$e^{2\pi i n x}$ 。Dirichlet算子是对称的,亏指标为$(1,1)$ ——它有一个一参数的自共轭延拓族,每个对应于不同的相位条件$f(1) = e^{i\theta}f(0)$ 。最小算子在$C_c^\infty$ 上是对称的,亏指标也是$(1,1)$ ,原因相同。

这种对定义域的敏感性不是缺陷——而是物理。边界条件编码了物理设置:周期边界意味着区间是一个圆,Dirichlet条件表示硬壁,相位条件$f(1) = e^{i\theta}f(0)$ 描述了一个带有磁通量$\theta$ 的粒子。不同的物理需要不同的定义域,谱理论忠实地反映了这一点。

$\mathbb{R}$ (整个实轴)上,情况更简单:$-id/dx$$H^1(\mathbb{R})$ 为定义域本质上是自共轭的(其闭包是自共轭的),并且其唯一的自共轭实现具有谱$\mathbb{R}$ (全连续,没有特征值)。“特征函数”$e^{ikx}$ 不在$L^2$ 中——它们是分布意义下的广义特征函数。无穷远处没有边界消除了有界区间中的模糊性。

一个容易让学生困惑的细节:Laplacian $-d^2/dx^2$$[0,1]$ 上的亏指标为$(2,2)$ (不像一阶算子是$(1,1)$ ),因为ODE $-f'' = \pm if$$[0,1]$ 上有两个线性独立的$L^2$ 解。自共轭延拓形成一个四实参数族,对应于最一般的“边界条件”连接$f(0), f'(0), f(1), f'(1)$ 。Dirichlet条件($f(0)=f(1)=0$ )、Neumann条件($f'(0)=f'(1)=0$ )和周期条件($f(0)=f(1), f'(0)=f'(1)$ )各自从这个族中选出一个延拓。从这个角度看,Sturm-Liouville问题的谱理论就是二阶微分算子在区间上的自共轭延拓分类。

物理原则:指定自共轭延拓所需的边界条件数量等于亏指标(对于正则问题,等于微分方程的阶数)。一阶算子需要一个条件;二阶算子需要两个。微分方程的阶数、边界条件的数量和亏指标之间的这种对应关系是该理论中最清晰的结构结果之一。

Laplacian 的谱与特征函数

闭算子、可闭性与图#

无界算子的基本正则性条件是闭性图 $G(T) = \{(x, Tx) : x \in D(T)\} \subset H \times H$ 是乘积希尔伯特空间的一个线性子空间。如果 $G(T)$$H \times H$ 中是闭的,那么算子 $T$ 就是闭的等价地说,如果 $x_n \in D(T)$$x_n \to x$ ,且 $Tx_n \to y$ ,那么 $x \in D(T)$$Tx = y$

闭算子:图在乘积拓扑中是闭的

闭性比有界性弱,但比单纯的线性强得多。对于有界算子,闭图像定理(第6篇文章)说明闭等于有界。对于无界算子,闭性是一个实质性的条件,它使得谱理论成为可能:闭算子的预解集和谱是明确定义的。

闭性为什么重要。如果 $T$ 是闭的,预解式 $(\lambda I - T)^{-1}$ (当它作为有界算子存在时)会自动由闭图像定理保证是有界的。这使得谱 $\sigma(T) = \{\lambda : (\lambda I - T)^{-1} \text{ 不存在为 } H \text{ 上的有界算子}\}$ 成为 $\mathbb{C}$ 的一个闭子集,并具有合理的性质。如果没有闭性,预解式可能存在但无界,使谱理论变得毫无意义。

一个算子是可闭的,如果它的图的闭包 $\overline{G(T)}$ 本身也是一个算子的图(即,$(0, y) \in \overline{G(T)}$ 意味着 $y = 0$ )。对称算子总是可闭的。闭包 $\bar{T}$ 的定义域是 $D(\bar T) = \{x : \exists (x_n) \subset D(T), x_n \to x, Tx_n \text{ 收敛}\}$ ,并且 $\bar T x = \lim Tx_n$

举个例子。导数算子 $d/dx$$C^1[0,1] \subset L^2[0,1]$ 上是可闭的。它的闭包的定义域是 $H^1[0,1]$ (Sobolev 空间,绝对连续函数且其导数平方可积)。序列 $f_n(x) = x^{1/2 + 1/n}$ 属于 $C^1$ ,且 $f_n \to x^{1/2}$$L^2$ 中收敛,$f_n' = (1/2 + 1/n)x^{-1/2+1/n} \to \frac{1}{2}x^{-1/2}$$L^2$ 中也收敛。因为 $x^{1/2} \in H^1$ (它是绝对连续的且其导数平方可积),所以极限是一致的:闭包将 $x^{1/2}$ 接受为其定义域的一部分。像 $|x-1/2|^{1/4}$ 这样的函数(在 $L^2$ 中但其导数不在 $L^2$ 中)不在闭包的定义域内。

图范数 $\|x\|_T = (\|x\|^2 + \|Tx\|^2)^{1/2}$ 使 $D(T)$ 成为一个希尔伯特空间(当 $T$ 是闭的时)。闭图像定理可以重新表述为:在巴拿赫空间之间的闭算子,如果其定义域等于整个源空间,则必须是有界的。反过来说:无界算子必须有适当的定义域。图范数是 $D(T)$ 上的自然拓扑,许多偏微分方程中的估计(先验估计、正则性定理)都是关于算子在图范数下的有界性的陈述。

闭性的一个关键应用是预解恒等式 $(\lambda - T)^{-1} - (\mu - T)^{-1} = (\mu - \lambda)(\lambda - T)^{-1}(\mu - T)^{-1}$ 对闭算子成立,并给出了预解式 $\lambda \mapsto (\lambda - T)^{-1}$ 作为 $B(H)$ 值函数在预解集上的解析性。这种解析性是 Dunford-Taylor 泛函演算(预解式的围道积分)的基础,并将算子理论与复分析联系起来。

伴随算子与自伴随性:微妙的区别#

对于一个定义域为 $D(T)$ 的稠密定义算子 $T$ ,其伴随算子 $T^*$ 定义在域 $D(T^*) = \{y \in H : x \mapsto \langle Tx, y\rangle \text{ 在 } D(T) \text{ 上有界}\}$ 上。对于这样的 $y$ ,有界线性泛函 $x \mapsto \langle Tx, y\rangle$ 可以通过 Riesz 表示定理扩展到整个 $H$ ,从而得到唯一的 $T^*y$ ,使得 $\langle Tx, y\rangle = \langle x, T^*y\rangle$ 对所有 $x \in D(T)$ 成立。

对称 vs 自伴:定义域的区别

关键点是:$D(T^*)$ 可能比 $D(T)$ 大或小。稠密定义算子的伴随算子总是闭的(它的图像是 $G(T)$$H \times H$ 中旋转后的正交补)。但是 $D(T^*)$ 依赖于 $T$ 及其定义域的方式并不直观。

  • 如果 $D(T) \subseteq D(T^*)$ 且对所有 $x \in D(T)$$T^*x = Tx$ ,那么 $T$ 是对称的等价地,$\langle Tx, y\rangle = \langle x, Ty\rangle$ 对所有 $x, y \in D(T)$ 成立。
  • 如果 $T = T^*$ ,即 $D(T) = D(T^*)$ 且对所有 $x \in D(T)$$Tx = T^*x$ ,那么 $T$ 是自伴随的
  • 如果 $T$ 的闭包 $\bar T$ 是自伴随的,那么 $T$ 是本质自伴随的

对称性比自伴随性弱。对称算子满足 $D(T) \subseteq D(T^*)$ ;自伴随性要求两者相等。$D(T)$$D(T^*)$ 之间的差距由亏指标 $n_\pm = \dim\ker(T^* \mp iI)$ 来衡量。von Neumann 定理指出:一个闭对称算子有自伴随延拓当且仅当 $n_+ = n_-$ ,并且当 $n_+ = n_- = n$ 时,这些延拓由酉群 $U(n)$ 参数化。

这种区别在物理上很重要。对称算子可以生成收缩半群,但不一定生成酉群。自伴随性(通过 Stone 定理)等价于生成一个一参数酉群——即在量子力学中定义一致的时间演化。仅仅是对称但不是自伴随的可观测量不能生成良好的动力学。数学上的区别直接关系到物理内容:只有自伴随算子才对应真正的物理可观测量。

Laplacian 作为 L^2 上的无界自伴随算子

实例:半直线上的动量算子:$P = -id/dx$ 作用在 $L^2[0,\infty)$ 上,定义域为 $D(P) = \{f \in H^1(0,\infty) : f(0) = 0\}$ 。这是对称的:分部积分给出 $\langle Pf, g\rangle - \langle f, Pg\rangle = i[\bar g f]_0^\infty = -if(0)\overline{g(0)} = 0$ (利用 $f(0) = 0$ 和无穷远处的衰减)。计算 $P^*$$g \in D(P^*)$ 的条件是 $f \mapsto \langle Pf, g\rangle = \langle -if', g\rangle$$\{f \in H^1 : f(0)=0\}$ 上有界。分部积分给出 $\langle -if', g\rangle = \langle f, -ig'\rangle$$g \in H^1$ 成立(不需要 $g$ 的边界条件,因为 $f(0) = 0$ 消除了边界项)。因此 $D(P^*) = H^1(0,\infty)$ ——严格大于 $D(P)$ 。这个算子是对称的但不是自伴随的。

亏指标:$\ker(P^* - i) = \{g : -ig' = ig, g \in L^2(0,\infty)\} = \{ce^{-x}\}$ ,是一维的。$\ker(P^* + i) = \{g : -ig' = -ig, g \in L^2(0,\infty)\} = \{ce^x\}$ ,但 $e^x \notin L^2(0,\infty)$ ,所以这是 $\{0\}$ 。亏指标 $(n_+, n_-) = (1, 0)$ 。由于 $n_+ \neq n_-$ ,von Neumann 定理表明:不存在自伴随延拓。半直线上的动量算子没有自伴随实现。

物理解释:一个被无限墙限制在 $[0, \infty)$ 区间的粒子不能有一个良好定义的动量可观测量。它可以向右移动,但墙壁阻止了向左移动,打破了正负动量之间的对称性。亏指标的不对称性($n_+ = 1$$n_- = 0$ )反映了这种物理不对称性。这不是数学上的病态——而是正确功能分析编码的物理现象。

相比之下,在全直线 $\mathbb{R}$ 上,$H^1(\mathbb{R})$ 上的动量 $-id/dx$ 有亏指标 $(0,0)$ (既 $e^x$ 也不 $e^{-x}$ 属于 $L^2(\mathbb{R})$ ),因此它是本质自伴随的。粒子可以在两个方向自由移动,动量是一个真实的可观测量。定义域(全直线 vs 半直线)编码了物理设置,而亏指标检测该设置是否允许一致的可观测量。

无界自伴算子的谱定理#

谱定理可以扩展到无界自伴算子,其表述与有界算子基本相同,但谱可能延伸到无穷远。

核心结论是:对于希尔伯特空间 $H$ 上的每个自伴算子 $T$ (可能是无界的),存在唯一的投影值测度 $E$$(\mathbb{R}, \mathcal{B})$ 上,使得 $T = \int_{-\infty}^{\infty} \lambda\,dE(\lambda)$ 。这意味着对 $x \in D(T)$ 和所有 $y \in H$$\langle Tx, y\rangle = \int \lambda\,d\langle E(\lambda)x, y\rangle$ 。算子 $T$ 的定义域是 $D(T) = \{x : \int \lambda^2\,d\|E(\lambda)x\|^2 < \infty\}$

函数演算也扩展了:对于任何 Borel 函数 $f: \mathbb{R} \to \mathbb{C}$ ,定义 $f(T) = \int f(\lambda)\,dE(\lambda)$ ,其定义域为 $\{x : \int |f(\lambda)|^2\,d\|E(\lambda)x\|^2 < \infty\}$ 。这样可以得到 $e^{itT}$ (实数 $t$ 时为酉算子,生成 Stone 定理)、$(T - \lambda)^{-1}$ (预解式)、$|T| = \int|\lambda|\,dE(\lambda)$ 等。

$\mathbb{R}^n$ 上的拉普拉斯算子#

算子 $-\Delta$$H^2(\mathbb{R}^n)$ 为定义域是自伴的。通过傅里叶变换,$\widehat{(-\Delta f)}(\xi) = |\xi|^2\hat f(\xi)$ ,所以 $-\Delta$$L^2(\mathbb{R}^n)$ 上乘以 $|\xi|^2$ 的算子在酉等价。谱为 $[0,\infty)$ ,全是连续谱(没有特征值——没有 $L^2$ 函数满足 $|\xi|^2\hat f = \lambda\hat f$ 对于固定的 $\lambda$ )。谱测度是 $E(B)f = \mathcal{F}^{-1}(\mathbf{1}_{|\xi|^2 \in B}\hat f)$

这种酉等价(通过傅里叶变换)是谱定理应用的一个典型例子。抽象陈述“存在一个酉算子 $U$ 和一个测度空间 $(\Sigma, \mu)$ 使得 $UTU^{-1}$ 是乘以可测函数”变得具体:$U = \mathcal{F}$ (傅里叶变换),函数是 $|\xi|^2$ 。函数演算也非常明确:$f(-\Delta)\psi = \mathcal{F}^{-1}(f(|\xi|^2)\hat\psi)$ 。例如,热半群 $e^{t\Delta}\psi = \mathcal{F}^{-1}(e^{-t|\xi|^2}\hat\psi)$ ——与高斯核 $(4\pi t)^{-n/2}e^{-|x|^2/(4t)}$ 卷积。

极小极大原理#

$$\lambda_k = \min_{\substack{V \subset D(T) \\ \dim V = k}} \max_{\substack{x \in V \\ \|x\| = 1}} \langle Tx, x\rangle.$$

这是 Courant-Fischer 极小极大表征。它是量子化学中变分方法的基础(选择有限维子空间给出特征值的上界),也是区域单调性结果的基础(扩大区域——放松边界条件——只能减少特征值)。对于有界区域 $\Omega$ 上的狄利克雷拉普拉斯算子,极小极大原理给出 Weyl 渐近公式:$\lambda_k \sim C_n(|\Omega|^{-2/n})k^{2/n}$$k \to \infty$ 时,将特征值的增长率与区域体积联系起来。

本质谱 vs 离散谱

谱分为离散$\sigma_d$ = 孤立特征值且重数有限)和本质$\sigma_{ess}$ = 其余部分)。Weyl 定理:对于紧算子 $K$$\sigma_{ess}(T + K) = \sigma_{ess}(T)$ 。对于薛定谔算子 $-\Delta + V$$V \to 0$ 在无穷远处,$\sigma_{ess}(-\Delta + V) = [0,\infty) = \sigma(-\Delta)$ (本质谱不受衰减势的影响)。低于零的离散特征值对应于束缚态——氢原子有 $\sigma_d = \{-1/(4n^2) : n \geq 1\}$$\sigma_{ess} = [0,\infty)$

自伴性判据:加藤-雷利希定理与弗里德里希斯延拓#

证明自伴性是大多数应用中的难点。一旦确立了自伴性,谱定理和斯通定理就会自动适用。这里有三个标准工具:

亏指数决定自伴扩张

(a) 加藤-雷利希定理。如果 $T_0$ 是自伴的,并且 $V$ 是对称的,满足 $D(V) \supseteq D(T_0)$ 以及 $\|Vf\| \leq a\|T_0 f\| + b\|f\|$ 对于某个 $a < 1$ 成立,那么 $T_0 + V$$D(T_0)$ 上是自伴的。这是处理薛定谔算子的主要工具:取 $T_0 = -\Delta$$V$ 为势能,条件要求 $V$ 在某种精确意义上被动能“支配”。

举个例子:氢原子。$H = -\Delta - 1/|x|$$L^2(\mathbb{R}^3)$ 上。哈代不等式 $\int|f|^2/|x|^2 \leq 4\int|\nabla f|^2$ 给出 $\|f/|x|\|_{L^2} \leq 2\|\nabla f\|_{L^2}$ ,而索伯列夫不等式则将 $\|\nabla f\|_{L^2}$ 限制在 $\epsilon\|\Delta f\|_{L^2} + C_\epsilon\|f\|_{L^2}$ 内。结合两者:$\|Vf\| = \|f/|x|\| \leq \epsilon\|\Delta f\| + C\|f\|$ 对任意 $\epsilon > 0$ 成立。取 $\epsilon < 1$ 满足加藤-雷利希假设。结论是:$-\Delta - 1/|x|$$H^2(\mathbb{R}^3)$ 上是自伴的。

从这个单一的自伴性结果中,可以推导出所有氢原子物理。谱定理给出谱:离散部分 $\sigma_d = \{-1/(4n^2) : n = 1, 2, 3, \ldots\}$ (波尔能级,退化度为 $n^2$ )和本质谱 $\sigma_{ess} = [0, \infty)$ (散射态)。斯通定理给出时间演化 $e^{-iHt}$ ——量子电子在库仑场中的动力学。散射理论(在电离阈值 $E = 0$ 以上)描述了平面波如何被势能偏转。这一切都来自一个加藤-雷利希估计。

这种方法可以扩展到多电子原子和分子。氦原子的哈密顿量是 $H = -\Delta_1 - \Delta_2 - 2/|x_1| - 2/|x_2| + 1/|x_1 - x_2|$$L^2(\mathbb{R}^6)$ 上。加藤-雷利希定理仍然适用(每个库仑项都满足相对界条件),从而在 $H^2(\mathbb{R}^6)$ 上得到自伴性。对于更重的原子,证明会更技术化(需要处理多体库仑奇异性),但原则不变:证明扰动 $V$ 相对于动能 $-\Delta$ 是小的,得出自伴性,然后应用谱理论。

(b) 弗里德里希斯延拓。对于半有界对称算子($\langle Tx, x\rangle \geq c\|x\|^2$ 对所有 $x \in D(T)$ 成立),存在一个规范的自伴延拓——弗里德里希斯延拓——它保留了下界。通过在“能量范数” $\|x\|_T = (\langle Tx, x\rangle + (1-c)\|f\|^2)^{1/2}$ 下完成 $D(T)$ 并通过紧包含将其识别为 $H$ 的子空间来构造。弗里德里希斯延拓是唯一一个其定义域包含在形式域内的自伴延拓。它是变分问题中最自然的延拓。

举个例子:狄利克雷拉普拉斯算子。设 $\Omega \subset \mathbb{R}^n$ 是一个有界开集。定义 $T_0 = -\Delta$$C_c^\infty(\Omega) \subset L^2(\Omega)$ 上。这是对称且正的:$\langle -\Delta f, f\rangle = \|\nabla f\|^2 \geq 0$ 。能量形式是 $Q(f, g) = \int_\Omega \nabla f \cdot \nabla g$ ,并且 $C_c^\infty(\Omega)$ 在范数 $(\|\nabla f\|^2 + \|f\|^2)^{1/2}$ 下的完备化是 $H^1_0(\Omega)$ ——具有零边界值的索伯列夫空间。弗里德里希斯延拓是 $-\Delta_D$ (狄利克雷拉普拉斯算子),其定义域为 $D(-\Delta_D) = H^2(\Omega) \cap H^1_0(\Omega)$ :在迹意义上在 $\partial\Omega$ 上消失的 $H^2$ 函数。这就是出现在狄利克雷问题 $-\Delta u = f$$u|_{\partial\Omega} = 0$ 中的算子,其特征值 $0 < \lambda_1 < \lambda_2 \leq \cdots$ 是固定边界的膜的共振频率。弗里德里希斯构造从 $-\Delta|_{C_c^\infty}$ 的许多可能自伴延拓中选择了“狄利克雷”延拓,而且它不需要任何正则性理论——只需要能量形式 $\int|\nabla f|^2$

(c) 通过亏缺指标的本质自伴性。对称算子本质上是自伴的当且仅当 $n_+ = n_- = 0$ ,即 $T^* \pm i$ 有平凡核。等价地,$\text{Range}(T \pm i)$$H$ 中稠密。对于 $-\Delta$$C_c^\infty(\mathbb{R}^n)$ 上,本质自伴性可以从 $-\Delta$ 的椭圆性和光滑函数在图范数下的稠密性得出。

量子力学中的位置、动量和哈密顿算子

量子力学:典型算子#

量子力学中的算子及其定义域:

本质谱与离散谱

位置 $X$$L^2(\mathbb{R})$ 上:$(Xf)(x) = xf(x)$ ,定义域 $D(X) = \{f : xf \in L^2\}$ 。自共轭(乘以实值函数)。谱为 $\mathbb{R}$ ,全部连续。谱测度:$E(B)f = \mathbf{1}_B f$

动量 $P = -id/dx$$L^2(\mathbb{R})$ 上:定义域 $H^1(\mathbb{R})$ 。自共轭(通过傅里叶变换与 $X$ 单位等价:$\mathcal{F}P\mathcal{F}^{-1} = M_\xi$ )。谱为 $\mathbb{R}$ ,全部连续。

正则对易关系 $[X, P] = iI$ 在包含施瓦茨空间 $\mathcal{S}(\mathbb{R})$ 的共同定义域上。这不能由有界算子实现(取迹得到 $0 = \text{tr}([X,P]) = i\,\text{tr}(I)$ ,这是不可能的)。无界性对于量子力学是必不可少的。

Weyl 形式 的对易关系避免了定义域问题:$e^{isX}e^{itP} = e^{-ist}e^{itP}e^{isX}$ 对所有 $s, t \in \mathbb{R}$ 成立。这是一个在有界(酉)算子之间的关系,处处定义良好。Stone-von Neumann 定理 表明,在可分希尔伯特空间上,Weyl 关系的不可约表示最多只有一个——即 $L^2(\mathbb{R})$ 上的薛定谔表示。这个唯一性定理解释了为什么量子力学具有唯一的运动学结构(每个希尔伯特空间中的一自由度都与 $L^2(\mathbb{R})$ 上的位置和动量单位等价)。对于无限多个自由度(量子场论),Stone-von Neumann 定理失效——存在不等价的表示,选择一个表示是物理的一部分(不同的真空态,不同的超选择扇区)。

谐振子 $H = -\frac{1}{2}d^2/dx^2 + \frac{1}{2}x^2$$L^2(\mathbb{R})$ 上:自共轭,定义域 $\{f \in H^2 : x^2f \in L^2\}$ 。纯点谱 $\{n + 1/2 : n = 0, 1, 2, \ldots\}$ ,特征函数是厄米特函数 $\phi_n(x) = c_n H_n(x)e^{-x^2/2}$ ,形成 $L^2(\mathbb{R})$ 的正交基。这是最干净的非平凡量子系统——显式可对角化且谱离散。通过创建和湮灭算子 $a^\pm = (X \mp iP)/\sqrt{2}$ 满足 $[a^-, a^+] = I$ 的代数方法,可以得到 $H = a^+a^- + 1/2$ ,并且谱从阶梯结构得出:$a^+$ 将特征值提高 1,$a^-$ 降低特征值,$a^-\phi_0 = 0$ 确定了基态 $\phi_0(x) = \pi^{-1/4}e^{-x^2/2}$ 。整个谱理论归结为代数。这种代数方法推广到量子场论(Fock 空间),其中每个模式的创建和湮灭算子构造希尔伯特空间本身。

Stone 定理 连接自共轭性和动力学:稠密定义的算子 $T$ 生成强连续的一参数酉群 $U(t) = e^{itT}$ 当且仅当 $T$ 是自共轭的。物理内容:自共轭算子生成对称性(哈密顿量的时间演化,动量的空间平移,角动量的旋转)。数学陈述“可观测量是自共轭的”等价于“可观测量生成连续对称性”。这是 Noether 定理的算子理论内容。

如果使用仅仅是对称(不是自共轭)的算子,$e^{itT}$ 不能定义为所有 $t$ 的酉算子——群性质失效。自共轭性正是时间演化一致性的条件。对称算子可能生成一参数半群(向前演化时收缩),但不能生成群(没有时间反演)。这种过去和未来的不对称性在物理上有意义:耗散系统(热方程,阻尼振荡器)由非自共轭算子生成,其时间演化是不可逆的。

Stone 定理的一个方向的证明直接来自谱定理:给定自共轭 $T$ 和谱测度 $E$ ,定义 $U(t) = e^{itT} = \int e^{it\lambda}\,dE(\lambda)$ 。由于 $|e^{it\lambda}| = 1$ ,这是酉算子。强连续性从控制收敛定理得出:$\|U(t)x - x\|^2 = \int|e^{it\lambda} - 1|^2\,d\|E(\lambda)x\|^2 \to 0$$t \to 0$ 。相反方向(每个强连续酉群都有自共轭生成元)是更深层次的部分,需要从群重构生成元:$Tx = \lim_{t\to 0} (U(t)x - x)/(it)$

实例:自由薛定谔方程$i\partial_t\psi = -\frac{1}{2}\Delta\psi$$L^2(\mathbb{R}^3)$ 上。哈密顿量 $H = -\frac{1}{2}\Delta$$H^2(\mathbb{R}^3)$ 上自共轭。Stone 定理给出解 $\psi(t) = e^{-iHt}\psi_0 = e^{it\Delta/2}\psi_0$ 。在傅里叶空间中:$\hat\psi(t, \xi) = e^{-it|\xi|^2/2}\hat\psi_0(\xi)$ ——每个频率以依赖于 $|\xi|$ 的二次相位传播。色散性质(高频传播更快)导致波包扩散,并且衰减 $\|\psi(t)\|_{L^\infty} \leq C t^{-3/2}\|\psi_0\|_{L^1}$ 从驻相法得出。所有这些都由应用于单个自共轭算子的谱定理组织起来。

常见困惑、Trotter公式和数值方法#

值得明确指出的五个常见困惑:

  1. “对称意味着自共轭。” 错误。差距(亏指标)可能很大。Hellinger-Toeplitz定理:一个处处定义的对称算子是有界的。因此,真正无界的对称算子必然定义在适当的子空间上,而自共轭是一个真正的附加条件。

  2. “对称算子的闭包是自共轭的。” 错误。闭包是对称的,但仍然可能有非零的亏指标。当存在自共轭延拓时,它们会将定义域扩展到闭包之外。

  3. “所有自共轭延拓具有相同的谱。” 错误。不同的边界条件给出不同的谱:Dirichlet Laplacian在$[0,1]$ 上的特征值是$n^2\pi^2$ ;Neumann Laplacian的特征值包括$n=0$ 的情况,即$n^2\pi^2$ ;周期条件给出$(2\pi n)^2$ 且重数为二。

  4. “伴随 = Hermitian共轭。” 对于矩阵来说,是的。对于无界算子,伴随涉及仔细的定义域指定。“形式伴随”(分部积分)仅在特定定义域上等于算子论中的伴随。

  5. “每个稠密定义的对称算子都有自共轭延拓。” 错误:亏指标必须满足$n_+ = n_-$ 。半直线动量算子(前面的例子)有$(n_+, n_-) = (1, 0)$ ,没有自共轭延拓。物理教训是,并不是每一个写下来的“可观测量”都是真正的可观测量。如果你提出一个量(如半直线上的动量),发现亏指标不相等,这意味着这个量没有良好的概率分布——它不是一个真正的量子可观测量。

另一个值得提及的困惑:“$T$ 限制在不变子空间上的谱等于$\sigma(T)$ 的限制。” 这对于无界算子严重失效。如果$M$$T$ 的闭不变子空间,$T|_M$ 表示限制(定义域为$D(T) \cap M$ ),那么$\sigma(T|_M)$ 可能比$\sigma(T) \cap$ “相关部分”大得多。限制的谱可能会增加点(由于投影有效施加的边界条件),即使全算子在该处有间隙。这是定义域敏感性的又一表现:限制定义域可以创建新的谱。

$$e^{it(A+B)} = \lim_{n \to \infty}\left(e^{itA/n} e^{itB/n}\right)^n$$

在强算子拓扑中成立。对于量子力学,其中$A = -\Delta/2$ (动能)和$B = V$ (势能),这将时间演化分解成交替的“自由传播”和“势能踢”。物理图像是:在一个短时间间隔$\Delta t = t/n$ 内,粒子首先自由传播(根据自由Schrödinger方程扩散),然后从势能中获得相位踢。当$n \to \infty$ 时,交替变得连续,恢复精确演化。

计算上,这是分裂步法(也称为分裂算子或Strang分裂)——模拟Schrödinger方程的工作马算法。自由传播步骤在Fourier空间中是对角的($e^{-it|\xi|^2/(2n)}$ ),势能步骤在物理空间中是对角的($e^{-itV(x)/n}$ )。使用FFT在这两者之间交替,每一步的时间复杂度为$O(N\log N)$ ,基本Lie-Trotter分裂的误差为$O((\Delta t)^2)$ ,对称Strang分裂的误差为$O((\Delta t)^3)$ 。这种方法用于光纤模拟、玻色-爱因斯坦凝聚动力学、量子计算模拟以及无数其他应用。数学内容纯粹是Trotter公式加上$A + B$ 的自共轭性。

数值谱计算:实际上,通过取有限区间和离散化来计算无界自共轭算子的谱。在$[-L, L]$ 区间内的$N$ 个网格点上,$-\Delta$ 变成一个$N \times N$ 的三对角矩阵(二阶差分矩阵),其特征值近似连续算子的低能量谱。收敛理论有两个部分:(1) 谱逼近定理 保证离散化算子的特征值在$N \to \infty$$L \to \infty$ 时收敛到连续算子的特征值(对于离散谱),(2) Weyl准则 描述本质谱:$\lambda \in \sigma_{ess}(T)$ 当且仅当存在Weyl序列$(x_n)$ ,使得$\|x_n\| = 1$$x_n \rightharpoonup 0$ 弱收敛,且$(T-\lambda)x_n \to 0$ 。数值上,本质谱表现为随着离散化细化而不收敛到孤立点的密集特征值簇。

一个严重的陷阱是谱污染在本质谱的间隙中出现虚假特征值,这些特征值不会收敛到任何真实特征值。当用于离散化的有限维子空间(Galerkin投影)不尊重算子结构时,这种情况会发生。对于Dirac算子(其本质谱为$(-\infty, -mc^2] \cup [mc^2, \infty)$ ,在$(-mc^2, mc^2)$ 间隙中包含离散特征值),朴素的有限元离散化会在整个间隙中产生虚假特征值。解决方法包括平衡基(选择尊重Dirac算子块结构的试验函数)和二次方法(通过在网格上找到$\det(T_N - \lambda)$ 穿过零点的位置来计算$\sigma(T)$ ,而不是直接对角化投影矩阵)。

对于有下界的自共轭算子,min-max原理保证Galerkin特征值总是真实特征值的上界——Rayleigh-Ritz方法不会低估。这就是为什么变分计算基态如此稳健:投影到任何有限基都会给出$\lambda_1$ 的上界,增加基只会改进估计。收敛速度取决于基如何捕捉真实特征函数:多项式基给出代数收敛(对于$C^k$ 特征函数,速率为$O(N^{-k})$ ),而谱方法(Fourier, Hermite)对于解析特征函数给出指数收敛。谐振子特征函数是高斯函数与整函数的乘积,Hermite谱方法对其近似呈指数快速——这是自共轭性保证的光滑性的实际结果。

数值算例:$-d^2/dx^2$$[0, 1]$ 上的谱依赖于定义域#

定义域不是装饰——它精确决定谱。同一个微分表达式 $-d^2/dx^2$ 在不同边界条件下给出不同算子,谱完全不同。

Dirichlet 边界($f(0) = f(1) = 0$ ): 特征函数 $\sin(n\pi x)$ ,特征值 $\lambda_n = (n\pi)^2$$n = 1, 2, \ldots$ 。具体数:$\lambda_1 = \pi^2 \approx 9.870$$\lambda_2 = 4\pi^2 \approx 39.48$$\lambda_{10} = 100\pi^2 \approx 986.96$ 。基态非零,Weyl 渐近 $\lambda_n \sim n^2 \pi^2$

Neumann 边界($f'(0) = f'(1) = 0$ ): 特征函数 $\cos(n\pi x)$ ,特征值 $\lambda_n = (n\pi)^2$$n = 0, 1, 2, \ldots$ (包含 $\lambda_0 = 0$ 对应常函数)。具体数:$\lambda_0 = 0$$\lambda_1 = \pi^2$$\lambda_2 = 4\pi^2$ 。区别:Neumann 谱多了一个 $\lambda_0 = 0$ (常函数是基态)。

周期边界($f(0) = f(1)$$f'(0) = f'(1)$ ): 特征值 $\lambda_n = (2 n \pi)^2$$n = 0, 1, 2, \ldots$ ,每个 $n \geq 1$ 重数 $2$ (成对正余弦)。$\lambda_0 = 0, \lambda_1 = 4\pi^2 \approx 39.48, \lambda_2 = 16\pi^2 \approx 157.9$

Dirichlet-Neumann 混合($f(0) = 0, f'(1) = 0$ ): 特征函数 $\sin((n - 1/2)\pi x)$$\lambda_n = ((n - 1/2)\pi)^2$$\lambda_1 = \pi^2/4 \approx 2.467$$\lambda_2 = 9\pi^2/4 \approx 22.21$ 。基态半频,谱整体往低拉。

对照: 物理意义是不同的——Dirichlet 是"两端固定的弦",Neumann 是"两端自由的弦",周期是"圆周上的弦",混合是"一端固定一端自由"。同一个微分表达式 $-d^2/dx^2$ 配上不同 $\mathcal{D}(T)$ 给出四个不同的自共轭算子,谱可以差一个零、可以差因子四、可以差重数。在量子力学这转化为不同的物理系统:盒子里的粒子(Dirichlet)、绝热盒(Neumann)、环(周期)、半盒(混合),能级表全不同。

反例:对称不蕴含自伴#

最常见的混淆:把对称(symmetric, $T \subset T^*$ )当成自伴(self-adjoint, $T = T^*$ )。这两条对有界算子等价,对无界算子完全不等价。

经典反例:$T = -i \, d/dx$$\mathcal{D}(T) = C_c^\infty(0, 1)$ (紧支撑光滑函数)上。$T$ 对称:$\langle Tf, g\rangle = -i \int_0^1 f' \bar g\,dx = i \int_0^1 f \bar g'\,dx = \langle f, Tg\rangle$ ,分部积分边界项消失因为 $f, g$ 紧支撑。

$T$ 自伴:伴随的定义域 $\mathcal{D}(T^*) = H^1(0, 1) = \{f \in L^2 : f' \in L^2\}$$\mathcal{D}(T) = C_c^\infty$ 大得多($T^*$ 的元素不需要在边界为零)。所以 $T \subsetneq T^*$ ,严格包含但不等。

亏指标: 计算 $\dim \ker(T^* \mp i)$ 。求 $-i f' \pm i f = 0$$f' = \mp f$ ,解 $f = c e^{\mp x}$ 。两个亏空间各是 $1$ 维,亏指标 $(1, 1)$ 。亏指标相等说明 $T$ 有自伴延拓,但不唯一——一族延拓由相位 $\theta \in [0, 2\pi)$ 参数化,对应不同的边界条件 $f(1) = e^{i\theta} f(0)$$\theta = 0$ 是周期边界,$\theta = \pi$ 是反周期,每个不同 $\theta$ 给出谱不同的自伴算子。

陷阱: 写"$T$ 自伴"前要确认 $\mathcal{D}(T) = \mathcal{D}(T^*)$ 。验证 $T$ 对称很容易(一次分部积分),验证 $T$ 自伴要算亏指标或验证 $\mathcal{D}(T^*)$ 不超过 $\mathcal{D}(T)$ 。把这两件事混淆,谱定理就会被错误地引用——对称非自伴的算子可能没有谱测度,可能没有完整的函数演算。

常见陷阱:Hamilton 算子的本质自伴性需要单独验证#

物理书常常默默假设"$H = -\Delta + V$ 自伴"。事实上自伴性是一条非平凡的性质,依赖 $V$ 的具体形式。

具体陷阱:对 $V(x) = -1/|x|$ (氢原子),$H = -\Delta - 1/|x|$$\mathbb{R}^3$本质自伴$C_c^\infty(\mathbb{R}^3 \setminus \{0\})$ (Kato-Rellich 定理给出,因为 $1/|x|$$\mathbb{R}^3$ 上是 $\Delta$ -紧的)。但在 $\mathbb{R}$ 上同一公式 $H = -d^2/dx^2 - 1/|x|$ 不自伴——奇点 $x = 0$ 在一维需要单独处理(亏指标可能非零)。

第二个陷阱:势 $V$$L^p$ 类条件决定自伴性。Kato-Rellich 要求 $V = V_1 + V_2$ 其中 $V_1 \in L^2_{loc}$$\Delta$ -相对有界(界 $< 1$ ),$V_2 \in L^\infty$ 。验证一个具体的 $V$ 满足这条件需要算 Sobolev 嵌入,不是几何直觉能猜的。本系列后面提到"$H$ 自伴"时都暗示了这条验证已经做过,但读者引用时要核对——错过这一步就跌进了"对称非自伴"的反例库。

下一步#

下一篇文章将探讨无界自伴算子在动力学中的应用,构建它们生成的一参数半群。Hille-Yosida 定理刻画了强连续压缩半群的生成元——这些算子 $A$ 使得 $e^{tA}$$t \geq 0$ 时存在且为有界算子,并满足 $\|e^{tA}\| \leq 1$ 。对于酉情况(自伴生成元),这是 Stone 定理;对于耗散情况(增生算子,不一定自伴),则需要完整的 Hille-Yosida 机制。Lumer-Phillips 定理给出了一个优雅的重新表述:一个稠密定义的闭算子 $A$ 生成压缩半群当且仅当 $A$$A^*$ 都是耗散的(即对所有 $x \in D(A)$$\text{Re}\langle Ax, x\rangle \leq 0$ )。

动画:热半群对初始数据的光滑化

这些半群可以解决演化偏微分方程的初值问题:热方程 $\partial_t u = \Delta u$ (由 $\Delta$ 生成,耗散),波动方程 $\partial_t^2 u = \Delta u$ (简化为一阶系统,由能量空间上的斜自伴算子生成),薛定谔方程 $i\partial_t u = Hu$ (由 $-iH$ 生成,酉),以及 Fokker-Planck 方程 $\partial_t \rho = \nabla\cdot(D\nabla\rho - b\rho)$ (由非自伴二阶算子生成)。这个框架将时间依赖的偏微分方程转化为固定 Hilbert 空间上的算子理论——这是我们已经建立的谱理论的自然延续。关键在于,解 PDE 的时间演化等价于指数化一个(通常是无界的)算子,而这种指数化产生良好解的条件正是 Hille-Yosida 定理验证的条件(闭性、定义域的稠密性、耗散性)。


这是 泛函分析 系列文章的第 9 部分(共 12 篇)。

本系列

泛函分析 12 篇

  1. 01 泛函分析(一):度量空间 —— 距离、收敛与完备性
  2. 02 泛函分析(二):赋范空间与Banach空间
  3. 03 泛函分析(三):Hilbert 空间 —— 无限维空间中的几何
  4. 04 泛函分析(四):对偶空间与 Hahn-Banach 定理 —— 线性泛函的驯服
  5. 05 泛函分析(五):弱拓扑和弱*拓扑 —— 当范数收敛太强时
  6. 06 泛函分析(六):有界线性算子与三大定理
  7. 07 泛函分析(七):紧算子——通往有限维的桥梁
  8. 08 泛函分析(八):谱理论 —— 分解算子
  9. 09 泛函分析(九):无界算子 —— 当有界性失效时 当前
  10. 10 泛函分析(十):算子半群 — 无限维空间中的演化方程
  11. 11 泛函分析(十一):分布与Sobolev空间 — 广义解
  12. 12 泛函分析(十二):泛函分析在行动 —— 偏微分方程和量子力学

读有所得?

GitHub 关注我 → 新文周更

GitHub