格林函数与势论

当时的困境

到 1820 年代，静电学已经有了它的方程。

拉普拉斯为引力问题写下了 $\nabla^2 \phi = 0$，泊松 1813 年把它推广到有源的情形 $\nabla^2 \phi = -4\pi\rho$。这些方程同时适用于引力和静电，形式完全一样。

方程有了，解不出来。

困难在边界。真空中一个孤立点电荷的势是教科书上的一行公式。但真实的问题几乎总是带边界的：一个电荷放在一块接地的金属板附近，金属表面会感应出电荷，感应电荷又改变原来的场，而感应电荷的分布本身正是要求的未知量。

这类问题当时只能一个一个硬解，每换一种几何形状就得从头再来。泊松花了很大力气算球形导体的情形，那已经是当时能处理的最复杂的边界之一。

还有一个更基础的问题：那个函数还没有名字。 拉普拉斯和泊松都在用它——它的梯度给出力，它满足那条方程——但在文献里它只是”那个函数”。一个没有名字的概念，很难成为思考的对象。

洞见时刻

1828 年，诺丁汉一个三十五岁的磨坊工人自费印了一本小册子，一共只征到 51 位订阅者。

他叫乔治·格林，一生只上过约一年半的学，数学是在诺丁汉订阅图书馆自学的——为了读拉普拉斯的原文，他自学了法文。

《论数学分析在电磁理论中的应用》做了三件事。

第一，给那个函数命名。 他叫它 potential function——势函数。这个词沿用至今。

第二，给出一条把体积积分换成边界积分的恒等式（今天叫格林恒等式）：

$\int_V (\phi \nabla^2 \psi - \psi \nabla^2 \phi)\, dV = \oint_S (\phi \nabla\psi - \psi \nabla\phi) \cdot d\mathbf{S}$

物理上的意思很朴素：要知道一个区域内部的场，只需要知道边界上发生了什么。

第三，也是真正的核心——一个求解策略。

与其去解那个带任意源、任意边界的难题，不如先解一个最简单的情形：一个单位点源放在区域内某处、边界条件相同时，场是什么样？

$\nabla^2 G(\mathbf{r}, \mathbf{r}') = \delta(\mathbf{r} - \mathbf{r}')$

把这个解出来，得到一个双变量的函数 $G$——它回答的是”源在 $\mathbf{r}'$ 时，$\mathbf{r}$ 处的响应是多少”。

然后，因为方程是线性的，任意的源分布都可以看成无数点源的叠加，解就是一个积分：

$\phi(\mathbf{r}) = \int G(\mathbf{r}, \mathbf{r}')\, \rho(\mathbf{r}')\, dV'$

难题被拆成了一次性的工作加上一次积分。

这个策略的力量在于它把”解一个具体问题”变成了”解一次基本问题，然后无限次复用”。同一个边界几何，无论源怎么变，$G$ 都不用重算。

$G$ 后来被称为格林函数。

（他在同一本书里还用了后来叫”镜像法”的技巧：接地平面对电荷的全部影响，等价于在平面另一侧放一个符号相反的虚拟电荷。汤姆孙 1848 年独立重新发现并系统化了它。）

涟漪效应

这本书差点就消失了。

格林 1841 年去世时，他的《论文》在学界几乎无人知晓。1845 年，二十一岁的威廉·汤姆孙——后来的开尔文勋爵——从导师那里听说有这本册子，找了很久，最后在离开剑桥的前一天借到三份。他带去巴黎给刘维尔和斯图姆看；据他回忆，刘维尔读完激动得难以自持。

在汤姆孙的推动下，《论文》1850–1854 年在克雷尔的杂志上分三部分重印，这才真正进入学界。从出版到被人读懂，隔了二十二年，而作者已经去世九年。

此后它成了物理学最通用的工具之一，而且远远超出了静电学。

关键在于线性。 只要方程是线性的，格林的策略就成立——而物理学里绝大多数基本方程都是线性的，或者可以在某个近似下线性化。

于是同一个思路反复出现，只是换了名字：

• 电磁学：天线的辐射、波导的模式、散射问题。
• 量子力学：定态散射理论的李普曼-施温格方程就是一个格林函数积分方程。
• 量子场论：费曼传播子就是格林函数。费曼图里每一条内线，画的都是一个格林函数——粒子从一处传到另一处的振幅。整个微扰论就是”点源解 + 叠加”这个策略在场论里的展开。
• 多体理论：凝聚态用的单粒子格林函数（含时的、有限温度的松原格林函数），是计算材料电子结构的主力工具。GW 近似、动力学平均场理论，名字里的 G 就是格林函数。
• 弹性与地震：地震波从震源传到台站的响应函数就是格林函数。今天的震源反演与层析成像，做的正是从记录反推震源和介质。
• 工程：边界元法之所以能只在物体表面划网格（而不是整个空间），靠的就是格林恒等式把体积积分化成了面积分。

格林本人研究的是静电学。他大概想不到，同一个式子会用来算电子在晶体里的传播、地震波穿过地幔、以及两个夸克交换胶子。

如果没有它

电磁仿真。天线的方向图、手机内部的电磁兼容、雷达散射截面、芯片封装的信号完整性——工业界用的矩量法与边界元法，内核就是格林函数。

量子场论。没有传播子，就没有费曼图，就没有能算的微扰论。QED 那个精确到小数点后十二位的电子反常磁矩，是一层层格林函数积出来的。

材料计算。预测一种新材料是导体还是绝缘体、算出它的光吸收谱，今天的标准方法（GW、DMFT）都建立在格林函数形式之上。

地球内部的图像。我们知道地核是液态的、地幔里有下沉的板块残骸，靠的是全球地震台网的记录加上格林函数反演。

逆问题的通用语言。医学的超声与光学层析、无损探伤、油气勘探——凡是”从边界上的测量反推内部结构”的问题，格林函数都是那个把两者联系起来的桥。

一本卖出 51 份、几乎无人读懂的小册子，最后成了物理学解方程的默认办法。