相空间概念与刘维尔定理
在经典力学中，我们通常关注单个质点的运动轨迹。然而，当研究包含 $N$ 个粒子的复杂系统（如气体或液体）时，试图追踪每一个粒子的位置和动量在计算上是极其困难且不切实际的。为了从整体上描述系统的演化，统计力学引入了“相空间”这一几何概念，并通过“刘维尔定理”揭示了系统在相空间中演化的深刻规律。
在经典力学框架下，一个系统的完整状态由所有粒子的位置坐标和动量坐标共同决定。
1. 微观状态与坐标表示
假设一个系统由 $N$ 个粒子组成，每个粒子在三维空间中运动。为了完全描述该系统在某一时刻的状态，我们需要知道：
- 每个粒子的位置坐标：$\mathbf{q}_1, \mathbf{q}_2, \dots, \mathbf{q}_N$（共 $3N$ 个分量）
- 每个粒子的动量坐标：$\mathbf{p}_1, \mathbf{p}_2, \dots, \mathbf{p}_N$（共 $3N$ 个分量）
将这 $6N$ 个变量组合在一起，构成了一个高维空间中的一个点。这个点被称为系统的微观状态（Microstate）。
2. 相空间的定义
所谓的相空间（Phase Space），就是由所有可能的微观状态构成的数学空间。对于一个 $N$ 粒子的系统，其相空间是一个 $6N$ 维的流形。
- 维度：$6N$ 维。
- 点：相空间中的每一个点 $\Gamma = (\mathbf{q}_1, \dots, \mathbf{q}_N, \mathbf{p}_1, \dots, \mathbf{p}_N)$ 都唯一对应系统的一个瞬时状态。
- 轨迹：随着时间的推移，由于系统遵循哈密顿力学规律，相空间中的点会沿着特定的路径运动，这条路径被称为相轨迹（Phase Trajectory）。
通过将复杂的粒子动力学问题转化为相空间中点的运动问题，我们可以利用几何工具来研究系统的统计性质。
相空间中的动力学演化
系统的演化由哈密顿方程（Hamilton's Equations）严格支配：
$$\dot{q}_i = \frac{\partial H}{\partial p_i}, \quad \dot{p}_i = -\frac{\partial H}{\partial q_i}$$
其中 $H(\mathbf{q}, \mathbf{p}, t)$ 是系统的哈密顿量（通常为总能量）。
在相空间中，这些方程定义了一个“速度场”。每一个点 $\Gamma$ 都有一个确定的速度向量 $\mathbf{v} = (\dot{\mathbf{q}}, \dot{\mathbf{p}})$。这意味着，如果我们有一群具有不同初始条件的系统（即一个“系综”），这些系统在相空间中就像是一团流动的“流体”。
刘维尔定理 (Liouville's Theorem)
刘维尔定理是统计力学中最核心的定理之一，它描述了相空间中“概率密度”随时间演化的特性。
1. 概率密度函数
在统计力学中，由于我们无法精确获知系统的确切状态，我们使用一个概率密度函数 $\rho(\mathbf{q}, \mathbf{p}, t)$ 来描述系统处于相空间某区域的可能性。$\rho$ 的物理意义是：在相空间体积元 $d\Gamma$ 内找到系统的概率为 $\rho d\Gamma$。
2. 定理的内容
刘维尔定理指出：在相空间中，随着系统沿着相轨迹运动，概率密度 $\rho$ 保持不变。
用数学语言表达，如果我们将 $\rho$ 视为随时间变化的标量场，那么其全导数（或称物质导数）为零：
$$\frac{d\rho}{dt} = \frac{\partial \rho}{\partial t} + \sum_{i=1}^{3N} \left( \frac{\partial \rho}{\partial q_i} \dot{q}_i + \frac{\partial \rho}{\partial p_i} \dot{p}_i \right) = 0$$
3. 物理直观理解：不可压缩性
一个更直观的理解方式是：相空间中的“概率流体”是不可压缩的。
想象在相空间中有一个由许多点组成的区域（代表一组可能的初始状态）。随着时间的演化，这个区域可能会发生剧烈的变形、拉伸或扭曲，但这个区域在相空间中所占据的体积始终保持不变。
我们可以通过连续性方程来证明这一点。在相空间中，概率守恒要求满足：
$$\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$$
其中 $\nabla \cdot$ 是相空间中的散度。展开散度项：
$$\frac{\partial \rho}{\partial t} + \mathbf{v} \cdot \nabla \rho + \rho (\nabla \cdot \mathbf{v}) = 0$$
注意到 $\frac{\partial \rho}{\partial t} + \mathbf{v} \cdot \nabla \rho = \frac{d\rho}{dt}$。因此，只要证明相空间速度场的散度 $\nabla \cdot \mathbf{v} = 0$，刘维尔定理就成立。
利用哈密顿方程计算散度：
$$\nabla \cdot \mathbf{v} = \sum_{i=1}^{3N} \left( \frac{\partial \dot{q}_i}{\partial q_i} + \frac{\partial \dot{p}i}{\partial p_i} \right) = \sum{i=1}^{3N} \left( \frac{\partial^2 H}{\partial q_i \partial p_i} - \frac{\partial^2 H}{\partial p_i \partial q_i} \right) = 0$$
由于二阶偏导数满足交换律，散度恒为零。这证明了相空间流体的不可压缩性。
示例与应用
示例：单摆的相空间
考虑一个简谐振子（小角度摆动）。其坐标为 $q$（角度），动量为 $p$（动量）。
- 其哈密顿量为 $H = \frac{p^2}{2m} + \frac{1}{2}kq^2$。
- 在相空间 $(q, p)$ 中，能量等值线表现为一系列椭圆。
- 随着时间演化，相空间中的一个点会沿着椭圆轨迹运动。
- 如果你在相空间中画一个小的圆圈，随着摆动，这个圆圈会沿着椭圆路径移动，虽然它可能会因为非线性效应变成长条状，但其面积（即相空间体积）始终保持恒定。
统计力学的意义
刘维尔定理为系综理论（Ensemble Theory）奠定了基础：
- 平衡态的定义：如果一个系统的概率分布 $\rho$ 不随时间变化（即 $\frac{\partial \rho}{\partial t} = 0$），那么根据刘维尔定理，该分布必须满足 $\sum (\frac{\partial \rho}{\partial q_i} \dot{q}_i + \frac{\partial \rho}{\partial p_i} \dot{p}_i) = 0$。这为寻找平衡态分布（如正则系综中的玻尔兹曼分布）提供了约束条件。
- 信息守恒：相空间体积不变意味着系统演化过程是可逆的，微观层面的信息在动力学演化过程中既不会增加也不会减少。
- 熵的讨论：虽然相空间体积不变，但在宏观观测中，系统往往表现出向平衡态趋近的过程。这引出了关于“粗粒化”（Coarse-graining）和统计熵增加的深刻讨论，即宏观上的熵增实际上是由于相空间体积虽然守恒，但分布变得极其破碎、复杂，导致宏观观测无法分辨细微结构的结果。