探索经典力学、量子力学与统计物理中的核心纽带 —— 位力定理 (Virial Theorem)。从牛顿引力到星系动力学，从原子结构到量子色动力学，理解能量在动能与势能之间的深刻平衡。
在物理学中，位力定理（Virial Theorem，又译维里定理）是一个普遍成立的定理，它描述了在一个稳定的、自引力束缚系统或受幂律势场束缚的系统中，平均动能（Kinetic Energy）与平均势能（Potential Energy）之间的关系。该定理由鲁道夫·克劳修斯（Rudolf Clausius）于1870年首次提出，最初用于气体动力学理论，后来被广泛应用于天体物理学、量子力学和统计力学领域。
对于一个受幂律势 V(r) ∝ rⁿ 束缚的系统，位力定理指出：
2⟨T⟩ = n⟨V⟩
其中 ⟨T⟩ 是时间平均动能，⟨V⟩ 是时间平均势能。
在引力系统中，势能遵循平方反比律，即 n = -1。此时定理简化为：
2⟨T⟩ = -⟨V⟩
这意味着系统的总能量 E = ⟨T⟩ + ⟨V⟩ = -⟨T⟩ = ½⟨V⟩。负号表示系统处于束缚态。
在量子力学中，位力定理以算符期望值的形式存在。对于哈密顿量 Ĥ = T̂ + V̂，若势能是坐标的 k 次齐次函数，则：
2⟨T⟩ = k⟨V⟩
这一形式在计算原子和分子基态能量时极为重要。
为了深入理解 位力定理，我们需要从不同的物理框架出发进行推导。以下是三种最常见的证明路径：经典力学、统计力学和量子力学。
考虑一个由 N 个粒子组成的系统，第 i 个粒子的质量为 mᵢ，位置为 rᵢ，受力为 Fᵢ。定义位力（Virial）G 为：
G = Σᵢ (pᵢ · rᵢ)
其中 pᵢ 是动量。对时间求导：
dG/dt = Σᵢ (dpᵢ/dt · rᵢ + pᵢ · drᵢ/dt)
根据牛顿第二定律 dpᵢ/dt = Fᵢ，且 pᵢ/mᵢ = vᵢ = drᵢ/dt，故 pᵢ · drᵢ/dt = mᵢvᵢ² = 2Tᵢ（Tᵢ 为动能）。因此：
dG/dt = Σᵢ (Fᵢ · rᵢ) + 2T
对于有界运动（束缚系统），G 是有界的。在一个长时间内取平均，dG/dt 的平均值为零（因为 G 不会无限增长）。因此：
2⟨T⟩ = -⟨Σᵢ (Fᵢ · rᵢ)⟩
如果力来自势能 V，且 V 是 k 次齐次函数（即 V(λr) = λᵏV(r)），根据欧拉齐次函数定理，Σᵢ (rᵢ · ∇ᵢV) = kV。由于 F = -∇V，则 Σᵢ (Fᵢ · rᵢ) = -kV。代入上式得：
2⟨T⟩ = k⟨V⟩
这就完成了经典力学中 位力定理 的证明。
在量子力学中，我们考虑哈密顿量 Ĥ = T̂ + V̂。对于定态（能量本征态），任何算符 A 的期望值不随时间变化，即 d⟨Â⟩/dt = 0。根据海森堡运动方程：
d⟨Ĝ⟩/dt = (i/ħ)⟨[Ĥ, Ĝ]⟩ = 0
选择 Ĝ = r · p（位置与动量的点积算符）。我们需要计算对易子 [Ĥ, r·p]。由于 [T̂, r·p] 和 [V̂, r·p] 分别贡献动能和势能项，经过详细的算符运算（利用 [xᵢ, pⱼ] = iħδᵢⱼ），可以得到：
[T̂, r·p] = -2iħT̂
对于势能项，如果 V(r) 是 k 次齐次函数，则：
[V̂, r·p] = iħkV̂
将对易子代入运动方程：
(i/ħ)⟨-2iħT̂ + iħkV̂⟩ = 0 ⇒ -2⟨T̂⟩ + k⟨V̂⟩ = 0
即：
2⟨T⟩ = k⟨V⟩
这证明了在量子束缚态中，位力定理 同样成立。
在统计力学中，我们可以利用相空间积分来证明 位力定理。考虑正则系综，系统处于能量 E 的概率分布为 ρ ∝ exp(-βH)。我们需要计算积分：
⟨xᵢ ∂H/∂xⱼ⟩ = ∫ xᵢ (∂H/∂xⱼ) e^(-βH) dΓ
利用分部积分法（假设边界项为零，因为玻尔兹曼因子在无穷远处衰减极快）：
∫ xᵢ ∂(e^(-βH))/∂xⱼ dΓ = - ∫ (∂xᵢ/∂xⱼ) e^(-βH) dΓ
右边积分结果为 -δᵢⱼ Z（Z 为配分函数），归一化后得到：
⟨xᵢ ∂H/∂xⱼ⟩ = δᵢⱼ k_B T
对于动能 T = Σ p²/2m，取 xᵢ = pᵢ，则 ∂H/∂pᵢ = pᵢ/m。令 i=j 并求和：
⟨pᵢ · (pᵢ/m)⟩ = 3k_B T (对于三维粒子)
由于 p²/2m = T，故 2⟨T⟩ = 3Nk_B T。结合势能项的类似推导，即可得到广义的 位力定理 关系。
位力定理 不仅是理论推导的工具，更是解决实际问题的一把钥匙。以下是其在不同学科中的具体应用案例。
这是 位力定理 最著名的应用之一。1933年，弗里茨·兹威基（Fritz Zwicky）研究后发座星系团（Coma Cluster）时，测量了其中星系的径向速度弥散度 σ。根据位力定理：
M ≈ (R σ²) / G
其中 M 是星系团总质量，R 是星系团的特征半径，G 是引力常数。计算结果显示，所需的质量比可见物质（星系和气体）的质量大了约 400 倍。这一巨大差异直接导致了暗物质假说的诞生。
在恒星内部，位力定理 解释了恒星的热平衡机制。对于理想气体恒星，动能对应于热运动能量（内能 U），势能对应于引力势能 Ω。由 2T + Ω = 0 可得：
总能量 E = T + Ω = -T = -½Ω。
这意味着当恒星因辐射损失能量（E 减小，变得更负）时，其动能 T 必须增加。换句话说，恒星在收缩时会变热。这一机制是恒星主序前阶段（如赫比格-罗天体）收缩发光的主要能源，也是触发核聚变的必要条件。
在计算化学中，位力定理 是检验波函数质量的重要准则。对于原子和分子，如果势能是库仑势（k=-1），则必须满足 2⟨T⟩ = -⟨V⟩。如果在哈特里-福克（Hartree-Fock）计算中此关系不满足，说明基组（basis set）不完备或计算未收敛。它提供了独立于能量数值之外的一个验证标准。
在受控核聚变装置（如托卡马克）中，位力定理 用于分析磁约束等离子体的稳定性。等离子体压力与磁压强之间的平衡关系可以通过位力原理进行推导，帮助物理学家预测等离子体是否会发生磁流体不稳定性（MHD instabilities）。
位力定理 的发展跨越了两个世纪，见证了物理学从经典到量子的范式转移。
克劳修斯在研究气体动力学理论时，首次引入了“位力”（Virial，源自拉丁语 vis，意为力）的概念，并证明了气体压强与分子动能及分子间作用力位力之间的关系。这是 位力定理 的起源。
弗里茨·兹威基将 位力定理 应用于后发座星系团，发现了质量-光度异常，间接预言了暗物质的存在，这是该定理在天体物理学中的里程碑式应用。
随着海森堡和薛定谔建立量子力学，位力定理 被推广到量子领域。泡利（Pauli）和费曼（Feynman）等人利用量子位力定理分析了原子结构和分子键合性质。
在计算机模拟大规模结构形成（N-body simulations）中，位力定理 被用于初始化模型和验证数值算法的准确性。它也是理解星系团合并过程中能量重新分布的关键工具。
以下是关于 位力定理 证明及应用中最常见的问题解答。
不是。位力定理主要适用于有界系统（bounded systems）或稳定系统。对于无界系统（如自由粒子散射）或非稳态系统（如正在剧烈坍缩或膨胀的系统），定理的标准形式可能不直接适用，或者需要引入额外的时间依赖项。此外，系统必须满足势能是坐标的齐次函数这一条件。
在引力系统中，n=-1，故 2⟨T⟩ = -⟨V⟩。总能量 E = ⟨T⟩ + ⟨V⟩ = -½⟨V⟩ + ⟨V⟩ = ½⟨V⟩。由于引力势能 V 定义为负值（无穷远处为零），故 ⟨V⟩ < 0，导致 E < 0。负的总能量是系统处于束缚态的标志，意味着需要输入能量才能使系统解体。
有用。在 QCD 中，胶子和夸克的相互作用势能近似为线性势（V ∝ r，即 n=1）加库仑势。位力定理帮助物理学家理解强子（如质子、中子）的质量来源。值得注意的是，强子的绝大部分质量并非来自夸克的静止质量，而是来自夸克的动能和胶子场的势能，这正是位力定理在微观尺度上的体现。
在分子动力学（MD）模拟中，可以通过长时间积分粒子的轨迹，计算平均动能 ⟨T⟩ 和平均势能 ⟨V⟩，然后检查 2⟨T⟩/⟨V⟩ 是否等于势能函数的齐次次数 n。如果偏差较大，通常意味着模拟步长过大或积分算法不稳定。