一次统计力学之旅(二):从无知中涌现的秩序

系列第(二)篇,共两篇。第(一)篇埋下一条公理——孤立系统以相等概率拜访每一个可达到的状态——并从中长出整座热力学的骨架:熵 \(\sigma = \ln g\)、温度 \(1/\tau = \partial\sigma/\partial U\)、压强、化学势,还有那些全都流自 \(d\sigma = 0\) 的热力学恒等式。我们停在一句别扭的自白上:整台机器都靠知道多重度 \(g\) 运转,而清点 \(10^{23}\) 个粒子的微观态根本没指望。这一篇就来补上这个窟窿。

我们就从那句自白接着往下走。我们花了四章把热力学正着建起来,从微观计数一路向上——它在概念上美得很。麻烦在于,概念逻辑和计算逻辑指向相反的两个方向。

  • 概念上: \(g \to \sigma \to \tau \to \text{平衡}\)
  • 计算上: \(g\) 几乎永远是一场噩梦——你试试当真去数 \(10^{23}\) 个粒子的排列组合。

于是物理学家凿出了一条”反向通道”。不从 \(g\) 出发,而是写下能量模型,算出一个叫配分函数 \(Z\) 的量,读出自由能 \(F\)然后再下降到熵 \(\sigma\) 和多重度 \(g\)。我们会遇见这条路上的两位明星——玻尔兹曼分布(固定 \(N\))和吉布斯分布(变量 \(N\))。引路的问题不再是”到底有多少个状态?“,而变成了远为好对付的一个:把系统丢进一个巨大的热库——找到它处于某一特定状态的概率是多少?

第五章:热物理学的罗塞塔石碑

不过在倒着跑任何东西之前,先停一停,把我们在第(一)篇里测绘过的地形通览一遍,因为它会归拢成某种很美的东西。四章下来,我们从三个视角审视了系统——微观计数(熵 \(\sigma\))、能量守恒(内能 \(U\))、实用主义(自由能 \(F\))。这有点”盲人摸象”的味道:三段描述,摸的是三个不同的部位。把这三个势 \(\sigma, U, F\) 对齐,列出它们的偏导,一张惊人对称的表格就出现了:

物理量熵表示 σ(U,V,N)内能表示 U(σ,V,N)自由能表示 F(τ,V,N)
温度 \(\tau\)\(\frac{1}{\tau}=\left(\frac{\partial \sigma}{\partial U}\right)_{V, N}\)\(\tau=\left(\frac{\partial U}{\partial \sigma}\right)_{V, N}\)\(\sigma=-\left(\frac{\partial F}{\partial \tau}\right)_{N,V}\)
压强 \(P\)\(\frac{P}{\tau}=\left(\frac{\partial \sigma}{\partial V}\right)_{U, N}\)\(-P=\left(\frac{\partial U}{\partial V}\right)_{\sigma, N}\)\(-P=\left(\frac{\partial F}{\partial V}\right)_{\tau, N}\)
化学势 \(\mu\)\(-\frac{\mu}{\tau}=\left(\frac{\partial \sigma}{\partial N}\right)_{U, V}\)\(\mu=\left(\frac{\partial U}{\partial N}\right)_{\sigma, V}\)\(\mu=\left(\frac{\partial F}{\partial N}\right)_{\tau, V}\)

这块罗塞塔石碑里,有三层意思要读出来。

对角线上的”定义”

每一格在数学上都成立,但物理学家倚重从左上到右下的对角线来确立定义

  • 第一行第一列: \(\frac{1}{\tau} = \frac{\partial \sigma}{\partial U}\)——温度的统计力学定义。
  • 第二行第二列: \(-P = \frac{\partial U}{\partial V}\)——压强的力学定义(一个广义力)。
  • 第三行第三列: \(\mu = \frac{\partial F}{\partial N}\)——化学势的实用定义(在恒温 \(T\) 下增加一个粒子所付的自由能代价)。

“唯一的真理”与冗余

这 9 个格子你都得背下来吗?谢天谢地,不用。后两列严格说是冗余的。所有物理都住在第一列,那个熵列里。

  • 一旦你拥有 \(\sigma = \sigma(U, V, N)\),你就拥有了整个宇宙。
  • 第二列(\(U\))只是它的反函数。
  • 第三列(\(F\))只是它的勒让德变换。

物理学家把三列都留着,不是为了多出什么物理,而是图坐标上的方便——现实里,把 \(\tau\) 固定住,远比把 \(\sigma\) 固定住容易。

数学桥梁

怎么从第一列跨到第二列?纯粹是微积分,别无其他。

桥梁 1:倒数关系。 对比第一行:\(\frac{1}{\tau} = (\frac{\partial \sigma}{\partial U})\)\(\tau = (\frac{\partial U}{\partial \sigma})\)。这无非是反函数法则 \(\frac{d y}{d x} = \frac{1}{\frac{d x}{d y}}\),证明两种温度定义彼此相符。

桥梁 2:循环关系。 对比第三行:\(-\frac{\mu}{\tau} = (\frac{\partial \sigma}{\partial N})_U\)\(\mu = (\frac{\partial U}{\partial N})_\sigma\)。那个别扭的 \(-\frac{\mu}{\tau}\) 是哪来的?三项乘积法则:

$$\left(\frac{\partial y}{\partial x}\right)_z = - \frac{\left(\frac{\partial z}{\partial x}\right)_y}{\left(\frac{\partial z}{\partial y}\right)_x}$$

\(y=\sigma, x=N, z=U\)

$$\left(\frac{\partial \sigma}{\partial N}\right)_{U} = - \frac{\left(\frac{\partial U}{\partial N}\right)_\sigma}{\left(\frac{\partial U}{\partial \sigma}\right)_N}$$

代入第二列的定义(分子 \(\mu\),分母 \(\tau\)),便得 \(\left(\frac{\partial \sigma}{\partial N}\right)_{U} = -\frac{\mu}{\tau}\)。那个负号和对 \(\tau\) 的除法,是数学上的必然,不是新的物理定律。

第六章:大逆转 —— 从配分函数到吉布斯求和

现在我们可以把引擎倒着开了。回想那个计划:写能量模型,算 \(Z\),得 \(F = -\tau\ln Z\),推导一切。我们来把它搭起来。

玻尔兹曼分布

把一个系统 \(S\) 放进和一个庞大热库 \(R\) 的接触里。它们交换能量,停在共同的 \(\tau\) 上。妙处在这个设问方式:我们把 \(S\)一个精确微观态钉死(于是它自己的多重度为 1),然后问这个态有多大概率。既然 \(S\) 的状态固定了,所有的自由都落在热库身上——概率完全由热库有多少种方式容纳剩下的能量来支配。于是概率之比就是热库多重度之比:

$$\frac{P(E_1)}{P(E_2)} = \frac{g_R(U_{total} - E_1)}{g_R(U_{total} - E_2)}$$

\(g = e^\sigma\),对热库的熵做泰勒展开(\(\sigma_R(U - E) \approx \sigma_R(U) - E/\tau\))。常数项在分子分母间抵消,幸存下来的,是那个著名的玻尔兹曼因子

$$P(E) \propto e^{-E / \tau}$$

一段信息熵的旁白: 如果你换个问法,要在已知平均能量 \(U\) 的前提下,求那个最无偏的概率分布——即最大化香农熵 \(\sigma_I = - \sum p_n \ln p_n\) 的那个——拉格朗日乘数法会把 \(p_n = C e^{-E_n/\tau}\) 原样还给你。玻尔兹曼分布,在一个精确的意义上,是大自然最诚实的猜测:和我们所知相容的、最平摊的那个分布。

配分函数:统计力学的”圣杯”

要把正比关系 \(P(E_s) \propto e^{-E_s / \tau}\) 变成等式,我们需要一个归一化的因子——配分函数 \(Z\)

$$Z \equiv \sum_{\text{所有状态 } s} e^{-E_s / \tau}$$

于是 \(P(s) = \frac{1}{Z} e^{-E_s / \tau}\)。为什么把它封作”圣杯”?写 \(P(s) = e^{(F - E_s)/\tau}\),要求概率之和为 1(即 \(e^{F/\tau} Z = 1\))。微观求和与宏观势之间的桥梁就掉了出来:

$$F = - \tau \ln Z$$

就是我们许下的那个计算大回旋,整条反向通道,四步走完:

  1. 写出微观模型——列出能级 \(E_n\)
  2. 求和成 \(Z\)——通常是几何级数或高斯积分。
  3. 读出 \(F\),由 \(F = -\tau \ln Z\)
  4. 得到其余一切——\(\sigma\)\(P\)\(\mu\)\(U\)——靠对 \(F\) 求导。

我们从头到尾都没直接去数 \(g\)。我们只算一个和,热力学就从它那里铺展开来。

吉布斯求和:当粒子流动时

现在让系统不光交换能量、也交换粒子——巨正则的设定。这时热库的熵 \(\sigma_R\) 要对能量和粒子数两者的变化作出响应,于是泰勒展开多出第二项:

$$\Delta \sigma_R = \frac{1}{\tau} (-E) + \left(-\frac{\mu}{\tau}\right) (-N) = \frac{\mu N - E}{\tau}$$

这给出吉布斯分布

$$P(N, E) \propto e^{(\mu N - E) / \tau}$$

其归一化因子是巨配分函数,也叫吉布斯求和 \(\mathcal{Z}\)

$$\mathcal{Z} \equiv \sum_{N=0}^{\infty} \sum_{s(N)} e^{(\mu N - E_s) / \tau}$$

引入绝对活性 \(\lambda \equiv e^{\mu/\tau}\),把它收拾干净:

$$P(N, E) = \frac{1}{\mathcal{Z}} \lambda^N e^{-E / \tau}$$

一个悖论? 在玻尔兹曼那里,\(e^{-E/\tau}\) 说能量越高概率越低——合情合理。可在吉布斯这里,\(e^{\mu N / \tau}\) 似乎在说粒子越多概率越高。系统岂不要撑破肚皮、爆炸了?不会——因为对经典理想气体,化学势 \(\mu\) 是个很大的数,于是每多一个粒子是受罚、而非受赏。危机解除。

第七章:量子统计 —— 当粒子失去姓名

我们终于要直面那个绕不开的难题:量子力学。 一路走来,我们都在不动声色地把粒子当成贴了标签的小台球——粒子 A 在这儿,粒子 B 在那儿。可全同的量子粒子根本没有标签。我们必须面对不可区分性

它是谁?不——它是什么

在量子世界里,同类粒子在根本上是可以互换的。交换其中两个,概率密度 \(|\Psi|^2\) 不能改变,这就逼出 \(\Psi(1, 2) = \pm \Psi(2, 1)\)。那个不起眼的 \(\pm\) 把宇宙劈成两半:

  • 玻色子(+): 爱热闹,乐于挤进同一个状态(光子、氦-4)。
  • 费米子(−): 反社交,受泡利不相容原理约束,绝不两个挤在一个状态里(电子、质子)。

而这改变了计数,计数才是全部要害。把 2 个粒子放进 4 个能级:

  • 经典: \(4^2 = 16\) 种方式。
  • 玻色子: 10 种方式(抱团被增强)。
  • 费米子: 6 种方式(共享一个状态被彻底禁止)。

既然 \(g\) 本身变了,下游的 \(\sigma\) 和每一项宏观性质也跟着变。量子不可区分性不是个脚注——它改写了我们整套框架赖以为食的那些计数。

从粒子视角到轨道视角

驯服这件事的窍门,是换个问法。别再问”粒子 A 在哪儿?“——粒子根本没名字可追。改问:“能量为 \(\varepsilon\) 的能级上坐了几个粒子?“由于一个能级的占据数可以变动,这恰恰是个粒子数涨落的情形——正是吉布斯求和当初为之而生的那种。

费米–狄拉克分布与玻色–爱因斯坦分布

费米子。 泡利只允许占据数 \(n = 0\)\(1\),所以单个能级的吉布斯求和只有两项:

$$\mathcal{Z} = 1 + \lambda e^{-\varepsilon/\tau}$$

平均占据数算出来是

$$\langle n \rangle = \frac{1}{e^{(\varepsilon-\mu)/\tau} + 1}$$

分母里那个 \(+1\),是费米子”社交焦虑”的数学指纹——无论你怎么逼,它都把 \(\langle n \rangle\) 封顶在 1。

玻色子。 这下 \(n\) 可以是 \(0, 1, \dots, \infty\),吉布斯求和是一条几何级数:

$$\mathcal{Z} = \frac{1}{1 - \lambda e^{-\varepsilon/\tau}}$$

给出

$$\langle n \rangle = \frac{1}{e^{(\varepsilon-\mu)/\tau} - 1}$$

孤零零一个符号的翻转——\(-1\) 而非 \(+1\)——后果却惊天动地:当 \(\varepsilon \approx \mu\) 时,分母趋于零,占据数发散。这场失控的堆积,就是玻色–爱因斯坦凝聚。

经典极限:殊途同归

那为什么日常生活从不向我们揭示这道费米/玻色的大裂缝?因为我们深居经典区间——高温或低密度,\((\varepsilon - \mu)/\tau \gg 1\)。在那里指数项把 \(\pm 1\) 压成了一个舍入误差:

$$\frac{1}{e^{(\varepsilon-\mu)/\tau} \pm 1} \approx e^{\mu/\tau} e^{-\varepsilon/\tau}$$

——而这正是玻尔兹曼分布。量子与经典两幅图景,恰恰在该融合的地方严丝合缝地融合了。可一旦把气体冷到量子浓度 \(n_Q\) 以下,那个 \(\pm 1\) 就轰然复活:费米子堆叠成一片费米海(撑住白矮星、抵抗引力的简并压),而玻色子则一齐坍缩成超流体。同样这两个符号,在室温下隐身不见,到了极端处却成了奇异物质的建筑师。

结语:从无知中涌现的秩序

我们的旅程在此告一段落,那就回头看看我们究竟建起了什么——因为它的形状,本身就是要点所在。

  1. 我们从承认无知出发:等概率原理,那句”我们不能也不会去追踪每一个粒子”的自白。
  2. 我们经由计数,把”可能性”变成了可量化的东西。
  3. 通过最大化这个熵,我们看着一座座宏观岛屿——温度、压强、化学势——从混沌的微观海洋里升起。
  4. 我们学会以自由能指引航向,以配分函数为罗盘来计算。
  5. 最后,注入量子本性,同一套框架预言了物质在最极端条件下所取的状态。

注意我们几乎从没做过的一件事:点名某种物质。我们从没算过铜的比热,也没算过氧分子的转动。这恰恰是统计力学的魅力——它是一套普适的语法,而不是一个特定的故事。 它说的是的逻辑,从不说的身份。气体里的分子、腔内的光子、恒星心脏里的核物质:把哈密顿量喂进去,宏观行为就自己组装起来。

于是我们绕回了第一页那杯静止的水——只是如今,它的平静读起来已是另一番意味。那份宁静并非风暴的缺席;它是风暴的统计签名,是 \(10^{23}\) 种盲目而等可能的可能性合谋酿出的、不可避免的秩序。这或许是物理学能献上的最深的一句诗:上帝不需要掷骰子——混沌本身就是秩序的源泉。