该笔记主要参考热力学部分——《热学热力学与统计物理（上册）》（曹烈兆周子舫），统计物理部分——《热学热力学与统计物理（下册）》（周子舫曹烈兆），以及《热力学与统计物理》（汪志诚）。
与Gemini Pro 🤖共同完成

第1章热力学的基本规律#

1.1 热力学状态与第零定律#

1.1.1 热力学平衡状态#

要建立热力学这座大厦，第一步是确定我们要研究的对象——平衡态。宏观体系包含极其巨大的粒子数，试图追踪每一个粒子的力学运动是绝望的。热力学的聪明之处在于，它只关注系统在不受外界干扰、经过足够长时间演化后，最终“稳定”下来的状态。这里的“稳定”绝不是微观上的死寂，而是宏观量不再随时间变化，而微观粒子依然在剧烈且无序地运动，即所谓 “宏观静止，微观涨落”。
从非平衡走向平衡所需的时间尺度，我们称为弛豫时间 ( $\tau$ )。这个概念在物理图像中极其关键。在实际的物理过程里，所有的测量和操作所经历的时间，都必须与弛豫时间进行对比。只有当外界条件变化的时间尺度( $t$ )远大于 $\tau$ 时（ $t \gg \tau$ ），系统才有时间对外界变化做出响应并随时重新达到平衡，我们才能把整个过程视为一连串平衡态的连缀（即准静态过程），此时定义出来的热力学量才是确切且有意义的。

1.1.2 热力学第零定律与温度#

既然定义了单个系统的平衡态，我们自然会进一步考察多个系统之间的相互作用。如果把两个独立的系统进行热接触，它们最终会达到一个共同的新平衡态。那么，如何在接触前，就预判它们之间会不会发生热量交换？这就需要为“热平衡”寻找一个普适的宏观判据。
为了从逻辑上严密地建立这个判据，物理学引入了热力学第零定律：
“如果两个热力学系统中的每一个都与第三个热力学系统处于热平衡，那么它们也必定互相处于热平衡。”
从日常经验上看，这似乎是一句理所当然的废话，但在物理逻辑上，它建立了一个严格的数学“等价关系”。这种等价关系的物理本质在于，它强行要求必须存在一个共有的状态函数——温度 ( $T$ )。
第零定律告诉我们，温度不仅仅是感知上的“冷热程度”，更是决定两个系统接触时是否会发生热交换的唯一宏观判据。有了这个定律兜底，我们就可以理直气壮地选取某个特定系统（温度计）的特定物理属性（如水银的体积）作为标准，建立起经验温标。

1.1.3 物态方程与热力学系数#

通过第零定律确立了温度 $T$ 后，加上我们在力学和几何上熟知的压强 $p$ 和体积 $V$ 等宏观参量，我们就有了一套描述系统的语言。对于一个确定的均匀物质系统，其实际的经验告诉我们，只要给定其中几个独立的参量，它的平衡态也就随之唯一确定了。由此引入物态方程（通常写为 $f(p, V, T) = 0$ ）对物质特定平衡状态进行宏观描述。它告诉我们，在这个特定的平衡态下，这几个宏观量是如何具体联系在一起的。
然而在真实的物理世界里，除了结构简单的气体，对于大多数固体或液体，我们几乎不可能从理论上写出一个完美的、全局适用的物态方程解析式。那在实验室里，引入三个极其重要的宏观实验测量量（热力学系数）：
1. 体胀系数 (Volume Expansivity) $\alpha = \frac{1}{V} \left( \frac{\partial V}{\partial T} \right)_p$
  在压强保持不变（等压过程）的条件下，温度每升高 1K 时，系统体积的相对增加量。它直接反映了物质“热胀冷缩”的剧烈程度。
2. 等温压缩系数 (Isothermal Compressibility) $\kappa_T = -\frac{1}{V} \left( \frac{\partial V}{\partial p} \right)_T$
  在温度保持不变（等温过程）的条件下，压强每增加 1Pa 时，系统体积的相对缩小量。
3. 压强系数 (Pressure Coefficient) $\beta = \frac{1}{p} \left( \frac{\partial p}{\partial T} \right)_V$
  在体积保持不变（等容过程）的条件下，温度每升高 1K 时，系统内部压强的相对增加量。它反映了如果把物质锁死在一个刚性容器里加热，其内部压强随温度攀升的敏锐程度。
这三个系数构成了联系物质宏观状态变化的“铁三角”，通过数学上的偏导数循环关系，它们满足严格的制约方程： $\alpha = p \beta \kappa_T$ (推到利用了 $\left( \frac{\partial V}{\partial p} \right)_{T} \left( \frac{\partial p}{\partial T} \right)_{V} \left( \frac{\partial T}{\partial V} \right)_{p} = -1$ )。在实际研究中，只要测出其中两个，就能直接算出第三个。在不知道复杂物质确切物态方程的情况下，只要我们在实验室里测出了该物质的这三个表观系数，就能在局部范围内完全掌握它的状态演化规律。
虽然复杂物质极度依赖这些实验系数，但对于最简单的一类物质——气体，我们确实可以建立直观的宏观物态方程模型。最简单的玩具模型是理想气体物态方程（ $pV = nRT$ ），它在极高温度和极低压强下成立，代价是直接抹除了微观细节：把粒子视为没有体积的几何点，并切断了所有相互作用。
理想气体物态方程的实验推导
理想气体的宏观物态方程是直接建立在三个控制变量的经典物理实验之上：
玻意耳定律（等温实验）：保持气体温度 $T$ 和物质的量 $n$ 不变，压缩或膨胀气体。实验观测发现，气体的体积与压强严格成反比： $V \propto \frac{1}{p}$ 。
查理-盖·吕萨克定律（等压实验）：保持气体压强 $p$ 和物质的量 $n$ 不变，给气体加热。实验观测发现，气体的体积与绝对温度成正比： $V \propto T$ 。
阿伏伽德罗定律（粒子数实验）：在相同的温度 $T$ 和压强 $p$ 下，任何气体的体积都直接正比于其包含的分子总数（即物质的量 $n$ ）： $V \propto n$ 。
既然气体的体积 $V$ 既正比于温度 $T$ 和物质的量 $n$ ，又反比于压强 $p$ ，我们在数学上可以引入一个普适的比例常数——理想气体常数 $R$ （由实验精确标定得出）,直接将这三个独立的比例关系合并为一个统一的表达式：
$pV = nRT$
但只要稍微靠近真实的物理世界，比如引入范德瓦尔斯方程： $(p + \frac{an^2}{V^2})(V - nb) = nRT$ ，物理图像立马就丰满了。等式左边多出来的这两个修正项，本质上是对微观粒子间相互作用的“粗视化”补偿：
内压项修正 ( $+\frac{an^2}{V^2}$ )：源于分子间的长程引力。靠近容器壁的分子被内部的分子往回拉，导致它们撞击器壁的实际压强变小了。因此，我们在宏观测量的压强 $p$ 上，必须把它损失的这部分给“加回来”，才能描述系统内部真实的受力状态。
共体积修正 ( $-nb$ )：源于分子间的短程斥力。粒子本身占据了一定的体积，存在不可穿透的硬核，导致气体分子真正能自由活动的空间其实比容器的宏观总体积 $V$ 要小，所以要减去一个“禁带”体积。
范德瓦尔斯方程修正项推导
设理想气体状态方程为 $p_{id} V_{id} = nRT$ 。实际气体在微观上具有分子本身体积（表现为短程斥力）和分子间相互作用力（表现为长程引力）。
1. 体积修正（共体积项）
假设气体分子为直径为 $d$ 的刚性球，单个分子的真实体积为：
$v_0 = \frac{4}{3}\pi \left(\frac{d}{2}\right)^3 = \frac{1}{6}\pi d^3$
当两个分子发生碰撞时，它们的中心距离最小为 $d$ 。即以任意一个分子中心为球心，半径为 $d$ 的球状空间内，不能有其他分子中心进入。两个分子共同的“排除体积”（Excluded Volume）为：
$V_{excl} = \frac{4}{3}\pi d^3 = 8 \left(\frac{1}{6}\pi d^3\right) = 8v_0$
分摊到平均每个分子上的排除体积为：
$v_b = \frac{1}{2} V_{excl} = 4v_0$
对于 $1 \text{ mol}$ 气体（ $N_A$ 个分子），其共体积常量 $b$ 为：
$b = N_A v_b = 4 N_A v_0$
对于 $n \text{ mol}$ 气体，分子实际可自由运动的有效体积 $V_{id}$ 为系统总体积 $V$ 减去所有分子的排除体积：
$V_{id} = V - nb$
2. 压强修正（内压项）
处于容器内部的分子受周围分子引力对称，合力为零。但处于器壁表面厚度为分子引力有效作用半径的表面层内的分子，受到内部气体的净引力（指向容器内部）。导致实际气体对器壁的碰撞压强 $p$ 相比理想压强 $p_{id}$ 减小了一个“内压” $\Delta p$ ：
$p = p_{id} - \Delta p \implies p_{id} = p + \Delta p$
内压 $\Delta p$ 的大小取决于两个因素：
1. 表面层内受到内部引力拉扯的分子数密度 $\rho_{surf}$ ；
2. 内部对表面层分子施加引力的分子数密度 $\rho_{int}$ 。
由于在平衡态下，各处分子数密度是均匀的，即 $\rho_{surf} \propto \frac{n}{V}$ ， $\rho_{int} \propto \frac{n}{V}$ 。因此，内压 $\Delta p$ 正比于密度的平方：
$\Delta p \propto \rho_{surf} \cdot \rho_{int} \propto \left(\frac{n}{V}\right) \cdot \left(\frac{n}{V}\right) = \left(\frac{n}{V}\right)^2$
引入引力比例常数 $a$ ，得到内压项：
$\Delta p = a \frac{n^2}{V^2}$ $p_{id} = p + a \frac{n^2}{V^2}$
将 $p_{id}$ 和 $V_{id}$ 的修正表达式代入理想气体状态方程 $p_{id} V_{id} = nRT$ ，得到范德瓦尔斯方程：
$\left( p + a \frac{n^2}{V^2} \right) (V - nb) = nRT$

1.2 热力学第一定律与能量转化#

1.2.1 功与内能#

要改变一个热力学系统的状态，宏观上只有两种手段：做功（如压缩气体、搅拌液体）和传热（如放在火上加热）。功 ( $W$ ) 和热量 ( $Q$ ) 都是伴随“过程”而生的，它们不仅取决于系统的初末状态，还严重依赖于系统经历的具体路径。在数学上，它们不是全微分（通常用 $\delta W$ 和 $\delta Q$ 表示）。
在热力学中，功不仅局限于机械体积功（ $p dV$ ）。任何宏观做功过程都可以写成广义力 ( $Y$ ) 与广义位移 ( $dx$ ) 的乘积： $\delta W = Y dx$ 。我们需要先计算外界对系统做的功 $đW_{ex}$ ，则系统对外的元功 $\delta W = - đW_{ex}$ 。接下来以电磁介质为例，展示元功的推导过程。在处理电磁介质时，最核心的物理直觉是：必须把“建立真空电磁场消耗的能量”和“真正改变介质内部状态（极化/磁化）所消耗的能量”剥离开来。 只有后者才属于电/磁介质这个热力学系统的“内能”变化。
1. 电介质的极化功
物理模型：考虑一个充满各向同性电介质的平行板电容器。极板面积为 $A$ ，间距为 $d$ ，介质体积 $V_0 = Ad$ 。
外界总电功：当外界电源给极板微小充电 $dq$ 时，极板间电势差为 $U$ 。外界做功为：
$đW_{ex} = U dq$
引入场量：根据电磁学，电势差 $U = Ed$ ，极板自由电荷 $q = DA$ （ $D$ 为电位移）。对其微分有 $dq = A dD$ 。代入上式：
$đW_{ex} = (Ed)(A dD) = V_0 E dD$
剥离真空场能量：介质中的电位移本构关系为 $D = \varepsilon_0 E + P$ （ $P$ 为极化强度）。代入微分式：
$đW_{ex} = V_0 E d(\varepsilon_0 E) + V_0 E dP = d\left(\frac{1}{2}\varepsilon_0 E^2 V_0\right) + V_0 E dP$
等式右边第一项 $d(\frac{1}{2}\varepsilon_0 E^2 V_0)$ 是即使电容器里是真空也必须提供的能量，它储存在真空电场中，与电介质材料的自身热力学状态无关。等式右边第二项 $V_0 E dP$ 才是外界电场真正为了抵抗分子内部束缚力、拉扯电荷中心产生极化，从而对电介质材料本身做的功。令系统总电偶极矩 $p = V_0 P$ ，则外界对介质做的极化功为 $E dp$ 。因此，电介质系统对外所做的元功表达式为：
$\delta W = - E dp$
1. 磁介质的磁化功
物理模型：考虑一个内部充满磁介质的细长螺线管。截面积为 $A$ ，长度为 $l$ ，体积 $V_0 = Al$ ，总匝数为 $N$ 。
外界总电磁功：当线圈电流 $I$ 发生变化时，会引起磁通量 $\Phi = NBA$ 的变化，从而产生感应电动势 $\mathcal{E} = -d\Phi/dt$ 。外界电源为了维持电流，必须克服感应电动势做功：
$đW_{ex} = -\mathcal{E} I dt = I d\Phi = I N A dB$
引入场量：根据安培环路定理，磁场强度 $H = NI / l$ ，即 $IN = Hl$ 。代入上式：
$đW_{ex} = (Hl) A dB = V_0 H dB$
剥离真空场能量：介质中的磁感应强度本构关系为 $B = \mu_0 (H + M)$ （ $M$ 为磁化强度）。代入微分式：
$đW_{ex} = V_0 H d(\mu_0 H) + V_0 H d(\mu_0 M) = d\left(\frac{1}{2}\mu_0 H^2 V_0\right) + \mu_0 V_0 H dM$
同样地，第一项 $d(\frac{1}{2}\mu_0 H^2 V_0)$ 是建立真空磁场所需的能量。第二项 $\mu_0 V_0 H dM$ 才是外界磁场为了使介质内部微观磁矩发生偏转定向，从而对磁介质材料本身做的功。令系统总磁矩 $m = V_0 M$ ，则外界对磁介质做的磁化功为 $\mu_0 H dm$ 。因此，磁介质系统对外所做的元功表达式为：
$\delta W = - \mu_0 H dm$
但物理学家在实验（焦耳的绝热做功实验）中发现：在一个绝热系统中（没有热量交换， $\delta Q=0$ ），无论你用什么方式对系统做机械功，只要系统的初末状态确定了，外界对系统所做的总功总是完全相等的。
这个实验事实意味着，虽然功 $W$ 本身是依赖路径的，但在“绝热”这个特定限制下，它变成了一个只由初末状态决定的量。在物理逻辑上，这就强制要求系统内部必须存在一个状态函数，它的变化量刚好等于绝热过程中的功。将这个状态函数定义为内能 ( $U$ )，其是系统内在的、由微观粒子无规则运动动能和分子间相互作用势能总和所构成的宏观状态函数。

1.2.2 热力学第一定律#

在真实的物理世界中，系统不仅可以与外界交换机械功，还可以通过温差直接交换热量，热量的进入显然也会引起系统状态的变化。宏观物理学给出了一个极其坚定的论断——无论是通过做功传递的宏观机械能，还是通过传热传递的微观分子热运动能量，它们在改变系统内能上是完全等效的，且总能量必须守恒。这就是热力学第一定律的物理内核，是能量守恒定律在涉及热现象的宏观过程中的直接表述。
用数学语言来刻画这个物理图景，就是：系统吸收的热量 ( $\delta Q$ )，一部分用于增加自身的内能 ( $dU$ )，另一部分用于对外做膨胀功 ( $\delta W$ )，即：
$\delta Q = dU + \delta W$
在这个微分方程中，左边的传热和右边的做功都是依赖路径的“过程量”，但神奇的是，当它们按照这个守恒法则组合在一起时，其差值 ( $dU$ ) 却是一个与路径无关的态函数。

1.2.3 焓 ( $H$ ) 与热容#

如果把系统锁死在一个刚性容器里（等容过程， $dV = 0$ ），系统对外不做体积功，吸收的热量就完完全全变成了内能的增加： $\delta Q_V = dU$ 。但在实验室里，甚至在大自然中，绝大多数的物理和化学过程并不是等容的，而是等压的（比如在敞口烧杯中发生的反应，系统始终承受着恒定的大气压）。系统在吸热的同时往往伴随着体积的膨胀，这意味着它必须推开周围的空气做功（ $\delta W = p dV$ ）。
把这个恒定压强下的做功代入第一定律，我们会得到 $\delta Q_p = dU + p dV$ 。因为压强 $p$ 是常数，它可以被毫无阻碍地塞进微分号里，变成：
$\delta Q_p = d(U + pV)$
物理学家为了处理这类极其常见的等压过程，直接把括号里这个由态函数组合而成的新态函数 $(U + pV)$ 定义为焓 ( $H$ )。它的出现完全是为了物理分析上的便利：在等压过程中，系统吸收的热量直接等于其焓的增量 ( $\delta Q_p = dH$ )。
为了描述系统吸热量与温度变化的关系，我们定义了两个宏观可测的热容：
定容热容 ( $C_V$ )： $C_V = \left(\frac{\partial U}{\partial T}\right)_V$
定压热容 ( $C_p$ )： $C_p = \left(\frac{\partial H}{\partial T}\right)_p$
等压升温时，系统吸收的热量不仅要用于增加内能（升温），还要分出一部分能量去抵抗外压做膨胀功，这导致宏观上 $C_p$ 必然大于 $C_V$ 。
理想气体内能表达式的推导
在焦耳自由膨胀实验中，理想气体向绝热的真空容器中膨胀。因为绝热（ $\delta Q = 0$ ）且向真空膨胀不克服任何阻力（不对外做功， $\delta W = 0$ ），根据热力学第一定律，系统的内能在膨胀前后保持恒定（ $dU = 0$ ）。同时，实验精确观测到气体的温度也未发生改变（ $dT = 0$ ）。既然体积 $V$ 经历了剧烈增加，而内能 $U$ 和温度 $T$ 却都不变，这在物理上直接导出了一个坚实的结论（焦耳定律）：理想气体的内能与体积完全无关，它仅仅是温度的单变量函数，即 $U = U(T)$ 。
既然明确了 $U$ 只是 $T$ 的函数，我们在定义定容热容 $C_V = \left(\frac{\partial U}{\partial T}\right)_V$ 时所用的偏导数，就可以理直气壮地退化为全导数形式：
$C_V = \frac{dU}{dT} \implies dU = C_V dT$
对该微分方程进行积分，即可直接得到理想气体的宏观内能表达式：
$U(T) = \int C_V dT + U_0$
（其中 $U_0$ 为积分常数，代表参考态下的内能。由于理想气体的定容热容 $C_V$ 在相当宽的温度范围内可视为常数，该表达式在实际计算中通常直接简化为 $U = C_V T + U_0$ ）。

1.2.4 理想气体的绝热过程#

绝热过程的定义是系统与外界无热交换，即 $\delta Q = 0$ 。根据热力学第一定律 $\delta Q = dU + \delta W$ ，可直接得到微积分方程：
$dU = -\delta W \implies C_V dT = -p dV$
气体对外膨胀做功（ $p dV > 0$ ）的能量，完全由消耗系统自身的内能（ $C_V dT < 0$ ，温度降低）来支付。将理想气体状态方程 $pV = nRT$ 两边取全微分：
$p dV + V dp = nR dT$
考虑等压升温比等容升温多吸收的热量，刚好等于气体对外做的膨胀功： $C_p dT = C_V dT + p dV$ ，由理想气体状态方程推得等压膨胀功为 $p dV = nR dT$ 。代入上式并约掉 $dT$ ，得 $C_p - C_V = nR$ ，利用此关系代入上式消去 $nR$ ：
$p dV + V dp = (C_p - C_V) dT$
将上述 $dT$ 的表达式代回第一定律的方程 $C_V dT = -p dV$ 中并整理等式两边：
$C_V V dp + C_p p dV = 0$
为了求解这个微分方程，将两边同除以 $C_V p V$ ，并引入绝热指数 $\gamma = \frac{C_p}{C_V}$ ：
$\frac{dp}{p} + \gamma \frac{dV}{V} = 0$
对等式两边直接积分：
$\ln p + \gamma \ln V = \text{常数}$
最终得到绝热过程的核心特征方程（泊松方程）：
$p V^\gamma = \text{常数}$
对特征方程求导可得绝热线的斜率为 $\frac{dp}{dV} = -\gamma \frac{p}{V}$ ，而等温线的斜率为 $-\frac{p}{V}$ 。因为 $\gamma > 1$ ，绝热线比等温线更陡峭。在绝热膨胀中，压强的下降由双重因素叠加导致——不仅是因为体积变大（空间变得稀疏），更是因为内能转化为功导致了温度的剧烈下降（分子热运动动能锐减）。

1.3 第二定律与不可逆性#

热力学第一定律确立了能量在转化过程中的守恒关系，但并未对过程的演化方向做出任何限制。在实际的宏观物理世界中，能量的转化和传递具有明确的方向性，第二定律正是对这种宏观不可逆性的严格物理界定。

1.3.1 过程的方向性与第二定律的表述#

自然界中的自发宏观过程均具有单向性。例如，机械功可以自发且无代价地完全转化为热（如摩擦生热），而热量却不能在不产生其他宏观影响的条件下完全转化为功；同样，热量必然自发地从高温物体传向低温物体，而无法自动逆向传递。为了对这种方向性进行确切的物理限定，热力学第二定律给出了两种经典的表述：
1. 开尔文表述（针对热功转换）：不可能从单一热源吸取热量，使之完全变为有用功而不产生其他影响。
  实质：确立了热机效率的上限必然小于 100%。热机在持续输出机械功的过程中，必须向低温热源排放废热。
2. 克劳修斯表述（针对热量传递）：不可能把热量从低温物体传到高温物体而不产生其他影响。
  实质：指明了热量传递的自发方向。制冷机若要实现热量的逆向传递，必须依赖外界对其做功。
这两种表述在物理上是完全等价的。若其中之一不成立，便可以推出另一个也不成立。它们从不同视角共同确立了宏观热力学过程的不可逆性。
1. 若开尔文表述不成立，则必然导致克劳修斯表述不成立
假设存在一台违反开尔文表述的“完美热机”（记为机 A）。该热机能够从单一的高温热源（温度为 $T_1$ ）吸收热量 $Q$ ，并将其完全转化为宏观机械功 $W$ （即 $W=Q$ ），且不对外界产生任何其他影响。
现引入一台遵循物理规律的普通制冷机（记为机 B），使其在相同的高温热源 $T_1$ 和低温热源 $T_2$ 之间工作。利用机 A 输出的功 $W$ 来直接驱动机 B。机 B 消耗功 $W$ ，从低温热源吸取热量 $Q_2$ ，并向高温热源排放热量 $Q_1'$ 。根据第一定律的能量守恒， $Q_1' = W + Q_2 = Q + Q_2$ 。
将机 A 与机 B 视为一个闭合的联合系统。分析该系统完成一个循环后的宏观净效应：
外界净做功：联合系统对外输出功 $W$ （机 A），同时内部消耗功 $W$ （机 B），与外界的净功交换为零。 低温热源 $T_2$ ：被机 B 吸取了热量 $Q_2$ 。 高温热源 $T_1$ ：向机 A 输出了热量 $Q$ ，同时接收了机 B 排放的热量 $Q + Q_2$ 。其净效应为接收了热量 $Q_2$ 。
综合来看，该联合系统的唯一最终效果是：自发地将热量 $Q_2$ 从低温热源传到了高温热源，而没有消耗外界的功，也没有引起其他任何状态变化。这直接违反了克劳修斯表述。
2. 若克劳修斯表述不成立，则必然导致开尔文表述不成立
假设存在一台违反克劳修斯表述的“完美制冷机”（记为机 C）。该制冷机能够在不消耗任何外界功的条件下，自动从低温热源 $T_2$ 吸取热量 $Q_2$ ，并将其全部排入高温热源 $T_1$ 。
现引入一台遵循物理规律的普通热机（记为机 D），使其在相同的两热源间工作。设定机 D 的规模，使其在每个循环中恰好向低温热源 $T_2$ 排放热量 $Q_2$ 。同时，机 D 从高温热源吸取热量 $Q_1$ ，并对外做功 $W$ 。根据第一定律， $W = Q_1 - Q_2$ 。
将机 C 与机 D 构成一个联合系统。分析其宏观净效应：
低温热源 $T_2$ ：机 D 向其排放热量 $Q_2$ ，机 C 恰好从其吸取热量 $Q_2$ 。该热源的净热量交换为零，状态完全未变。 高温热源 $T_1$ ：向机 D 输出了热量 $Q_1$ ，同时接收了机 C 排放的热量 $Q_2$ 。其净效应为向系统输出了热量 $Q_1 - Q_2$ 。 外界净做功：系统对外输出净功 $W = Q_1 - Q_2$ 。
综合来看，该联合系统的唯一最终效果是：从单一热源 $T_1$ 吸收了热量 $(Q_1 - Q_2)$ ，并将其完全转化为有用的机械功 $W$ ，且未对低温热源或外界环境产生任何其他影响。这直接违反了开尔文表述。

1.3.2 卡诺循环与卡诺定理#

在第二定律确立了热功转化的方向性约束后，确定热机在给定高低温热源之间工作的理论效率上限便成为核心问题。为此，卡诺引入了由两个准静态等温过程和两个准静态绝热过程构成的理想可逆循环——卡诺循环。
基于卡诺循环模型及第二定律的开尔文表述，通过逻辑反证法可以严格推导出卡诺定理：
1. 在相同的高温热源和相同的低温热源之间工作的一切实际（不可逆）热机，其效率都不可能大于可逆热机。
2. 在相同的高温热源和低温热源之间工作的一切可逆热机，其效率都相等，且与工作物质的具体性质无关。
卡诺定理为热功转换设定了不可逾越的理论极限。其第二条推论表明，可逆热机的效率 $\eta$ 仅由高、低温热源的状态决定，彻底摆脱了具体热机结构与工作介质属性的制约。

1.3.3 热力学绝对温标的建立#

卡诺定理中“可逆热机效率与工作物质无关”的结论，为建立完全独立于测温物质的温度标尺提供了严格的理论基础。
根据卡诺定理，在温度为 $t_1$ 和 $t_2$ 的两个热源间工作的可逆热机，其吸热量 $Q_1$ 与放热量 $Q_2$ 之比 $\frac{Q_1}{Q_2}$ 必然是一个不依赖于具体物质的普适函数： $\frac{Q_1}{Q_2} = f(t_1, t_2)$ 。
开尔文利用这一普适关系，直接规定两个热源的绝对温度（记为 $T_1$ 和 $T_2$ ）之比，等于在此两热源间工作的可逆卡诺热机吸收与放出的热量之比，即：
$\frac{T_1}{T_2} = \frac{Q_1}{Q_2}$
由此定义的温标称为热力学绝对温标（或开氏温标）。它建立在热力学第二定律的坚实基础之上，彻底排除了经验温标对特定测温物质（如水银、酒精等）物理属性的依赖。进一步地，若将理想气体作为工作物质代入卡诺循环进行推导，可证明热力学绝对温标与理想气体温标在数值上是严格等价的。

1.4 熵与热力学基本方程#

卡诺定理不仅确立了不依赖具体物质的绝对温标，其更深刻的数学推论在于揭示了一个全新的状态函数。本节将从热量传递的普遍规律中严格导出“熵”，并建立宏观热力学的基本微分方程。

1.4.1 克劳修斯等式与熵函数的导出#

在可逆卡诺循环中，工作物质从高温热源 $T_1$ 吸热 $Q_1$ ，向低温热源 $T_2$ 放热 $Q_2$ 。根据绝对温标的定义，有 $\frac{Q_1}{Q_2} = \frac{T_1}{T_2}$ 。若引入热力学通用的符号法则（系统吸热为正，放热为负），将放热量记为 $-Q_2$ ，则上式可改写为：
$\frac{Q_1}{T_1} + \frac{Q_2}{T_2} = 0$
对于任意复杂的闭合可逆循环，可通过引入无数条绝热线，将其分割为无数个微小的可逆卡诺循环。对所有微小循环求和取极限，可得到克劳修斯等式：
$\oint \frac{\delta Q_{\text{rev}}}{T} = 0$
这表明沿任意可逆闭合路径的积分为零。在数学上，这严格证明了从平衡态 A 到平衡态 B 的积分 $\int_A^B \frac{\delta Q_{\text{rev}}}{T}$ 其数值完全独立于具体的积分路径，而仅仅取决于系统的初末状态。既然该积分值只与状态有关，物理学便由此严格定义出一个新的状态函数——熵 ( $S$ )。其微元定义式为：
$dS = \frac{\delta Q_{\text{rev}}}{T}$
熵的出现，使得过程量热量 $\delta Q$ 能够通过积分因子 $\frac{1}{T}$ 转化为全微分 $dS$ 。

1.4.2 热力学基本方程#

将热力学第一定律 $\delta Q = dU + \delta W$ 与熵的微元定义相结合。对于一个只做可逆体积功的简单可压缩系统，有 $\delta W = p dV$ 且 $\delta Q = T dS$ 。代入第一定律，即得到热力学基本微分方程：
$dU = T dS - p dV$
该方程虽然是借助“可逆过程”的条件（ $\delta Q = T dS$ 和 $\delta W = p dV$ ）推导出来的，但由于方程中包含的物理量（ $U, T, S, p, V$ ）全部都是状态函数，状态函数的微小变化仅取决于相邻的初末平衡态。因此，该基本方程适用于简单系统相邻两个平衡态之间的任何过程（包括不可逆过程）。它是连接热力学第一定律与第二定律的核心枢纽。
为了将抽象的熵函数具体化，我们以理想气体为例，推导其宏观熵的解析表达式。这不仅是热力学计算的基础，也能直观地展示态函数如何通过状态参量来完全确定。
1. 以 $(T, V)$ 为独立变量的推导
从热力学基本微分方程出发：
$T dS = dU + p dV \implies dS = \frac{dU}{T} + \frac{p}{T} dV$
将理想气体的两个核心属性代入该式：
- 焦耳定律：内能仅为温度的函数， $dU = C_V dT$
- 物态方程： $pV = nRT \implies \frac{p}{T} = \frac{nR}{V}$
代入后得到熵的微元微分式：
$dS = C_V \frac{dT}{T} + nR \frac{dV}{V}$
假定在所讨论的温度范围内，定容热容 $C_V$ 为常数，对等式两边直接积分，得到以 $T, V$ 为自变量的熵函数表达式：
$S(T, V) = C_V \ln T + nR \ln V + S_0$
2. 以 $(T, p)$ 为独立变量的推导
在许多实际过程中（如等压过程），选取温度和压强作为独立变量更为方便。我们借助焓的微元定义 $dH = dU + p dV + V dp$ ，将基本方程改写为：
$T dS = dH - V dp \implies dS = \frac{dH}{T} - \frac{V}{T} dp$
同理，代入理想气体的焓变性质（ $dH = C_p dT$ ）与物态方程（ $\frac{V}{T} = \frac{nR}{p}$ ）：
$dS = C_p \frac{dT}{T} - nR \frac{dp}{p}$
在定压热容 $C_p$ 为常数的假设下，积分可得：
$S(T, p) = C_p \ln T - nR \ln p + S'_0$
式中的 $S_0$ 和 $S'_0$ 是积分常数。在经典宏观热力学框架内，熵只能确定到相差一个任意常数，我们只能计算过程前后的熵变 ( $\Delta S$ )。熵的绝对数值零点，必须等到引入热力学第三定律和量子统计物理后才能严格标定。上述简洁的对数表达式严格依赖于“热容是常数”这一近似。如果温度变化跨度极大（激发了分子的振动自由度），热容将成为温度的函数 $C_V(T)$ ，此时必须退回保留积分号的原始形式： $S = \int \frac{C_V(T)}{T} dT + nR \ln V + S_0$ 。
在 $S(T, V)$ 表达式中，体积项前是 正号 ( $+$ )。物理图像：体积膨胀，分子活动空间增大，系统微观状态数增多（混乱度增加），导致熵增。在 $S(T, p)$ 表达式中，压强项前是 负号 ( $-$ )。物理图像：压强增大意味着系统被压缩，分子活动空间受限，混乱度降低，导致熵减。
由于熵是态函数，不论理想气体经历的是可逆过程还是极度复杂的不可逆过程（如绝热自由膨胀），计算其初末态的熵变 $\Delta S$ 时，只需将初态和末态的参量代入上述公式直接相减即可，完全不需要考虑中间的实际路径。

1.4.3 克劳修斯不等式与熵增加原理#

对于包含不可逆过程的实际循环，根据卡诺定理，其效率必定小于相同热源间的可逆卡诺循环。经过严格的数学推演，可将克劳修斯等式推广为包含不可逆过程的克劳修斯不等式：
$\oint \frac{\delta Q}{T} \le 0$
（等号适用于全可逆循环，小于号适用于存在不可逆性的循环）。
若系统经任意过程（可逆或不可逆）从初态 A 演化至终态 B，结合熵的定义，可由克劳修斯不等式导出：
$S_B - S_A \ge \int_A^B \frac{\delta Q}{T}$
其微元形式为：
$dS \ge \frac{\delta Q}{T}$
熵增加原理：
将上述微元关系应用于孤立系统（系统与外界无物质和能量交换，即必定有 $\delta Q = 0$ ），不等式直接化简为：
$dS \ge 0$
孤立系统中发生的一切自发过程，必然朝着熵增加的方向进行；当系统的熵达到特定宏观约束下的极大值时，过程停止，系统进入稳定的热力学平衡态。熵增加原理为自然界宏观演化的方向性和限度，提供了绝对且严密的数学判据。

1.5 辅助特性函数的初步引入#

熵增加原理（ $dS \ge 0$ ）为宏观过程的方向性提供了判据，但其严格的适用前提是孤立系统（即系统的内能 $U$ 和体积 $V$ 保持不变）。在实际的物理和化学过程中，系统往往并非孤立，而是处于等温、等容或等压等特定约束条件下。为了在这些常见约束下寻找系统平衡的判据，需要引入新的热力学状态函数。

1.5.1 自由能 ( $F$ ) 与等温等容系统的判据#

考察一个与恒温热源（温度为 $T$ ）接触的系统。根据克劳修斯不等式，系统在微小过程中吸收的热量 $\delta Q$ 满足：
$\delta Q \le T dS$
结合热力学第一定律 $\delta Q = dU + \delta W$ （其中 $\delta W$ 为系统对外做功），可得：
$dU + \delta W \le T dS$
由于热源温度 $T$ 恒定（等温过程）， $T dS$ 可改写为 $d(TS)$ 。移项整理得到：
$d(U - TS) \le -\delta W$
为此，定义一个新的状态函数——自由能（Helmholtz Free Energy, $F$ ）：
$F \equiv U - TS$
代入上式，得到等温过程的核心不等式：
$dF \le -\delta W$
在等温过程中，系统对外所做的功，永远不可能超过其自由能的减少量。自由能的减少量是在该条件下系统做功的理论上限（可逆过程取等号）。自由能 $F$ 即为内能中能够“自由”转化为对外做功的那部分能量。
等温等容判据：若系统进一步被限制在等容条件下（ $dV = 0$ ），且无其他形式的功，则 $\delta W = 0$ 。上述不等式化简为：
$dF \le 0$
在等温等容条件下，系统内部发生的一切不可逆过程，必然朝着自由能减少的方向进行。当自由能达到极小值时，系统进入稳定的平衡态。

1.5.2 吉布斯函数 ( $G$ ) 与等温等压系统的判据#

实验室中最常见的条件是系统同时与恒温热源（温度 $T$ ）和恒压环境（大气压 $p$ ）接触。在此条件下，系统对外的体积功为 $\delta W = p dV$ 。代入等温过程的不等式 $d(U - TS) \le -\delta W$ ，有：
$d(U - TS) \le -p dV$
由于环境压强 $p$ 恒定， $p dV$ 可改写为 $d(pV)$ 。移项合并所有全微分项，得到：
$d(U - TS + pV) \le 0$
为此，定义另一个新的状态函数——吉布斯函数（Gibbs Free Energy, $G$ ）：
$G \equiv U - TS + pV = H - TS$
代入上式，得到等温等压过程的核心不等式：
$dG \le 0$
等温等压判据：在等温等压条件下，系统自发进行的过程必然导致吉布斯函数减小；当吉布斯函数达到极小值时，系统处于热力学平衡态。

1.5.3 判据的绝对对应关系#

孤立系统（恒 $U, V$ ） $\rightarrow$ 熵 $S$ 极大
等温等容系统（恒 $T, V$ ） $\rightarrow$ 自由能 $F$ 极小
等温等压系统（恒 $T, p$ ） $\rightarrow$ 吉布斯函数 $G$ 极小
切忌在条件不匹配时乱用判据。例如，在绝热膨胀过程中强行使用吉布斯函数极小值作为判据是严重的物理错误。判据的有效性严格绑定于其推导时所引入的外部约束。

第2章：热力学函数与应用#

在第一章中，我们确立了热力学基本方程 $dU = TdS - pdV$ 。但该方程在实际应用中存在一个缺陷：它包含了内能 $U$ 和熵 $S$ 这样无法在实验中直接测量的“绝对量”。为了将理论推导与实验观测（如压强 $p$ 、体积 $V$ 、温度 $T$ 及热容）对接，我们必须引入一套严格的宏观数学转换工具。

2.1 热力学基本函数与特性的数学推演#

2.1.1 勒让德变换与四个基本热力学函数#

热力学基本方程 $dU = TdS - pdV$ 表明，内能 $U$ 的“自然变量”是 $(S, V)$ 。但在实验室中，熵 $S$ 是极其难以控制的变量，我们更习惯控制温度 $T$ 和压强 $p$ 。
为了实现独立变量的转换而不丢失系统的热力学信息，数学上采用勒让德变换（Legendre Transformation）。其核心思想是：若要把函数 $f(x)$ 的自变量转换为其导数 $y = \frac{\partial f}{\partial x}$ ，只需构造新函数 $g = f - xy$ 即可。利用该变换，我们将包含热量交换的项 $TS$ 与包含做功的项 $pV$ 进行组合，严格定义出宏观热力学的四个基本函数（及其全微分式）：
1. 内能 (Internal Energy) $U(S, V)$
  $dU = TdS - pdV$
2. 焓 (Enthalpy) $H(S, p) \equiv U + pV$
  $dH = TdS + Vdp$
3. 自由能 (Helmholtz Free Energy) $F(T, V) \equiv U - TS$
  $dF = -SdT - pdV$
4. 吉布斯函数 (Gibbs Free Energy) $G(T, p) \equiv U - TS + pV$
  $dG = -SdT + Vdp$

2.1.2 麦克斯韦关系#

四个基本函数的全微分不仅确立了自然变量，还直接导出了热力学中最强大的数学等式。对于任意状态函数（全微分） $dz = M dx + N dy$ ，数学上必然满足二阶混合偏导数相等： $\left(\frac{\partial M}{\partial y}\right)_x = \left(\frac{\partial N}{\partial x}\right)_y$ 。将此纯数学定理直接应用于上述四个基本函数的微分表达式，立即得到四大麦克斯韦关系式：
1. 由 $dU$ ： $\left(\frac{\partial T}{\partial V}\right)_S = -\left(\frac{\partial p}{\partial S}\right)_V$
2. 由 $dH$ ： $\left(\frac{\partial T}{\partial p}\right)_S = \left(\frac{\partial V}{\partial S}\right)_p$
3. 由 $dF$ ： $\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V$
4. 由 $dG$ ： $\left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p$
麦克斯韦关系是热力学理论走向实验的桥梁。特别是第3和第4个关系式，它们成功地将包含熵 $S$ 的不可测偏导数，完美地转化为了由 $p, V, T$ 构成的完全可测偏导数。
雅可比行列式恒等式在热力学中，雅可比行列式有一个极为方便的性质，它可以将任何偏导数转化为行列式的商：
$\left( \frac{\partial u}{\partial x} \right)_y = \frac{\partial(u, y)}{\partial(x, y)}$
利用这一工具，我们可以得到一个极其重要的恒等式：
$\frac{\partial(T, S)}{\partial(p, V)} = 1$
证明：
直接取由自由能 $dF$ 导出的麦克斯韦关系：
$\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V$
利用雅可比行列式的链式乘法法则，将等式两边同乘 $\frac{\partial(T, V)}{\partial(p, V)}$ ，即得证：
$\frac{\partial(T, S)}{\partial(p, V)} = 1$
该恒等式表明，从温熵坐标 $(T, S)$ 变换到压容坐标 $(p, V)$ 时，其雅可比行列式（即变换的模）严格等于 1。这不仅意味着在平衡态热力学的可逆过程中相空间的“面积”是守恒的，更证明了 $(T, S)$ 和 $(p, V)$ 作为两对共轭变量，在热力学结构中具有完全对等的物理地位。

2.1.3 特性函数及其应用#

如果把某一个热力学函数表示为特定的一组独立变量的函数，使得不用借助于任何其他的函数（如由实验测定的物态方程或热容量），只要对该函数进行偏导数运算，就可以求得均匀系的全部热力学函数，从而把系统的全部热力学性质完全确定下来，那么这个函数就被称为在该组特定独立变量下的特性函数。常见的热力学特性函数有：
内能的特性函数： $U(S, V)$
焓的特性函数： $H(S, p)$
自由能的特性函数： $F(T, V)$
吉布斯函数的特性函数： $G(T, p)$
1. 以 $(T, V)$ 为独立变量的表述（自由能表示）
在等温过程中，最方便的特性函数是自由能 $F(T, V)$ 。一旦获知 $F(T, V)$ 的解析式，系统的所有性质可按如下逻辑导出：
基本参数：
根据 $dF = -SdT - pdV$ ，直接求一阶偏导：
熵： $S = -\left( \frac{\partial F}{\partial T} \right)_V$
压强： $p = -\left( \frac{\partial F}{\partial V} \right)_T$ （此即系统的状态方程）
能量参数：
内能： $U = F + TS = F - T\left( \frac{\partial F}{\partial T} \right)_V$
焓： $H = U + pV = F - T\left( \frac{\partial F}{\partial T} \right)_V - V\left( \frac{\partial F}{\partial V} \right)_T$
特性参量：
吉布斯函数： $G = F + pV = F - V\left( \frac{\partial F}{\partial V} \right)_T$
定容热容： $C_V = \left( \frac{\partial U}{\partial T} \right)_V = -T\left( \frac{\partial^2 F}{\partial T^2} \right)_V$
1. 以 $(U, V)$ 为独立变量的表述（熵表示）
根据热力学基本方程的另一种形式 $dS = \frac{1}{T}dU + \frac{p}{T}dV$ ，可知在以 $(U, V)$ 为独立变量时，熵 $S(U, V)$ 是该系统的特性函数（属于熵表示，Entropy Representation）。
已知 $S(U, V)$ 后，推导逻辑如下：
基本参数：
温度：根据 $\frac{1}{T} = \left( \frac{\partial S}{\partial U} \right)_V$ ，可得 $T = \left[ \left( \frac{\partial S}{\partial U} \right)_V \right]^{-1}$
压强：根据 $\frac{p}{T} = \left( \frac{\partial S}{\partial V} \right)_U$ ，可得 $p = T \left( \frac{\partial S}{\partial V} \right)_U = \frac{(\partial S / \partial V)_U}{(\partial S / \partial U)_V}$
其他热力学函数：
利用已求得的 $T$ 和 $p$ ，可以直接代入定义式求出其他态函数：
$F = U - TS$
$H = U + pV$
$G = U - TS + pV$
特性函数的“全能性”严格依赖于变量的选择。例如，虽然 $U$ 是状态函数，但如果将其表达为 $U(T, V)$ ，它就失去了特性函数的地位。因为从 $U(T, V)$ 出发，我们只能得到热容信息，却无法唯一地推导出系统的状态方程（压强 $p$ ）。
假设已知系统的内能函数为 $U(T, V)$ ，试着去求该系统的状态方程 $p(T, V)$ 。
借助能态方程（由麦克斯韦关系导出）：
$\left(\frac{\partial U}{\partial V}\right)_T = T\left(\frac{\partial p}{\partial T}\right)_V - p$
将上式整理为关于未知量 $p$ 的偏微分方程：
$\left(\frac{\partial p}{\partial T}\right)_V - \frac{1}{T}p = \frac{1}{T}\left(\frac{\partial U}{\partial V}\right)_T$
等式两边同乘积分因子 $\frac{1}{T^2}$ ，左边凑成全微分形式：
$\left[ \frac{\partial}{\partial T} \left( \frac{p}{T} \right) \right]_V = \frac{1}{T^2}\left(\frac{\partial U}{\partial V}\right)_T$
在保持体积 $V$ 不变的条件下，对温度 $T$ 进行偏积分：
$\frac{p}{T} = \int \frac{1}{T^2}\left(\frac{\partial U}{\partial V}\right)_T dT + f(V)$
得到压强的最终表达式：
$p(T, V) = T \int \frac{1}{T^2}\left(\frac{\partial U}{\partial V}\right)_T dT + T \cdot f(V)$
由于我们进行的是偏积分，产生的积分常数不再是一个具体的数值，而是一个完全未知的、关于体积 $V$ 的任意函数 $f(V)$ 。这表明，即使我们完全掌握了 $U(T, V)$ 的解析式，也无法唯一确定压强 $p$ 的表达式（即无法给出明确的物态方程）。因此， $U(T, V)$ 无法仅通过偏导数就给出系统的全部信息，它失去了特性函数的地位。

2.1.4 广义热容关系#

广义功与变量的物理对应
为了展现热力学框架的极致普适性，我们推导任意共轭参量 $(x, y)$ 下的热容关系。在普遍的热力学系统中，系统在微小过程中的“做功”项（ $\delta W$ ）总可以抽象为两个共轭变量的乘积：
$\delta W = y dx$
$x$ （广义位移）：通常代表系统的某种广延量（如体积 $V$ 、总电偶极矩 $p_e$ 、总磁矩 $m$ ）。
$y$ （广义力）：通常代表促使系统状态改变的强度量（如压强 $p$ 、外加电场 $E$ 、外加磁场 $H$ ）。
(注：以最常见的流体膨胀功 $p dV$ 为例，我们通常直接映射 $x \leftrightarrow V$ ， $y \leftrightarrow p$ )。
基于此，设系统保持广义位移 $x$ 恒定时的热容为 $C_x$ ，保持广义力 $y$ 恒定时的热容为 $C_y$ 。其热力学严格定义为：
$C_x = T\left(\frac{\partial S}{\partial T}\right)_x, \quad C_y = T\left(\frac{\partial S}{\partial T}\right)_y$
1. 广义热容之差 ( $C_y - C_x$ )
选取温度 $T$ 和广义位移 $x$ 作为独立变量，将熵写为 $S(T, x)$ ，取全微分：
$dS = \left(\frac{\partial S}{\partial T}\right)_x dT + \left(\frac{\partial S}{\partial x}\right)_T dx$
等式两边同除以 $dT$ ，并在保持广义力 $y$ 恒定的条件下取偏导数：
$\left(\frac{\partial S}{\partial T}\right)_y = \left(\frac{\partial S}{\partial T}\right)_x + \left(\frac{\partial S}{\partial x}\right)_T \left(\frac{\partial x}{\partial T}\right)_y$
等式两边同乘温度 $T$ ，代入热容定义，直接得到广义热容差的普遍公式：
$C_y - C_x = T \left(\frac{\partial S}{\partial x}\right)_T \left(\frac{\partial x}{\partial T}\right)_y$
将 $V, p$ 代入上式，并调用麦克斯韦关系 $\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V$ ，公式立刻退化为我们熟知的形式：
$C_p - C_V = T \left(\frac{\partial p}{\partial T}\right)_V \left(\frac{\partial V}{\partial T}\right)_p = \frac{T V \alpha^2}{\kappa_T}$
$C_y$ 之所以不等于 $C_x$ ，是因为在等广义力 $y$ （如等压）升温时，系统必须改变其广义位移 $x$ （如体积膨胀）。这一变化不仅涉及克服内部相互作用势能的改变，还需要对外输出广义功（ $\delta W = y dx$ ），从而必须吸收额外的热量。
2. 广义热容之比 ( $C_y / C_x$ )
利用微积分中的循环关系定理 $\left(\frac{\partial A}{\partial B}\right)_C \left(\frac{\partial B}{\partial C}\right)_A \left(\frac{\partial C}{\partial A}\right)_B = -1$ ，将偏导数展开：
对于 $(S, T, y)$ 系统： $\left(\frac{\partial S}{\partial T}\right)_y = - \left(\frac{\partial S}{\partial y}\right)_T \left(\frac{\partial y}{\partial T}\right)_S$
对于 $(S, T, x)$ 系统： $\left(\frac{\partial S}{\partial T}\right)_x = - \left(\frac{\partial S}{\partial x}\right)_T \left(\frac{\partial x}{\partial T}\right)_S$
将两者相比，并利用链式法则化简 $\frac{(\partial S / \partial y)_T}{(\partial S / \partial x)_T} = \left(\frac{\partial x}{\partial y}\right)_T$ 以及 $\frac{(\partial y / \partial T)_S}{(\partial x / \partial T)_S} = \left(\frac{\partial y}{\partial x}\right)_S$ ：
$\frac{C_y}{C_x} = \frac{(\partial S / \partial T)_y}{(\partial S / \partial T)_x} = \left(\frac{\partial x}{\partial y}\right)_T \left(\frac{\partial y}{\partial x}\right)_S$
整理得到广义热容比的普遍公式：
$\frac{C_y}{C_x} = \frac{\left(\frac{\partial y}{\partial x}\right)_S}{\left(\frac{\partial y}{\partial x}\right)_T}$
将 $V, p$ 代入上式：
$\frac{C_p}{C_V} = \frac{\left(\frac{\partial p}{\partial V}\right)_S}{\left(\frac{\partial p}{\partial V}\right)_T}$
由于等温压缩系数 $\kappa_T = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_T$ ，绝热压缩系数 $\kappa_S = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_S$ ，该比值直接映射为力学参量之比：
$\frac{C_p}{C_V} = \frac{\kappa_T}{\kappa_S}$
热学响应之比（热容比）与力学响应之比（压缩系数比）存在着完美的对应关系。由于必然有 $C_p > C_V$ ，因此在相同的体积压缩量下，绝热压缩引起的压强攀升一定比等温压缩更剧烈（ $\kappa_T > \kappa_S$ ）。

2.2 获取低温的物理过程与原理#

2.2.1 节流膨胀与焦耳-汤姆孙效应（等焓制冷）#

气体在持续的压力差下，通过多孔塞（或节流阀）从高压侧 ( $p_1$ ) 缓慢渗漏到低压侧 ( $p_2$ ) 的过程，称为节流膨胀。假设多孔塞及管道完全绝热（ $Q = 0$ ）。
分析特定质量气体在通过节流阀前后的能量转化：
外界对气体做推挤功： $W_1 = p_1 V_1$
气体膨胀对外界做功： $W_2 = -p_2 V_2$
根据第一定律 $\Delta U = W_1 + W_2$ ，有 $U_2 - U_1 = p_1 V_1 - p_2 V_2$ 。
移项整理得： $U_1 + p_1 V_1 = U_2 + p_2 V_2$ 。
这证明了节流膨胀前后的初末态满足 $H_1 = H_2$ 。在热力学处理中，该过程被抽象为一个等焓过程。
为了定量评估节流后的温度变化，定义焦耳-汤姆孙系数 $\mu_{JT} = \left(\frac{\partial T}{\partial p}\right)_H$ 。我们利用上一节的雅可比行列式性质进行极简推导：
$\mu_{JT} = \left(\frac{\partial T}{\partial p}\right)_H = \frac{\partial(T, H)}{\partial(p, H)} = \frac{\frac{\partial(T, H)}{\partial(T, p)}}{\frac{\partial(p, H)}{\partial(T, p)}} = -\frac{\left(\frac{\partial H}{\partial p}\right)_T}{\left(\frac{\partial H}{\partial T}\right)_p}$
分母即为定压热容 $C_p$ 。对于分子，将基本方程 $dH = TdS + Vdp$ 在等温条件下对压强求导：
$\left(\frac{\partial H}{\partial p}\right)_T = T\left(\frac{\partial S}{\partial p}\right)_T + V$
调用由吉布斯函数 $dG$ 导出的麦克斯韦关系 $\left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p$ ：
$\left(\frac{\partial H}{\partial p}\right)_T = -T\left(\frac{\partial V}{\partial T}\right)_p + V$
代回原式，得到焦耳-汤姆孙系数的普适解析式：
$\mu_{JT} = \frac{1}{C_p} \left[ T \left(\frac{\partial V}{\partial T}\right)_p - V \right] = \frac{V}{C_p} (T\alpha - 1)$
（其中 $\alpha = \frac{1}{V}\left(\frac{\partial V}{\partial T}\right)_p$ 为体胀系数）。
1. 理想气体无效应：对于理想气体 $V = nRT/p$ ，必然有 $T\alpha = 1$ ，导致 $\mu_{JT} = 0$ 。理想气体内能仅依赖于温度，等焓膨胀既不改变内能，也不改变温度。
2. 真实气体的制冷与转化：真实气体的 $\mu_{JT}$ 受分子间相互作用主导。
  当 $\mu_{JT} > 0$ 时，压强下降 ( $dp<0$ ) 必然导致温度下降 ( $dT<0$ )，实现制冷。（微观机制：分子间距变大，克服分子间吸引力做功消耗了分子动能）。
  当 $\mu_{JT} < 0$ 时，压强下降导致温度上升，即致热。（微观机制：气压极高时斥力占优，膨胀时斥力做正功转化为动能）。
  令 $\mu_{JT} = 0$ 得到的 $T-p$ 关系即为转化曲线。在工程液化中，必须先将气体预冷至转化曲线内部，方可利用节流阀实现持续降温。

2.2.2 绝热可逆膨胀（等熵制冷）#

若气体在绝热气缸中缓慢推动无摩擦活塞对外做功，过程严格可逆且绝热（ $\delta Q_{\text{rev}} = 0$ ），这是一个等熵过程。
定义等熵膨胀降温系数 $\mu_S = \left(\frac{\partial T}{\partial p}\right)_S$ 。同样利用雅可比行列式：
$\mu_S = \frac{\partial(T, S)}{\partial(p, S)} = \frac{\frac{\partial(T, S)}{\partial(T, p)}}{\frac{\partial(p, S)}{\partial(T, p)}} = -\frac{\left(\frac{\partial S}{\partial p}\right)_T}{\left(\frac{\partial S}{\partial T}\right)_p}$
分母为 $\frac{C_p}{T}$ 。分子再次调用麦克斯韦关系 $\left(\frac{\partial S}{\partial p}\right)_T = -\left(\frac{\partial V}{\partial T}\right)_p = -V\alpha$ ：
$\mu_S = \frac{T V \alpha}{C_p}$
比较两种制冷机制的系数差值：
$\mu_S - \mu_{JT} = \frac{TV\alpha}{C_p} - \frac{V}{C_p}(T\alpha - 1) = \frac{V}{C_p} > 0$
由于 $\mu_S > \mu_{JT}$ ，在经历相同的压强降幅（ $\Delta p < 0$ ）时，绝热可逆膨胀的降温幅度永远大于节流膨胀。从能量守恒的角度看这是必然的：节流膨胀气体仅克服内聚力做功；而绝热膨胀气体不仅要克服内聚力，还要向外输出宏观机械功，导致其内能消耗更为剧烈。

2.2.3 磁制冷：顺磁体的绝热去磁#

气体液化技术在 1K 左右遭遇物理极限（氦气的蒸气压变得极低）。要进入 mK（毫开尔文）级别的极低温领域，体积膨胀已失效，物理学转而利用固体内部的电子磁矩系统。
设系统总磁矩为 $m$ ，外加磁场强度为 $H$ 。磁场对顺磁体系统做功为 $\delta W = H dm$ （注：当系统磁矩增加时，外界磁场对其做正功）。结合第一定律 $\delta Q = dU - \delta W$ ，热力学基本方程改写为：
$dU = TdS + Hdm$
为了利用麦克斯韦关系，定义磁介质的“吉布斯函数” $G = U - TS - Hm$ 。取全微分并代入基本方程：
$dG = -SdT - mdH$
由此立刻得到对应的麦克斯韦关系： $\left(\frac{\partial S}{\partial H}\right)_T = \left(\frac{\partial m}{\partial T}\right)_H$ 。
- 绝热去磁的核心是考察等熵条件下减弱磁场引起的温度变化，即求 $\left(\frac{\partial T}{\partial H}\right)_S$ 。
利用雅可比行列式：
$\left(\frac{\partial T}{\partial H}\right)_S = -\frac{\left(\frac{\partial S}{\partial H}\right)_T}{\left(\frac{\partial S}{\partial T}\right)_H} = -\frac{T}{C_H} \left(\frac{\partial m}{\partial T}\right)_H$
根据顺磁体的居里定律，在恒定磁场下，温度升高导致热运动加剧，磁化强度必然减弱，即 $\left(\frac{\partial m}{\partial T}\right)_H < 0$ 。
代入上式，得出绝对结论： $\left(\frac{\partial T}{\partial H}\right)_S > 0$ 。绝热状态下，减小磁场（ $dH < 0$ ），必然导致温度下降（ $dT < 0$ ）。

2.3 热力学第三定律与负温度#

2.3.1 能斯特定理与绝对零度不可达#

能斯特在总结大量低温化学反应和相变实验后发现，随着温度趋近于绝对零度，系统不仅熵的绝对值在减小，等温过程中的熵变（ $\Delta S$ ）也严格趋于零。其热力学数学表述为：对于任何处于内部平衡态的凝聚系，当温度趋于绝对零度时，系统经历任何可逆等温过程的熵变趋于零。即：
$\lim_{T \to 0} \Delta S = 0$
这就是热力学第三定律的能斯特表述。普朗克随后将其进一步强化为绝对熵的标定：当 $T \to 0$ 时，所有凝聚系的绝对熵趋于零（ $\lim_{T \to 0} S = 0$ ）。在微观统计物理中，这对应着系统基态的非简并性（即微观状态数 $\Omega = 1$ ， $S = k_B \ln \Omega = 0$ ）。
能斯特定理在逻辑上直接宣判了绝对零度是宏观热力学过程的禁区。
回顾上一节的磁制冷循环：制冷依赖于“高场等温磁化（熵减小）”和“绝热去磁（等熵降温）”的交替。如果能斯特定理成立，当系统逼近 0 K 时，任意两个不同外场参数（如 $H_1$ 和 $H_2$ ）下的熵-温度曲线（ $S-T$ 曲线）在 $T=0$ 处必然收敛于同一点（因为等温熵变 $\Delta S \to 0$ ）。无论在何种高场下开始绝热退磁（沿水平等熵线向左移动），由于 $S-T$ 曲线在底部必定汇聚并变平，等熵线永远无法在有限次循环内与 $T=0$ 的纵轴相交。
这就引出了第三定律的等价表述：不可能通过有限的手续将任何系统冷却到绝对零度。
第三定律不仅限制了制冷机的下限，还给所有物质的低温热力学响应系数下了数学铁律。结合第二章的麦克斯韦关系，当 $T \to 0$ 时：
1. 热容趋于零：由于热容定义为 $C_x = T(\partial S / \partial T)_x$ 。在 $T \to 0$ 时，熵的绝对值趋于常数（0），其对温度的导数不发散，因此必然有 $\lim_{T \to 0} C_V = 0$ 且 $\lim_{T \to 0} C_p = 0$ 。这意味着在极低温下，极其微小的热量扰动就能引起巨大的温度剧变。
2. 体胀系数趋于零：由麦克斯韦关系 $(\partial S/\partial p)_T = -V\alpha$ 。因为 $T \to 0$ 时等温熵变 $\Delta S \to 0$ ，即熵对压强完全不敏感（ $(\partial S/\partial p)_T \to 0$ ），必然有 $\lim_{T \to 0} \alpha = 0$ 。物质在绝对零度附近将完全丧失热胀冷缩的力学性质。

2.3.2 负温度的统计物理本质#

在热力学中，温度的最严密定义是 $\frac{1}{T} = \left(\frac{\partial S}{\partial U}\right)_{V, N}$ 。要使温度 $T$ 出现负值，也就是说，随着系统总能量 $E$ 的增加，系统可占据的微观状态数必须减少。负温度系统必须满足一个严苛的先决条件：系统的总能量必须有严格的上限，否则能量越高系统越混乱，无法达到负温度状态。
对于处于热力学平衡的系统，处于能量为 $\varepsilon_i$ 的微观状态上的粒子数 $n_i$ ，服从玻尔兹曼分布律：
$n_i \propto \exp\left(-\frac{\varepsilon_i}{k_B T}\right)$
或者考察高能级（ $\varepsilon_2$ ）与低能级（ $\varepsilon_1$ ）的粒子数之比：
$\frac{n_2}{n_1} = \exp\left(-\frac{\varepsilon_2 - \varepsilon_1}{k_B T}\right)$
这个简单的指数公式，决定了不同温度下达到平衡时系统的总内能（ $U = \sum n_i \varepsilon_i$ ）大小。
设想一个总能量受限的双能级自旋系统（低能级 $\varepsilon_1$ ，高能级 $\varepsilon_2$ ）：
1. 常规状态（正温度 $T > 0$ 且能量较低）：
  由于 $T > 0$ ，指数项必定小于 1。因此 $n_2 < n_1$ 。高能级上的粒子数永远少于低能级。系统极其“守规矩”，绝大多数粒子沉寂在低能态。此时总内能 $U$ 处于较低水平。如果系统继续吸收热量，粒子逐渐向高能级跃迁，混乱度（熵 $S$ ）增加，表现为 $T$ 逐渐升高。
2. 极限混乱状态（ $T \to \pm\infty$ 且能量居中）：
  当系统吸收了足够多的能量，使得 $n_2 = n_1$ （高低能级粒子数完全相等）时，系统达到了绝对的最混乱状态，熵 $S$ 达到抛物线的顶点。
  代入玻尔兹曼公式： $\frac{n_2}{n_1} = 1 \implies \exp\left(-\frac{\Delta \varepsilon}{k_B T}\right) = 1$ 。这要求分母 $T \to \pm\infty$ 。
3. 粒子数反转与负温度（ $T < 0$ 且能量极高）：
  如果系统继续被强行注入能量，高能级上的粒子数将超过低能级，即 $n_2 > n_1$ 。这在物理上被称为粒子数反转（Population Inversion）（这也是激光产生的原理）。
  回到玻尔兹曼公式，由于 $\frac{n_2}{n_1} > 1$ ，指数必然大于 0。由于 $\Delta \varepsilon$ 和 $k_B$ 都是正数，唯一的数学可能是：温度 $T$ 必须变为负值！
  此时，随着能量 $U$ 继续增加（更多粒子跑到高能级），系统反而变得越来越有序，混乱度 $S$ 开始下降。这就完美呼应了宏观的热力学定义 $\frac{1}{T} = \frac{\partial S}{\partial U} < 0$ 。
通过玻尔兹曼分布，我们可以极其直观地看到：
1. 任何正温度系统，无论温度多高，其高能级的粒子数总是小于低能级粒子数。
2. 任何负温度系统，其微观前提就是粒子数反转，即高能级粒子数大于低能级粒子数。
由于内能是所有粒子能量的总和（ $U = \sum n_i \varepsilon_i$ ），一个大量粒子占据高能级的系统，其总内能绝对、严格地大于任何由低能级主导的正常系统。
既然负温度系统蕴含极高的内能，当它与普通正温度系统接触时，负温度系统必然充当“高温热源”向外释放热量（这也满足第二定律使联合系统总熵增加的要求）。因此，负温度并非比绝对零度更冷，而是跨越了正无穷后，比任何正温度都更热的极端状态。真实的宏观温度高低演化排序（从极冷到极热、内能从低到高）应严格表述为：
$+0 \text{ K} \rightarrow +300 \text{ K} \rightarrow +\infty \text{ K} \equiv -\infty \text{ K} \rightarrow -300 \text{ K} \rightarrow -0 \text{ K}$

第3章相变#

在实际的热力学过程中，物质经常以多个相（如固、液、气相）共存的形式出现，并在特定条件下发生相与相之间的转变。本章将基于热力学基本方程和平衡判据，讨论复相平衡的普遍规律及相变的宏观物理特征。由于相变通常在等温等压条件下进行，吉布斯函数极小值判据是分析相变问题的基础。

3.1 热力学平衡与稳定性判据#

在讨论相变与复相平衡之前，需确立系统处于平衡状态的热力学判据。这涉及热力学势的一阶变分（平衡条件）和二阶变分（稳定性条件）。特别需要指出的是，当选取不同类型的热力学参量（广延量与强度量）作为独立变量时，系统稳定性所对应的数学特性存在显著差异。
考虑置于等温等压环境中的封闭系统，其状态由温度 $T$ 、压强 $p$ 以及系统内部参量 $\xi$ （例如某相的物质的量、化学反应进度或序参量）共同描述。根据热力学第二定律，系统在此约束下演化时，吉布斯函数 $G(T, p; \xi)$ 趋于减小。
1. 平衡条件（一阶变分）：系统达到平衡时，吉布斯函数对内部参量的一阶偏导数为零：
$\left(\frac{\partial G}{\partial \xi}\right)_{T,p} = 0$
1. 稳定条件（二阶变分）：仅一阶变分为零不足以保证平衡态的稳定性。为了使系统能够抵抗微小涨落，平衡态对应的必须是局部极小值。因此， $G$ 作为内部参量 $\xi$ 的函数，在平衡点处应满足凸性条件：
$\left(\frac{\partial^2 G}{\partial \xi^2}\right)_{T,p} > 0$
当系统处于平衡态时，内部参量 $\xi$ 不再是独立变量，而是环境参量（如 $T, p$ ）的函数。此时，热力学势仅依赖于其各自的自然变量。系统稳定性的微观涨落要求，在不同自然变量空间中表现为不同符号的二阶偏导数。
1. 广延量空间的凸性
对于以内能 $U(S, V)$ 为代表的势函数，其自然变量均为广延量。热力学稳定性要求总内能对任何局部的熵或体积涨落具有恢复力，这在数学上等价于要求 $U$ 是其自变量 $S$ 和 $V$ 的凸函数。其二阶变分要求 $d^2U = dTdS - dpdV > 0$ 。由此导出系统的基本响应函数必须为正：
热容要求： $\left(\frac{\partial^2 U}{\partial S^2}\right)_V = \frac{T}{C_V} > 0 \implies C_V > 0$
压缩系数要求： $\left(\frac{\partial^2 U}{\partial V^2}\right)_S = \frac{1}{V\kappa_S} > 0 \implies \kappa_S > 0$
1. 为考察热学稳定性，施加等容约束（ $dV = 0$ ），不等式简化为：
$dTdS > 0$
由等容热容定义 $\delta Q_V = TdS = C_V dT$ ，可导出熵变 $dS = \frac{C_V}{T} dT$ 。代入上式有：
$\frac{C_V}{T} (dT)^2 > 0$
由于绝对温度 $T > 0$ 且涨落项 $(dT)^2 > 0$ 恒成立，为使该不等式严格大于零，必须得出：
$C_V > 0$
1. 为考察纯力学稳定性，施加等温约束（ $dT = 0$ ），不等式简化为：
$-dpdV > 0$
由等温压缩系数定义 $\kappa_T = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_T$ ，可导出体积变分 $dV = -V\kappa_T dp$ 。代入上式有：
$V\kappa_T (dp)^2 > 0$
由于系统宏观体积 $V > 0$ 且压强涨落项 $(dp)^2 > 0$ 恒成立，为使该不等式严格大于零，必然得出：
$\kappa_T > 0$
1. 强度量空间的凹性
当通过勒让德变换引入强度量作为独立变量时，热力学势的凹凸性会发生相应改变。以吉布斯函数 $G(T, p) = U - TS + pV$ 为例，对其自然变量 $T$ 和 $p$ 求二阶偏导数：
对温度的二阶偏导：
$\left(\frac{\partial^2 G}{\partial T^2}\right)_p = -\left(\frac{\partial S}{\partial T}\right)_p = -\frac{C_p}{T}$
由于热力学稳定性要求定压热容 $C_p > 0$ ，因此必然得到 $\left(\frac{\partial^2 G}{\partial T^2}\right)_p < 0$ 。
对压强的二阶偏导：
$\left(\frac{\partial^2 G}{\partial p^2}\right)_T = \left(\frac{\partial V}{\partial p}\right)_T = -V\kappa_T$
同理，由于等温压缩系数 $\kappa_T > 0$ ，必然得到 $\left(\frac{\partial^2 G}{\partial p^2}\right)_T < 0$ 。
系统热力学稳定性所施加的物理约束表明：平衡态的热力学势对其广延参量呈现凸性（二阶导数为正），而对其强度参量呈现凹性（二阶导数为负）。因此，尽管系统在趋向平衡时吉布斯函数随内部参量 $\xi$ 取极小值，但平衡态的吉布斯函数 $G_{eq}(T, p)$ 作为一个多变量曲面，在 $(T, p)$ 坐标系中是一个凹函数。这是勒让德变换的必然代数结果。

3.2 复相平衡条件与相图#

设一封闭系统内存在同一物质的两个相（记为相 1 和相 2），系统处于恒温（ $T$ ）和恒压（ $p$ ）环境中。根据热力学第二定律，在等温等压条件下，系统内部发生的任何自发过程必然导致总吉布斯函数 $G$ 减小；当系统达到热力学平衡时，总吉布斯函数取极小值，即：
$dG = 0$
系统的总吉布斯函数为各相吉布斯函数之和（ $G = G_1 + G_2$ ）。由于相变过程中伴随着物质在两相间的转移，需要引入广义热力学基本方程：
$dG = -S dT + V dp + \mu_1 dn_1 + \mu_2 dn_2$
其中， $n_i$ 为相 $i$ 的物质的量， $\mu_i = \left(\frac{\partial G_i}{\partial n_i}\right)_{T,p}$ 为相 $i$ 的化学势（对于纯物质，即为摩尔吉布斯函数）。
在等温等压条件下（ $dT = 0, dp = 0$ ），系统的状态演化需满足：
$dG = \mu_1 dn_1 + \mu_2 dn_2 \le 0$
由于封闭系统总物质的量守恒，必有 $dn_1 = -dn_2$ 。代入上式得到：
$(\mu_2 - \mu_1) dn_2 \le 0$
根据该不等式，可得出以下结论：
1. 自发过程的方向：若 $\mu_1 > \mu_2$ ，则必定有 $dn_2 > 0$ 。即物质总是自发地从化学势较高的相转移到化学势较低的相。
2. 平衡条件：当相变停止、两相达到共存的平衡态时， $dG = 0$ 。由于 $dn_2$ 是相互独立的微小变分，要使 $(\mu_2 - \mu_1) dn_2 = 0$ 恒成立，必然要求括号内为零。
综合热学和力学条件，两相平衡必须同时满足以下三个等式：
1. 热平衡条件： $T_1 = T_2$
2. 力学平衡条件： $p_1 = p_2$ （忽略气液等界面的表面张力影响）
3. 相变平衡条件： $\mu_1(T, p) = \mu_2(T, p)$
相变平衡条件 $\mu_1(T, p) = \mu_2(T, p)$ 是决定系统状态的核心方程。
两相共存线：对于确定的物质， $\mu_1$ 和 $\mu_2$ 是 $T$ 和 $p$ 的已知函数。方程 $\mu_1(T, p) = \mu_2(T, p)$ 在 $p-T$ 平面上确定了一条曲线。这条曲线即为两相共存线（如汽化线、熔解线）。在曲线上的任意状态点，两相可以平衡共存；在曲线两侧的区域，化学势较低的单一相成为唯一的稳定相。
三相点：若系统中有三个相平衡共存，必须同时满足：
$\mu_1(T, p) = \mu_2(T, p) = \mu_3(T, p)$
这构成了一个包含两个独立方程的方程组，对应两个未知数 $(T, p)$ 。一般情况下，该方程组只有唯一解。因此，在 $p-T$ 相图上，三条两相共存线必然交汇于一个确定的点，即三相点。在三相点处，系统的温度和压强被完全确定，系统自由度为零。

3.3 一阶相变与克拉珀龙方程（Clapeyron Equation）#

在相平衡点（两相共存线上），相 1 和相 2 的化学势必须严格相等： $\mu_1(T, p) = \mu_2(T, p)$ 。即化学势在相变过程中是连续的。
但是，化学势的一阶偏导数通常并不连续。根据热力学基本方程 $d\mu = -s dT + v dp$ （其中 $s$ 和 $v$ 分别为摩尔熵和摩尔体积），我们有：
$\left(\frac{\partial \mu}{\partial T}\right)_p = -s, \quad \left(\frac{\partial \mu}{\partial p}\right)_T = v$
如果在相变发生时，系统状态跨越共存线，虽然 $\mu$ 连续，但其一阶偏导数发生跃变，即：
$s_1 \neq s_2, \quad v_1 \neq v_2$
这种相变被称为一阶相变。此时有：
摩尔熵的突变（ $\Delta s = s_2 - s_1 \neq 0$ ）意味着相变过程中系统需要吸收或放出热量。在等温等压下，这部分热量被称为摩尔相变潜热 $L = T(s_2 - s_1)$ 。
摩尔体积的突变（ $\Delta v = v_2 - v_1 \neq 0$ ）意味着相变伴随着宏观的体积膨胀或收缩。
固-液熔化、液-气汽化、固-气升华以及不同固体结晶形态之间的转变，均属于典型的一阶相变。
在 $p-T$ 相图上，两相共存线不仅是一条几何曲线，它的斜率 $\frac{dp}{dT}$ 蕴含着深层的热力学必然性。假设系统沿两相共存线发生状态改变，两相必须分别保持化学势平衡：
$\mu_1(T, p) = \mu_2(T, p)$
在两相共存线上取微分有：
$d\mu_1 = d\mu_2$
代入化学势的全微分表达式 $d\mu = -s dT + v dp$ ：
$-s_1 dT + v_1 dp = -s_2 dT + v_2 dp$
重新排列各项，分离出 $p$ 和 $T$ 的微分：
$(v_2 - v_1) dp = (s_2 - s_1) dT$ $\frac{dp}{dT} = \frac{s_2 - s_1}{v_2 - v_1}$
利用相变潜热定义 $L = T(s_2 - s_1)$ 替换掉极难直接测量的摩尔熵差，我们得到克拉珀龙方程：
$\frac{dp}{dT} = \frac{L}{T(v_2 - v_1)} = \frac{L}{T \Delta v}$
克拉珀龙方程是一阶相变中最核心的数学方程，它将相图上的几何特征（共存线斜率）与物质的宏观热力学响应（潜热、体积差）进行了完美的绝对绑定。
1. 正常物质的熔化（如二氧化碳、绝大多数金属）：
  熔化过程吸收热量，因此潜热 $L > 0$ ；温度恒定为正 $T > 0$ ；固态熔化为液态时体积膨胀，因此 $\Delta v = v_l - v_s > 0$ 。
  代入克拉珀龙方程可得 $\frac{dp}{dT} > 0$ 。
  推论：正常物质的固-液共存线在相图上具有正斜率。外界压强增大，会导致熔点升高。
2. 水/冰体系的反常膨胀（复冰现象）：
  冰熔化成水同样吸热（ $L > 0$ ），但由于冰的特殊氢键晶体结构，固态的摩尔体积竟然大于液态的摩尔体积（即冰浮在水面上），因此 $\Delta v = v_l - v_s < 0$ 。
  代入方程可得绝对结论： $\frac{dp}{dT} < 0$ 。
  推论：冰的固-液共存线斜率为负。这意味着外界施加的压强越大，冰的熔点反而越低。这正是溜冰鞋能够滑行的热力学基础（冰刀下方的极高局部压强使得冰在冰点以下局部熔化为水，形成润滑层）。

3.4 范德瓦尔斯等温线与麦克斯韦等面积法则#

实际的气液一阶相变在等温条件下是一个等压过程（ $p-v$ 图上为水平线段）。范德瓦尔斯方程引入了分子体积和分子间吸引力，试图给出气液两相的统一描述，但在临界温度以下其解析形式与物理稳定性条件存在冲突。
范德瓦尔斯方程为：
$\left(p + \frac{a}{v^2}\right)(v - b) = RT$
将压强 $p$ 表示为摩尔体积 $v$ 的函数，在 $p-v$ 图上绘制等温线。
1. 当 $T > T_c$ 时， $\left(\frac{\partial p}{\partial v}\right)_T < 0$ ，等温线单调下降，系统仅存在单一气相。
2. 当 $T < T_c$ 时，等温线呈现非单调的极值结构。在极小值和极大值之间的区段，满足：
  $\left(\frac{\partial p}{\partial v}\right)_T > 0$
根据热力学稳定性判据，系统的等温压缩系数 $\kappa_T = -\frac{1}{v}\left(\frac{\partial v}{\partial p}\right)_T$ 。在上述区段内 $\kappa_T < 0$ ，系统处于力学不稳定状态，微小的密度涨落将导致系统自发坍缩或膨胀。因此，该波浪线区段在自然界中不能稳定存在。
为了修正非物理的振荡，确定两相共存的实际压强 $p_0$ ，必须引入复相平衡的化学势判据。
设饱和液相状态为 $A(v_l, p_0)$ ，饱和气相状态为 $C(v_g, p_0)$ 。两相平衡共存的必要条件是化学势相等：
$\mu(T, p_0)_A = \mu(T, p_0)_C$
利用化学势的全微分 $d\mu = -s dT + v dp$ 。在等温条件下（ $dT = 0$ ），沿理论上的范德瓦尔斯等温线（包括不稳定的振荡段）从态 $A$ 积分至态 $C$ ：
$\mu_C - \mu_A = \int_A^C d\mu = \int_A^C v dp$
代入平衡条件 $\mu_C = \mu_A$ ，得出：
$\int_A^C v dp = 0$
几何意义：在 $p-v$ 图上，实际相变的水平直线 $p = p_0$ 必须切割范德瓦尔斯等温线，使得水平线与理论曲线之间围成的上方区域面积（液相侧）与下方区域面积（气相侧）严格相等。此即为确定相变压强的麦克斯韦等面积法则。
随着温度升高，等温线上的极大值与极小值相互逼近。当 $T = T_c$ 时，两极值点重合，波浪结构消失，液相和气相的摩尔体积差 $\Delta v \to 0$ 。
临界点是临界等温线上的拐点，且该点切线水平。其数学条件为压强对体积的一阶和二阶偏导数同时为零：
$\left(\frac{\partial p}{\partial v}\right)_{T=T_c} = 0, \quad \left(\frac{\partial^2 p}{\partial v^2}\right)_{T=T_c} = 0$
将范德瓦尔斯方程代入上述条件，解得临界常数：
$v_c = 3b, \quad p_c = \frac{a}{27b^2}, \quad T_c = \frac{8a}{27Rb}$
在临界点处，由于 $\left(\frac{\partial p}{\partial v}\right)_T = 0$ ，等温压缩系数发散： $\kappa_T \to \infty$ 。这意味着流体在临界点的微小压强扰动会导致巨大的密度涨落，一阶相变的宏观边界在此终结。

3.5 连续相变与临界现象#

3.5.1 埃伦费斯特相变分类与二阶相变方程#

为了将不伴随潜热和体积突变的相变纳入严格的热力学框架，埃伦费斯特提出按化学势（或摩尔吉布斯函数 $g$ ）偏导数出现不连续性的最低阶次来对相变进行分类。
设系统在温度 $T$ 和压强 $p$ 下发生相变（由相 1 转变为相 2），两相平衡的核心条件始终是摩尔吉布斯函数连续：
$g_1(T, p) = g_2(T, p)$
1. 一阶相变： $g$ 连续，但 $g$ 的一阶偏导数不连续。
  由于 $\left(\frac{\partial g}{\partial T}\right)_p = -s$ ， $\left(\frac{\partial g}{\partial p}\right)_T = v$ 。
  一阶偏导数跃变意味着存在摩尔熵差 $\Delta s \neq 0$ （即潜热 $L = T\Delta s \neq 0$ ）和摩尔体积差 $\Delta v \neq 0$ 。
2. 二阶相变： $g$ 连续，且 $g$ 的一阶偏导数也连续，但 $g$ 的二阶偏导数发生有限的跃变。
  一阶导数连续意味着：
  $s_1 = s_2 \implies \Delta s = 0 \quad (\text{无相变潜热})$ $v_1 = v_2 \implies \Delta v = 0 \quad (\text{无体积跃变})$
  二阶导数跃变对应着系统的三个宏观响应函数发生阶跃：
  定压摩尔热容： $c_p = -T\left(\frac{\partial^2 g}{\partial T^2}\right)_p \implies \Delta c_p \neq 0$
  等温压缩系数： $\kappa_T = -\frac{1}{v}\left(\frac{\partial^2 g}{\partial p^2}\right)_T \implies \Delta \kappa_T \neq 0$
  体胀系数： $\alpha = \frac{1}{v}\frac{\partial^2 g}{\partial T \partial p} \implies \Delta \alpha \neq 0$
对于一阶相变，我们利用 $g_1 = g_2$ 推导出了决定相界斜率的克拉珀龙方程。
对于二阶相变，由于一阶偏导数在相界线上处处相等，我们可以利用 $s_1(T, p) = s_2(T, p)$ 和 $v_1(T, p) = v_2(T, p)$ 建立新的相界斜率方程。
设系统沿二阶相变共存线发生微小位移 $(dT, dp)$ ，两相的熵和体积必须保持同步变化：
$ds_1 = ds_2 \implies d(\Delta s) = 0$ $dv_1 = dv_2 \implies d(\Delta v) = 0$
1. 由熵的连续性推导第一方程：
  将 $\Delta s$ 视为 $(T, p)$ 的函数，取全微分：
  $d(\Delta s) = \left(\frac{\partial \Delta s}{\partial T}\right)_p dT + \left(\frac{\partial \Delta s}{\partial p}\right)_T dp = 0$
  利用热力学关系 $\left(\frac{\partial s}{\partial T}\right)_p = \frac{c_p}{T}$ ，以及麦克斯韦关系 $\left(\frac{\partial s}{\partial p}\right)_T = -\left(\frac{\partial v}{\partial T}\right)_p = -v\alpha$ 。代入上式：
  $\frac{\Delta c_p}{T} dT - v \Delta \alpha dp = 0$
  移项整理，得到埃伦费斯特第一方程：
  $\frac{dp}{dT} = \frac{\Delta c_p}{T v \Delta \alpha}$
2. 由体积的连续性推导第二方程：
  同理，对 $\Delta v$ 取全微分：
  $d(\Delta v) = \left(\frac{\partial \Delta v}{\partial T}\right)_p dT + \left(\frac{\partial \Delta v}{\partial p}\right)_T dp = 0$
  利用定义 $\left(\frac{\partial v}{\partial T}\right)_p = v\alpha$ 和 $\left(\frac{\partial v}{\partial p}\right)_T = -v\kappa_T$ 。代入上式：
  $v \Delta \alpha dT - v \Delta \kappa_T dp = 0$
  移项整理，得到埃伦费斯特第二方程：
  $\frac{dp}{dT} = \frac{\Delta \alpha}{\Delta \kappa_T}$
埃伦费斯特方程（ $\frac{dp}{dT} = \frac{\Delta c_p}{T v \Delta \alpha} = \frac{\Delta \alpha}{\Delta \kappa_T}$ ）在数学上结构严密，它预言了在二阶相变点，热容和压缩系数只是发生一个有限高度的阶跃。
然而，当物理学家对超导转变、液氦的正常流体-超流体转变（ $\lambda$ 相变）以及铁磁体的居里点进行高精度热力学测量时，发现了致命的矛盾：在逼近真实相变温度 $T_c$ 时， $c_p$ 、 $\alpha$ 和 $\kappa_T$ 并非发生有限的阶跃，而是趋于无穷大（发散）。在数学上，如果 $\Delta c_p \to \infty$ 且 $\Delta \alpha \to \infty$ ，埃伦费斯特方程将变成 $\frac{\infty}{\infty}$ 的不定型，完全失去了描述相界行为的能力。这一物理危机表明：埃伦费斯特的纯宏观泰勒展开式分类法，在面临临界点附近巨大的微观涨落时彻底失效。物理学必须抛弃“有限阶跃”的图像，转而寻找能够处理“无穷大发散”的新理论。

3.5.2 朗道二级相变理论#

为了从微观对称性机制上解释二阶相变，朗道提出了一套基于自由能极值与稳定性分析的唯象理论。相变过程伴随着物质宏观对称性的改变。高温相通常对称性较高（无序态），低温相通常对称性较低（有序态）。引入序参量 $\eta$ 定量描述系统偏离高对称态的程度。规定在高温无序相中 $\eta = 0$ ，在低温有序相中 $\eta \neq 0$ 。
在相变点 $T_c$ 附近， $\eta$ 是一个极小量。假设系统的吉布斯函数 $G(T, p, \eta)$ 在 $T_c$ 附近解析，可按 $\eta$ 的幂级数展开。由于系统在高对称相具有反演对称性（ $\eta$ 与 $-\eta$ 物理状态等价），展开式中奇数次幂项系数必为零：
$G(T, p, \eta) = G_0(T, p) + A(T, p)\eta^2 + B(T, p)\eta^4 + \cdots$
为了保证系统在任意温度下处于热力学稳定状态（即 $\eta \to \infty$ 时自由能 $G$ 不至于发散至负无穷而导致系统崩溃），必须有最高次项系数 $B(T, p) > 0$ 。
系统在热力学平衡时，自由能必须取极小值。求一阶和二阶偏导数：
平衡条件： $\frac{\partial G}{\partial \eta} = 2\eta (A + 2B\eta^2) = 0$
稳定条件： $\frac{\partial^2 G}{\partial \eta^2} = 2A + 12B\eta^2 > 0$
考察高对称相（ $\eta = 0$ ）的稳定性。将 $\eta = 0$ 代入二阶导数得： $\left.\frac{\partial^2 G}{\partial \eta^2}\right|_{\eta=0} = 2A(T, p)$ 。
1. 在高温区 ( $T > T_c$ )：实验事实表明，系统平衡于高对称相。即方程必须仅有唯一的实数解 $\eta = 0$ 。
  将 $\eta = 0$ 代入稳定条件，要求该极小值点满足 $\left.\frac{\partial^2 G}{\partial \eta^2}\right|_{\eta=0} = 2A > 0$ 。因此严格要求： $A(T, p) > 0$ 。
2. 在低温区 ( $T < T_c$ )：系统发生对称性破缺，平衡于低对称相。即平衡方程除了 $\eta=0$ 之外，必须存在非零的实数解 $\eta \neq 0$ 。
  由平衡条件得出 $A + 2B\eta^2 = 0$ ，解得：
  $\eta^2 = -\frac{A}{2B}$
  由于物理上要求序参量 $\eta$ 必须为实数（即 $\eta^2 > 0$ ），且前提条件已定 $B > 0$ ，因此要使该非零解存在，代数上严格要求：
  (注：此时将非零解代入稳定条件，得到二阶导数为 $2A + 12B(-\frac{A}{2B}) = -4A$ 。既然 $A < 0$ ，则 $-4A > 0$ ，证明该非零解确实是一个稳定的极小值点。)
3. 在临界点 ( $T = T_c$ )：作为高温相与低温相的边界，系数 $A(T, p)$ 必须连续地从正值过渡到负值。
  因此，发生二级相变的绝对必要条件为： $A(T_c, p) = 0$ 。
既然确定了 $A(T_c, p) = 0$ ，在非常逼近相变点的温度区间内，将系数 $A(T, p)$ 围绕 $T_c$ 进行泰勒展开，保留最低阶非零项（线性项）：
$A(T, p) \approx \left( \frac{\partial A}{\partial T} \right)_{T_c} (T - T_c)$
令 $a = \left( \frac{\partial A}{\partial T} \right)_{T_c}$ ，即得到：
$A(T, p) = a(T - T_c)$
同时，当 $T > T_c$ 时 $A > 0$ ，必然要求系数 $a > 0$ 。
将 $A = a(T - T_c)$ 代回平衡条件 $2\eta (A + 2B\eta^2) = 0$ 进行求解：
1. 高温相 ( $T \ge T_c$ )： $A \ge 0$ ，方程仅有唯一的实数解 $\eta = 0$ 。
2. 低温相 ( $T < T_c$ )： $A < 0$ ，除了不稳定的 $\eta=0$ 外，方程出现两个对称的非零实数解（代表极小值点）：
$\eta = \pm \sqrt{\frac{-A}{2B}} = \pm \sqrt{\frac{a(T_c - T)}{2B}}$
这在数学上完美预言了序参量随温度下降呈现指数为 $\beta = 1/2$ 的连续增长。
将低温相的非零解 $\eta^2 = \frac{a(T_c - T)}{2B}$ 代回自由能表达式：
$G(T < T_c) = G_0 - \frac{a^2(T_c - T)^2}{4B}$
求熵 $s = -\left(\frac{\partial G}{\partial T}\right)_p$ ：
$s(T < T_c) = s_0 - \frac{a^2(T_c - T)}{2B}$
在 $T = T_c$ 时， $\Delta s = 0$ ，符合二级相变无潜热特征。
求定压热容 $c_p = T\left(\frac{\partial s}{\partial T}\right)_p$ ：
$c_p(T < T_c) = c_{p0} + \frac{a^2 T}{2B}$
而高温相 $c_p(T > T_c) = c_{p0}$ 。因此在相变点，热容发生大小为 $\Delta c_p = \frac{a^2 T_c}{2B}$ 的有限跃变。朗道理论在平均场近似下，从微观对称性出发，极其自洽地推导出了埃伦费斯特二阶相变的全部宏观特征。

3.5.3 涨落的影响与临界指数#

朗道理论本质上是一种“平均场近似”，它假设系统内的序参量是均匀分布的，忽略了热力学涨落（即局部状态偏离宏观平均值的现象）。
通常情况下涨落的宏观影响可以忽略。但当系统极其逼近相变点 $T_c$ 时，决定系统恢复力的系数 $A(T) \to 0$ 。此时，微小的局部热涨落会被无限放大并达到宏观尺度。
这种剧烈的涨落彻底破坏了平均场假设。因此，朗道理论预测的热容等参量的“有限跃变”，在真实的实验中变为了趋于无穷大的“发散”。
朗道理论的泰勒展开在临界区失效，物理学转而采用幂律关系来描述发散行为。
定义无量纲的约化温度 $t = \frac{T - T_c}{T_c}$ 。在临界区（ $|t| \ll 1$ ），系统的宏观热力学量服从如下渐近关系：
序参量 $\eta$ ： $\eta \propto (-t)^\beta \quad (T < T_c)$
定压热容 $c_p$ ： $c_p \propto |t|^{-\alpha}$
等温压缩系数 $\kappa_T$ ： $\kappa_T \propto |t|^{-\gamma}$
这些指数（ $\alpha, \beta, \gamma$ ）即为临界指数。真实系统的实验测定值偏离了朗道理论的预言（例如实验中 $\alpha > 0$ 表现为热容发散，而朗道预言 $\alpha = 0$ 即有限跃变）。
实验表明，临界指数具有深刻的普适性：具有完全不同微观分子间作用力的系统（如气液系统与铁磁自旋系统），只要它们的空间维数和序参量对称性相同，其测得的临界指数就完全一致。
此外，各临界指数之间并非相互独立，而是满足严格的热力学代数约束（称为标度律），例如：
$\alpha + 2\beta + \gamma = 2$
这表明在连续相变的临界点附近，主导系统宏观奇异行为的不再是微观相互作用的细节，而是系统的全局对称结构。

第4章多元系复相平衡与化学平衡#

在前三章中，我们处理的主要是处于封闭状态的单组分系统。本章将热力学框架扩展到包含多个组分，且组分的物质的量可以发生变化（由于相变或化学反应）的系统。

4.1 多元系的热力学基本方程与吉布斯-杜亥姆方程#

设系统由 $k$ 种不同的化学组分构成，各组分的物质的量分别为 $n_1, n_2, \dots, n_k$ 。
当系统的组分发生微小变化时，由于引入或移出粒子会改变系统的总能量，热力学基本方程必须引入与粒子数相关的化学功项 $\sum_{i=1}^k \mu_i dn_i$ 。
四个核心热力学函数的全微分方程推广为：
$dU = TdS - pdV + \sum_{i=1}^k \mu_i dn_i$ $dH = TdS + Vdp + \sum_{i=1}^k \mu_i dn_i$ $dF = -SdT - pdV + \sum_{i=1}^k \mu_i dn_i$ $dG = -SdT + Vdp + \sum_{i=1}^k \mu_i dn_i$
由上述方程可得出组分 $i$ 的化学势 $\mu_i$ 的严密数学定义，最常用的是基于吉布斯函数的定义：
$\mu_i = \left(\frac{\partial G}{\partial n_i}\right)_{T, p, n_{j \neq i}}$
物理意义：在等温、等压且其他组分物质的量不变的条件下，增加一摩尔组分 $i$ 所引起的系统总吉布斯函数的增量。
温度 $T$ 、压强 $p$ 是强度量，而吉布斯函数 $G$ 和物质的量 $n_i$ 是广延量。
在等温等压下，如果将系统中所有组分的物质的量同时放大 $\lambda$ 倍，系统的总吉布斯函数也必将放大 $\lambda$ 倍。这在数学上表明 $G$ 是关于各组分 $n_i$ 的一阶齐次函数：
$G(T, p, \lambda n_1, \dots, \lambda n_k) = \lambda G(T, p, n_1, \dots, n_k)$
根据数学上的欧拉齐次函数定理，一阶齐次函数满足：
$G = \sum_{i=1}^k n_i \left(\frac{\partial G}{\partial n_i}\right)_{T,p,n_{j \neq i}}$
代入化学势的定义，得到多元系吉布斯函数的绝对积分表达式：
$G = \sum_{i=1}^k n_i \mu_i$
这表明，多元系的吉布斯函数严格等于各组分化学势与其物质的量的加权和。
对上述积分形式 $G = \sum n_i \mu_i$ 取全微分：
$dG = \sum_{i=1}^k \mu_i dn_i + \sum_{i=1}^k n_i d\mu_i$
将其与广义热力学基本方程 $dG = -SdT + Vdp + \sum \mu_i dn_i$ 联立，两式相减，消去 $dG$ 和 $\sum \mu_i dn_i$ 项，得到：
$\sum_{i=1}^k n_i d\mu_i = -SdT + Vdp$
此即为吉布斯-杜亥姆方程。
推论与物理意义：
1. 系统的强度状态参量（ $T, p, \mu_1, \dots, \mu_k$ ）并不是完全独立的。它们的变化受到吉布斯-杜亥姆方程的严格约束。
2. 在等温等压过程中（ $dT=0, dp=0$ ），方程退化为：
  $\sum_{i=1}^k n_i d\mu_i = 0$
  或者除以总物质的量 $n = \sum n_i$ ，使用摩尔分数 $x_i = n_i/n$ 表达：
  $\sum_{i=1}^k x_i d\mu_i = 0$
  这表明在恒温恒压下的混合物中，如果某一成分的化学势升高，必然导致其他某些成分的化学势降低，以维持总体的热力学制约关系。
对于单组分理想气体，在等温条件下 $d\mu = v dp = \frac{RT}{p} dp$ 。积分可得单组分气体的化学势：
$\mu(T, p) = \mu^0(T) + RT \ln p$
其中 $\mu^0(T)$ 是 $p=1$ 标准大气压下的标准化学势。
对于多元理想气体混合物，根据道尔顿分压定律，组分 $i$ 在混合物中的行为如同它单独占据整个容器时的状态。因此，组分 $i$ 的化学势仅取决于温度及其分压 $p_i$ ：
$\mu_i(T, p) = \mu_i^0(T) + RT \ln p_i$
又因为分压 $p_i = x_i p$ （ $x_i$ 为该组分的摩尔分数， $p$ 为系统总压），可将其改写为：
$\mu_i(T, p) = \mu_i^0(T) + RT \ln p + RT \ln x_i$
该方程将化学势清晰地分离为了温度依赖项、总压依赖项以及组分浓度（摩尔分数）依赖项。

4.2 独立组分数与吉布斯相律 (Gibbs Phase Rule)#

在明确了多元系的热力学基本方程后，我们需要解决一个宏观热力学中最普遍的拓扑问题：对于一个包含多个组分和多个相的系统，它究竟拥有多少个可以独立改变的强度状态参量（如温度、压强、浓度）？这个问题的最终答案即为吉布斯相律。
设系统由 $c$ 个独立组分构成，共存在 $\phi$ 个相（记为 $\alpha, \beta, \dots, \phi$ ）。当系统达到热力学平衡时，必须同时满足以下三个条件：
1. 热平衡：所有相的温度相等。
$T^{(\alpha)} = T^{(\beta)} = \dots = T^{(\phi)} = T$
1. 力学平衡：忽略表面张力和重力场，所有相的压强相等。
$p^{(\alpha)} = p^{(\beta)} = \dots = p^{(\phi)} = p$
1. 相变与化学平衡：每一种组分在各个相中的化学势必须严格相等。对于任意组分 $i$ ( $i = 1, 2, \dots, c$ )：
$\mu_i^{(\alpha)}(T, p, x_1^{(\alpha)}, \dots) = \mu_i^{(\beta)}(T, p, x_1^{(\beta)}, \dots) = \dots = \mu_i^{(\phi)}(T, p, x_1^{(\phi)}, \dots)$
（其中 $x_i^{(\alpha)}$ 为组分 $i$ 在 $\alpha$ 相中的摩尔分数）。
在实际化学系统中，总物种数（记为 $K$ ）不一定等于独立决定系统组成所需的变量数。如果系统中各物种之间发生化学反应，或存在其他特定的比例约束，组分数会随之减少。
独立组分数 $c$ 的计算公式为：`
$c = K - R - R'$
$K$ ：系统中的总化学物种数。
$R$ ：物种之间能够达到化学平衡的独立化学反应数目。
$R'$ ：同一相中各物种浓度之间的其他独立限制条件（例如电中性条件，或由初始投料比决定的特定浓度比例）。
自由度 ( $f$ ) 定义为在保持系统相数 $\phi$ 不变的前提下，系统可以独立改变的强度参量（ $T, p$ 及各相浓度）的数目。
根据代数方程理论，自由度等于系统中总强度参量的数目减去这些参量必须满足的独立约束方程数目。
1. 总强度参量数目：
  整个系统有一个共同的温度 $T$ 和压强 $p$ （共2个）。
  对于每一个相，其组成由 $c$ 个组分的摩尔分数决定。由于满足归一化条件 $\sum_{i=1}^c x_i = 1$ ，每个相有 $c - 1$ 个独立的摩尔分数。
  因此， $\phi$ 个相的总强度参量数目为： $2 + \phi(c - 1)$ 。
2. 独立约束方程数目：
  复相平衡条件要求每一组分在各相中的化学势相等。对于单一组分 $i$ ，连等式 $\mu_i^{(\alpha)} = \mu_i^{(\beta)} = \dots = \mu_i^{(\phi)}$ 包含了 $\phi - 1$ 个独立的方程式。
  系统共有 $c$ 个独立组分，因此化学势等式构成的总约束方程数目为： $c(\phi - 1)$ 。
3. 计算自由度 $f$ ：
  将上述两项相减：
  $f = [2 + \phi(c - 1)] - c(\phi - 1)$ $f = 2 + \phi c - \phi - c\phi + c$ $f = c - \phi + 2$
这就是极其著名的吉布斯相律 (Gibbs Phase Rule)。式中的“2”代表温度和压强这两个外界物理变量。
由于自由度 $f$ 必须是非负整数（ $f \ge 0$ ），相律对系统可能存在的相数给出了绝对的上限：
$\phi \le c + 2$
单组分系统 ( $c=1$ )：
$f = 3 - \phi$ 。
由于 $f \ge 0$ ，最大相数 $\phi_{\max} = 3$ 。当存在三个相（固、液、气）时， $f = 0$ 。这证明了单组分系统的三相点在 $p-T$ 图上是一个唯一确定的无自由度点，不可能出现四相共存。
二元系统 ( $c=2$ )：
$f = 4 - \phi$ 。
最大相数 $\phi_{\max} = 4$ 。如果保持压强恒定（去除了1个自由度），则条件自由度 $f^* = c - \phi + 1 = 3 - \phi$ 。在恒压下的二元合金相图中，最多只能观察到三相共存（如共晶反应点）。

4.3 化学平衡条件与质量作用定律#

在多元相平衡中，我们假设物质在各相间转移但不发生种类的改变。本节将吉布斯极小值判据应用于发生化学反应的系统，推导出决定反应进行程度的普适法则。
设一封闭系统中发生如下一般形式的化学反应：
$\nu_1 A_1 + \nu_2 A_2 + \dots + \nu_k A_k = 0$
简写为 $\sum_{i=1}^k \nu_i A_i = 0$ 。其中 $A_i$ 代表化学物种， $\nu_i$ 为化学计量数（规定：生成物 $\nu_i > 0$ ，反应物 $\nu_i < 0$ ）。
引入反应进度（Extent of Reaction, $\lambda$ ）。当反应发生极小的推进 $d\lambda$ 时，各组分物质的量的变化必定与其化学计量数成严格的比例关系：
$dn_i = \nu_i d\lambda$
在等温等压条件下，系统自发演化的方向必须满足吉布斯函数减小（ $dG \le 0$ ）。代入多元系全微分方程：
$dG = \sum_{i=1}^k \mu_i dn_i = \left( \sum_{i=1}^k \nu_i \mu_i \right) d\lambda \le 0$
1. 反应自发方向：
  若 $\sum \nu_i \mu_i < 0$ ，必然要求 $d\lambda > 0$ ，反应自发向正方向进行。
  若 $\sum \nu_i \mu_i > 0$ ，必然要求 $d\lambda < 0$ ，反应自发向逆方向进行。
2. 化学平衡的绝对条件：当系统达到彻底的热力学平衡时，总吉布斯函数处于极小值，微小变分 $dG = 0$ 。由于 $d\lambda$ 可正可负，括号内的项必须严格为零：
  $\sum_{i=1}^k \nu_i \mu_i = 0$
这就是推导一切化学平衡规律的根本方程。它表明在平衡态下，生成物化学势的加权和等于反应物化学势的加权和。
将理想气体混合物组分的化学势公式 $\mu_i(T, p) = \mu_i^0(T) + RT \ln p_i$ 代入平衡条件 $\sum \nu_i \mu_i = 0$ 中：
$\sum_{i=1}^k \nu_i \left[ \mu_i^0(T) + RT \ln p_i \right] = 0$
展开并重新组合各项：
$RT \sum_{i=1}^k \nu_i \ln p_i = -\sum_{i=1}^k \nu_i \mu_i^0(T)$
利用对数运算性质 $\sum \nu_i \ln p_i = \ln \left( \prod p_i^{\nu_i} \right)$ ，并定义标准摩尔反应吉布斯函数 $\Delta G^0(T) = \sum \nu_i \mu_i^0(T)$ ，上式变为：
$RT \ln \left( \prod_{i=1}^k p_i^{\nu_i} \right) = -\Delta G^0(T)$
对两边取指数，我们得到了极其重要的质量作用定律（Law of Mass Action）：
$\prod_{i=1}^k p_i^{\nu_i} = K_p(T)$
其中， $K_p(T) = \exp\left(-\frac{\Delta G^0(T)}{RT}\right)$ 定义为压强平衡常数。
物理推论：
1. 质量作用定律表明，在达到化学平衡时，各参与气体分压的幂乘积是一个绝对的常数。
2. 因为 $\Delta G^0(T)$ 仅仅是温度 $T$ 的函数，所以 $K_p$ 仅仅是温度的函数，它与系统的总压强、反应容器体积以及反应物的初始投料比例毫无关系。
在实际工程计算中，我们更关心物质的浓度（摩尔分数 $x_i$ ）。利用分压与总压的关系 $p_i = x_i p$ ，代入质量作用定律：
$\prod_{i=1}^k (x_i p)^{\nu_i} = \left( \prod_{i=1}^k x_i^{\nu_i} \right) p^{\sum \nu_i} = K_p(T)$
定义摩尔分数平衡常数 $K_x = \prod_{i=1}^k x_i^{\nu_i}$ ，反应前后的气体摩尔数变化为 $\Delta \nu = \sum \nu_i$ 。移项可得：
$K_x = K_p(T) \cdot p^{-\Delta \nu}$
物理推论：
1. 虽然 $K_p$ 仅依赖于温度，但摩尔分数平衡常数 $K_x$ 却依赖于系统的总压强 $p$ 。
2. 当反应是一个气体分子数增加的过程（ $\Delta \nu > 0$ ）时，增大总压 $p$ 会导致 $p^{-\Delta \nu}$ 减小，从而使 $K_x$ 减小，意味着平衡向逆方向（体积缩小方向）移动。这就是我们在下一节探讨勒夏特列原理的数学基础。

4.4 反应热与化学平衡的移动#

在前一节中，我们导出了决定化学平衡态的核心常数 $K_p(T)$ 。本节将利用经典热力学基本方程，定量考察外界条件（温度和压强）的变化如何导致平衡态的移动，并严格推导出勒夏特列原理的数学表达。
在等温等压且不作非体积功的条件下，系统发生化学反应所吸收的热量 $Q_p$ 等于系统的焓变 $\Delta H$ 。定义标准摩尔反应焓（反应热）：
$\Delta H^0(T) = \sum_{i=1}^k \nu_i H_i^0(T)$
为了寻找平衡常数 $K_p$ 与温度的关系，我们需要借助吉布斯函数 $G$ 与焓 $H$ 之间的基本联系。由 $G = H - TS$ ，两边同除以 $T$ 并对 $T$ 求恒压偏导：
$\left[ \frac{\partial (G/T)}{\partial T} \right]_p = \frac{1}{T}\left(\frac{\partial G}{\partial T}\right)_p - \frac{G}{T^2}$
代入 $\left(\frac{\partial G}{\partial T}\right)_p = -S$ ，得到：
$\left[ \frac{\partial (G/T)}{\partial T} \right]_p = -\frac{S}{T} - \frac{H - TS}{T^2} = -\frac{H}{T^2}$
这就是著名的吉布斯-亥姆霍兹方程。将其应用于标准态下的化学反应，即有：
$\frac{d}{dT} \left( \frac{\Delta G^0}{T} \right) = -\frac{\Delta H^0}{T^2}$
（注：因为标准态仅为温度的函数，偏微商变为全微商）。
根据质量作用定律，压强平衡常数 $K_p$ 满足：
$\Delta G^0 = -RT \ln K_p \implies \frac{\Delta G^0}{T} = -R \ln K_p$
将其代入上述吉布斯-亥姆霍兹方程中：
$-R \frac{d \ln K_p}{dT} = -\frac{\Delta H^0}{T^2}$
整理得到决定化学平衡随温度移动的根本法则——范霍夫等压方程：
$\frac{d \ln K_p}{dT} = \frac{\Delta H^0}{RT^2}$
勒夏特列原理的热力学证明
(1) 温度对平衡移动的影响
由范霍夫等压方程可知， $K_p$ 随温度的变化率完全由反应热 $\Delta H^0$ 的符号决定（由于 $RT^2 > 0$ 恒成立）：
- 吸热反应 ( $\Delta H^0 > 0$ )：
$\frac{d \ln K_p}{dT} > 0$ 。升高温度 ( $dT > 0$ )，必导致 $\ln K_p$ 增大，平衡常数变大，平衡向**正方向（产物方向）**移动。
- 放热反应 ( $\Delta H^0 < 0$ )：
$\frac{d \ln K_p}{dT} < 0$ 。升高温度 ( $dT > 0$ )，必导致 $\ln K_p$ 减小，平衡常数变小，平衡向**逆方向（反应物方向）**移动。
(2) 压强对平衡移动的影响
在 4.3 节中我们导出了摩尔分数平衡常数 $K_x$ 与总压 $p$ 的关系： $K_x = K_p(T) p^{-\Delta \nu}$ ，其中 $\Delta \nu = \sum \nu_i$ 是反应前后气体分子总数的相对增量。
对上式两边取对数，并在等温条件下对压强 $p$ 求导：
$\ln K_x = \ln K_p(T) - \Delta \nu \ln p$ $\left( \frac{\partial \ln K_x}{\partial p} \right)_T = -\frac{\Delta \nu}{p}$
由于 $p > 0$ 恒成立，摩尔分数组成的移动方向完全由 $\Delta \nu$ 的符号决定：
- 分子数增加的反应 ( $\Delta \nu > 0$ )：
$\left( \frac{\partial \ln K_x}{\partial p} \right)_T < 0$ 。增大系统总压 ( $dp > 0$ )，导致 $K_x$ 减小，平衡向**逆方向（分子数减少的方向）**移动。
- 分子数减少的反应 ( $\Delta \nu < 0$ )：
$\left( \frac{\partial \ln K_x}{\partial p} \right)_T > 0$ 。增大系统总压 ( $dp > 0$ )，导致 $K_x$ 增大，平衡向**正方向（也是分子数减少的方向）**移动。
- 分子数不变的反应 ( $\Delta \nu = 0$ )：
压强变化不影响 $K_x$ ，平衡不发生移动。

第5章近独立子系组成的统计理论#

本章从微观粒子的力学性质出发，引入统计物理的基本假设，建立系统宏观状态与其微观状态之间的定量联系。我们首先对微观态进行量子力学界定，并确立等概率原理作为一切统计推导的逻辑起点。

5.1 统计物理的基本假设与微观态的量子力学描述#

5.1.1 相空间与微观态的量子化描述#

在经典力学中，单粒子的微观运动状态由广义坐标 $\mathbf{q} = (q_1, q_2, q_3)$ 和广义动量 $\mathbf{p} = (p_1, p_2, p_3)$ 共同确定。这 6 个变量构成的 6 维空间称为单粒子相空间（ $\mu$ 空间）。经典粒子的状态对应于 $\mu$ 空间中的一个几何点。
然而，根据量子力学的海森堡不确定性原理，粒子在某一方向上的坐标与动量不能同时被无限精确地确定，其不确定度满足：
$\Delta q_i \Delta p_i \ge \frac{h}{4\pi}$
（在玻尔-索末菲量子化条件或半经典极限下，通常取体积元尺度为普朗克常数 $h$ ）。
因此，单粒子微观态不再是相空间中的一个数学点，而是一个体积为 $h^3$ 的相格（Phase Cell）。在体积元 $d^3\mathbf{q} d^3\mathbf{p}$ 内，系统可能存在的微观量子态数目为：
$d\Omega = \frac{d^3\mathbf{q} d^3\mathbf{p}}{h^3}$
一维相空间量子化体积 $\Delta q \Delta p = h$ 的推导
考虑一个质量为 $m$ 的自由粒子，被限制在长度为 $L$ 的一维无限深势阱中（坐标域 $0 \le q \le L$ ）。
根据量子力学的薛定谔方程，粒子波函数必须在边界处为零（即形成驻波），其德布罗意波长 $\lambda$ 必须满足驻波条件：
$L = n \frac{\lambda_n}{2} \quad (n = 1, 2, 3, \dots)$
结合德布罗意关系 $p = \frac{h}{\lambda}$ ，该粒子的动量本征值（绝对值）被严格量子化为：
$p_n = \frac{n h}{2L}$
现在切换到经典相空间（由坐标 $q$ 和动量 $p$ 构成的二维空间）。
对于处于第 $n$ 个量子态的粒子，其动量大小为 $p_n$ 。在经典图景中，该粒子在势阱中以恒定速率来回碰撞，因此其动量在 $+p_n$ 和 $-p_n$ 之间交替，坐标 $q$ 则在 $0$ 到 $L$ 之间均匀取值。
该粒子在相空间中的运动轨迹是一个面积被封闭的矩形线圈。根据玻尔-索末菲量子化条件或经典作用量积分，第 $n$ 个量子态所包围的相空间总面积 $A_n$ 为：
$A_n = \oint p \, dq = \int_0^L (+p_n) \, dq + \int_L^0 (-p_n) \, dq = 2 p_n L$
将量子化动量 $p_n = \frac{n h}{2L}$ 代入相空间面积公式中：
$A_n = 2 \left( \frac{n h}{2L} \right) L = n h$
这表明，第 $n$ 个量子态所包围的相空间面积严格等于 $n$ 倍的普朗克常数。
进一步计算两个相邻量子态（第 $n$ 态与第 $n-1$ 态）在相空间中所夹的环形面积（即一个微观量子态所独占的相空间面积）：
$\Delta A = A_n - A_{n-1} = n h - (n-1) h = h$
无论势阱的长度 $L$ 如何变化，每一个独立的一维微观量子态在相空间中必然占据绝对恒定的面积 $h$ 。
将此一维模型直接推广到三维独立自由粒子，坐标与动量分量相互正交，每个维度的相面积元均为 $h$ 。因此，三维单粒子相空间（ $\mu$ 空间）中，一个微观量子态所对应的相体积基元严格等于：
$\Delta \Omega = (\Delta q_x \Delta p_x) (\Delta q_y \Delta p_y) (\Delta q_z \Delta p_z) = h^3$
这在数学上赋予了连续的经典相体积积分 $\frac{1}{h^3} \int d^3\mathbf{q} d^3\mathbf{p}$ 以严格的离散量子态计数物理意义。
系统允许的单粒子能量通常是量子化的离散值，记为能级 $\epsilon_1, \epsilon_2, \dots, \epsilon_l, \dots$ 。当存在多个不同的量子态具有相同的能量 $\epsilon_l$ 时，该能级被称为是简并的，其包含的量子态数目定义为简并度 $g_l$ 。

5.1.2 近独立子系模型与全同粒子分类#

实际的宏观系统由大量粒子构成。若系统极度稀薄或粒子间相互作用能远小于粒子的平均动能，系统可以被抽象为近独立子系（Nearly Independent Subsystems）。
在此模型中，考虑孤立系统，系统总能量 $E$ 严格等于各个子系（粒子）的能量之和。若第 $l$ 个能级 $\epsilon_l$ 上分布有 $a_l$ 个粒子，则系统总能量为：
$E = \sum_l a_l \epsilon_l$
同时满足总粒子数守恒：
$N = \sum_l a_l$
根据量子力学的全同性原理，系统中的微观粒子按其不可分辨性及自旋特征严格分为三类：
1. 定域子系（Localized Subsystems）：虽然粒子全同，但由于被固定在晶格等空间特定格点上，粒子可以通过其空间位置被绝对分辨（如爱因斯坦固体模型中的晶格振子）。
2. 非定域玻色子（Non-localized Bosons）：自旋为整数的粒子（如光子、声子、氦-4原子）。粒子在空间中自由运动，绝对不可分辨。其多体波函数是对称的，同一个微观量子态可以容纳任意数量的玻色子。
3. 非定域费米子（Non-localized Fermions）：自旋为半奇数的粒子（如电子、质子）。粒子同样绝对不可分辨。其多体波函数是反对称的，严格服从泡利不相容原理（Pauli Exclusion Principle），即一个微观量子态内最多只能容纳一个费米子（ $a_l \le g_l$ ）。

5.1.3 等概率原理（Principle of Equal A Priori Probabilities）#

统计物理学必须在动力学规律与概率论之间建立公理联系。对于一个与外界无任何能量与物质交换的孤立系统（微正则系综，宏观约束为 $N, V, E$ ），由于不受外界干扰，力学规律并不偏袒任何一个符合宏观守恒条件的微观状态。
由此确立统计物理学最根本的假设——等概率原理：
对于处于平衡态的孤立系统，系统出现在任何一个满足宏观约束条件的微观状态的概率是绝对相等的。
基于等概率原理，系统在宏观上必然表现出拥有微观状态数最多的那个分布状态。这表明，求解宏观热力学平衡态的问题，在数学上已被完全转化为 计算在给定 $(N, V, E)$ 约束下特定分布所包含的微观状态数（热力学概率 $\Omega$ ） 的组合数学问题。

5.2 近独立子系微观状态数 $\Omega$ 的计算#

对于包含 $N$ 个近独立子系的宏观平衡态系统，设其处于特定的能量分布 $\{a_1, a_2, \dots, a_l, \dots\}$ （即第 $l$ 能级 $\epsilon_l$ 上的占据粒子数为 $a_l$ ，对应的能级简并度为 $g_l$ ）。系统对应于该分布的总微观状态数 $\Omega$ 的代数计算，严格依赖于微观粒子的可分辨性及泡利不相容原理的物理约束。

5.2.1 麦克斯韦-玻尔兹曼（MB）统计：定域子系#

定域子系在空间上是绝对可分辨的。计算系统处于给定分布 $\{a_l\}$ 的微观状态数分为两个独立步骤：
1. 粒子分组组合：将 $N$ 个可分辨粒子划分为若干组，使得第 $l$ 能级确切分配到 $a_l$ 个粒子。该过程的组合数可通过连乘形式严格展开：
$C_N^{a_1} C_{N-a_1}^{a_2} C_{N-a_1-a_2}^{a_3} \cdots = \frac{N!}{a_1!(N-a_1)!} \frac{(N-a_1)!}{a_2!(N-a_1-a_2)!} \cdots = \frac{N!}{\prod_l a_l!}$
1. 能级内部态分配：对于分配至第 $l$ 能级的 $a_l$ 个可分辨粒子，每个粒子均可独立占据该能级所属的 $g_l$ 个量子态之一。因此，该能级内部的微观占据方式总数为 $g_l^{a_l}$ 。
各能级的量子态占据过程相互独立，系统的总微观状态数为粒子分组组合数与各能级内部占据方式连乘积的乘积：
$\Omega_{MB} = \left( \frac{N!}{\prod_l a_l!} \right) \prod_l g_l^{a_l} = N! \prod_l \frac{g_l^{a_l}}{a_l!}$

5.2.2 费米-狄拉克（FD）统计：非定域费米子#

费米子为全同不可分辨粒子，严格服从泡利不相容原理，每个量子态至多容纳一个粒子，即必须满足数学约束 $a_l \le g_l$ 。
1. 单能级状态数：对于第 $l$ 能级，由于粒子不可分辨，微观状态的区别仅在于 $g_l$ 个量子态中有哪 $a_l$ 个量子态被占据。互换处于不同量子态的粒子不产生新的微观态。因此，该能级的微观状态数等于从 $g_l$ 中选取 $a_l$ 的组合数：
$\Omega_{l, FD} = C_{g_l}^{a_l} = \frac{g_l!}{a_l!(g_l - a_l)!}$
1. 总热力学概率：由于各能级量子态的占据相互独立，整个系统的总微观状态数为各单能级微观状态数的连乘积：
$\Omega_{FD} = \prod_l \Omega_{l, FD} = \prod_l C_{g_l}^{a_l} = \prod_l \frac{g_l!}{a_l!(g_l - a_l)!}$

5.2.3 玻色-爱因斯坦（BE）统计：非定域玻色子#

玻色子为全同不可分辨粒子，且量子态的占据数不受任何限制。
1. 单能级状态数：对于第 $l$ 能级，将 $a_l$ 个不可分辨粒子分配至 $g_l$ 个量子态的数学模型，等价于将 $a_l$ 个全同元素与 $g_l - 1$ 个全同隔板进行一维线形全排列。总排列位置数为 $a_l + g_l - 1$ ，从中选取 $a_l$ 个位置放置粒子的组合数为：
$\Omega_{l, BE} = C_{a_l + g_l - 1}^{a_l} = C_{a_l + g_l - 1}^{g_l - 1} = \frac{(a_l + g_l - 1)!}{a_l!(g_l - 1)!}$
1. 总热力学概率：各能级的状态分配相互独立，系统对应的总微观状态数为各单能级微观状态数的连乘积：
$\Omega_{BE} = \prod_l \Omega_{l, BE} = \prod_l C_{a_l + g_l - 1}^{a_l} = \prod_l \frac{(a_l + g_l - 1)!}{a_l!(g_l - 1)!}$

5.3 极值条件与三种最概然分布的导出#

在热力学极限（ $N \to \infty, V \to \infty$ ）下，系统微观状态数 $\Omega$ 在某一个特定的分布 $\{a_l^*\}$ 处会表现出极其尖锐的极大值。该分布所包含的微观状态数占据了系统所有可能状态数的绝对主导地位，因此被称为最概然分布。求系统平衡态的统计规律，在数学上即转化为在宏观约束条件下求解 $\ln \Omega$ 的极大值问题。

5.3.1 斯特林公式与连续化近似#

为了应用微积分求解泛函极值，必须将离散的阶乘连续化。当系统中各能级上的粒子数 $a_l \gg 1$ 且简并度 $g_l \gg 1$ 时，可引入斯特林公式（Stirling’s approximation）：
$\ln(x!) \approx x \ln x - x$
基于此，我们分别对三种统计体系的 $\ln \Omega$ 进行严格展开与一阶变分求导。

5.3.2 三种统计系统的泛函一阶变分#

麦克斯韦-玻尔兹曼 (MB) 系统：
利用斯特林公式展开 $\ln \Omega_{MB}$ ：
$\ln \Omega_{MB} = \ln N! + \sum_l (a_l \ln g_l - \ln a_l!) \approx N \ln N - N + \sum_l \left[ a_l \ln g_l - (a_l \ln a_l - a_l) \right]$
对粒子数分布 $a_l$ 发生微小变分 $\delta a_l$ 时，求一阶变分（由于 $N$ 守恒， $\delta N = 0$ ）：
$\delta \ln \Omega_{MB} = \sum_l \left[ \ln g_l - \frac{\partial}{\partial a_l}(a_l \ln a_l - a_l) \right] \delta a_l = \sum_l (\ln g_l - \ln a_l - 1 + 1) \delta a_l = \sum_l \ln \left( \frac{g_l}{a_l} \right) \delta a_l$
费米-狄拉克 (FD) 系统：
$\ln \Omega_{FD} = \sum_l \left[ \ln g_l! - \ln a_l! - \ln(g_l - a_l)! \right]$ $\approx \sum_l \left[ g_l \ln g_l - g_l - (a_l \ln a_l - a_l) - ((g_l - a_l) \ln(g_l - a_l) - (g_l - a_l)) \right]$
化简后对其求一阶变分（注意仅 $a_l$ 为变量， $g_l$ 为常数）：
$\delta \ln \Omega_{FD} = \sum_l \left[ - (\ln a_l + 1) - \frac{\partial}{\partial a_l} ((g_l - a_l) \ln(g_l - a_l)) \right] \delta a_l$ $= \sum_l \left[ - \ln a_l - 1 - ( - \ln(g_l - a_l) - 1 ) \right] \delta a_l = \sum_l \ln \left( \frac{g_l - a_l}{a_l} \right) \delta a_l = \sum_l \ln \left( \frac{g_l}{a_l} - 1 \right) \delta a_l$
玻色-爱因斯坦 (BE) 系统：
由于 $a_l \gg 1, g_l \gg 1$ ，可近似取 $a_l + g_l - 1 \approx a_l + g_l$ 以及 $g_l - 1 \approx g_l$ ：
$\ln \Omega_{BE} \approx \sum_l \left[ (a_l + g_l) \ln(a_l + g_l) - (a_l + g_l) - (a_l \ln a_l - a_l) - (g_l \ln g_l - g_l) \right]$
对其求一阶变分：
$\delta \ln \Omega_{BE} = \sum_l \left[ \ln(a_l + g_l) + 1 - (\ln a_l + 1) \right] \delta a_l = \sum_l \ln \left( \frac{a_l + g_l}{a_l} \right) \delta a_l = \sum_l \ln \left( \frac{g_l}{a_l} + 1 \right) \delta a_l$

5.3.3 统一变分表达式与拉格朗日乘子法#

观察上述三个变分结果，我们可以通过引入一个符号参量 $\delta$ ，将其严格统一为唯一的代数形式：
$\delta \ln \Omega = \sum_l \ln \left( \frac{g_l}{a_l} - \delta \right) \delta a_l$
其中：
$\delta = 0$ 对应 MB 统计
$\delta = 1$ 对应 FD 统计
$\delta = -1$ 对应 BE 统计
对于处于热力学平衡的孤立系统，必须同时满足两个绝对的宏观守恒约束：
1. 总粒子数守恒：
  $\sum_l \delta a_l = 0$
2. 总能量守恒：
  $\sum_l \epsilon_l \delta a_l = 0$
为了在受约束条件下求解 $\delta \ln \Omega = 0$ ，应用拉格朗日乘子法（Lagrange multipliers）。分别引入两个待定常数乘子 $-\alpha$ 和 $-\beta$ ，乘以两个守恒约束方程并加到变分方程中，构造无条件极值方程：
$\sum_l \left[ \ln \left( \frac{g_l}{a_l} - \delta \right) - \alpha - \beta \epsilon_l \right] \delta a_l = 0$
由于引入了拉格朗日乘子，上述方程中的各个 $\delta a_l$ 现已在数学上转化为完全独立的变分量。为了保证该连加等式对任意的 $\delta a_l$ 组合均严格恒等于零，方程各项方括号内的系数必须分别恒等于零：
$\ln \left( \frac{g_l}{a_l} - \delta \right) - \alpha - \beta \epsilon_l = 0$ $\frac{g_l}{a_l} - \delta = e^{\alpha + \beta \epsilon_l}$
通过代数移项，我们得到了主宰微观统计体系的最核心、最普适的统一最概然分布律：
$a_l = \frac{g_l}{e^{\alpha + \beta \epsilon_l} + \delta} \quad (\delta = 0, 1, -1)$

5.4 拉格朗日乘子的热力学意义与配分函数#

为了确定最概然分布中引入的拉格朗日乘子 $\alpha$ 和 $\beta$ 的物理意义，我们从宏观系统的内能全微分出发，通过与热力学基本方程的严格比对，自然导出宏观热力学量并证实玻尔兹曼关系。

5.4.1 内能的全微分与热功分解#

系统的内能 $U$ 为处于各能级上的粒子能量总和：
$U = \sum_l a_l \epsilon_l$
考虑系统经历一个极微小的可逆热力学过程，内能的全微分包含两部分：
$dU = \sum_l \epsilon_l da_l + \sum_l a_l d\epsilon_l$
这两项具有明确的宏观热力学对应：
1. 第一项 $\sum \epsilon_l da_l$ ：在外部约束（能级 $\epsilon_l$ ）固定的条件下，粒子在各能级间重新分布引起的能量变化。这完全对应于系统与外界交换的热量，即可逆过程中的微小热量 $\delta Q = \sum \epsilon_l da_l$ 。
2. 第二项 $\sum a_l d\epsilon_l$ ：对于近独立子系，能级 $\epsilon_l$ 的数值仅取决于系统受到的外部约束（如系统的体积 $V$ ）。在粒子分布 $\{a_l\}$ 不变的条件下，改变体积导致能级移动所引起的能力变化，即为系统对外界做功。因此 $- \sum a_l d\epsilon_l = p dV$ 。

5.4.2 利用微观变分条件替换热量项#

当系统处于平衡态时，由于系统绝大部分时间停留在最概然分布，微观状态数 $\Omega$ 可用最概然分布的状态数代替。对 $\ln \Omega$ 求全微分：
$d(\ln \Omega) = \sum_l \frac{\partial \ln \Omega}{\partial a_l} da_l$
将 5.3 节导出的极值条件 $\frac{\partial \ln \Omega}{\partial a_l} = \alpha + \beta \epsilon_l$ 代入上式：
$d(\ln \Omega) = \sum_l (\alpha + \beta \epsilon_l) da_l = \alpha \sum_l da_l + \beta \sum_l \epsilon_l da_l$
由于 $\sum da_l = dN$ ，且 $\sum \epsilon_l da_l = \delta Q$ ，可将热量项解出：
$\delta Q = \frac{1}{\beta} d(\ln \Omega) - \frac{\alpha}{\beta} dN$

5.4.3 与热力学基本方程的比对与玻尔兹曼关系的自然导出#

将 $\delta Q$ 和功的表达式代回内能全微分 $dU$ 中：
$dU = \frac{1}{\beta} d(\ln \Omega) - \frac{\alpha}{\beta} dN - p dV$
将此式与宏观开系的热力学基本方程进行对比：
$dU = T dS + \mu dN - p dV$
通过各项系数的对应关系，我们得到以下严格结论：
推导温度与 $\beta$ 的关系：对于封闭系统（ $dN = 0$ ），对比得出 $T dS = \frac{1}{\beta} d(\ln \Omega)$ ，即：
$dS = \frac{1}{\beta T} d(\ln \Omega)$
因为热力学熵 $S$ 是广延量，而对于两个相互独立的系统，总概率 $\Omega = \Omega_1 \Omega_2$ 使得 $\ln \Omega$ 具有严格的加和性（广延性）。因此，比例系数 $\frac{1}{\beta T}$ 必须是一个与具体系统无关的普适常数，定义为玻尔兹曼常数 $k_B$ 。
由此得出乘子 $\beta$ 的物理意义：
$\beta = \frac{1}{k_B T}$
同时对其积分，自然导出了著名的玻尔兹曼关系（略去积分常数）：
$S = k_B \ln \Omega$
推导化学势与 $\alpha$ 的关系：对比 $dN$ 的系数：
$\mu = -\frac{\alpha}{\beta} = -\alpha k_B T \implies \alpha = -\frac{\mu}{k_B T}$

5.4.4 麦克斯韦-玻尔兹曼分布与单粒子配分函数 $Z_1$ #

将 $\alpha, \beta$ 代回 5.3 节的分布律。对于遵守 MB 统计的系统（ $\delta = 0$ ），其各能级平均粒子数为：
$a_l = g_l e^{\mu/k_B T} e^{-\epsilon_l/k_B T}$
为了消去化学势 $\mu$ ，利用总粒子数归一化条件 $N = \sum a_l = e^{\mu/k_B T} \sum g_l e^{-\epsilon_l/k_B T}$ 。
定义单粒子配分函数（Partition Function） $Z_1$ 为：
$Z_1 = \sum_l g_l e^{-\epsilon_l/k_B T}$
由此得到 $e^{\mu/k_B T} = \frac{N}{Z_1}$ ，将其代回分布公式，得到以配分函数表达的分布律：
$a_l = N \frac{g_l e^{-\epsilon_l/k_B T}}{Z_1}$

5.4.5 宏观热力学量的统计计算框架#

确立了 $Z_1$ 后，系统的宏观热力学量均可通过对 $\ln Z_1$ 及其参量（ $T, V$ ）求偏导数直接导出。
内能 $U$ ：
$U = \sum_l a_l \epsilon_l = \frac{N}{Z_1} \sum_l \epsilon_l g_l e^{-\epsilon_l/k_B T} = -N \frac{\partial \ln Z_1}{\partial \beta} = N k_B T^2 \left( \frac{\partial \ln Z_1}{\partial T} \right)_V$
亥姆霍兹自由能 $F$ ：
由热力学关系 $\mu = \left(\frac{\partial F}{\partial N}\right)_{T,V}$ ，将 $e^{\mu/k_B T} = \frac{N}{Z_1}$ 取对数得 $\mu = -k_B T \ln\left(\frac{Z_1}{N}\right)$ 。对该式积分（对于定域子系），可得：
$F = -N k_B T \ln Z_1$
压强 $p$ （物态方程）：
$p = -\left(\frac{\partial F}{\partial V}\right)_T = N k_B T \left( \frac{\partial \ln Z_1}{\partial V} \right)_T$

5.5 经典极限条件#

在 5.3 节导出的统一分布律中，单粒子能级的平均占据数（分布函数）表示为：
$f_l = \frac{a_l}{g_l} = \frac{1}{e^{\alpha + \beta \epsilon_l} + \delta} \quad (\delta = 0, 1, -1)$
当系统处于特定宏观状态，使得所有允许能级上的平均占据数均远小于 1（即 $f_l \ll 1$ ）时，系统满足经典极限条件，在统计物理中亦称非简并条件（Non-degeneracy Condition）。
由于能级能量 $\epsilon_l \ge 0$ 且 $\beta = \frac{1}{k_B T} > 0$ ，要使得分母中的指数项远大于 1 对应任意能级严格恒成立，其数学充要条件为：
$e^\alpha \gg 1 \quad \text{或等价表示为} \quad -\frac{\mu}{k_B T} \gg 1$
在此极限下，分母中的量子修正项 $\delta = \pm 1$ 的代数权重可以忽略（即量子统计的零阶近似）。费米分布与玻色分布在数学上同时收敛于经典的麦克斯韦-玻尔兹曼分布形式：
$a_l \approx g_l e^{-\alpha} e^{-\beta \epsilon_l}$
利用总粒子数守恒条件 $N = \sum a_l$ ，在经典极限下，总粒子数与单粒子配分函数 $Z_1 = \sum g_l e^{-\beta \epsilon_l}$ 满足：
$N \approx e^{-\alpha} \sum_l g_l e^{-\beta \epsilon_l} = e^{-\alpha} Z_1$
由此解得经典极限下逸度 $z$ 与拉格朗日乘数 $\alpha$ 的一阶表达式：
$z = e^{\mu/k_B T} = e^{-\alpha} = \frac{N}{Z_1}$
将 5.6 节在连续能量近似下导出的单原子理想气体平动配分函数 $Z_1 = V \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2}$ 代入上式，得到以宏观热力学量（温度 $T$ 、体积 $V$ 、粒子数 $N$ ）表征的经典极限显式判据：
$e^\alpha = \frac{V}{N} \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2} \gg 1$
引入分子的平均平均间距 $d = \left(\frac{V}{N}\right)^{1/3} = n^{-1/3}$ ，以及表征微观粒子波动性空间扩展尺度的热德布罗意波长（Thermal De Broglie Wavelength）：
$\lambda_{dB} = \frac{h}{\sqrt{2\pi m k_B T}}$
则上述经典极限判据可严格改写为无量纲的几何比值形式：
$\left( \frac{d}{\lambda_{dB}} \right)^3 \gg 1 \implies d \gg \lambda_{dB}$
该结果表明，经典极限条件的物理实质在于：气体的平均粒子间距远大于其微观热德布罗意波长。此时，粒子波包的空间重叠区域趋于零，由全同性原理引发的量子统计关联（玻色聚集或费米排斥）在宏观统计平均中不再显现。

5.6 经典理想气体的统计分布与双原子分子内部自由度#

在弱简并极限下，系统的统计行为遵循经典的麦克斯韦-玻尔兹曼分布。对于宏观容器中的平动自由粒子，其相邻能级间的能量差极小，离散的求和运算在数学上可严格转化为对连续相空间的积分。对于多原子分子（如双原子分子），除了质心的平动之外，还必须引入分子内部的转动与振动自由度。

5.6.1 连续能量积分与平动配分函数#

对于处于体积 $V$ 内的单原子理想气体，粒子的平动能仅取决于其动量，色散关系为 $\epsilon = \frac{p^2}{2m} = \frac{p_x^2 + p_y^2 + p_z^2}{2m}$ 。
在单粒子相空间（ $\mu$ 空间）中，处于微小相体积元 $d^3\mathbf{q} d^3\mathbf{p}$ 内的微观量子态数目为 $\frac{d^3\mathbf{q} d^3\mathbf{p}}{h^3}$ 。将 5.4 节中配分函数的离散求和 $\sum_l g_l e^{-\beta \epsilon_l}$ 转化为积分：
$Z_t = \iint \frac{1}{h^3} e^{-\beta \frac{p^2}{2m}} d^3\mathbf{q} d^3\mathbf{p} = \frac{1}{h^3} \left( \iiint_V d^3\mathbf{q} \right) \left( \iiint_{-\infty}^{\infty} e^{-\frac{p_x^2 + p_y^2 + p_z^2}{2m k_B T}} dp_x dp_y dp_z \right)$
由于空间积分直接给出系统体积 $V$ ，动量积分可分离为三个独立的一维高斯积分（利用公式 $\int_{-\infty}^{\infty} e^{-ax^2} dx = \sqrt{\frac{\pi}{a}}$ ）：
$Z_t = \frac{V}{h^3} \left( \int_{-\infty}^{\infty} e^{-\frac{p_x^2}{2m k_B T}} dp_x \right) \left( \int_{-\infty}^{\infty} e^{-\frac{p_y^2}{2m k_B T}} dp_y \right) \left( \int_{-\infty}^{\infty} e^{-\frac{p_z^2}{2m k_B T}} dp_z \right) = V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2}$

5.6.2 麦克斯韦速度与速率分布律#

系统处于动量区间 $\mathbf{p}$ 到 $\mathbf{p} + d\mathbf{p}$ （即相体积元 $V d^3\mathbf{p}$ ）内的粒子数 $dN_{\mathbf{p}}$ 可由 MB 分布律直接给出：
$dN_{\mathbf{p}} = N \left( \frac{1}{2\pi m k_B T} \right)^{3/2} e^{-\frac{p_x^2 + p_y^2 + p_z^2}{2m k_B T}} dp_x dp_y dp_z$
引入动量与速度的关系 $\mathbf{p} = m\mathbf{v}$ ，则 $d p_x = m d v_x$ ，代换可得处于速度区间 $\mathbf{v}$ 到 $\mathbf{v}+d\mathbf{v}$ 的粒子数（麦克斯韦速度分布律）：
$dN_{\mathbf{v}} = N \left( \frac{m}{2\pi k_B T} \right)^{3/2} e^{-\frac{m(v_x^2 + v_y^2 + v_z^2)}{2k_B T}} dv_x dv_y dv_z$
将速度空间的直角坐标积分变换为球坐标系，体积元 $dv_x dv_y dv_z$ 对应于球壳体积 $4\pi v^2 dv$ ，得到速率概率密度函数 $f(v) = \frac{1}{N} \frac{dN_v}{dv}$ （麦克斯韦速率分布律）：
$f(v) = 4\pi \left( \frac{m}{2\pi k_B T} \right)^{3/2} v^2 e^{-\frac{mv^2}{2k_B T}}$

5.6.3 多自由度解耦与双原子分子平动#

对于构造更复杂的理想双原子分子系统，分子的总能量不仅包含其质心的平动能 $\epsilon_t$ ，还包含分子绕质心转动的能量 $\epsilon_r$ 以及两原子沿轴线振动的能量 $\epsilon_v$ 。在非极端条件下，这些自由度可视为相互解耦的独立变量。因此，双原子分子的单粒子哈密顿量可表达为加和形式：
$\epsilon = \epsilon_t + \epsilon_r + \epsilon_v$
由于能量表达式满足代数可加性，单粒子配分函数 $Z_1$ 能够严格地分解为各项独立自由度配分函数的乘积：
$Z_1 = \sum e^{-\beta(\epsilon_t + \epsilon_r + \epsilon_v)} = Z_t \cdot Z_v \cdot Z_r$
平动部分的对比：
设双原子分子由质量分别为 $m_1$ 和 $m_2$ 的原子构成，其总质量为 $m = m_1 + m_2$ 。双原子分子的平动完全由其质心运动描述，其质心色散关系同样为 $\epsilon_t = \frac{p^2}{2m}$ 。因此，在连续能量近似下，双原子分子的平动配分函数在数学形式上与单原子分子完全一致：
$Z_t = V \left( \frac{2\pi (m_1+m_2) k_B T}{h^2} \right)^{3/2}$

5.6.4 谐振子模型与振动配分函数#

双原子分子两核间的相对位移在微振动近似下可视为一维理想量子谐振子。设其特征振动频率为 $\nu$ 。
量子谐振子的本征能级为：
$\epsilon_n = \left(n + \frac{1}{2}\right)h\nu \quad (n = 0, 1, 2, \dots)$
由于一维谐振子能级是非简并的，其简并度 $g_n = 1$ 。构建振动配分函数 $Z_v$ 表现为严格的等比级数形式：
$Z_v = \sum_{n=0}^\infty e^{-\beta (n + 1/2)h\nu} = e^{-\frac{h\nu}{2k_B T}} \sum_{n=0}^\infty \left( e^{-\frac{h\nu}{k_B T}} \right)^n$
由于公比 $e^{-\frac{h\nu}{k_B T}} < 1$ ，应用等比级数求和公式 $\sum_{n=0}^\infty q^n = \frac{1}{1-q}$ ，直接导出振动配分函数的精确解析式：
$Z_v = \frac{e^{-\frac{h\nu}{2k_B T}}}{1 - e^{-\frac{h\nu}{k_B T}}}$
引入特征振动温度 $\Theta_v = \frac{h\nu}{k_B}$ 。分析其在不同温度区间下的渐近行为：
1. 低温极限（ $T \ll \Theta_v$ ）：指数项 $e^{-\frac{\Theta_v}{T}} \to 0$ ，分母趋于 1，配分函数收敛于 $Z_v \approx e^{-\frac{h\nu}{2k_B T}}$ 。此时系统完全冻结在零点振动基态，对热容等宏观物理量无微商贡献。
2. 高温经典极限（ $T \gg \Theta_v$ ）：将分母的指数项进行一阶泰勒展开（ $e^{-\frac{h\nu}{k_B T}} \approx 1 - \frac{h\nu}{k_B T}$ ），同时分子近似为 1，振配分函数自动退化为：
$Z_v \approx \frac{1}{1 - (1 - \frac{h\nu}{k_B T})} = \frac{k_B T}{h\nu} = \frac{T}{\Theta_v}$
由此，理想双原子分子的总单粒子配分函数在经典高温极限下可完整表述为：
$Z_1 = \left[ V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2} \right] \cdot \left[ \frac{2I k_B T}{\hbar^2} \right] \cdot \left[ \frac{k_B T}{h\nu} \right]$
对该式应用 5.4 节的热力学偏导数框架，即可严密解析出双原子气体包含全自由度贡献的内能及定容热容。

5.6.5 刚性转子模型与转动配分函数#

双原子分子的转动通常抽象为绕质心旋转的刚性转子模型。设分子的折合质量为 $\mu = \frac{m_1 m_2}{m_1 + m_2}$ ，两核间距为 $R_0$ ，则绕垂直于核轴线的转动惯量为 $I = \mu R_0^2$ 。
根据量子力学，角动量量子化导致转动能级呈现离散分布：
$\epsilon_J = \frac{\hbar^2}{2I} J(J+1) \quad (J = 0, 1, 2, \dots)$
每个能级对应的空间取向简并度为：
$g_J = 2J + 1$
由此构建转动配分函数 $Z_r$ 的级数形式：
$Z_r = \sum_{J=0}^\infty (2J+1) e^{-\frac{\hbar^2 J(J+1)}{2I k_B T}}$
引入特征转动温度 $\Theta_r = \frac{\hbar^2}{2I k_B}$ 。在常规及高温极限下（即 $T \gg \Theta_r$ ），转动能级间隔远小于热运动能量 $k_B T$ ，离散求和可转换为连续积分：
$Z_r \approx \int_0^\infty (2J+1) e^{-\frac{\Theta_r}{T} J(J+1)} dJ$
令变元 $x = J(J+1)$ ，则 $dx = (2J+1)dJ$ ，积分化简为：
$Z_r = \int_0^\infty e^{-\frac{\Theta_r}{T} x} dx = \frac{T}{\Theta_r} = \frac{2I k_B T}{\hbar^2}$
(注：上述推导适用于异核双原子分子。若为全同同核双原子分子，由于空间旋转 $180^\circ$ 后的状态在拓扑上不可分辨，需引入对称性因子 $\sigma = 2$ ，配分函数相应修正为 $Z_r = \frac{T}{2\Theta_r}$ )。

5.7 理想气体熵的吉布斯佯谬与全同性修正#

本节探讨全同粒子的不可分辨性对系统宏观广延性质的决定性影响，并展示如何通过量子修正消除经典统计力学在计算混合熵时出现的非广延性缺陷。

5.7.1 经典理想气体熵与吉布斯佯谬的数学展现#

若不考虑量子全同性，将气体分子视为经典力学中可分辨的独立粒子，则含有 $N$ 个粒子的理想气体系统总配分函数直接写为 $Z_N = (Z_1)^N$ 。将 5.6 节导出的单粒子平动配分函数 $Z_1 = V \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2} = V f(T)$ 代入，系统的亥姆霍兹自由能为：
$F = -k_B T \ln Z_N = -N k_B T \ln [V f(T)]$
根据热力学关系 $S = -\left(\frac{\partial F}{\partial T}\right)_{V,N}$ ，推导出未修正的经典理想气体熵表达式：
$S = N k_B \ln V + N k_B \left[ \ln f(T) + T \frac{f'(T)}{f(T)} \right]$
现考虑吉布斯佯谬的物理模型：在一个由隔板分为两个对称区域（体积均为 $V$ 、各含有 $N$ 个同温同压的同种气体分子）的孤立容器中，抽去隔板发生扩散过程。
混合前：两部分气体的总初熵为各项之和：
$S_i = 2 \cdot \left\{ N k_B \ln V + N k_B \left[ \ln f(T) + T \frac{f'(T)}{f(T)} \right] \right\}$
混合后：若粒子被视为可分辨的，抽去隔板后总粒子数变为 $2N$ ，气体可到达的总体积变为 $2V$ ，系统终熵为：
$S_f = (2N) k_B \ln (2V) + (2N) k_B \left[ \ln f(T) + T \frac{f'(T)}{f(T)} \right]$
计算混合熵变 $\Delta S$ ：
$\Delta S = S_f - S_i = 2N k_B \ln (2V) - 2N k_B \ln V = 2N k_B \ln 2 > 0$
物理佯谬分析：根据热力学理论，相同温度和压强的同种气体相互混合是一个完全不改变系统宏观状态的可逆物理过程，其实际宏观熵变 $\Delta S$ 必须严格为零。经典统计推导出的非零混合熵（ $2N k_B \ln 2$ ）破坏了热力学熵的广延性约束（即体积和粒子数同时加倍时，熵未能表现出严格的加倍特征），此矛盾结果即为吉布斯佯谬（Gibbs Paradox）。

5.7.2 全同性因子的引入与萨库尔-特特罗德方程#

吉布斯佯谬的微观根源在于经典力学允许跟踪粒子轨迹，从而错误地赋予了全同粒子以可分辨性。量子力学全同性原理表明，同种微观粒子在本质上是绝对不可分辨的。因此，当 $N$ 个全同粒子在弱简并条件下分布于不同的量子态时，任何在粒子之间进行的置换操作均不产生新的微观多体物理状态。经典配分函数 $Z_N = (Z_1)^N$ 将系统状态数重复计算了 $N!$ 次。
为了恢复统计理论的自洽性，必须引入全同性因子 $\frac{1}{N!}$ 对系统的总配分函数进行修剪：
$Z_N = \frac{1}{N!} (Z_1)^N$
利用斯特林近似 $\ln N! \approx N \ln N - N$ ，重新计算系统的亥姆霍兹自由能：
$F = -k_B T \ln \left[ \frac{1}{N!} (Z_1)^N \right] \approx -k_B T \left[ N \ln Z_1 - (N \ln N - N) \right] = -N k_B T \left[ \ln \left( \frac{Z_1}{N} \right) + 1 \right]$
将 $Z_1 = V \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2}$ 代入上式：
$F = -N k_B T \left[ \ln \left( \frac{V}{N} \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2} \right) + 1 \right]$
依据 $S = -\left(\frac{\partial F}{\partial T}\right)_{V,N}$ ，对上式关于温度 $T$ 求偏导数（注意 $\frac{V}{N}$ 项不随温度变化）：
$S = N k_B \left[ \ln \left( \frac{V}{N} \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2} \right) + 1 \right] + N k_B T \cdot \frac{\partial}{\partial T} \left[ \frac{3}{2} \ln T \right] = N k_B \left[ \ln \left( \frac{V}{N} \left(\frac{2\pi m k_B T}{h^2}\right)^{3/2} \right) + \frac{5}{2} \right]$
此式即为严格符合热力学广延性的单原子理想气体绝对熵公式——萨库尔-特特罗德方程（Sackur-Tetrode Equation）。
佯谬的代数消除检验：
利用修正后的方程重新审视同种气体的混合过程：
混合前：系统总熵为两部分之和：
$S_i = 2 \cdot N k_B \left[ \ln \left( \frac{V}{N} f(T) \right) + \frac{5}{2} \right]$
混合后：由于对数函数内的自变量表现为特异体积 $\frac{\text{总体积}}{\text{总粒子数}} = \frac{2V}{2N} = \frac{V}{N}$ ，该比值在等密度混合过程中维持严格恒定：
$S_f = (2N) k_B \left[ \ln \left( \frac{2V}{2N} f(T) \right) + \frac{5}{2} \right] = 2N k_B \left[ \ln \left( \frac{V}{N} f(T) \right) + \frac{5}{2} \right]$
最终熵变：
$\Delta S = S_f - S_i = 0$

5.8 弱简并理想玻色气体和费米气体#

在弱简并条件（ $e^\alpha \gg 1$ ）下，系统相空间中的量子态极为充裕，能级平均占据数远小于 $1$ 。但若系统并未绝对逼近经典极限（即处于弱简并与中度简并的过渡区），占据数分布函数中分母上的量子特征项 $\delta = \pm 1$ 将对系统的宏观热力学量产生不可忽略的一阶微扰。

5.8.1 分布函数的一阶渐近展开#

引入参量逸度（Fugacity） $z = e^{-\alpha} = e^{\mu/k_B T}$ 。在弱简并条件下，化学势 $\mu$ 表现为绝对值极大的负数，因此逸度是一个极小量，即 $0 < z \ll 1$ 。
将统一的最概然分布律用逸度 $z$ 表示，并在 $z \ll 1$ 附近将其对 $z e^{-\beta \epsilon_l}$ 进行一阶泰勒展开：
$a_l = \frac{g_l}{z^{-1}e^{\beta \epsilon_l} + \delta} = g_l z e^{-\beta \epsilon_l} \frac{1}{1 + \delta z e^{-\beta \epsilon_l}} \approx g_l z e^{-\beta \epsilon_l} \left( 1 - \delta z e^{-\beta \epsilon_l} \right)$
展开式中的第一项即为经典的麦克斯韦-玻尔兹曼分布，而第二项则代表了由于全同粒子统计特性（ $\delta = -1$ 对应玻色子， $\delta = +1$ 对应费米子）引入的量子修正一阶项。

5.8.2 粒子数与内能的连续积分及一阶修正展开#

利用三维自由粒子态密度 $D(\epsilon) = \frac{2\pi V}{h^3} (2m)^{3/2} \epsilon^{1/2}$ ，将离散求和转化为连续积分，计算包含一阶修正的总粒子数 $N$ ：
$N = \int_0^\infty a(\epsilon) D(\epsilon) d\epsilon \approx \frac{2\pi V}{h^3} (2m)^{3/2} \int_0^\infty \epsilon^{1/2} \left[ z e^{-\beta \epsilon} - \delta z^2 e^{-2\beta \epsilon} \right] d\epsilon$
利用伽马积分 $\int_0^\infty x^{1/2} e^{-ax} dx = \frac{\sqrt{\pi}}{2} a^{-3/2}$ ，计算可得：
$N \approx \frac{2\pi V}{h^3} (2m)^{3/2} \left[ z \frac{\sqrt{\pi}}{2} \beta^{-3/2} - \delta z^2 \frac{\sqrt{\pi}}{2} (2\beta)^{-3/2} \right]$
引入系统的热德布罗意波长 $\lambda = \frac{h}{\sqrt{2\pi m k_B T}}$ ，化简上式得：
$N = \frac{V}{\lambda^3} \left( z - \frac{\delta}{2^{3/2}} z^2 \right) \quad \implies \quad n \lambda^3 = z - \frac{\delta}{2^{3/2}} z^2$
（此处 $n = N/V$ 为粒子数密度）。由于在弱简并下 $z \ll 1$ ，可利用迭代法将逸度 $z$ 用宏观可测的参量 $n$ 和 $T$ 反向表达。零阶近似为 $z \approx n \lambda^3$ ；将其代入一阶项中，得到反解关系：
$z \approx n \lambda^3 \left( 1 + \frac{\delta}{2^{3/2}} n \lambda^3 \right)$
同理，计算系统的总内能 $U$ 。内能积分核包含 $\epsilon \cdot \epsilon^{1/2} = \epsilon^{3/2}$ ：
$U = \int_0^\infty \epsilon a(\epsilon) D(\epsilon) d\epsilon \approx \frac{2\pi V}{h^3} (2m)^{3/2} \int_0^\infty \epsilon^{3/2} \left[ z e^{-\beta \epsilon} - \delta z^2 e^{-2\beta \epsilon} \right] d\epsilon$
利用积分 $\int_0^\infty x^{3/2} e^{-ax} dx = \frac{3\sqrt{\pi}}{4} a^{-5/2}$ ：
$U = \frac{2\pi V}{h^3} (2m)^{3/2} \left[ z \frac{3\sqrt{\pi}}{4} \beta^{-5/2} - \delta z^2 \frac{3\sqrt{\pi}}{4} (2\beta)^{-5/2} \right] = \frac{3}{2} k_B T \frac{V}{\lambda^3} \left( z - \frac{\delta}{2^{5/2}} z^2 \right)$
为了将此式与经典统计学结果进行直接对比，我们将前面求得的 $z$ 的反解表达式代入 $U$ 中，并在展开时仅保留至 $(n \lambda^3)^2$ 的一阶修正项（即取 $z^2 \approx (n\lambda^3)^2$ ）：
$U \approx \frac{3}{2} k_B T \frac{V}{\lambda^3} \left[ n \lambda^3 \left( 1 + \frac{\delta}{2^{3/2}} n \lambda^3 \right) - \frac{\delta}{2^{5/2}} (n \lambda^3)^2 \right]$
提取公因子 $n \lambda^3$ 并利用恒等式 $\frac{V}{\lambda^3} n \lambda^3 = N$ ，对方括号内合并同类项（注意 $\frac{1}{2^{3/2}} - \frac{1}{2^{5/2}} = \frac{1}{2\sqrt{2}} - \frac{1}{4\sqrt{2}} = \frac{1}{2^{5/2}}$ ）：
$U \approx \frac{3}{2} N k_B T \left[ 1 + \left( \frac{1}{2^{3/2}} - \frac{1}{2^{5/2}} \right) \delta n \lambda^3 \right] = \frac{3}{2} N k_B T \left[ 1 + \frac{\delta}{2^{5/2}} n \lambda^3 \right]$
推导证明：弱简并气体的内能严格等于经典理想气体平动内能（ $\frac{3}{2}N k_B T$ ）附加一个由量子全同性引出的微小修正项。

5.8.3 物态方程的量子修正与统计等效相互作用#

对于非相对论性理想气体系统（无论经典还是量子分布），其宏观压强与平动内能的严密体积微商关系均满足 $p = \frac{2}{3} \frac{U}{V}$ 。将 5.8.2 节化简完毕的内能表达式直接代入此关系中，立即得到包含一阶量子修正的物态方程：
$p = \frac{2}{3V} \cdot \frac{3}{2} N k_B T \left[ 1 + \frac{\delta}{2^{5/2}} n \lambda^3 \right]$
整理为标准的状态方程形式：
$pV \approx N k_B T \left[ 1 + \frac{\delta}{2^{5/2}} n \lambda^3 \right] = N k_B T \left[ 1 + \delta \frac{N}{V} \frac{h^3}{16 (\pi m k_B T)^{3/2}} \right]$
对比经典理想气体状态方程 $pV = N k_B T$ ，方括号中的修正项完全起源于全同粒子的量子统计关联：
1. 对于理想玻色气体（ $\delta = -1$ ）：修正项为负，实际压强小于相同条件下的经典理想气体。这表明由于玻色子倾向于共同占据同一量子态的对称波函数特征，宏观上表现出一种“统计等效吸引力”，使得气体分子在热力学效应上显得彼此聚拢。
2. 对于理想费米气体（ $\delta = +1$ ）：修正项为正，实际压强大于经典的理想气体。这直接体现了泡利不相容原理的反对称波函数限制，排斥多个粒子占据同一相体积元，在宏观上表现出一种“统计等效排斥力”。

5.9 理想玻色气体的强简并条件与玻色-爱因斯坦凝聚#

当系统处于低温和高密度状态时，粒子的热德布罗意波长变得与粒子间平均间距相当，波函数发生显著的宏观交叠。此时，弱简并条件失效，系统进入强简并状态。本节专门探讨理想玻色气体在强简并极限下展现出的宏观量子特征。

5.9.1 玻色系统的化学势限制#

在平衡状态下，理想玻色气体在能级 $\epsilon_l$ 上的平均占据数为：
$a_l = \frac{g_l}{e^{(\epsilon_l - \mu)/k_B T} - 1}$
玻色子不受泡利不相容原理限制，但作为真实的物理系统，任何能级上的粒子数都必须保证为非负实数，即对于所有可能的 $\epsilon_l$ ，必须满足：
$a_l \ge 0 \implies e^{(\epsilon_l - \mu)/k_B T} > 1$
该不等式等价于 $\mu < \epsilon_l$ 。为了使该条件对系统内所有的能级均严格成立，化学势 $\mu$ 必须小于系统的最低能级（基态能量 $\epsilon_0$ ）。若取基态能量为能量零点（ $\epsilon_0 = 0$ ），则玻色气体的化学势受到绝对的代数限制：
$\mu \le 0$
这一限制是玻色系统引发宏观凝聚相变的前提。

5.9.2 连续能量积分与粒子数的级数展开#

对于宏观体积 $V$ 中的三维自由粒子系统，引入态密度函数 $D(\epsilon) = \frac{2\pi V}{h^3} (2m)^{3/2} \epsilon^{1/2}$ 。考虑一般情况（化学势 $\mu < 0$ ），引入参量逸度（Fugacity） $z = e^{\mu/k_B T}$ ，由上述限制可知 $0 < z < 1$ 。
将激发态上的粒子数 $N_e$ 转化为连续积分：
$N_e = \int_0^\infty \frac{D(\epsilon)}{e^{(\epsilon - \mu)/k_B T} - 1} d\epsilon = \frac{2\pi V}{h^3} (2m)^{3/2} \int_0^\infty \frac{\epsilon^{1/2}}{z^{-1}e^{\epsilon/k_B T} - 1} d\epsilon$
引入无量纲变量 $x = \epsilon/k_B T$ ，由于 $z e^{-x} < 1$ ，可对被积函数的分母进行几何级数展开：
$\frac{1}{z^{-1}e^x - 1} = \frac{z e^{-x}}{1 - z e^{-x}} = \sum_{n=1}^\infty z^n e^{-nx}$
代入积分并交换求和与积分顺序：
$N_e = \frac{2\pi V}{h^3} (2m k_B T)^{3/2} \sum_{n=1}^\infty z^n \int_0^\infty x^{1/2} e^{-nx} dx$
利用伽马函数积分 $\int_0^\infty x^{1/2} e^{-nx} dx = \frac{\Gamma(3/2)}{n^{3/2}} = \frac{\sqrt{\pi}}{2 n^{3/2}}$ ，化简可得化学势不为零时，激发态粒子数的严格级数解析式：
$N_e = V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2} \sum_{n=1}^\infty \frac{z^n}{n^{3/2}}$
定义玻色-爱因斯坦积分函数 $g_{3/2}(z) = \sum_{n=1}^\infty \frac{z^n}{n^{3/2}}$ ，则系统总粒子数方程可精确表述为：
$N = N_0 + N_e = N_0 + V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2} g_{3/2}(z)$

5.9.3 粒子数容纳上限与玻色-爱因斯坦凝聚 (BEC)#

随着温度 $T$ 的降低，因子 $(k_B T)^{3/2}$ 减小；为了维持总粒子数 $N$ 恒定，逸度 $z$ 必须单调递增（即 $\mu$ 趋近于 $0$ ）。当 $\mu \to 0$ 时， $z \to 1$ 。
此时，激发态所能容纳的粒子数达到绝对上限 $N_{\max}(T)$ ：
$N_{\max}(T) = V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2} g_{3/2}(1)$
数学上，级数 $g_{3/2}(1) = \sum_{n=1}^\infty \frac{1}{n^{3/2}}$ 正好对应于黎曼 $\zeta$ 函数（Riemann zeta function）在 $s = 3/2$ 时的取值，即 $\zeta(3/2)$ 。因此，连续激发态的最大容纳量为：
$N_{\max}(T) = V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2} \zeta\left(\frac{3}{2}\right)$
相变判据：当温度降至某一临界值 $T_c$ 时，激发态的粒子容纳上限恰好等于系统总粒子数 $N$ ，即 $N_{\max}(T_c) = N$ 。解得临界温度解析式：
$T_c = \frac{h^2}{2\pi m k_B} \left[ \frac{N}{V \zeta(3/2)} \right]^{2/3}$
当系统温度继续降低（ $T < T_c$ ）时，由于 $\mu$ 无法突破 $0$ 的上限，连续谱的激发态无法再容纳系统中全部的宏观粒子。
多余的粒子必须被迫“坍缩”至被连续积分所忽略的单粒子极低能态（动量 $\mathbf{p}=0$ 的基态）。基态粒子数 $N_0$ 将占据宏观比例：
$N_0 = N - N_e(T) = N - V \left( \frac{2\pi m k_B T}{h^2} \right)^{3/2} \zeta\left(\frac{3}{2}\right)$
利用临界温度 $T_c$ 的定义代换，可得到极其优美的凝聚体分数方程：
$\frac{N_0}{N} = 1 - \left( \frac{T}{T_c} \right)^{3/2}$
数值结果与物理结论：
代入黎曼 $\zeta$ 函数的具体数值 $\zeta(3/2) \approx 2.612$ ，临界温度最终确认为：
$T_c \approx \frac{h^2}{2\pi m k_B} \left( \frac{n}{2.612} \right)^{2/3}$
当 $T < T_c$ 时，系统发生玻色-爱因斯坦凝聚（BEC）。宏观数量级的玻色子聚集在基态，形成巨型相干的物质波。该相变完全由粒子的不可分辨性及玻色统计规律所驱动，不存在经典意义上的空间气液相界面，是一种纯粹在动量空间内发生的连续量子相变。

5.10 理想费米气体的强简并条件与低温微扰展开#

当系统的温度远低于其特征简并温度（即 $T \ll T_F$ ）时，费米-狄拉克分布函数的分母中指数项不再能被视作大标量，系统进入强简并状态。由于泡利不相容原理的刚性限制，全同费米子体系在低温下表现出与经典和玻色系统截然不同的力学与热力学本征行为。

5.10.1 绝对零度下的费米海与费米能量 ( $T = 0\text{ K}$ )#

在绝对零度下，系统处于能量最低的量子基态。为了满足泡利不相容原理（每个微观量子态至多容纳一个费米子）， $N$ 个全同费米子不得不自底向上依次填满单粒子能级。
此时，费米-狄拉克分布函数退化为一个严格的阶跃函数：
$f(\epsilon) = \frac{1}{e^{(\epsilon - \mu)/k_B T} + 1} \xrightarrow{T \to 0} \Theta(E_F - \epsilon) = \begin{cases} 1, & \epsilon \le E_F \\ 0, & \epsilon > E_F \end{cases}$
式中 $E_F = \mu(0)$ 定义为绝对零度时的化学势，称为费米能量（Fermi Energy）。系统所有被占据的量子态在动量空间中构成的边界曲面称为费米面（Fermi Surface）。
引入粒子的内秉自旋简并度 $g_s = 2S + 1$ （对于电子， $S=1/2$ ， $g_s = 2$ ），三维平动自由费米子的态密度函数为：
$D(\epsilon) = g_s \frac{2\pi V}{h^3} (2m)^{3/2} \epsilon^{1/2}$
利用总粒子数守恒条件，通过连续谱积分严格确定费米能量 $E_F$ ：
$N = \int_0^{E_F} D(\epsilon) d\epsilon = g_s \frac{2\pi V}{h^3} (2m)^{3/2} \int_0^{E_F} \epsilon^{1/2} d\epsilon = g_s \frac{2\pi V}{h^3} (2m)^{3/2} \cdot \frac{2}{3} E_F^{3/2} = \frac{4\pi g_s V}{3h^3} (2m E_F)^{3/2}$
由此解得费米能量以宏观数密度 $n = N/V$ 表达的解析式：
$E_F = \frac{h^2}{2m} \left( \frac{3N}{4\pi g_s V} \right)^{2/3} = \frac{h^2}{2m} \left( \frac{3n}{4\pi g_s} \right)^{2/3}$
定义费米温度为 $T_F = E_F / k_B$ 。
进一步计算绝对零度下费米子体系的总基态内能 $U_0$ ：
$U_0 = \int_0^{E_F} \epsilon D(\epsilon) d\epsilon = g_s \frac{2\pi V}{h^3} (2m)^{3/2} \int_0^{E_F} \epsilon^{3/2} d\epsilon = g_s \frac{2\pi V}{h^3} (2m)^{3/2} \cdot \frac{2}{5} E_F^{5/2}$
将前面求得的 $N$ 的表达式代入上式中进行消项化简：
$U_0 = \left[ \frac{4\pi g_s V}{3h^3} (2m E_F)^{3/2} \right] \cdot \frac{3}{2 E_F^{3/2}} \cdot \frac{2}{5} E_F^{5/2} = \frac{3}{5} N E_F$
利用压强与内能的普适关系 $p = \frac{2}{3}\frac{U}{V}$ ，导出绝对零度下的零点简并压（Degeneracy Pressure）：
$p_0 = \frac{2}{3} \frac{U_0}{V} = \frac{2}{5} n E_F$
该代数结果表明，由于泡利不相容原理限制，即使在绝对零度下，费米气体也具有宏观的非零动能与极其巨大的压强，这构成了解释白矮星力学支撑及金属电子不坍缩的微观物理基础。

5.10.2 费米-狄拉克积分与索末菲展开原理 ( $T \ll T_F$ )#

当系统温度略微升高但仍处于强简并区间（ $0 < T \ll T_F$ ）时，费米面附近的阶跃分布函数发生微弱的平滑模糊。计算热力学量需要求解形如 $I = \int_0^\infty H(\epsilon) f(\epsilon) d\epsilon$ 的积分（其中 $H(\epsilon)$ 为关于能量的光滑函数）。为此，必须引入索末菲展开（Sommerfeld Expansion）进行低温级数微扰处理。
定义原函数 $G(\epsilon) = \int_0^\epsilon H(\epsilon') d\epsilon'$ ，显然有 $G(0) = 0$ 且 $\frac{dG}{d\epsilon} = H(\epsilon)$ 。对积分 $I$ 应用分部积分法：
$I = \int_0^\infty f(\epsilon) dG(\epsilon) = \left[ F(\epsilon) G(\epsilon) \right]_0^\infty - \int_0^\infty G(\epsilon) \frac{df}{d\epsilon} d\epsilon = \int_0^\infty G(\epsilon) \left( -\frac{df}{d\epsilon} \right) d\epsilon$
微商项 $-\frac{df}{d\epsilon} = \frac{\beta e^{\beta(\epsilon-\mu)}}{(e^{\beta(\epsilon-\mu)}+1)^2}$ 具有强烈的局域拓扑特征：当 $T \to 0$ 时，它在 $\epsilon = \mu$ 处收敛为狄拉克 $\delta$ 函数。因此，可在 $\epsilon = \mu$ 附近将 $G(\epsilon)$ 严格展开为泰勒级数：
$G(\epsilon) = G(\mu) + G'(\mu)(\epsilon - \mu) + \frac{1}{2} G''(\mu)(\epsilon - \mu)^2 + \frac{1}{6} G'''(\mu)(\epsilon - \mu)^3 + \dots$
将级数项代回积分链条中，并引入无量纲代换变量 $x = \beta(\epsilon - \mu) = \frac{\epsilon - \mu}{k_B T}$ ：
$I = \int_{-\frac{\mu}{k_B T}}^\infty \left[ G(\mu) + G'(\mu)(k_B T)x + \frac{1}{2}G''(\mu)(k_B T)^2 x^2 + \dots \right] \frac{e^x}{(e^x+1)^2} dx$
由于强简并下 $\frac{\mu}{k_B T} \gg 1$ ，可将积分下限严格延展至 $-\infty$ 。分别计算各项定积分：
1. 零阶项： $G(\mu) \int_{-\infty}^\infty \frac{e^x}{(e^x+1)^2} dx = G(\mu) \left[ -\frac{1}{e^x+1} \right]_{-\infty}^\infty = G(\mu) \cdot 1 = \int_0^\mu H(\epsilon) d\epsilon$
2. 一阶项（奇函数积分）：由于积分区间对称且被积函数为奇函数，该项严格为零：
$G'(\mu)(k_B T) \int_{-\infty}^\infty x \frac{e^x}{(e^x+1)^2} dx = 0$
1. 二阶项：利用留数定理或标准积分公式计算偶函数积分 $\int_{-\infty}^\infty \frac{x^2 e^x}{(e^x+1)^2} dx = \frac{\pi^2}{3}$ ：
$\frac{1}{2} G''(\mu)(k_B T)^2 \int_{-\infty}^\infty x^2 \frac{e^x}{(e^x+1)^2} dx = \frac{1}{2} H'(\mu) (k_B T)^2 \left( \frac{\pi^2}{3} \right) = \frac{\pi^2}{6} (k_B T)^2 H'(\mu)$
截断至二阶微扰项，索末菲展开的核心近似公式确立为：
$\int_0^\infty H(\epsilon) f(\epsilon) d\epsilon \approx \int_0^\mu H(\epsilon) d\epsilon + \frac{\pi^2}{6} (k_B T)^2 H'(\mu)$

5.10.3 低温下化学势与内能的解析修正#

利用 5.10.2 节的索末菲展开公式，可以精确计算封闭费米系统在有限低温下的热力学参量修正。
化学势 $\mu(T)$ 的温度依赖性：
令 $H(\epsilon) = D(\epsilon) = C \epsilon^{1/2}$ （其中 $C = g_s \frac{2\pi V}{h^3} (2m)^{3/2}$ ）。代入总粒子数公式：
$N = \int_0^\infty D(\epsilon) f(\epsilon) d\epsilon \approx \int_0^\mu C \epsilon^{1/2} d\epsilon + \frac{\pi^2}{6} (k_B T)^2 \left.\frac{d(C \epsilon^{1/2})}{d\epsilon}\right|_\mu = \frac{2}{3}C \mu^{3/2} + \frac{\pi^2}{12} C (k_B T)^2 \mu^{-1/2}$
我们已知在绝对零度下 $N = \frac{2}{3}C E_F^{3/2}$ 。建立等式关系并提取公因子：
$\frac{2}{3}C E_F^{3/2} \approx \frac{2}{3}C \mu^{3/2} \left[ 1 + \frac{\pi^2}{8} \left( \frac{k_B T}{\mu} \right)^2 \right] \implies E_F \approx \mu \left[ 1 + \frac{\pi^2}{8} \left( \frac{k_B T}{\mu} \right)^2 \right]^{2/3}$
由于 $\frac{k_B T}{\mu} \ll 1$ ，应用二项式展开 $(1+x)^{2/3} \approx 1 + \frac{2}{3}x$ ，进一步反解出化学势 $\mu$ ：
$\mu(T) \approx E_F \left[ 1 + \frac{\pi^2}{12} \left( \frac{k_B T}{E_F} \right)^2 \right]^{-1} \approx E_F \left[ 1 - \frac{\pi^2}{12} \left( \frac{k_B T}{E_F} \right)^2 \right]$
代数结果表明，随着温度升高，强简并费米气体的化学势以温度的二次方规律缓慢下降。
总内能 $U(T)$ 的微扰计算：
令 $H(\epsilon) = \epsilon D(\epsilon) = C \epsilon^{3/2}$ 。代入内能积分式中应用索末菲展开：
$U(T) = \int_0^\infty \epsilon D(\epsilon) f(\epsilon) d\epsilon \approx \int_0^\mu C \epsilon^{3/2} d\epsilon + \frac{\pi^2}{6} (k_B T)^2 \left.\frac{d(C \epsilon^{3/2})}{d\epsilon}\right|_\mu = \frac{2}{5}C \mu^{5/2} + \frac{\pi^2}{4} C (k_B T)^2 \mu^{1/2}$
将前面得到的化学势修正关系 $\mu \approx E_F [ 1 - \frac{\pi^2}{12} ( \frac{k_B T}{E_F} )^2 ]$ 代入上式的各项中，并严格保留至 $T^2$ 项：
$\mu^{5/2} \approx E_F^{5/2} \left[ 1 - \frac{\pi^2}{12} \left( \frac{k_B T}{E_F} \right)^2 \right]^{5/2} \approx E_F^{5/2} \left[ 1 - \frac{5\pi^2}{24} \left( \frac{k_B T}{E_F} \right)^2 \right]$
代回内能展开式中，保留至二阶项，合并整理：
$U(T) \approx \frac{2}{5}C E_F^{5/2} \left[ 1 - \frac{5\pi^2}{24} \left( \frac{k_B T}{E_F} \right)^2 \right] + \frac{\pi^2}{4} C (k_B T)^2 E_F^{1/2} = \frac{2}{5}C E_F^{5/2} + C E_F^{5/2} (k_B T / E_F)^2 \left[ -\frac{\pi^2}{12} + \frac{\pi^2}{4} \right]$
由于 $-\frac{\pi^2}{12} + \frac{\pi^2}{4} = \frac{\pi^2}{6}$ ，且由 5.10.1 节知 $U_0 = \frac{2}{5}C E_F^{5/2} = \frac{3}{5}NE_F$ ，上式可重写为直接与经典能量标度对应的形式：
$U(T) \approx U_0 + \frac{\pi^2}{6} C E_F^{1/2} (k_B T)^2 = U_0 + \frac{\pi^2}{6} \left( \frac{3N}{2 E_F^{3/2}} \right) E_F^{1/2} (k_B T)^2 = U_0 + \frac{\pi^2}{4} N \frac{(k_B T)^2}{E_F}$
将其完全表示为基态内能 $U_0$ 的修正倍数形式：
$U(T) \approx \frac{3}{5}N E_F \left[ 1 + \frac{5\pi^2}{12} \left( \frac{k_B T}{E_F} \right)^2 \right]$

5.10.4 简并费米子气体热容的温度依赖性#

基于包含二阶量子修正的系统总内能表达式，可以通过对温度的一阶微商直接导出强简并理想费米气体的定容热容 $C_V$ ：
$C_V = \left( \frac{\partial U}{\partial T} \right)_V = \frac{\partial}{\partial T} \left[ U_0 + \frac{\pi^2}{4} N \frac{k_B^2 T^2}{E_F} \right] = \frac{\pi^2}{2} N k_B \left( \frac{k_B T}{E_F} \right)$
引入费米温度 $T_F = E_F / k_B$ ，热容方程写为：
$C_V = \frac{\pi^2}{2} N k_B \frac{T}{T_F}$
物理机制解析：
根据经典能量均分定理，完全自由的可分辨电子气体的定容热容应当恒等于定值 $\frac{3}{2}N k_B$ ，这导致了经典金属热容理论在历史上与实验数据的严重偏离。
由统计力学导出的上式严格表明，强简并费米气体的热容不是常数，而是与绝对温度 $T$ 呈严格的线性一阶正比关系。
这一现象的物理本质在于泡利不相容原理造成的量子冻结：在满足 $T \ll T_F$ 的条件下，深埋在费米海内部的绝大多数费米子由于其周围的量子态已被完全占满，无法接受热激发。真正能够通过热涨落跃迁至激发态的，仅限于费米面附近厚度约为 $k_B T$ 狭窄能量区间内的少数费米子。该部分粒子数占总人数的比例约为 $T/T_F$ 。因此，只有这部分粒子能够贡献热容，使得宏观热容受到经典值附加一阶线性系数（ $\sim T/T_F$ ）的强烈压制。

5.11 玻色费米统计应用例#

5.11.1 晶格振动与固体热容的量子统计理论（声子气体）#

将近独立子系模型应用于晶格振动系统。晶体由 $N$ 个原子组成，每个原子在三维空间中进行微振动，系统总共具有 $3N$ 个振动自由度。在简谐近似下，晶格的微观热运动可等效为 $3N$ 个独立的线性量子谐振子。
爱因斯坦模型（Einstein Model）
爱因斯坦模型做出了最严格的零阶近似：假设晶体中所有的 $3N$ 个独立谐振子均以同一恒定频率 $\omega_E$ 振动。
单个量子谐振子的能级为 $\epsilon_n = \left(n + \frac{1}{2}\right)\hbar\omega_E$ 。根据 5.6 节推导的单谐振子配分函数 $Z_v$ （将 $h\nu$ 替换为 $\hbar\omega_E$ ），单振子平均能量为：
$\langle \epsilon \rangle = -\frac{\partial \ln Z_v}{\partial \beta} = \frac{1}{2}\hbar\omega_E + \frac{\hbar\omega_E}{e^{\beta\hbar\omega_E} - 1}$
由于 $3N$ 个振子相互独立，系统总内能 $U$ 为：
$U = 3N \langle \epsilon \rangle = \frac{3}{2}N\hbar\omega_E + \frac{3N\hbar\omega_E}{e^{\hbar\omega_E/k_B T} - 1}$
对温度求导得到固体定容热容 $C_V$ ：
$C_V = \left( \frac{\partial U}{\partial T} \right)_V = 3N k_B \left( \frac{\hbar\omega_E}{k_B T} \right)^2 \frac{e^{\hbar\omega_E/k_B T}}{(e^{\hbar\omega_E/k_B T} - 1)^2}$
定义爱因斯坦温度 $\Theta_E = \hbar\omega_E / k_B$ ：
1. 高温极限（ $T \gg \Theta_E$ ）：展开指数项 $e^{\Theta_E/T} \approx 1 + \Theta_E/T$ ，得到 $C_V \approx 3N k_B$ ，严格退化为经典的杜隆-珀蒂定律（Dulong-Petit law）。
2. 低温极限（ $T \ll \Theta_E$ ）：指数项主导分母，得到 $C_V \approx 3N k_B (\Theta_E/T)^2 e^{-\Theta_E/T}$ 。热容随温度呈指数衰减。
模型局限：实验表明，绝缘体在极低温下的真实热容呈 $T^3$ 规律衰减，而非指数衰减。爱因斯坦模型失效的根源在于其忽略了晶体原子的长波集体耦合振动模式。
德拜模型（Debye Model）与声子气体
德拜模型将晶体视为连续的弹性介质，将 $3N$ 个独立振子修正为具有不同频率的连续格波。由于格波能量是量子化的，能量量子 $\hbar\omega$ 被称为声子（Phonon）。声子服从玻色-爱因斯坦统计，且由于声子数不守恒（类似于光子），其化学势 $\mu \equiv 0$ 。
在连续介质近似下，格波色散关系呈线性： $\omega = v q$ （ $v$ 为声速， $q$ 为波矢）。三维空间中的波矢态密度为 $\frac{V}{(2\pi)^3} 4\pi q^2 dq$ 。由于弹性波包含 1 个纵波分支和 2 个横波分支，设平均声速为 $\bar{v}$ ，声子频率的态密度 $D(\omega)$ 为：
$D(\omega) d\omega = 3 \frac{V}{2\pi^2 \bar{v}^3} \omega^2 d\omega$
晶体具有有限的晶格常数，因而波长存在下限，频率存在上限（德拜截止频率 $\omega_D$ ）。截止频率由系统总自由度必须严格等于 $3N$ 的代数条件确定：
$\int_0^{\omega_D} D(\omega) d\omega = \int_0^{\omega_D} \frac{3V}{2\pi^2 \bar{v}^3} \omega^2 d\omega = \frac{V}{2\pi^2 \bar{v}^3} \omega_D^3 = 3N \implies \omega_D = \bar{v} \left( \frac{6\pi^2 N}{V} \right)^{1/3}$
系统总内能为各频率格波能量之和：
$U = U_0 + \int_0^{\omega_D} \frac{\hbar\omega}{e^{\hbar\omega/k_B T} - 1} D(\omega) d\omega = U_0 + \frac{3V \hbar}{2\pi^2 \bar{v}^3} \int_0^{\omega_D} \frac{\omega^3}{e^{\hbar\omega/k_B T} - 1} d\omega$
引入德拜温度 $\Theta_D = \hbar\omega_D / k_B$ ，并作无量纲代换 $x = \hbar\omega/k_B T$ ：
$U = U_0 + 9N k_B T \left( \frac{T}{\Theta_D} \right)^3 \int_0^{\Theta_D/T} \frac{x^3}{e^x - 1} dx$
求导得到德拜热容公式：
$C_V = 9N k_B \left( \frac{T}{\Theta_D} \right)^3 \int_0^{\Theta_D/T} \frac{x^4 e^x}{(e^x - 1)^2} dx$
低温极限展开（ $T \ll \Theta_D$ ）：
当温度趋于绝对零度时，积分上限 $\Theta_D/T \to \infty$ 。积分项收敛为常数 $\int_0^\infty \frac{x^4 e^x}{(e^x - 1)^2} dx = \frac{4\pi^4}{15}$ 。代入原式化简：
$C_V \approx \frac{12\pi^4}{5} N k_B \left( \frac{T}{\Theta_D} \right)^3$
此结果即为著名的德拜 $T^3$ 定律，与晶体在极低温下的实验测量结果严格吻合。

5.11.2 强简并费米气体的应用——金属泡利顺磁性#

考察外磁场对强简并费米气体（如金属中的自由电子气）宏观磁化强度的影响。电子具有自旋磁矩 $\mu_B$ （玻尔磁子， $\mu_B = \frac{e\hbar}{2m_e}$ ）。
磁场中的能带劈裂与粒子分布
在外加均匀弱磁场 $\mathbf{B}$ 中，电子自旋与磁场耦合，产生塞曼效应（Zeeman Effect）。自旋平行于磁场的电子（自旋向上， $+$ ）能量降低，自旋反平行于磁场的电子（自旋向下， $-$ ）能量升高。单电子的色散关系修正为：
$\epsilon_{\pm} = \epsilon_0 \mp \mu_B B$
其中 $\epsilon_0 = \frac{p^2}{2m_e}$ 为无磁场时的平动动能。
由于自旋简并被解除，自旋向上与向下的态密度函数分离，各自为总态密度（无自旋简并）的一半，且在能量坐标上发生反向平移：
$D_{\pm}(\epsilon) = \frac{1}{2} D_0(\epsilon \pm \mu_B B) \approx \frac{1}{2} D_0(\epsilon) \pm \frac{1}{2} \mu_B B D_0'(\epsilon)$
（其中 $D_0(\epsilon) = \frac{4\pi V}{h^3} (2m_e)^{3/2} \epsilon^{1/2}$ 为零磁场且包含自旋简并度的费米气体态密度）。
磁化强度的代数推导
在平衡态下，自旋向上和向下的电子气体共用同一个费米能级 $E_F$ （由系统总粒子数 $N = N_+ + N_-$ 的守恒严格约束）。
在绝对零度及低温强简并条件（ $T \ll T_F$ ）下，系统处于阶跃分布。自旋向上的粒子数 $N_+$ 与自旋向下的粒子数 $N_-$ 分别为连续谱积分：
$N_{\pm} = \int_0^{E_F} D_{\pm}(\epsilon) d\epsilon = \frac{1}{2} \int_0^{E_F} D_0(\epsilon \pm \mu_B B) d\epsilon$
由于外加磁场通常是弱场，满足磁耦合能量远小于费米能量（ $\mu_B B \ll E_F$ ），可以通过变量代换和泰勒级数将积分边界统一化：
$N_{\pm} = \frac{1}{2} \int_{\pm \mu_B B}^{E_F \pm \mu_B B} D_0(x) dx \approx \frac{1}{2} \int_0^{E_F} D_0(x) dx \pm \frac{1}{2} \int_{E_F}^{E_F \pm \mu_B B} D_0(x) dx$
由于下限处的态密度 $D_0(0) = 0$ ，下限的微扰不产生影响。第二项积分在边界狭窄区间内可用矩形面积近似：
$N_{\pm} \approx \frac{N}{2} \pm \frac{1}{2} \mu_B B D_0(E_F)$
系统宏观磁化强度 $M$ 等于单位体积内净磁矩的总和：
$M = \mu_B (N_+ - N_-) = \mu_B \left[ \left( \frac{N}{2} + \frac{1}{2} \mu_B B D_0(E_F) \right) - \left( \frac{N}{2} - \frac{1}{2} \mu_B B D_0(E_F) \right) \right] = \mu_B^2 B D_0(E_F)$
泡利磁化率的确定
由于零场时费米面处的态密度可表示为 $D_0(E_F) = \frac{3N}{2E_F}$ ，将其代入磁化强度公式：
$M = \mu_B^2 B \frac{3N}{2E_F}$
宏观体积磁化率 $\chi_P$ （单位体积内的磁化强度对外磁场的响应）定义为 $\chi_P = \frac{M}{V B}$ ，由此导出泡利顺磁磁化率（Pauli Paramagnetic Susceptibility）的解析解：
$\chi_P = \frac{3n \mu_B^2}{2E_F}$
根据经典统计物理推导的居里定律（Curie’s Law），系统的顺磁磁化率应与绝对温度成反比（ $\chi \propto 1/T$ ）。然而，上述量子统计推导严密证明：强简并费米气体的顺磁磁化率是一个极其微弱的正值，且在 $T \ll T_F$ 极限下与温度完全无关。
其物理机制在于泡利不相容原理的刚性限制。外磁场只能迫使处于费米面附近深度为 $\mu_B B$ 的极少量电子发生自旋翻转（粒子数比例约为 $\mu_B B / E_F$ ），而费米海深处的绝大多数电子处于双占据状态，它们的自旋响应被泡利不相容原理绝对锁死，从而导致宏观顺磁响应受到强烈的量子压制。

第6章系综理论#

本章将抛弃“近独立子系”模型的局限性。对于内部存在任意复杂相互作用的真实宏观系统，单粒子的 $\mu$ 空间无法再提供自洽的描述。我们必须将整个由 $N$ 个粒子构成的宏观系统作为一个单一的力学体系，在更高维度的相空间中研究其统计规律。这是现代统计物理学中最具普适性和严格性的理论框架。

6.1 $\Gamma$ 空间与刘维尔定理#

6.1.1 $\Gamma$ 空间的概念建立#

对于一个由 $N$ 个微观粒子构成的系统，若每个粒子具有 $3$ 个空间自由度，则整个系统的力学状态需要由 $s = 3N$ 个广义坐标 $q = (q_1, q_2, \dots, q_s)$ 和 $s$ 个广义动量 $p = (p_1, p_2, \dots, p_s)$ 共同完全确定。
这 $2s = 6N$ 个变量所张成的多维正交空间被称为系统的相空间，即 $\Gamma$ 空间（ $\Gamma$ -space）。
在 $\Gamma$ 空间中：
1. 相点（Representative Point）：在任意给定时刻 $t$ ，整个宏观系统所有的微观运动状态被 $\Gamma$ 空间中的一个唯一确定的几何点 $(q, p)$ 所代表。
2. 相轨线（Phase Trajectory）：随着时间演化，系统状态的改变在 $\Gamma$ 空间中描绘出一条连续的轨迹。其动力学演化严格服从经典力学的哈密顿正则方程：
$\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i} \quad (i = 1, 2, \dots, s)$
其中 $H(q, p)$ 为系统的总哈密顿量。只要哈密顿量不显含时间，相轨线在 $\Gamma$ 空间中绝对不会自身相交（除非在奇点）。

6.1.2 相体积与系综密度函数#

为了运用概率论，吉布斯（J. W. Gibbs）引入了系综（Ensemble）的概念。系综是指在相同宏观条件下（如相同的 $N, V, T$ 或 $E$ ），包含着大量在微观上互不干涉、由系统所有可能微观状态构成的假想系统集合。
在 $\Gamma$ 空间中，系综表现为海量相点构成的“流体”。定义相体积元为 $d\Gamma = dq_1 \dots dq_s dp_1 \dots dp_s$ 。
引入系综密度函数 $\rho(q, p, t)$ ，它表示在时刻 $t$ ，位于相点 $(q, p)$ 附近单位相体积内的系统数目（或概率密度）。
对于包含 $\mathcal{N}$ 个系统的系综，位于 $d\Gamma$ 内的系统数为 $d\mathcal{N} = \mathcal{N} \rho(q, p, t) d\Gamma$ 。
分布函数满足全空间归一化条件：
$\int \rho(q, p, t) d\Gamma = 1$

6.1.3 刘维尔定理 (Liouville’s Theorem) 的严格推导#

系综相点在 $\Gamma$ 空间中的运动不会凭空产生或消失，必须严格服从类似于流体力学的连续性方程：
$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$
在 $6N$ 维 $\Gamma$ 空间中，散度算符展开为对所有广义坐标和动量的微商，速度矢量 $\mathbf{v}$ 的分量为 $(\dot{q}_1, \dots, \dot{q}_s, \dot{p}_1, \dots, \dot{p}_s)$ 。连续性方程写为：
$\frac{\partial \rho}{\partial t} + \sum_{i=1}^s \left[ \frac{\partial}{\partial q_i}(\rho \dot{q}_i) + \frac{\partial}{\partial p_i}(\rho \dot{p}_i) \right] = 0$
利用乘积法则展开括号内的项：
$\frac{\partial \rho}{\partial t} + \sum_{i=1}^s \left( \frac{\partial \rho}{\partial q_i}\dot{q}_i + \frac{\partial \rho}{\partial p_i}\dot{p}_i \right) + \rho \sum_{i=1}^s \left( \frac{\partial \dot{q}_i}{\partial q_i} + \frac{\partial \dot{p}_i}{\partial p_i} \right) = 0$
应用哈密顿正则方程，考察速度场的散度项（即最后一项的括号内部分）：
$\frac{\partial \dot{q}_i}{\partial q_i} + \frac{\partial \dot{p}_i}{\partial p_i} = \frac{\partial}{\partial q_i}\left( \frac{\partial H}{\partial p_i} \right) + \frac{\partial}{\partial p_i}\left( -\frac{\partial H}{\partial q_i} \right) = \frac{\partial^2 H}{\partial q_i \partial p_i} - \frac{\partial^2 H}{\partial p_i \partial q_i} \equiv 0$
这在物理上意味着 $\Gamma$ 空间中的“相点流体”是绝对不可压缩的。
将散度为零代回连续性方程，利用泊松括号（Poisson Bracket） $[\rho, H]$ 的定义，得到描述系综密度演化的方程：
$\frac{\partial \rho}{\partial t} + \sum_{i=1}^s \left( \frac{\partial \rho}{\partial q_i} \frac{\partial H}{\partial p_i} - \frac{\partial \rho}{\partial p_i} \frac{\partial H}{\partial q_i} \right) = \frac{\partial \rho}{\partial t} + [\rho, H] = 0$
根据多元函数的全导数定义，上式直接等价于：
$\frac{d\rho}{dt} = 0$
刘维尔定理结论：沿着系统在 $\Gamma$ 空间中运动的相轨线，系综密度函数 $\rho$ 是一个绝对常数。任意一块包含固定数目相点的相体积 $d\Gamma$ ，在随时间演化的过程中，其形状可能发生剧烈扭曲，但其包含的超体积绝对大小保持不变（即相体积守恒）。
对于处于宏观热力学平衡态的系统，其统计性质不随时间显含变化，即要求 $\frac{\partial \rho}{\partial t} = 0$ 。根据演化方程，平衡态系综必然满足：
$[\rho, H] = 0$
这意味着在平衡统计物理中，分布密度函数 $\rho$ 必须是系统运动积分（Constants of motion）的函数。通常由于宏观守恒量的限制， $\rho$ 仅能是系统总哈密顿量 $H(q, p)$ 的函数，即 $\rho = \rho(H)$ 。

6.1.4 密度矩阵与量子刘维尔方程的推广#

经典力学中的相空间点和连续分布在量子力学中是不适用的。在量子统计中，系统状态由希尔伯特空间中的多体波函数或态矢量 $|\psi\rangle$ 描述，系综分布被推广为密度矩阵算符（Density Operator） $\hat{\rho}$ 。
若系综以概率 $w_k$ 处于正交归一的纯态 $|\psi_k\rangle$ ，密度算符定义为：
$\hat{\rho} = \sum_k w_k |\psi_k\rangle \langle \psi_k|$
经典刘维尔方程在量子力学中直接映射为冯·诺依曼方程（von Neumann Equation）。利用薛定谔方程 $i\hbar \frac{\partial}{\partial t}|\psi\rangle = \hat{H}|\psi\rangle$ ，求导可得密度算符的演化律：
$i\hbar \frac{\partial \hat{\rho}}{\partial t} = [\hat{H}, \hat{\rho}]$
其中 $[\hat{H}, \hat{\rho}] = \hat{H}\hat{\rho} - \hat{\rho}\hat{H}$ 为对易子（Commutator），严格对应于经典的泊松括号 $i\hbar \{H, \rho\}$ 。
对于量子平衡态系统 $\frac{\partial \hat{\rho}}{\partial t} = 0$ ，必然要求 $[\hat{H}, \hat{\rho}] = 0$ 。这表明处于平衡态的量子系统，其密度算符 $\hat{\rho}$ 必须与哈密顿量算符 $\hat{H}$ 具有共同的本征态群，矩阵在能量表象中是对角化的。

6.2 统计物理的基本公设与系综分类#

在 6.1 节中，刘维尔定理证明了 $\Gamma$ 空间中代表系综的相流是不可压缩的，且平衡态密度函数 $\rho$ 必然是运动积分的函数。为了将这一纯粹的动力学演化与实际可观测的宏观热力学量联系起来，统计物理学必须确立其核心逻辑公设，并根据系统受到的宏观外部约束对系综进行严密的物理分类。

6.2.1 各态历经假说与等概率原理#

一个具有 $s = 3N$ 个自由度的宏观孤立系统，其总能量 $E$ 严格守恒。在动力学上，这限制了系统的代表相点 $(q, p)$ 只能在 $\Gamma$ 空间中由方程 $H(q, p) = E$ 所确定的 $(6N-1)$ 维能量超曲面（Energy Hypersurface）上运动。
为了将系统的长时间动力学行为与概率统计相联系，玻尔兹曼提出了各态历经假说（Ergodic Hypothesis）：
只要经过足够长的时间，孤立系统的代表相点将穿过（或无限逼近）能量超曲面上所有满足守恒律的微观相点。
由于相点在能量超曲面上是均匀“巡游”的，系统在其演化历史中，停留在超曲面上任意大小相等的微小相体积元 $d\Gamma$ 内的时间比例是相同的。将时间演化转化为系综空间分布，这直接导出了经典统计力学的核心基础——等概率原理（Principle of Equal A Priori Probabilities）：
对于一个处于平衡态的孤立系统，在满足系统宏观守恒条件（如能量 $E$ 、粒子数 $N$ 、体积 $V$ 恒定）的 $\Gamma$ 空间允许区域内，系综密度函数 $\rho(q, p)$ 处处相等（即常数）；而在该区域之外， $\rho(q, p)$ 严格为零。

6.2.2 宏观物理量的系综平均与时间平均#

在实验室中，任何宏观物理量（如压强 $p$ 、磁化强度 $M$ ）的测量过程都需要经历一段宏观时间 $\tau$ 。在此期间，微观相点在 $\Gamma$ 空间中飞速划过极其漫长的轨迹。因此，实验测得的物理量 $O_{obs}$ 在本质上是微观力学量 $O[q(t), p(t)]$ 在观测时间 $\tau$ 内的时间平均值：
$\overline{O} = \lim_{\tau \to \infty} \frac{1}{\tau} \int_0^\tau O[q(t), p(t)] dt$
然而，由于系统自由度巨大（ $N \sim 10^{23}$ ），求解 $6N$ 个哈密顿演化方程以求得时间平均值在解析上是绝对不可能的。
吉布斯系综理论的伟大之处在于，利用各态历经假说和等概率原理，将上述无法计算的“随时间积分”在数学上严格转换为“在 $\Gamma$ 空间中对所有可能状态的加权积分”。定义物理量 $O(q, p)$ 的系综平均值（Ensemble Average）为：
$\langle O \rangle = \int_{\Gamma} O(q, p) \rho(q, p) d\Gamma$
统计物理的基本公设断言，对于平衡系统，物理量的时间平均值严格等于系综平均值：
$O_{obs} = \overline{O} \equiv \langle O \rangle$
这一等式是整个现代宏观物理学能够通过微观第一性原理进行定量计算的绝对基石。

6.2.3 三大系综的物理界定#

自然界中的宏观系统根据其与外界的相互作用边界条件，受到不同的热力学约束。为了普适地处理各类边界条件，吉布斯构造了三种绝对的基础系综模型：
1. 微正则系综 (Microcanonical Ensemble)
  物理模型：代表完全孤立的系统，与外界既无能量交换，也无物质交换。
  宏观约束量：系统的粒子数 $N$ 、体积 $V$ 和总能量 $E$ 是严格固定的常数（考虑极其微弱的能量涨落，通常设能量在 $E$ 到 $E + \Delta E$ 的无限窄壳层内）。
  统计特征：在该能量壳层内， $\rho(q, p) = \text{const}$ ；壳层外 $\rho = 0$ 。这是最基础的系综，但数学积分极为困难。
2. 正则系综 (Canonical Ensemble)
  物理模型：代表被刚性透热壁包裹的封闭系统，它与一个温度为 $T$ 的绝对巨大的恒温热浴（Heat Bath）处于热接触平衡状态。
  宏观约束量：系统的粒子数 $N$ 、体积 $V$ 和绝对温度 $T$ 为给定常数。
  统计特征：系统的总能量 $E$ 不再守恒，而是在一个平均值周围发生热涨落。代表系统的相点分布于整个 $\Gamma$ 空间，不受能量超曲面的硬约束。
3. 巨正则系综 (Grand Canonical Ensemble)
  物理模型：代表被多孔透热壁包裹的开放系统，它与外界同时存在热量交换与粒子交换（与巨大的温度浴和粒子源相连）。
  宏观约束量：系统的绝对温度 $T$ 、体积 $V$ 和化学势 $\mu$ 为给定常数。
  统计特征：不仅系统的能量 $E$ 会发生涨落，系统内部的粒子总数 $N$ 也是一个动态涨落的随机变量。此时代表系统的相点存在于不同粒子数维度的 $\Gamma$ 空间的直和空间中。

6.3 微正则系综 (Microcanonical Ensemble) 与热力学基本方程#

在确立了等概率原理后，我们将这一公设首先应用于与外界绝对隔离的孤立系统。微正则系综是统计物理中最基本、最直接反映动力学守恒定律的系综模型，它是建立宏观绝对熵与微观动力学之间严格数学映射的起点。

6.3.1 微正则分布与相体积的几何界定#

对于一个由 $N$ 个粒子构成、体积为 $V$ 的绝对孤立系统，其总能量 $E$ 应当是一个严格的守恒量。然而，在真实的宏观物理测量中，绝对无限精确的能量隔离是不存在的，系统能量总是存在一个极其微小的观测不确定度或宏观涨落 $\Delta E$ （且 $\Delta E \ll E$ ）。
因此，系统的代表相点并不局限于一个无限薄的 $(6N-1)$ 维超曲面，而是被限制在由 $H(q, p) = E$ 和 $H(q, p) = E + \Delta E$ 两个超曲面所夹的高维能量壳层（Energy Shell）内。
根据各态历经假说与等概率原理，微正则系综的密度函数 $\rho(q, p)$ 在 $\Gamma$ 空间中表现为严格的“方波”截断分布：
$\rho(q, p) = \begin{cases} C, & E \le H(q, p) \le E + \Delta E \\ 0, & \text{其他区域} \end{cases}$
其中 $C$ 为归一化常数。系统在能量壳层内所占据的经典相空间超体积 $\Delta \Gamma$ 可通过对该壳层区域进行 $6N$ 重积分严格求得：
$\Delta \Gamma = \int_{E \le H \le E + \Delta E} d\Gamma$
由于等概率原理，系统出现在壳层内任意等体积相元中的概率完全相同。因此，这个宏观的相体积 $\Delta \Gamma$ 直接正比于系统可达的微观状态总数。

6.3.2 玻尔兹曼熵的绝对定义与状态数的计算方法#

为了将经典的连续相体积转化为有绝对物理意义的离散微观状态数 $\Omega$ ，需要引入量子力学的两项修正：相体积的量子化单元 $h^{3N}$ ，以及全同粒子的不可分辨性因子 $1/N!$ 。
据此，微正则系综的微观状态数 $\Omega(N, V, E)$ 首先被绝对定义为能量壳层内的积分：
$\Omega(N, V, E) = \frac{\Delta \Gamma}{N! h^{3N}} = \frac{1}{N! h^{3N}} \int_{E \le H \le E + \Delta E} d\Gamma$
微观状态数的实用计算方法：
在绝大多数高维物理模型中，直接对由两个等能面夹成的壳层进行多重几何积分是极其困难的。为了获得精确的代数解析解，必须引入能量面所包围的总相体积 $\Sigma(E)$ 的辅助函数法。
定义 $\Sigma(E)$ 为 $\Gamma$ 空间中所有满足 $H(q, p) \le E$ 的相点构成的超球体体积：
$\Sigma(E) = \int_{H(q, p) \le E} d\Gamma$
计算无约束超球体的体积 $\Sigma(E)$ 在代数上通常可以严格执行。由于能量涨落壳层厚度 $\Delta E$ 是宏观小量（ $\Delta E \ll E$ ），壳层的相体积 $\Delta \Gamma$ 可以在数学上转化为总相体积 $\Sigma(E)$ 对能量 $E$ 的一阶偏导数与 $\Delta E$ 的乘积：
$\Delta \Gamma = \Sigma(E + \Delta E) - \Sigma(E) \approx \frac{\partial \Sigma(E)}{\partial E} \Delta E$
将此微分关系代回绝对定义式中，我们得到了微正则系综中计算 $\Omega$ 的标准代数范式：
$\Omega(N, V, E) = \frac{1}{N! h^{3N}} \frac{\partial \Sigma(E)}{\partial E} \Delta E$
确立了可达微观状态数后，宏观热力学熵 $S$ 的绝对统计定义由玻尔兹曼关系直接给出：
$S(N, V, E) = k_B \ln \Omega(N, V, E)$
(注：在热力学极限 $N \to \infty$ 下，微观状态数随能量呈极度陡峭的指数级增长。数学上可以严格证明，无论能量壳层厚度 $\Delta E$ 取何种宏观小量，甚至直接用总相体积量度 $\frac{\Sigma(E)}{N!h^{3N}}$ 替代 $\Omega$ ，它们对 $\ln \Omega$ 的相对偏差仅为 $\mathcal{O}(\frac{\ln N}{N})$ 量级，在宏观极限下完全消隐。因此，宏观熵的定义与 $\Delta E$ 的具体数值绝对解耦。)

6.3.3 热力学基本方程的微观导出#

微正则系综的终极目的是在没有任何经验假设的前提下，纯粹通过对微观状态数函数的微积分操作推导出全部热力学定律。
已知系统熵是能量 $E$ 、体积 $V$ 和粒子数 $N$ 的多元函数 $S(E, V, N)$ 。对其取全微分：
$dS = \left( \frac{\partial S}{\partial E} \right)_{V,N} dE + \left( \frac{\partial S}{\partial V} \right)_{E,N} dV + \left( \frac{\partial S}{\partial N} \right)_{E,V} dN$
将此式与经典热力学基本方程 $dS = \frac{1}{T} dE + \frac{p}{T} dV - \frac{\mu}{T} dN$ 进行各项系数的绝对比对。这种等价性赋予了宏观热力学参量以极其冷峻的微观偏导数定义：
1. 绝对温度 $T$ ：温度在本质上是系统熵随能量增加的纯粹导数斜率：
$\frac{1}{T} = \left( \frac{\partial S}{\partial E} \right)_{V,N} = k_B \left( \frac{\partial \ln \Omega}{\partial E} \right)_{V,N}$
1. 宏观压强 $p$ ：压强本质上是系统微观状态数对外部空间自由度（体积）膨胀的响应率：
$p = T \left( \frac{\partial S}{\partial V} \right)_{E,N} = k_B T \left( \frac{\partial \ln \Omega}{\partial V} \right)_{E,N}$
1. 化学势 $\mu$ ：化学势衡量了引入额外粒子（增加内部自由度约束）时，系统相体积变化所必须付出的代数代价：
$\mu = -T \left( \frac{\partial S}{\partial N} \right)_{E,V} = -k_B T \left( \frac{\partial \ln \Omega}{\partial N} \right)_{E,V}$

6.4 正则系综 (Canonical Ensemble) 与能量涨落#

微正则系综虽然在逻辑上最为基础，但其要求系统绝对孤立、能量绝对守恒的约束在实际物理测量中极为苛刻，且 $\Gamma$ 空间中能量壳层积分的数学处理往往遇到不可逾越的代数障碍。为了获得具有强大可计算性的理论框架，我们必须引入正则系综，将系统状态的硬性守恒约束（ $E$ 恒定）松弛为平均守恒约束（与恒温热浴耦合）。

6.4.1 正则分布律的微正则导出#

考虑一个被刚性透热壁包裹的宏观系统 $A$ （具有确定的粒子数 $N$ 和体积 $V$ ），它与一个无限大的恒温热浴 $A'$ 处于热接触平衡。热浴极其庞大，满足自由度 $N' \gg N$ 且能量 $E' \gg E_A$ 。
复合系统 $A^{(0)} = A + A'$ 与外界完全隔离，构成一个微正则系综。其总能量 $E^{(0)} = E_i + E'$ 严格守恒（其中 $E_i$ 为系统 $A$ 处于某一特定微观态 $i$ 时的能量，相互作用能极小可忽略）。
根据等概率原理，复合系统 $A^{(0)}$ 处于任一微观态的概率均相等。因此，子系统 $A$ 处于特定微观态 $i$ 的概率 $P_i$ ，直接正比于此时热浴 $A'$ 所能容纳的微观状态总数 $\Omega'(E')$ ：
$P_i \propto \Omega'(E^{(0)} - E_i)$
由于 $E_i \ll E^{(0)}$ ，且 $\ln \Omega'$ 是能量的极其平缓的函数（广延量），可在 $E^{(0)}$ 处对 $\ln \Omega'(E^{(0)} - E_i)$ 进行泰勒展开：
$\ln \Omega'(E^{(0)} - E_i) \approx \ln \Omega'(E^{(0)}) - \left( \frac{\partial \ln \Omega'}{\partial E'} \right)_{E^{(0)}} E_i + \dots$
根据微正则系综中温度的绝对定义 $\frac{\partial \ln \Omega'}{\partial E'} = \frac{1}{k_B T}$ （其中 $T$ 为热浴温度，处于热平衡的子系统 $A$ 必然与之同温）：
$\ln \Omega'(E^{(0)} - E_i) \approx \text{const} - \frac{E_i}{k_B T}$
取指数并引入规范常数算符 $\beta = \frac{1}{k_B T}$ ，导出统计物理中最核心的玻尔兹曼分布律（正则分布）：
$P_i = \frac{1}{Z_N} e^{-\beta E_i}$

6.4.2 配分函数 $Z_N$ 与亥姆霍兹自由能 $F$ #

为了使概率分布满足归一化条件 $\sum_i P_i = 1$ ，归一化常数 $Z_N$ 被定义为系统配分函数（Partition Function）：
$Z_N = \sum_i e^{-\beta E_i}$
对于连续的经典 $\Gamma$ 空间，引入相体积量子化及全同性修正，配分函数表现为 $6N$ 重相空间积分：
$Z_N(T, V, N) = \frac{1}{N! h^{3N}} \int_{\Gamma} e^{-\beta H(q, p)} d\Gamma$
该积分为无约束全空间积分，在代数上（如高斯积分）极其容易解析。
为了确立正则系综与宏观热力学的桥梁，引入信息熵公式 $S = -k_B \sum_i P_i \ln P_i$ 。代入正则分布律 $P_i = \frac{e^{-\beta E_i}}{Z_N}$ ：
$S = -k_B \sum_i P_i (-\beta E_i - \ln Z_N) = k_B \beta \sum_i P_i E_i + k_B \ln Z_N \sum_i P_i$
由于系统平均内能 $U = \langle E \rangle = \sum_i P_i E_i$ ，上式化简为：
$S = \frac{U}{T} + k_B \ln Z_N \implies U - TS = -k_B T \ln Z_N$
根据热力学定义，亥姆霍兹自由能 $F = U - TS$ 。由此得到极其严密的宏观热力学映射方程：
$F = -k_B T \ln Z_N$
一旦通过微观哈密顿量算出了 $Z_N$ ，系统的全部热力学量均可通过对 $F$ 的偏导数直接导出（如 $S = -\left(\frac{\partial F}{\partial T}\right)_{V,N}$ , $p = -\left(\frac{\partial F}{\partial V}\right)_{T,N}$ , $\mu = \left(\frac{\partial F}{\partial N}\right)_{T,V}$ ）。

6.4.3 能量涨落定理与系综等价性#

在正则系综中，系统的能量不再是守恒量，而是在期望值 $\langle E \rangle$ 附近发生热涨落。
内能 $\langle E \rangle$ 的代数计算为：
$\langle E \rangle = \frac{1}{Z_N} \sum_i E_i e^{-\beta E_i} = -\frac{1}{Z_N} \frac{\partial Z_N}{\partial \beta} = -\frac{\partial \ln Z_N}{\partial \beta}$
衡量能量涨落幅度的绝对量为能量均方差 $\langle (\Delta E)^2 \rangle = \langle E^2 \rangle - \langle E \rangle^2$ 。对上述平均能量方程再次对 $\beta$ 求偏导：
$\frac{\partial \langle E \rangle}{\partial \beta} = -\frac{\partial}{\partial \beta} \left( \frac{1}{Z_N} \sum_i E_i e^{-\beta E_i} \right) = \frac{1}{Z_N^2} \frac{\partial Z_N}{\partial \beta} \sum_i E_i e^{-\beta E_i} - \frac{1}{Z_N} \sum_i (-E_i^2) e^{-\beta E_i}$ $= -\langle E \rangle^2 + \langle E^2 \rangle = \langle (\Delta E)^2 \rangle$
利用宏观热容定义 $C_V = \left(\frac{\partial \langle E \rangle}{\partial T}\right)_V = \left(\frac{\partial \langle E \rangle}{\partial \beta}\right)_V \frac{d\beta}{dT} = \left(\frac{\partial \langle E \rangle}{\partial \beta}\right)_V \left(-\frac{1}{k_B T^2}\right)$ ，导出能量涨落定理：
$\langle (\Delta E)^2 \rangle = k_B T^2 C_V$
由于系统的内能 $U \propto N$ ，热容 $C_V \propto N$ ，系统的相对涨落尺度为：
$\frac{\sqrt{\langle (\Delta E)^2 \rangle}}{\langle E \rangle} \propto \frac{\sqrt{N}}{N} = \frac{1}{\sqrt{N}}$
系综等价性结论：当系统处于热力学极限（ $N \sim 10^{23}$ ）时，相对能量涨落趋近于 $\sim 10^{-11}$ ，这远远超出了任何宏观实验的测量精度。
这意味着，在正则系综中，概率分布 $P(E)$ 表现为一个以 $\langle E \rangle$ 为中心、宽度极窄的狄拉克 $\delta$ 函数状尖峰。系统绝大部分时间严格停留在微正则系综所定义的能量壳层内。由此在数学上绝对证明了：微正则系综与正则系综在热力学极限下是完全等价的。

6.5 巨正则系综 (Grand Canonical Ensemble) 与粒子数涨落#

在微正则系综中，系统的能量 $E$ 和粒子数 $N$ 是固定不变的；在正则系综中，能量约束被松弛，系统与温度为 $T$ 的热浴进行能量交换。现在，我们进一步解除物质交换的限制。
考虑一个存在多孔透热壁的开放系统，它不仅与外界交换能量，还交换粒子。这种同时与恒温热浴和恒定粒子源（化学势浴）耦合的统计模型，即为巨正则系综。它是处理量子多体系统、相变涨落以及化学反应的强有力工具。

6.5.1 巨正则分布律的微正则导出#

设有一个被多孔透热壁包裹的小系统 $A$ （体积固定的开放系统），它被浸泡在一个宏观庞大的热浴兼粒子源 $A'$ 中。复合系统 $A^{(0)} = A + A'$ 是一个孤立的微正则系综，其总能量 $E^{(0)}$ 和总粒子数 $N^{(0)}$ 严格守恒。
此时，系统 $A$ 的微观状态由两个变量联合定义：它包含了 $N_i$ 个粒子，且处于由这 $N_i$ 个粒子构成的第 $i$ 号多体量子态（能量为 $E_{i, N_i}$ ）。
当系统 $A$ 处于特定的微观态 $(i, N_i)$ 时，热浴 $A'$ 的能量被扣除为 $E' = E^{(0)} - E_{i, N_i}$ ，粒子数被扣除为 $N' = N^{(0)} - N_i$ 。
根据等概率原理，系统 $A$ 处于该特定状态的概率 $P_{i, N_i}$ 正比于此时热浴所能容纳的微观状态数 $\Omega'(E', N')$ ：
$P_{i, N_i} \propto \Omega'(E^{(0)} - E_{i, N_i}, N^{(0)} - N_i)$
由于 $E_{i, N_i} \ll E^{(0)}$ 且 $N_i \ll N^{(0)}$ ，对广延量 $\ln \Omega'$ 在 $(E^{(0)}, N^{(0)})$ 处进行二元泰勒展开，并截断至一阶：
$\ln \Omega' \approx \ln \Omega'_0 - \left( \frac{\partial \ln \Omega'}{\partial E'} \right)_{V', N'} E_{i, N_i} - \left( \frac{\partial \ln \Omega'}{\partial N'} \right)_{E', V'} N_i$
利用微正则系综给出的热力学偏导数定义： $\frac{\partial \ln \Omega'}{\partial E'} = \frac{1}{k_B T}$ （其中 $T$ 为绝对温度）， $\frac{\partial \ln \Omega'}{\partial N'} = -\frac{\mu}{k_B T}$ 。代入展开式：
$\ln \Omega' \approx \text{const} - \frac{E_{i, N_i}}{k_B T} + \frac{\mu N_i}{k_B T}$
取指数，引入规范常数 $\beta = \frac{1}{k_B T}$ ，我们便得到了巨正则分布律（即系统具有 $N$ 个粒子且处于第 $i$ 态的概率）：
$P_{i, N} = \frac{1}{\Xi} e^{-\beta (E_{i,N} - \mu N)}$

6.5.2 巨配分函数 $\Xi$ 与巨热力学势 $J$ #

为了使概率分布满足归一化条件 $\sum_{N=0}^\infty \sum_i P_{i, N} = 1$ ，必须引入巨配分函数（Grand Partition Function） $\Xi$ ：
$\Xi(T, V, \mu) = \sum_{N=0}^\infty \sum_i e^{-\beta (E_{i,N} - \mu N)} = \sum_{N=0}^\infty \left[ e^{\beta \mu N} \left( \sum_i e^{-\beta E_{i,N}} \right) \right]$
括号内的求和正是系统在给定粒子数 $N$ 时的正则配分函数 $Z_N(T, V)$ 。引入逸度（Fugacity） $z = e^{\beta \mu}$ ，巨配分函数可表示为各阶正则配分函数的母函数（离散拉普拉斯变换）：
$\Xi(T, V, z) = \sum_{N=0}^\infty z^N Z_N(T, V)$
确立微观与宏观桥梁的步骤与正则系综类似。代入信息熵公式 $S = -k_B \sum_{N,i} P_{i,N} \ln P_{i,N}$ ：
$S = -k_B \sum_{N,i} P_{i,N} (-\beta E_{i,N} + \beta \mu N - \ln \Xi) = \frac{U}{T} - \frac{\mu \langle N \rangle}{T} + k_B \ln \Xi$
整理得到 $U - TS - \mu N = -k_B T \ln \Xi$ 。
根据经典热力学，系统内能的欧拉积分为 $U = TS - pV + \mu N$ ，即 $U - TS - \mu N = -pV$ 。
由此，巨配分函数直接给出了巨热力学势（Grand Potential） $J$ ，或者等效于系统的压强方程：
$J = -pV = -k_B T \ln \Xi$
巨热力学势 $J(T, V, \mu)$ 的全微分为 $dJ = -S dT - p dV - N d\mu$ 。只需通过对 $\ln \Xi$ 求解偏导数，即可获得系统的全套宏观参量：
平均粒子数： $N = \langle N \rangle = \frac{1}{\beta} \left( \frac{\partial \ln \Xi}{\partial \mu} \right)_{T, V} = z \left( \frac{\partial \ln \Xi}{\partial z} \right)_{T, V}$
宏观熵： $S = -\left( \frac{\partial J}{\partial T} \right)_{V, \mu} = k_B \ln \Xi + k_B T \left( \frac{\partial \ln \Xi}{\partial T} \right)_{V, \mu}$
状态方程： $p = -\frac{J}{V} = \frac{k_B T}{V} \ln \Xi$
这部分推导是连接纯统计涨落与宏观热力学响应的经典桥梁。为了满足简洁且不过于刻板的要求，我将这部分的数学推导重新梳理，剥离了冗余的修饰词，重点突出代数变形的逻辑链条。
你可以将以下内容直接替换进 6.5.3 节：

6.5.3 粒子数涨落与等温压缩率#

在巨正则系综中，不仅能量发生涨落，系统占据固定体积 $V$ 内的瞬时粒子数 $N$ 也是一个涨落的随机变量。
首先计算平均粒子数 $\langle N \rangle$ ：
$\langle N \rangle = \frac{1}{\beta} \frac{\partial \ln \Xi}{\partial \mu} = \frac{1}{\Xi} \sum_{N,i} N e^{-\beta(E_{i,N} - \mu N)}$
为了寻找均方涨落 $\langle (\Delta N)^2 \rangle = \langle N^2 \rangle - \langle N \rangle^2$ ，对上式继续关于化学势 $\mu$ 求偏导，并利用乘积法则展开：
$\frac{\partial \langle N \rangle}{\partial \mu} = -\frac{1}{\Xi^2} \left( \frac{\partial \Xi}{\partial \mu} \right) \sum_{N,i} N e^{-\beta(E_{i,N} - \mu N)} + \frac{1}{\Xi} \sum_{N,i} \beta N^2 e^{-\beta(E_{i,N} - \mu N)}$
由于 $\frac{1}{\Xi} \frac{\partial \Xi}{\partial \mu} = \beta \langle N \rangle$ ，上式可以自然化简为：
$\frac{\partial \langle N \rangle}{\partial \mu} = -\beta \langle N \rangle^2 + \beta \langle N^2 \rangle = \beta \langle (\Delta N)^2 \rangle$
引入 $\beta = 1/k_B T$ ，我们便得到了联系微观涨落与系统状态参量的基本公式：
$\langle (\Delta N)^2 \rangle = k_B T \left( \frac{\partial \langle N \rangle}{\partial \mu} \right)_{T, V}$
为了将上式中的纯统计微商 $\left( \frac{\partial \langle N \rangle}{\partial \mu} \right)_{T, V}$ 与可测量的宏观力学响应（等温压缩率 $\kappa_T$ ）对应，我们可以借助热力学恒等式进行变量代换。
引入粒子数密度 $n = \langle N \rangle / V$ 。在体积 $V$ 固定的前提下，偏导数可写为：
$\left( \frac{\partial \langle N \rangle}{\partial \mu} \right)_{T, V} = V \left( \frac{\partial n}{\partial \mu} \right)_T$
根据恒温条件下的吉布斯-杜亥姆方程 $d\mu = v dp = \frac{1}{n} dp$ ，可以得到化学势随密度变化的链式关系：
$\left( \frac{\partial \mu}{\partial n} \right)_T = \frac{1}{n} \left( \frac{\partial p}{\partial n} \right)_T \implies \left( \frac{\partial n}{\partial \mu} \right)_T = n \left( \frac{\partial n}{\partial p} \right)_T$
接下来，利用密度的微分定义 $dn = d(N/V) = -\frac{N}{V^2}dV = -\frac{n}{V}dV$ （在保持 $N$ 恒定考察体积压缩时），代入等温压缩率的定义式 $\kappa_T = -\frac{1}{V}\left(\frac{\partial V}{\partial p}\right)_T$ ：
$\left( \frac{\partial n}{\partial p} \right)_T = -\frac{n}{V} \left( \frac{\partial V}{\partial p} \right)_T = n \kappa_T$
将上述推导逐层代回，即可完成桥接：
$\left( \frac{\partial \langle N \rangle}{\partial \mu} \right)_{T, V} = V \cdot (n \cdot n \kappa_T) = \frac{\langle N \rangle^2}{V} \kappa_T$
将其代入第一步得出的涨落公式，便导出了经典的粒子数涨落定理：
$\langle (\Delta N)^2 \rangle = k_B T \frac{\langle N \rangle^2}{V} \kappa_T$
系统的相对粒子数涨落表现为：
$\frac{\sqrt{\langle (\Delta N)^2 \rangle}}{\langle N \rangle} = \sqrt{\frac{k_B T \kappa_T}{V}} \propto \frac{1}{\sqrt{\langle N \rangle}}$
在大数极限下（ $\langle N \rangle \to \infty$ ），只要等温压缩率 $\kappa_T$ 保持为有限值，粒子数的相对涨落将趋近于零。此时巨正则系综的分布概率 $P(N)$ 收缩为一个尖峰，巨正则系综在宏观上与正则/微正则系综等效。

6.6 系综理论对量子统计的严格重构#

在第五章中，我们曾使用最概然分布法（即微正则系综的近似化身），通过引入拉格朗日乘子和斯特林近似，通过组合数学推导出了玻色与费米分布律。
现在，利用巨正则系综的数学框架，我们将抛弃粒子数 $N$ 的硬性约束约束，直接从 $\Gamma$ 空间向独立量子态空间进行代数映射，量子统计分布律的推导将变得异常简洁和优美。

6.6.1 巨配分函数的能级解耦#

考虑一个由无相互作用的全同粒子组成的理想量子气体系统。设单粒子能级为 $\epsilon_l$ ，系统处于某一特定多体微观态时，各能级上的粒子占据数为 $\{a_1, a_2, \dots, a_l, \dots\}$ 。
系统的总能量和总粒子数可分别表示为占据数的线性叠加：
$E = \sum_l a_l \epsilon_l, \quad N = \sum_l a_l$
在巨正则系综中，系统的巨配分函数定义为对所有可能的微观量子态求和：
$\Xi = \sum_{\{a_l\}} e^{-\beta(E - \mu N)} = \sum_{a_1, a_2, \dots} e^{-\beta \sum_l a_l (\epsilon_l - \mu)}$
利用指数函数的乘法性质，将求和号内的指数展开为连乘积：
$\Xi = \sum_{a_1, a_2, \dots} \left( \prod_l e^{-\beta a_l (\epsilon_l - \mu)} \right)$
这是整个推导中最具决定性的一步： 如果在正则系综中，由于存在 $\sum_l a_l = N$ 的刚性约束，各个 $a_l$ 之间是相互牵制的，求和与连乘绝对无法交换顺序。但在巨正则系综中，总粒子数 $N$ 是浮动的，每个能级上的占据数 $a_l$ 都可以完全独立地遍历其所有可能取值！
因此，多重求和与连乘可以无障碍地交换顺序，配分函数实现了完美的数学解耦：
$\Xi = \prod_l \left( \sum_{a_l} e^{-\beta a_l (\epsilon_l - \mu)} \right) \equiv \prod_l \Xi_l$
其中， $\Xi_l = \sum_{a_l} e^{-\beta a_l (\epsilon_l - \mu)}$ 被定义为单能级巨配分函数。宏观系统的巨热力学势由此化为各独立能级贡献的简单加和：
$\ln \Xi = \sum_l \ln \Xi_l$

6.6.2 费米-狄拉克与玻色-爱因斯坦分布的导出#

对于任意能级 $l$ ，其平均占据数 $\langle a_l \rangle$ 可以通过对该能级的巨配分函数求偏导极其简便地获得：
$\langle a_l \rangle = \frac{\sum_{a_l} a_l e^{-\beta a_l (\epsilon_l - \mu)}}{\sum_{a_l} e^{-\beta a_l (\epsilon_l - \mu)}} = -\frac{1}{\beta} \frac{\partial \ln \Xi_l}{\partial \epsilon_l}$
接下来，只需代入全同粒子的自旋统计特征，即可瞬间解析出分布律。
费米-狄拉克 (FD) 统计：
费米子服从泡利不相容原理，每个能级上的占据数只能取 $a_l = 0$ 或 $1$ 。
单能级巨配分函数为一个仅含两项的极简和式：
$\Xi_l = \sum_{a_l=0}^1 e^{-\beta a_l (\epsilon_l - \mu)} = 1 + e^{-\beta (\epsilon_l - \mu)}$
代入平均占据数公式，直接求偏导：
$\langle a_l \rangle = -\frac{1}{\beta} \frac{\partial}{\partial \epsilon_l} \ln \left( 1 + e^{-\beta (\epsilon_l - \mu)} \right) = -\frac{1}{\beta} \frac{-\beta e^{-\beta (\epsilon_l - \mu)}}{1 + e^{-\beta (\epsilon_l - \mu)}}$
分子分母同乘 $e^{\beta (\epsilon_l - \mu)}$ ，顷刻间得出 FD 分布：
$\langle a_l \rangle = \frac{1}{e^{\beta (\epsilon_l - \mu)} + 1}$
玻色-爱因斯坦 (BE) 统计：
玻色子不受泡利原理限制，占据数可以取 $a_l = 0, 1, 2, \dots, \infty$ 。
只要满足 $\epsilon_l - \mu > 0$ ，其单能级巨配分函数就是一个标准的无穷等比级数（公比 $q = e^{-\beta (\epsilon_l - \mu)} < 1$ ）：
$\Xi_l = \sum_{a_l=0}^\infty \left[ e^{-\beta (\epsilon_l - \mu)} \right]^{a_l} = \frac{1}{1 - e^{-\beta (\epsilon_l - \mu)}}$
对其取对数 $\ln \Xi_l = -\ln \left( 1 - e^{-\beta (\epsilon_l - \mu)} \right)$ ，并求偏导：
$\langle a_l \rangle = -\frac{1}{\beta} \frac{\partial}{\partial \epsilon_l} \left[ -\ln \left( 1 - e^{-\beta (\epsilon_l - \mu)} \right) \right] = \frac{1}{\beta} \frac{\beta e^{-\beta (\epsilon_l - \mu)}}{1 - e^{-\beta (\epsilon_l - \mu)}}$
分子分母同乘 $e^{\beta (\epsilon_l - \mu)}$ ，顺理成章得出 BE 分布：
$\langle a_l \rangle = \frac{1}{e^{\beta (\epsilon_l - \mu)} - 1}$
回顾第五章中动用海量微观态组合排列、斯特林近似展开以及拉格朗日待定系数法所耗费的漫长篇幅，巨正则系综由于天然契合粒子数不守恒的量子微观态图像，使得推导过程宛如摧枯拉朽般自然、严谨。

6.6.3 理想量子气体普适物态方程的解析证明#

利用巨正则系综，我们还可以非常优美地证明非相对论理想气体的普适状态方程 $pV = \frac{2}{3}U$ 。
根据 6.5.2 节的巨热力学势关系 $pV = k_B T \ln \Xi = k_B T \sum_l \ln \Xi_l$ 。将 FD 和 BE 统计的单能级配分函数统一写为：
$pV = \pm k_B T \sum_l \ln \left( 1 \pm e^{-\beta (\epsilon_l - \mu)} \right)$
（式中上方符号“ $+$ ”对应费米子，下方符号“ $-$ ”对应玻色子）。
将离散求和转化为连续相空间积分，引入三维自由粒子的态密度 $D(\epsilon) = C \epsilon^{1/2}$ （其中 $C$ 包含自旋简并度及常数项）：
$pV = \pm k_B T \int_0^\infty D(\epsilon) \ln \left( 1 \pm e^{-\beta (\epsilon - \mu)} \right) d\epsilon = \pm k_B T \int_0^\infty C \epsilon^{1/2} \ln \left( 1 \pm e^{-\beta (\epsilon - \mu)} \right) d\epsilon$
对该定积分应用分部积分法。令：
$u = \ln \left( 1 \pm e^{-\beta (\epsilon - \mu)} \right)$ ，则 $du = \frac{\mp \beta e^{-\beta (\epsilon - \mu)}}{1 \pm e^{-\beta (\epsilon - \mu)}} d\epsilon = \mp \beta \frac{1}{e^{\beta (\epsilon - \mu)} \pm 1} d\epsilon$
$dv = C \epsilon^{1/2} d\epsilon$ ，则 $v = \frac{2}{3} C \epsilon^{3/2} = \frac{2}{3} \epsilon D(\epsilon)$
代入积分公式 $\int u dv = [uv]_0^\infty - \int v du$ ：
边界项 $[uv]_0^\infty$ 在上限 $\epsilon \to \infty$ 时指数衰减为零，在下限 $\epsilon = 0$ 时因 $\epsilon^{3/2} = 0$ 亦为零。因此仅剩下 $-\int v du$ 项：
$pV = \pm k_B T \int_0^\infty \left[ \frac{2}{3} \epsilon D(\epsilon) \right] \left[ \pm \beta \frac{1}{e^{\beta (\epsilon - \mu)} \pm 1} d\epsilon \right]$
注意到 $\pm \cdot \pm = +$ ，且 $k_B T \cdot \beta = 1$ ，常数项 $\frac{2}{3}$ 可提取至积分号外：
$pV = \frac{2}{3} \int_0^\infty \epsilon \left[ \frac{1}{e^{\beta (\epsilon - \mu)} \pm 1} \right] D(\epsilon) d\epsilon$
方括号中的项正是量子平均占据数 $\langle a(\epsilon) \rangle$ 。积分号内部的式子 $\int_0^\infty \epsilon \langle a(\epsilon) \rangle D(\epsilon) d\epsilon$ 恰好是系统总平动内能 $U$ 的标准定义式。
由此，我们完成了极其纯粹的解析证明：
$pV = \frac{2}{3} U$
该物态方程不依赖于任何关于经典或量子的特化近似，不依赖于温度或密度，它仅根源于三维非相对论粒子的色散关系（ $\epsilon \propto p^2$ ），是统计物理框架下极其坚固的普适定律。

第七章涨落理论#

在前几章的讨论中，我们主要关注了宏观热力学量在平衡态下的统计平均值。然而，微观粒子的无规则热运动会导致系统的实际状态在平均值附近发生微小的动态偏离，这种自发偏离被称为涨落（Fluctuations）。
涨落不仅是检验统计物理理论自洽性的重要部分，它更深刻地联系着系统自发偏离平衡的能力与其在外部扰动下作出响应的能力。本章将系统性地探讨涨落现象的统计与准热力学理论。

7.1 涨落的广义系综理论#

广义系综理论提供了一种自下而上的微观推导方法。通过将配分函数视为数学上的母函数（Generating Function），我们可以利用连续的代数微商操作，非常简洁地导出各物理量涨落的统计矩。

7.1.1 广义力与广义位移的统计框架#

为了建立普适的涨落公式，我们考察一个处于温度 $T$ 的平衡系统，并引入广义外力参量 $y_i$ （如外加磁场 $B$ 、压强 $p$ 或化学势 $\mu$ ）以及与之共轭的广义位移（或广义坐标） $x_i$ （如磁化强度 $M$ 、体积 $V$ 或粒子数 $N$ ）。
假设系统的总哈密顿量可以分离为一个无外场干扰的基础部分 $H_0$ ，以及一系列广义力与广义位移线性耦合的势能项：
$H(x, y) = H_0(x) - \sum_i y_i x_i$
在正则系综的框架下，系统处于某一微观状态的概率分布为：
$P(x) = \frac{1}{\mathcal{Z}} e^{-\beta \left( H_0(x) - \sum_i y_i x_i \right)}$
其中，归一化常数 $\mathcal{Z}(y_1, y_2, \dots)$ 被定义为系统的广义配分函数，它包含了所有广义力参量的信息：
$\mathcal{Z}(y) = \sum_{x} e^{-\beta H_0(x)} e^{\beta \sum_i y_i x_i}$

7.1.2 物理量的平均值与均方涨落#

利用广义配分函数对参量 $y_i$ 求对数偏导，可以直接得到广义位移 $x_i$ 的统计平均值 $\langle x_i \rangle$ ：
$\langle x_i \rangle = \frac{\sum_x x_i e^{-\beta H(x, y)}}{\sum_x e^{-\beta H(x, y)}} = \frac{1}{\beta \mathcal{Z}} \frac{\partial \mathcal{Z}}{\partial y_i} = \frac{1}{\beta} \frac{\partial \ln \mathcal{Z}}{\partial y_i}$
为了考察该物理量的涨落幅度，我们研究其偏离平均值的程度 $\Delta x_i = x_i - \langle x_i \rangle$ 。涨落的平均值显然为零（ $\langle \Delta x_i \rangle = 0$ ），因此衡量涨落尺度的核心量是均方涨落 $\langle (\Delta x_i)^2 \rangle = \langle x_i^2 \rangle - \langle x_i \rangle^2$ 。
对平均值公式再次关于共轭参量 $y_i$ 求偏导：
$\frac{\partial \langle x_i \rangle}{\partial y_i} = \frac{\partial}{\partial y_i} \left( \frac{\sum_x x_i e^{-\beta H}}{\mathcal{Z}} \right) = \frac{\sum_x \beta x_i^2 e^{-\beta H}}{\mathcal{Z}} - \frac{\sum_x x_i e^{-\beta H}}{\mathcal{Z}^2} \frac{\partial \mathcal{Z}}{\partial y_i}$
这可以化简为：
$\frac{\partial \langle x_i \rangle}{\partial y_i} = \beta \langle x_i^2 \rangle - \beta \langle x_i \rangle^2 = \beta \langle (\Delta x_i)^2 \rangle$
引入 $\beta = 1/k_B T$ ，我们得到了涨落的普适基本方程：
$\langle (\Delta x_i)^2 \rangle = k_B T \frac{\partial \langle x_i \rangle}{\partial y_i}$

7.1.3 涨落-耗散关系#

上述公式具有非常深刻的热力学意义。等式右侧的偏导数 $\chi_{ii} = \frac{\partial \langle x_i \rangle}{\partial y_i}$ 称为广义极化率（Generalized Susceptibility）或响应函数。它衡量了系统在受到微小外部广义力 $y_i$ 扰动时，广义位移 $\langle x_i \rangle$ 所作出的宏观线性响应程度。
这个普适关系表明：系统自发产生的均方涨落，直接正比于其对外部对应扰动的宏观响应能力。
在之前的章节中推导过的具体涨落定理，均可视为该普适框架的特例：
1. 若取 $y = -1/T$ ， $x = E$ （能量），广义极化率关联到定容热容 $C_V$ ，即可得到 $\langle (\Delta E)^2 \rangle = k_B T^2 C_V$ 。
2. 若取 $y = \mu$ ， $x = N$ （粒子数），响应函数关联到等温压缩率 $\kappa_T$ ，即可得到 $\langle (\Delta N)^2 \rangle = k_B T \frac{\langle N \rangle^2}{V} \kappa_T$ 。
3. 若考察顺磁体，取 $y = B$ （磁场）， $x = M$ （磁化强度），则磁化强度的涨落直接正比于等温磁化率 $\chi_T$ 。

7.1.4 多变量的交叉关联与响应矩阵的对称性#

真实系统中往往存在多个物理量共同涨落的情况。为了描述两个不同物理量 $x_i$ 与 $x_j$ 之间的统计相关性，可以定义它们的协方差（交叉关联函数） $\langle \Delta x_i \Delta x_j \rangle$ 。
类似于单一变量的推导，我们将 $\langle x_i \rangle$ 对另一个参量 $y_j$ 求偏导：
$\frac{\partial \langle x_i \rangle}{\partial y_j} = \beta \langle x_i x_j \rangle - \beta \langle x_i \rangle \langle x_j \rangle = \beta \langle \Delta x_i \Delta x_j \rangle$
由此得到交叉涨落公式：
$\langle \Delta x_i \Delta x_j \rangle = k_B T \frac{\partial \langle x_i \rangle}{\partial y_j}$
根据数学上的偏导数交换律，广义配分函数的二阶混合偏导数与求导顺序无关：
$\frac{\partial \langle x_i \rangle}{\partial y_j} = \frac{\partial}{\partial y_j} \left( \frac{1}{\beta} \frac{\partial \ln \mathcal{Z}}{\partial y_i} \right) = \frac{1}{\beta} \frac{\partial^2 \ln \mathcal{Z}}{\partial y_j \partial y_i} = \frac{\partial \langle x_j \rangle}{\partial y_i}$
这一关系式证明了宏观响应矩阵的对称性，即系统参量 $x_i$ 对力 $y_j$ 的响应率，等同于参量 $x_j$ 对力 $y_i$ 的响应率。这种由微观统计关联直接导致的宏观对称性，构成了不可逆过程热力学中著名的昂萨格倒易关系（Onsager Reciprocal Relations）的静态微观基础。

7.2 涨落的准热力学理论#

在 7.1 节中，我们利用微观的配分函数导出了涨落的一般统计规律。然而，爱因斯坦（A. Einstein）在 1910 年提出了一种极为深刻的逆向物理视角：他跳出了微观哈密顿量的繁琐计算，将宏观的经典热力学定律直接向下延伸，通过反演玻尔兹曼公式，建立了一套纯粹基于宏观热力学参量的唯象涨落理论，即准热力学理论。

7.2.1 玻尔兹曼原理的反演与爱因斯坦公式#

玻尔兹曼公式 $S = k_B \ln \Omega$ 建立了孤立系统宏观熵与微观状态数 $\Omega$ 的绝对对应关系。爱因斯坦将其逻辑反转：既然 $\Omega$ 正比于系统处于某宏观态的概率，那么系统自发处于熵为 $S$ 的状态的概率 $P$ 必然满足：
$P \propto e^{S / k_B}$
考察一个由研究对象（体积较小的局部子系统）和广阔环境（庞大的恒温恒压浴）共同构成的绝对孤立复合系统。由于环境极大，其温度 $T_0$ 和压强 $p_0$ 保持恒定。
当子系统内部发生微小的自发涨落时，整个复合系统偏离了绝对平衡态。复合系统的总熵变 $\Delta S_{total}$ 等于子系统熵变 $\Delta S$ 与环境熵变 $\Delta S_{env}$ 之和：
$\Delta S_{total} = \Delta S + \Delta S_{env}$
环境吸收的热量为 $\Delta Q_{env} = -\Delta U - p_0 \Delta V$ （其中 $\Delta U$ 和 $\Delta V$ 分别为子系统内能和体积的涨落增量）。由于环境足够大，其熵变可通过可逆热力学过程计算： $\Delta S_{env} = \frac{\Delta Q_{env}}{T_0} = -\frac{\Delta U + p_0 \Delta V}{T_0}$ 。
因此，总熵变为：
$\Delta S_{total} = \Delta S - \frac{\Delta U + p_0 \Delta V}{T_0} = -\frac{\Delta U - T_0 \Delta S + p_0 \Delta V}{T_0} = -\frac{\Delta W_{\min}}{T_0}$
其中 $\Delta W_{\min}$ 为可逆条件下使子系统产生该宏观偏离所必须对它做的最小功。
将其代回爱因斯坦概率反演公式，我们得到了描述子系统发生热力学涨落的爱因斯坦涨落公式：
$P \propto \exp\left( \frac{\Delta S_{total}}{k_B} \right) = \exp\left( -\frac{\Delta W_{\min}}{k_B T_0} \right)$

7.2.2 熵的二阶展开与高斯正态近似#

由于孤立复合系统在绝对平衡态时的总熵达到极大值， $\Delta S_{total}$ 泰勒展开的一阶变分为零。因此，决定涨落概率分布的只能是二阶微小量。
我们将最小功做二阶泰勒展开。对于纯粹的热力学 $p-V-T$ 系统，子系统的内能是独立变量 $S$ 和 $V$ 的函数 $U(S, V)$ 。对 $\Delta U$ 进行展开：
$\Delta U \approx \left( \frac{\partial U}{\partial S} \right)_V \Delta S + \left( \frac{\partial U}{\partial V} \right)_S \Delta V + \frac{1}{2} \left[ \frac{\partial^2 U}{\partial S^2} (\Delta S)^2 + 2 \frac{\partial^2 U}{\partial S \partial V} \Delta S \Delta V + \frac{\partial^2 U}{\partial V^2} (\Delta V)^2 \right]$
利用热力学基本关系 $\left(\frac{\partial U}{\partial S}\right)_V = T_0$ 且 $\left(\frac{\partial U}{\partial V}\right)_S = -p_0$ ，展开式中的一阶项恰好抵消了最小功定义式中的一阶项（即证明了平衡态下的一阶微商为零）。
对于二阶项（方括号内的部分），利用微积分法则：
$\Delta \left( \frac{\partial U}{\partial S} \right)_V = \Delta T, \quad \Delta \left( \frac{\partial U}{\partial V} \right)_S = -\Delta p$
可将其代数化简为极其优美的对称二次型：
$\Delta W_{\min} \approx \frac{1}{2} (\Delta T \Delta S - \Delta p \Delta V)$
将其代回爱因斯坦公式，得出描述多变量联合涨落的多元高斯分布律（正态分布）：
$P \propto \exp\left[ -\frac{1}{2 k_B T} (\Delta T \Delta S - \Delta p \Delta V) \right]$
(注：由于涨落处于极小区间，为了符号简洁，分母中的环境温度 $T_0$ 直接写作系统平均温度 $T$ 。)

7.2.3 热力学共轭变量的涨落与统计独立性#

上述多元高斯分布中包含了四个涨落变量（ $\Delta T, \Delta S, \Delta p, \Delta V$ ），但对于单组分系统，仅有两个变量是完全独立的。我们可以选择一对相互独立的参量（例如体积 $V$ 和温度 $T$ ）作为独立的随机基底，对分布函数进行变换。
将 $\Delta S$ 和 $\Delta p$ 按照独立变量 $T$ 和 $V$ 进行全微分展开：
$\Delta S = \left( \frac{\partial S}{\partial T} \right)_V \Delta T + \left( \frac{\partial S}{\partial V} \right)_T \Delta V = \frac{C_V}{T} \Delta T + \left( \frac{\partial p}{\partial T} \right)_V \Delta V$ $\Delta p = \left( \frac{\partial p}{\partial T} \right)_V \Delta T + \left( \frac{\partial p}{\partial V} \right)_T \Delta V$
(注：上式 $\Delta S$ 中使用了麦克斯韦关系式 $\left(\frac{\partial S}{\partial V}\right)_T = \left(\frac{\partial p}{\partial T}\right)_V$ )
将这两式代入二次型积分核 $(\Delta T \Delta S - \Delta p \Delta V)$ 中：
$\Delta T \Delta S - \Delta p \Delta V = \Delta T \left[ \frac{C_V}{T} \Delta T + \left( \frac{\partial p}{\partial T} \right)_V \Delta V \right] - \left[ \left( \frac{\partial p}{\partial T} \right)_V \Delta T + \left( \frac{\partial p}{\partial V} \right)_T \Delta V \right] \Delta V$
展开后可以发现，交叉项 $\left(\frac{\partial p}{\partial T}\right)_V \Delta T \Delta V$ 发生极其精确的正负抵消，结果完全解耦为单纯平方项的组合：
$\Delta T \Delta S - \Delta p \Delta V = \frac{C_V}{T} (\Delta T)^2 - \left( \frac{\partial p}{\partial V} \right)_T (\Delta V)^2$
联合分布概率随之严格分解为两个独立的高斯函数的乘积：
$P(\Delta T, \Delta V) \propto \exp\left[ -\frac{C_V}{2 k_B T^2} (\Delta T)^2 \right] \exp\left[ \frac{1}{2 k_B T} \left( \frac{\partial p}{\partial V} \right)_T (\Delta V)^2 \right]$
统计独立性与涨落的定量结论：
概率密度的因式分解从数学上确证了温度涨落 $\Delta T$ 与体积涨落 $\Delta V$ 之间不存在任何统计相关性。根据高斯分布的标准方差公式（ $\int x^2 e^{-ax^2} dx$ 对应的方差为 $1/2a$ ），我们直接获得了如下精确的热力学宏观偏离准则：
1. 统计独立性（交叉关联为零）：
$\langle \Delta T \Delta V \rangle = 0$
(同理可证，压强与熵的涨落也相互独立， $\langle \Delta p \Delta S \rangle = 0$ 。但 $T$ 与 $S$ 、 $p$ 与 $V$ 作为热力学共轭变量则是强烈相关的。)
1. 基础变量的均方涨落：
$\langle (\Delta T)^2 \rangle = \frac{k_B T^2}{C_V}$ $\langle (\Delta V)^2 \rangle = -k_B T \left( \frac{\partial V}{\partial p} \right)_T = k_B T V \kappa_T$
爱因斯坦的准热力学理论巧妙避开了配分函数求解的具体微观细节，以热力学定律的普适约束反推统计力学特性，其得出的结论与 7.1 节及更早期系综理论所求出的涨落方程绝对一致。这充分表明了吉布斯微观统计框架与经典热力学唯象宏观结构的完美自洽性。

第1章 热力学的基本规律#

1.1 热力学状态与第零定律#

1.1.1 热力学平衡状态#

1.1.2 热力学第零定律与温度#

1.1.3 物态方程与热力学系数#

1.2 热力学第一定律与能量转化#

1.2.1 功与内能#

1.2.2 热力学第一定律#

1.2.3 焓 (HHH) 与热容#

1.2.4 理想气体的绝热过程#

1.3 第二定律与不可逆性#

1.3.1 过程的方向性与第二定律的表述#

1.3.2 卡诺循环与卡诺定理#

1.3.3 热力学绝对温标的建立#

1.4 熵与热力学基本方程#

1.4.1 克劳修斯等式与熵函数的导出#

1.4.2 热力学基本方程#

1.4.3 克劳修斯不等式与熵增加原理#

1.5 辅助特性函数的初步引入#

1.5.1 自由能 (FFF) 与等温等容系统的判据#

1.5.2 吉布斯函数 (GGG) 与等温等压系统的判据#

1.5.3 判据的绝对对应关系#

第2章：热力学函数与应用#

2.1 热力学基本函数与特性的数学推演#

2.1.1 勒让德变换与四个基本热力学函数#

2.1.2 麦克斯韦关系#

2.1.3 特性函数及其应用#

2.1.4 广义热容关系#

2.2 获取低温的物理过程与原理#

2.2.1 节流膨胀与焦耳-汤姆孙效应（等焓制冷）#

2.2.2 绝热可逆膨胀（等熵制冷）#

2.2.3 磁制冷：顺磁体的绝热去磁#

2.3 热力学第三定律与负温度#

2.3.1 能斯特定理与绝对零度不可达#

2.3.2 负温度的统计物理本质#

第3章 相变#

3.1 热力学平衡与稳定性判据#

3.2 复相平衡条件与相图#

3.3 一阶相变与克拉珀龙方程（Clapeyron Equation）#

3.4 范德瓦尔斯等温线与麦克斯韦等面积法则#

3.5 连续相变与临界现象#

3.5.1 埃伦费斯特相变分类与二阶相变方程#

3.5.2 朗道二级相变理论#

3.5.3 涨落的影响与临界指数#

第4章 多元系复相平衡与化学平衡#

4.1 多元系的热力学基本方程与吉布斯-杜亥姆方程#

4.2 独立组分数与吉布斯相律 (Gibbs Phase Rule)#

4.3 化学平衡条件与质量作用定律#

4.4 反应热与化学平衡的移动#

第5章 近独立子系组成的统计理论#

5.1 统计物理的基本假设与微观态的量子力学描述#

5.1.1 相空间与微观态的量子化描述#

5.1.2 近独立子系模型与全同粒子分类#

5.1.3 等概率原理（Principle of Equal A Priori Probabilities）#

5.2 近独立子系微观状态数 Ω\OmegaΩ 的计算#

5.2.1 麦克斯韦-玻尔兹曼（MB）统计：定域子系#

5.2.2 费米-狄拉克（FD）统计：非定域费米子#

5.2.3 玻色-爱因斯坦（BE）统计：非定域玻色子#

5.3 极值条件与三种最概然分布的导出#

5.3.1 斯特林公式与连续化近似#

5.3.2 三种统计系统的泛函一阶变分#

5.3.3 统一变分表达式与拉格朗日乘子法#

5.4 拉格朗日乘子的热力学意义与配分函数#

5.4.1 内能的全微分与热功分解#

5.4.2 利用微观变分条件替换热量项#

5.4.3 与热力学基本方程的比对与玻尔兹曼关系的自然导出#

5.4.4 麦克斯韦-玻尔兹曼分布与单粒子配分函数 Z1Z_1Z1​#

5.4.5 宏观热力学量的统计计算框架#

5.5 经典极限条件#

5.6 经典理想气体的统计分布与双原子分子内部自由度#

5.6.1 连续能量积分与平动配分函数#

5.6.2 麦克斯韦速度与速率分布律#

5.6.3 多自由度解耦与双原子分子平动#

5.6.4 谐振子模型与振动配分函数#

5.6.5 刚性转子模型与转动配分函数#

5.7 理想气体熵的吉布斯佯谬与全同性修正#

5.7.1 经典理想气体熵与吉布斯佯谬的数学展现#

5.7.2 全同性因子的引入与萨库尔-特特罗德方程#

5.8 弱简并理想玻色气体和费米气体#

5.8.1 分布函数的一阶渐近展开#

第1章热力学的基本规律#

1.2.3 焓 ( $H$ ) 与热容#

1.5.1 自由能 ( $F$ ) 与等温等容系统的判据#

1.5.2 吉布斯函数 ( $G$ ) 与等温等压系统的判据#

第3章相变#

第4章多元系复相平衡与化学平衡#

第5章近独立子系组成的统计理论#

5.2 近独立子系微观状态数 $\Omega$ 的计算#

5.4.4 麦克斯韦-玻尔兹曼分布与单粒子配分函数 $Z_1$ #

5.10.1 绝对零度下的费米海与费米能量 ( $T = 0\text{ K}$ )#

5.10.2 费米-狄拉克积分与索末菲展开原理 ( $T \ll T_F$ )#

第6章系综理论#

6.1 $\Gamma$ 空间与刘维尔定理#

6.1.1 $\Gamma$ 空间的概念建立#

6.4.2 配分函数 $Z_N$ 与亥姆霍兹自由能 $F$ #

6.5.2 巨配分函数 $\Xi$ 与巨热力学势 $J$ #

第七章涨落理论#