刘维尔方程 Liouville Equation

刘维尔方程（Liouville equation）描述相空间中概率密度沿哈密顿流的演化，是经典统计力学的基石方程，也是连接确定性动力学与统计描述的桥梁。

方程表述

$N$ 个自由度的哈密顿系统，状态为相空间点 $(\mathbf{q},\mathbf{p})\in {\mathbb{R}}^{2N}$，概率密度 $\rho (\mathbf{q},\mathbf{p},t)$ 满足：

$$\frac{\partial \rho }{\partial t}+\sum _{i=1}^N\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)=0$$

用 Poisson 括号压缩为：

$$\frac{\partial \rho }{\partial t}+\{\rho ,H\}=0,\{A,B\}=\sum _i\left(\frac{\partial A}{\partial q_i}\frac{\partial B}{\partial p_i}-\frac{\partial A}{\partial p_i}\frac{\partial B}{\partial q_i}\right)$$

推导

相空间点遵从哈密顿方程：

$${\dot{q}}_i=\frac{\partial H}{\partial p_i},{\dot{p}}_i=-\frac{\partial H}{\partial q_i}$$

把相空间看作 $2N$ 维流体，密度守恒给出连续性方程：

$$\frac{\partial \rho }{\partial t}+{\nabla }_{q,p}\cdot (\rho \mathbf{v})=0,\mathbf{v}=(\dot{\mathbf{q}},\dot{\mathbf{p}})$$

关键在哈密顿流无散：

$${\nabla }_{q,p}\cdot \mathbf{v}=\sum _i\left(\frac{{\partial }^{2}H}{\partial q_i\partial p_i}-\frac{{\partial }^{2}H}{\partial p_i\partial q_i}\right)=0$$

于是 $\nabla \cdot (\rho \mathbf{v})=\mathbf{v}\cdot \nabla \rho$，连续性方程退化为刘维尔方程。这一步本质是 散度定理 语言下的相空间不可压缩性。

刘维尔定理

沿任意轨迹密度的物质导数为零：

$$\frac{d\rho }{dt}=\frac{\partial \rho }{\partial t}+\{\rho ,H\}=0$$

即密度沿轨迹守恒，等价表述为相体积守恒：哈密顿流把相空间区域映射到同体积的区域（不可压缩流）。这与耗散系统形成对照——吸引子上的体积收缩，不可能用哈密顿力学描述。

表述	含义
$d\rho /dt=0$	系综点随流携带各自的密度
相体积不变	辛结构保体积（Poincaré 不变量）
平衡解 $\{\rho ,H\}=0$	$\rho =f(H)$ 皆为定态，微正则系综取 $f$ 为常数

意义与应用

领域	角色
统计力学	系综演化的基本方程；平衡态统计的出发点
遍历理论	Poincaré 回复、混合性（mixing）等概念的载体
等离子体/天体力学	平均场近似下给出 Vlasov 方程（无碰撞 Boltzmann 方程）
量子力学对应	经典极限下 von Neumann 方程 $i\hslash {\partial }_t\hat{\rho }=[\hat{H},\hat{\rho }]$，Poisson 括号对应 $(i\hslash {)}^{-1}[\cdot ,\hat{H}]$
生成流模型	Continuous Normalizing Flow 的对数密度 ODE（瞬时 change-of-variables）即其推论

消歧

• 复分析的 Liouville 定理（有界整函数必为常数）是另一结果，仅共享人名，与方程无关。
• 几何中的二维 Liouville 方程 $\Delta u=K{e}^u$（曲面度量与曲率）是另一独立对象。

相关笔记

• 常微分方程 — 相空间轨迹的语言基础
• 散度定理 — 连续性方程的积分形式
• 变分法 — 哈密顿力学的变分起源