位力定理(,又稱维里定理、均功定理)是力學中描述穩定的多自由度孤立體系的總動能和總勢能時間平均之間的數學關係。考慮一個有N個質點的體系,其數學表達式爲:
:\langle T \rangle = - \frac 1 2 \sum_{k=1}^N \langle \boldsymbol F_k \cdot \boldsymbol r_k \rangle
其中:角括號表示對時間取平均,T是系统内部的总动能,\boldsymbol F_k是第k個質點所受的力,\boldsymbol r_k是第k個質點的位置向量;等式右邊稱作均位力積(,簡稱位力),反映體系內相互作用強度。英語virial一詞由德國物理學家魯道夫·克勞修斯於1870年根據拉丁語單詞vīs(意爲力、能量)命名。
特別地,若系統内任何粒子兩兩之間的力來自與粒子間距離r的n次冪成正比的勢能V(r) = \alpha r^n(其中\alpha ,n為常數),則定理簡化為:
:2 \langle T \rangle = n \langle V_\text{Total} \rangle
即:體系的總動能2倍等於總勢能的n倍。對於引力勢能,這裏的n = -1。
位力定理的一個意義在於,它允許計算平均總動能,即便是對於那些無法精確解的非常複雜的系統,例如在統計力學中考慮的那些;根據能量均分定理,該平均總動能與系統溫度有關。然而,維里定理不依賴於溫度的概念,甚至適用於不處於熱平衡的系統。維里定理已以各種方式推廣,特別是張量形式。
歷史
命題推導
簡單例子
考慮個質量相同的質點構成的孤立體系,它們受萬有引力相互作用。假設兩個質點分別以和的速度(大小均為,方向相反)圍繞共同質心做匀速圆周运动,半徑為,兩者分別受到作用力和(大小均爲,方向相反)。則體系的時間平均縂動能為:
:\langle T \rangle = \sum_{k=1}^N \frac12 m_k \left|\mathbf{v}_k \right|^2 = \frac12 m|\mathbf{v}_1|^2 + \frac12 m|\mathbf{v}_2|^2 = mv^2
以共同質心為原點,兩者的位置向量分別爲及(大小均爲常數)。引力方向朝向原點,與位置向量方向相反,故。又,向心力大小等於萬有引力大小:。代入得:
:-\frac12 \sum_{k=1}^N \bigl\langle \mathbf{F}_k \cdot \mathbf{r}_k \bigr\rangle = -\frac12(-Fr - Fr) = Fr = \frac{mv^2}{r} \cdot r = mv^2 = \langle T \rangle
一般推導
預先知識
對於 Virial theorem 的推導, 將需要用到齐次函数的如下性質,
既當 f(\vec{x}) 為 k 次 齊次函數時, 有:
\frac{d f(\alpha \vec{x})}{d(\alpha \vec{x})} \cdot \vec{x} = \frac{\partial}{\partial{\alpha}}[f(\alpha \vec{x})]
= \frac{\partial[\alpha^k f(\vec{x})]}{\partial \alpha} = k \alpha^{k-1} f(\vec{x})
對於a=1時有:
\vec{x} \vec{ \nabla}f(\vec{x} )= k f(\vec{x})
具體推導
注意到動能T是一個關於速度\vec{v}的2次齊次函數:
\vec{v}\cdot \frac{\partial T}{\partial \vec{v}} = 2T , 同時有 \frac{\partial{T}}{\partial\vec{v}} = \vec{p}, 從而得到
\ 2 T = \vec{v} \cdot \vec{p} = \frac{d}{d t}(\vec{p} \cdot \vec{x}) - \dot{\vec{p}} \cdot \vec{x}
計算上式對於時間\Delta t的平均:
_{\Delta t} = \frac{1}{\Delta t} \int_{0}^{\Delta t}(...)\text{( 等 式 右 邊 )}\ dt
我們關注\Delta t \rightarrow \infty 的情況, 假設系統的運動是有限的 ((\vec{p} \cdot \vec{x})不會有\infty出現的情況), 此時等式右邊的前半部分將趨近於0:
lim_{\Delta t \rightarrow \infty} \frac{1}{\Delta t} \int_{0}^{\Delta t}\frac{d}{dt}(\vec{p} \cdot \vec{x})dt
= lim_{\Delta t \rightarrow \infty} \frac{(\vec{p}\cdot \vec{x})|_{(\Delta t)} - (\vec{p} \cdot \vec{x})|_{(0)}}{\Delta t} \rightarrow 0
我們得到:
2_{t} = -_t
\dot{\vec{p}}可以通過系統的勢能 V(\vec{x}) 寫出: \dot{\vec{p}} = - \vec{\nabla}V(\vec{x}); 另外我們最終假設勢能 V(\vec{x}) 為,k次齊次函數, 並利用預先知識中a=1時的等式 就能夠得到位力定理:
2 _t = -_t = _t = k_t
與質點間勢能之關聯
對於冪定律力
關於時間平均
一般化
引入電磁場
相對論均匀系統
各學科中的應用
量子力學
狹義相對論
天體物理學
位力質量、位力半徑
太阳模型
\frac{dP}{dr} = -\frac{Gm(r)\rho(r)}{r^2} \Rightarrow dP=-\frac{Gm(r)\rho(r)}{r^{2}}dr.
考虑恒星的位力定理在天体物理学中的应用。如果将恒星的物质视为流体,则可以使用流体静力平衡条件来考虑恒星。这个假设条件允许将恒星的引力与其内部的压力建立关系,从而将引力势能与内能联系起来,即位力定理。基于引力势能 \Omega \propto r^{-1},我们期待内能与势能之间的关系为 E_{\text{int}} = -\frac{1}{2}\Omega。
下面是更详细的推导过程:将恒星视为正球体来简化推导过程。气压 P 是半径 r 的函数:
\frac{dP}{dr} = -\frac{Gm(r)\rho(r)}{r^2} \Rightarrow dP=-\frac{Gm(r)\rho(r)}{r^{2}}dr.
对于理想气体,内能 E_{\text{int}} 为:
E_{\text{int}} = \frac{3}{2}Nk_{B}T = \frac{3}{2}\frac{M}{m_g}k_{B}T,
其中 m_g 是粒子的平均质量,M 是恒星的质量,k_{B} 是玻尔兹曼常数,T 是温度。
考虑恒星的静力平衡条件,同时乘以体积 V=\frac{4\pi r^3}{3},并积分,得到:
\int_0^R Vdp = -\frac{1}{3}\int_0^R \frac{Gm(r)}{r}\rho(r)4\pi r^2 dr = -\frac{1}{3}\int_0^R \frac{Gm(r)}{r}dm(r).
右侧积分包含了恒星的重力势能 \Omega = -\int_0^R \frac{Gm(r)}{r}dm(r),所以我们可以得到:
\int_0^R Vdp = \frac{1}{3}\Omega.
对左侧积分使用分部积分可得:
\int_0^R VdP = PV|_0^R - \int_0^R PdV = -\int_0^R PdV,
其中 PV|_0^R = 0,因为恒星最外层压强为0,最内侧体积为0。对于理想气体,PV = Nk_{b}T,将其与理想气体内能公式结合,得到单位体积内的内能:
U = \frac{E_{\text{int}}}{V} = \frac{3}{2}P \Rightarrow P = \frac{2}{3} U,
将其应用到上面的积分,得到:
-\int_0^R PdV = -\frac{2}{3} \int_0^R UdV = -\frac{2}{3}E_{\text{int}}.
将两侧积分结果相等,得到:
-\frac{2}{3}E_{\text{int}} = \frac{1}{3}\Omega \Rightarrow E_{\text{int}} = -\frac{1}{2} \Omega.
这就是恒星在流体静力平衡下的位力定理。
通过这个公式,可以推算太阳的平均温度 \langle T_{\odot} \rangle 大约为 10^6 开尔文,对应的内能大约为 1 \text{ keV}。太阳的表面温度仅在 5774 开尔文左右,因此可以认为太阳的内部温度比表面高很多。由于电子的结合能仅为 1 \text{ eV},而太阳的平均内能远大于这个数值,因此可以认为太阳是离子气体。
統計物理
在統計物理中,:
:\mathbf{p}=\frac 2 V \left ( \frac 1 2 \left \langle \sum_{i=1}^N m_i \mathbf{v}_i \otimes \mathbf{v}_i \right \rangle + \frac 1 2 \left \langle \sum_{i
亦即體系壓強爲(與動能相關的)動理壓強和(與相互作用相關的)內壓強之和。上式中的第二項即爲均位力積相關項。
引用
评论 (0)