能量均分

在古典統計力學中，**能量均分定理**（Equipartition Theorem）是一種聯繫系統溫度及其平均能量的基本公式。能量均分定理又被稱作**能量均分定律**、**能量均分原理**、**能量均分**，或僅稱**均分**。能量均分的初始概念是熱平衡時能量被等量分到各種形式的運動中；例如，一個分子在平移運動時的平均動能應等於其做旋轉運動時的平均動能。

能量均分定理能夠作出定量預測。類似於均功定理，對於一個給定溫度的系統，利用均分定理，可以計算出系統的總平均動能及位能，從而得出系統的熱容。均分定理還能分別給出能量各個組分的平均值，如某特定粒子的動能又或是一個彈簧的位能。例如，它預測出在熱平衡時理想氣體中的每個粒子平均動能皆為 _(3/2)kBT_，其中 _kB_ 為波耳茲曼常數而 _T_ 為溫度。更普遍地，無論多複雜也好，它都能被應用於任何處於熱平衡的古典系統中。能量均分定理可用於推導古典理想氣體定律，以及固體比熱的杜隆 - 泊替定律。它亦能夠應用於預測恆星的性質，因為即使考慮相對論效應的影響，該定理依然成立。

儘管均分定理在一定條件下能夠對物理現象提供非常準確的預測，但是當量子效應變得顯著時（如在足夠低的溫度條件下），基於這一定理的預測就變得不準確。具體來說，當熱能 _kBT_ 比特定自由度下的量子能級間隔要小的時候，該自由度下的平均能量及熱容比均分定理預測的值要小。當熱能比能級間隔小得多時，這樣的一個自由度就說成是被「凍結」了。比方說，在低溫時很多種類的運動都被凍結，因此固體在低溫時的熱容會下降，而不像均分定理原測的一般保持恆定。對十九世紀的物理學家而言，這種熱容下降現象是表明古典物理學不再正確，而需要新的物理學的第一個徵兆。均分定理在預測電磁波的失敗（被稱為「紫外災變」）導致普朗克提出了光本身被量子化而成為光子，而這一革命性的理論對刺激量子力學及量子場論的發展起到了重要作用。

## 基本概念及簡易例子

應用波爾茲曼統計方法可以得到：氣體處於平衡態時，分子任何一個自由度的平均能量都相等，均為 kT/2, 這就是能量按自由度均分定理，簡稱能量均分定理。名字裏面的「均分」是指「攤分或類似於攤分」。能量均分定理的原始概念是，當系統達到熱平衡時，系統的總動能由各獨立分量所等分。均分定理也為這些能量做出量化的預測。例如它預測惰性氣體的每一個原子，當於溫度 _T_ 達至熱平衡時，會有平移平均動能 _(3/2)KBT_，其中 KB 為波茲曼常數。隨此引出的是，在等溫時氙的重原子速度會比氦的較輕原子要低。圖二顯示的是四種惰性氣體原子速度的馬克士威 - 波茲曼分布。

在這例子中，關鍵點是動能是速度的二次齊函數。均分定理顯示出於熱平衡時，任何在能量中只以二次出現的自由度（例如是一粒子的位置或速度的一個分量）有著等於 _½KBT_ 的平均能量，並因此向系統的熱容提供了 ½KB。這個結果有著許多的應用。

### 平移能量與理想氣體

一粒子質量為 _m_，速度為 _v_，其（牛頓力學）動能為：

　
: $H^{\mathrm {kin} }={\tfrac {1}{2}}m|\mathbf {v} |^{2}={\tfrac {1}{2}}m\left(v_{x}^{2}+v_{y}^{2}+v_{z}^{2}\right),$

其中 _vx_、_vy_ 及 _vz_ 是速度 _v_ 的直角坐標的分量。這裏，_H_ 是哈密頓量，由於哈密頓表述是均分定理一般形式的中心，故下文將以其作為能量的符號。

由於能量是速度各分量的二次方，均分這三分量得每分量在熱平衡時向平均動能提供 ½_kBT_。因此粒子的平均動能為 (3/2)_kB_T，跟上面惰性氣體的例子一樣。

更普遍地，理想氣體中的，總能量幾乎全為（平移）動能：假定粒子無內自由度且運動不受其他粒子影響。均分因此預測有 _N_ 個粒子的理想氣體有平均總能量 (3/2) _N kBT_。

而氣體的熱容則為 (3/2) _N kB_，因此這樣一莫耳氣體的熱容為 (3/2)_NAkB_=(3/2)_R_，其中 _NA_ 是亞佛加厥常數，而 _R_ 則是氣體常數。由於 _R_ ≈ 2 Cal/(mol・K)，均分預測理想氣體的莫耳比熱容約為 3 Cal/(mol・K)。這個預測已被實驗證實。

從平均動能可以求出氣體粒子的均方根速度 _vrms_：

　
: $v_{\mathrm {rms} }={\sqrt {\langle v^{2}\rangle }}={\sqrt {\frac {3k_{B}T}{m}}}={\sqrt {\frac {3RT}{M}}},$

其中 _M = NAm_ 是一莫耳氣體粒子的質量。這個結果對很多應用方面都有用處，例如逸散用的格銳目定律為鈾濃縮提供了一個方法。

### 旋轉能量與溶液中的分子滾翻

在另一個相近的例子中，有一粒子其主轉動慣量 _I1_、_I2_ 及 _I3_。它的旋轉能量是：

　
: $H^{\mathrm {rot} }={\tfrac {1}{2}}(I_{1}\omega _{1}^{2}+I_{2}\omega _{2}^{2}+I_{3}\omega _{3}^{2}),$

其中 _ω1_、_ω2_ 及 _ω3_ 是角速度的主分量。使用跟平移同一套的論證，均分意味著每個粒子的平均旋轉能量為 _(3/2)KBT_。同樣地，均分使計算出分子均方根角速度成為可能。

剛性粒子的滾翻 —— 即是分子於溶液中的隨機旋轉 —— 在核磁共振中觀測到弛緩中有著重要的角色，尤其是在蛋白質核磁共振及剩餘雙極耦合中。旋轉滲透可被其他生物物理探測法所觀測到，例如是螢光異向性、流動雙折射及介電質光譜學。

### 位能與諧波振蕩器

均分定理除可應用於動能外，還能被應用於位能計算：重要例子包括像彈簧這樣的諧波振蕩器，其二次位能為

　
: $H^{\mathrm {pot} }={\tfrac {1}{2}}aq^{2},\,$

其中常數 _a_ 描述彈簧的韌性，而 _q_ 則是由平衡導出的。假若這樣一個系統的質量為 _m_，那麼它的動能 _H_kin 為 ½_mv2_=_p_2/2_m_，其中 _v_ 及 _p_=_mv_ 代表振蕩器的速度和動量。聯合這些項可得總能量：

　
: $H=H^{\mathrm {kin} }+H^{\mathrm {pot} }={\frac {p^{2}}{2m}}+{\frac {1}{2}}aq^{2}.$

因此均分定理預測在熱平衡時，振蕩器有平均能量

　
: $\langle H\rangle =\langle H^{\mathrm {kin} }\rangle +\langle H^{\mathrm {pot} }\rangle ={\tfrac {1}{2}}k_{B}T+{\tfrac {1}{2}}k_{B}T=k_{B}T,$

其中角括號 $\left\langle \ldots \right\rangle $ 代表括號內的平均量 。

這個結果對任何種類的諧波振蕩器都是有效的，例如鐘擺，一個振動中的粒子或是被動的電子振蕩器。這樣的振蕩器在很多情況下都會出現；由均分可得，每個這樣的振蕩器都得到一個平均總能量 _kBT_ 並因此向系統熱容提供 _kB_。這個可以被用於導出熱雜音的公式 及固體莫耳比熱容的杜隆 - 泊替定律公式。後者在均分定理的歷史中尤其重要。

### 固體的比熱容

均分定理的一個重要應用是在於晶狀固體的比熱容。如此固體的每一個原子都能夠在三個獨立的方向下振蕩，因此該固體可以被視為一個擁有各自獨立的 _3N_ 個簡諧振子的系統，其中 _N_ 為晶格中的原子數。由於每一個諧振子都有平均能量 _kBT_，所以固體的平均總能量為 _3NkBT_，而比熱容則為 _3NkB_。如取 _N_ 為亞佛加厥常數 _NA_，並使用 _R = NAkB_ 這個聯繫氣體常數 _R_ 及波茲曼常數 _kB_ 的關係式，可得固體莫耳比熱容的杜隆 - 泊替定律的一個解釋，定律指出晶格中每莫耳的原子熱容為 _3R_ ≈ 6 cal/(mol·K)。

然而，由於量子效應的關係，這條定律在低溫時並不準確；這也不符合實驗導出的熱力學第三定律，第三定律指出莫耳比熱容於絕對零度時必為零。艾爾伯特・愛因斯坦（1907 年）及彼得・德拜（1911 年）在基礎上加入了量子效應，發展出一套更準確的理論。

固體中每個原子的振動不是獨立的，可以用一組組的耦合振子作為模型。如此振蕩子的模型可以被分解成簡振模，這跟鋼琴弦的振動模態及管風琴的共振模態是相近的。另一方面，均分定理被應用於這種系統時一般都會失敗，因為正常模態間是沒有能量交換的。在一個非常的情況下，模態獨立且它們的能量獨立地守恆。這個顯示出有某種的能量混合，正式叫做_遍歷性_，對於均分定理的成立是十分重要的。

### 粒子的澱積

位能並不一定跟位置成二次關係的。不過，均分定理指出若能量對自由度 _x_ 是 s 次齊次的（對一固定實數 _s_ 而言），則該部份於熱平衡時的平均能量為 _kBT/s_。

在重力下澱積的這個延伸有一個簡單的應用。例如在啤酒裏有時見到的薄霧能由一團團會散射光的蛋白質所組成。一段時間以後，這些蛋白質團因受重力影響而向下沉澱，使得近底下的部份比頂端的薄霧更多。不過，一個向相合方向作用的過程中，粒子也會向上滲透回到酒瓶的頂部。一達到平衡狀態時，就可以使用均分定理來斷定某一浮力質量 _mb_ 的蛋白質團的平均位置。對一支瓶高無限的啤酒而言，重力位能可由下式求出

　
: $H^{\mathrm {grav} }=m_{b}gz\,,$

其中 _z_ 為蛋白質團的高度，而 _g_ 則為重力加速度。由於 _s=1_ 的關係，蛋白質團的平均位能等於 _kBT_。因此，一個浮力質量為 10MDa（大體上為病毒的大小）的蛋白質團會於平衡狀態做出一股 2cm 高的薄霧。這樣一種往平衡的澱積由梅森 - 韋弗爾方程式所描述。

## 歷史

　
: : _本文使用非國際單位制單位的 Cal/(mol・K) 作為莫耳比熱容的單位，因為它能用單位數字提供更佳的準確度。 
如欲將該單位轉換成國際單位制的對應單位 J/(mol・K)，將原數值乘上 4.2J/Cal 可得約值。_

動能均分這個概念最早是在 1843 年，或更準確地說應是 1845 年，由約翰・詹姆斯・瓦塔斯頓提出的。於 1859 年，詹姆斯・克拉克・馬克士威主張氣體的動熱能由線性及旋轉能量所等量攤分 於 1876 年，路德維希・波茲曼因表明了平均能量是被一系統中各獨立分量所等分，而將原理進一步擴展。波茲曼亦應用了均分定理去為固體比熱容的杜隆 - 泊替定律提出了一個理論解釋。

能量均分定理的歷史與莫耳比熱容的歷史是密不可分的，兩者都是在十九世紀時被研究的。於 1819 年，法國物理學家皮埃爾・路易・杜隆（Pierre Louis Dulong）和阿列克西・泰雷茲・珀蒂（Alexis Thérèse Petit）發現了所有室溫下的固體比熱容幾乎都是相等的，約為 6 cal/(mol・K)。他們的定律曾在多年間被用作量度原子質量的一種技巧，然而，後來詹姆斯・杜瓦及海因里希・夫里德里希・韋伯的研究表明杜隆 - 泊替定律只於高溫時成立；在低溫時或像金剛石這種異常地硬的固體，比熱還要再低一點。

氣體比熱的實驗觀測也引起了對均分定理是否有效的質疑。定理預測簡單單原子氣體的莫耳比熱容應約為 3 cal/(mol・K)，而雙原子氣體則約為 7 cal/(mol・K)。實驗驗證了預測的前者，但卻發現雙原子氣體的典型莫耳比熱容約為 5 cal/(mol・K)，並於低溫時下跌到約 3 cal/(mol・K)。詹姆斯・克拉克・馬克士威於 1875 年指出實驗與均分定理的不合比這些數字暗示的要壞得多；由於原子有內部部份，熱能應該走向這些內部部分的運動，使得單原子及雙原子的比熱預測值比 3 cal/(mol・K) 7 cal/(mol・K) 要高得多。

第三個有關的不符之處是金屬的比熱。根據古典德魯德模型，金屬電子的舉止跟幾乎理想的氣體一樣，因此它們應該向 (3/2)_NekB_ 的熱容，其中 _Ne_ 為電子的數量。不過實驗指出電子對熱容的供給並不多：很多的金屬的莫耳比熱容與絕緣體幾乎一樣。

數個說明均分失敗原因的解釋被提出了。波茲曼辯護他的均分定理推導是正確的，但就提出氣體可能因為與以太交互作用而不處於熱平衡狀態。由於與實驗不符，克耳文勳爵提出均分定理的推導一定是不定確的，但卻說不出甚麼不正確。反而瑞立勳爵提出一個更徹底的看法，說均分定理及實驗時系統處於熱平衡的假設這兩者_都_正確；為使兩者相符，他指出需要一個能為均分定理提供「從破壞性的簡單中逃走的去路」的新原理。艾爾伯特・愛因斯坦就提供了這條去路，於 1907 年他表明了這些比熱的異數都是由量子效應引起的，尤其是固體彈性模態能量的量子化。愛因斯坦用了均分定理的失敗作為需要一個新物質量子理論的論據。瓦爾特・能斯特於 1910 年在低溫的比熱量度支持了愛因斯坦的理論，並引起了物理學家們對量子理論的廣泛承認。

## 能量均分定理的通用公式化

均分定理最通用的形式 明確表示在適當的假設下（下文會討論），對一個有哈密頓能量函數 _H_ 及自由度 _xn_ 的物理系統而言，以下均分公式於熱平衡時對任何值的指數 _m_ 及 _n_ 都有效：

這裏 _δmn_ 為克羅內克爾 δ，當 _m=n_ 時為一而其他情況為零。平均用的括號 $\left\langle \ldots \right\rangle $ 可以代表一個單系統長時間的平均，或，更多地是，相空間的系綜平均。定理內含的遍歷性假設意味著這兩個平均相符，而且都被用於複雜物理系統的內能估算。

通用均分定理在系統總能量恆定時於微正則系綜有效，在跟能交換能量的熱庫耦合時於正則系綜亦都有效。通式的推導在下文。

通式與以下兩式等價

1. 對所有 _n_， ${\Bigl \langle }x_{n}{\frac {\partial H}{\partial x_{n}}}{\Bigr \rangle }=k_{B}T$
2. 對所有 _m_≠_n_， ${\Bigl \langle }x_{m}{\frac {\partial H}{\partial x_{n}}}{\Bigr \rangle }=0$

若一自由度 _xn_ 只作為一個二次項 _anxn2_ 於哈密頓函數 _H_ 中出現，則可從第一式引出

　
: $k_{B}T={\Bigl \langle }x_{n}{\frac {\partial H}{\partial x_{n}}}{\Bigr \rangle }=2\langle a_{n}x_{n}^{2}\rangle ,$

即這自由度向平均能量 $\langle H\rangle $ 所提供的兩倍。因此二次能量系統的均分定理結果很容易就出來了。用一個相近的論點，將 2 換成 _s_，以 _anxns_ 的形式應用於能量中。

自由度 _xn_ 是系統相空間的坐標，因此一般被細分成廣義位置坐標 _qk_ 及廣義動量坐標 _pk_，其中 _pk_ 為 _qk_ 的共軛動量。此時，1 式對所有 _k_ 值而言

　
: ${\Bigl \langle }p_{k}{\frac {\partial H}{\partial p_{k}}}{\Bigr \rangle }={\Bigl \langle }q_{k}{\frac {\partial H}{\partial q_{k}}}{\Bigr \rangle }=k_{B}T.$

使用哈密頓力學的方程式，這些式子可被寫成

　
: ${\Bigl \langle }p_{k}{\frac {dq_{k}}{dt}}{\Bigr \rangle }=-{\Bigl \langle }q_{k}{\frac {dp_{k}}{dt}}{\Bigr \rangle }=k_{B}T.$

另外式 2 明確指出平均為

　
: ${\Bigl \langle }q_{j}{\frac {\partial H}{\partial q_{k}}}{\Bigr \rangle },\quad {\Bigl \langle }q_{j}{\frac {\partial H}{\partial p_{k}}}{\Bigr \rangle },\quad {\Bigl \langle }p_{j}{\frac {\partial H}{\partial p_{k}}}{\Bigr \rangle },\quad {\Bigl \langle }p_{j}{\frac {\partial H}{\partial q_{k}}}{\Bigr \rangle },\quad {\Bigl \langle }q_{k}{\frac {\partial H}{\partial p_{k}}}{\Bigr \rangle },$ 及 ${\Bigl \langle }p_{k}{\frac {\partial H}{\partial q_{k}}}{\Bigr \rangle }$

若 _j≠k_ 則皆為零。

### 與均功定理的關係

通用均分定理是均功定理的一個延伸。均功定理於 1870 年被提出，它明確指出

　
: ${\Bigl \langle }\sum _{k}q_{k}{\frac {\partial H}{\partial q_{k}}}{\Bigr \rangle }={\Bigl \langle }\sum _{k}p_{k}{\frac {\partial H}{\partial p_{k}}}{\Bigr \rangle }={\Bigl \langle }\sum _{k}p_{k}{\frac {dq_{k}}{dt}}{\Bigr \rangle }=-{\Bigl \langle }\sum _{k}q_{k}{\frac {dp_{k}}{dt}}{\Bigr \rangle },$

其中 _t_ 代表時間。兩個關鍵性的分別在於均功定理聯繫的是平均的_總和_而不是_個別_的平均，而且跟溫度 _T_ 沒有關係。另一個分別是均功定理的傳統推導使用時間的平均，而均分定理則用相空間的平均。

## 應用

### 理想氣體定律

理想氣體給能量均分定理提供一個重要應用。也為每粒子的平均動能提供了公式：

　
: ${\begin{aligned}\langle H^{\mathrm {kin} }\rangle &={\frac {1}{2m}}\langle p_{x}^{2}+p_{y}^{2}+p_{z}^{2}\rangle \\&={\frac {1}{2}}{\biggl (}{\Bigl \langle }p_{x}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{x}}}{\Bigr \rangle }+{\Bigl \langle }p_{y}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{y}}}{\Bigr \rangle }+{\Bigl \langle }p_{z}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{z}}}{\Bigr \rangle }{\biggr )}={\frac {3}{2}}k_{B}T\end{aligned}}$

能量均分定理可被用於從古典力學中導出理想氣體定律。若 **q**=(_qx_,_qy_,_qz_) 且 **p**=(_px_,_py_,_pz_) 代表氣體中一個粒子的位置向量及動量向量，而 **F** 為作用在該粒子上的淨力，則

　
: ${\begin{aligned}\langle \mathbf {q} \cdot \mathbf {F} \rangle &={\Bigl \langle }q_{x}{\frac {dp_{x}}{dt}}{\Bigr \rangle }+{\Bigl \langle }q_{y}{\frac {dp_{y}}{dt}}{\Bigr \rangle }+{\Bigl \langle }q_{z}{\frac {dp_{z}}{dt}}{\Bigr \rangle }\\&=-{\Bigl \langle }q_{x}{\frac {\partial H}{\partial q_{x}}}{\Bigr \rangle }-{\Bigl \langle }q_{y}{\frac {\partial H}{\partial q_{y}}}{\Bigr \rangle }-{\Bigl \langle }q_{z}{\frac {\partial H}{\partial q_{z}}}{\Bigr \rangle }=-3k_{B}T,\end{aligned}}$

其中第一條等式為牛頓第二定律，第二行用了哈密頓力學及均分公式。將系統中的 _N_ 個粒子的所有加起來得

　
: $3Nk_{B}T=-{\biggl \langle }\sum _{k=1}^{N}\mathbf {q} _{k}\cdot \mathbf {F} _{k}{\biggr \rangle }$

由牛頓第三定律及理想氣體假設得，系統的淨力是由它們的容器上的牆所施行的，而這個力可由氣體的氣壓 _P_ 所求得。因此

　
: $-{\biggl \langle }\sum _{k=1}^{N}\mathbf {q} _{k}\cdot \mathbf {F} _{k}{\biggr \rangle }=P\oint _{\mathrm {surface} }\mathbf {q} \cdot \mathbf {dS} ,$

其中'_dS_ 為沿著容器牆的無限小面積元件。由於位置向量 _q_ 散度為

　
: ${\boldsymbol {\nabla }}\cdot \mathbf {q} ={\frac {\partial q_{x}}{\partial q_{x}}}+{\frac {\partial q_{y}}{\partial q_{y}}}+{\frac {\partial q_{z}}{\partial q_{z}}}=3,$

而散度定理表明

　
: $P\oint _{\mathrm {surface} }\mathbf {q} \cdot \mathbf {dS} =P\int _{\mathrm {volume} }\left({\boldsymbol {\nabla }}\cdot \mathbf {q} \right)dV=3PV,$

其中 _dV_ 為容器內無限小的體積，而 _V_ 則為容器的總容量。

將這些等式結合得

　
: $3Nk_{B}T=-{\biggl \langle }\sum _{k=1}^{N}\mathbf {q} _{k}\cdot \mathbf {F} _{k}{\biggr \rangle }=3PV,$

即馬上導出 _N_ 個粒子理想氣體定律

　
: $PV=Nk_{B}T=nRT,\,$

其中 _n=N/NA_ 為氣體的莫耳數，而 _R=NAkB_ 則為氣體常數。

### 雙原子氣體

一雙原子氣體可用被一韌度為 _a_ 彈簧連接的兩質量 _m1_ 及 _m2_ 作為模型，稱為「剛性轉子 - 諧波振蕩器近似法」。此系統的古典能量為

　
: $H={\frac {\left|\mathbf {p} _{1}\right|^{2}}{2m_{1}}}+{\frac {\left|\mathbf {p} _{2}\right|^{2}}{2m_{2}}}+{\frac {1}{2}}aq^{2},$

其中 **p**1 與 **p**2 為兩原子的動量，而 _q_ 則為從平衡值得出的原子間距離的偏差。由於能量的每一自由度都為二次，並故此應向總平均能量提供 ½_kBT_，向熱容提供 ½_kB_。因此一有 _N_ 雙原子分子的氣體的熱容被預測為 _7N_·½_kB_：**p**_1_ 及 **p**_2_ 各提供三個自由度，_q_ 的延伸提供第七個。由此得一莫耳無其他自由度的雙原子氣體，其比熱容應為 (7/2)_NAkB_=(7/2)_R_，而因此莫耳比熱容應大概為 7 cal/(mol・K)。但是雙原子氣體摩熱比熱容的典型實驗結果約為 5 cal/(mol・K)，並於非常低溫時跌至 3 cal/(mol・K)。 這個均分預測與實驗莫耳熱容值的不符，不能被更複雜的分子模型所解釋，因為增加自由度只會_增加_預測比熱值，而不會減少。這個不符是需要物質量子理論的一個關鍵證據。

### 極度相對性理想氣體

均分定理被用於從牛頓力學中導出古典理想氣體定律。但是相對性效應在某些系統中變得顯著，例如是白矮星及中子星，而此時理想氣體定律一定要被稍作修改。均分定理為導出極度相對性理想氣體提供了方便的途徑。在這樣的個案中，一個單粒子的動能可由以下公式求得：

　
: $H^{\mathrm {kin} }\approx cp=c{\sqrt {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}}.$

以動量分量 _px_ 取 _H_ 的導數得公式：

　
: $p_{x}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{x}}}=c{\frac {p_{x}^{2}}{\sqrt {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}}}$

之後以同樣的方式處理分量 _py_ 及 _pz_。將三個分量全加起來得：

　
: ${\begin{aligned}\langle H^{\mathrm {kin} }\rangle &={\biggl \langle }c{\frac {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}{\sqrt {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}}}{\biggr \rangle }\\&={\Bigl \langle }p_{x}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{x}}}{\Bigr \rangle }+{\Bigl \langle }p_{y}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{y}}}{\Bigr \rangle }+{\Bigl \langle }p_{z}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{z}}}{\Bigr \rangle }\\&=3k_{B}T\end{aligned}}$

其中最後一條等式則是均分定理公式的結果。因此，一極度相對性理想氣體的平均總能量要比非相對性的大兩倍：當有 _N_ 個粒子時，總平均能量為 _3 N kBT_。

### 非理想氣體

在一理想氣體中，假設粒子只通過碰撞來交互作用。均分定理可以用於推導「非理想氣體」的能量和壓力，這氣體中的粒子也通過守恆力與其他粒子交互作用，守恆力的勢 _U_(_r_) 只受粒子間距離 _r_ 影響。要描述這個狀況，首先將注意力集中於一個粒子，之後把其餘的氣體當球狀對稱分佈看待。然後按慣例引入徑向分佈函數 _g(r)_，使得從某固定粒子於距離 _r_ 找到另一粒子的機率密度等於 4π_r2ρ g(r)_，其中 _ρ=N/V_ 為氣體的平均密度。可得某固定粒子跟其餘氣體的互相作用平均位能為

　
: $\langle h^{\mathrm {pot} }\rangle =\int _{0}^{\infty }4\pi r^{2}\rho U(r)g(r)\,dr$ 。

氣體的總平均位能因此為 $\langle H^{pot}\rangle ={\tfrac {1}{2}}N\langle h^{\mathrm {pot} }\rangle $ ，其中 _N_ 是氣體中的粒子數；在加起所有粒子數目的過程中會把每個粒子算兩次，所以需要用到 ½ 這個因子。

把動能及位能加起來，然後使用均分定理，可得_能量方程式_

　
: $H=\langle H^{\mathrm {kin} }\rangle +\langle H^{\mathrm {pot} }\rangle ={\frac {3}{2}}Nk_{B}T+2\pi N\rho \int _{0}^{\infty }r^{2}U(r)g(r)\,dr$ 。

用一套相近的論證下，可用於導出_壓力方程式_

　
: $3Nk_{B}T=3PV+2\pi N\rho \int _{0}^{\infty }r^{3}U'(r)g(r)\,dr$ 。

### 非諧振器

一非諧振器（跟簡單諧波振蕩器相對）的位能並不與延伸長度 _q_（量度平衡點偏差的廣義位置）成二次關係。這樣的振蕩器為均分定理提供了一個可作補充的觀點。簡單例子中用的是位能函數，其形式為

　
: $H^{\mathrm {pot} }=Cq^{s}$ ，

其中 _C_ 和 _S_ 為任意實數。這些個案中，均分定律預測

　
: $k_{B}T={\Bigl \langle }q{\frac {\partial H^{\mathrm {pot} }}{\partial q}}{\Bigr \rangle }=\langle q\cdot sCq^{s-1}\rangle =\langle sCq^{s}\rangle =s\langle H^{\mathrm {pot} }\rangle $ 。

因此，平均位能等於 _kBT/s_，而不像二次諧波振蕩器那樣等於 _kBT/2_（此時 _s_=2）。

更廣義地，一次元系統典型能量函數的有以延伸長度 _q_ 表示的泰勒級數：

　
: $H^{\mathrm {pot} }=\sum _{n=2}^{\infty }C_{n}q^{n}$

其中 _n_ 須為非負整數式子方能成立。沒有 _n_=1 的項是因為在平衡點上並無淨力作用，所以能量第一導數為零。毋需包括 _n_=0 的項是因為在約定俗成下平衡位置的能量可訂為零。此時，均分定律預測：

　
: $k_{B}T={\Bigl \langle }q{\frac {\partial H^{\mathrm {pot} }}{\partial q}}{\Bigr \rangle }=\sum _{n=2}^{\infty }\langle q\cdot nC_{n}q^{n-1}\rangle =\sum _{n=2}^{\infty }nC_{n}\langle q^{n}\rangle $ 。

相比以上引述的例子，均分定理公式

　
: $\langle H^{\mathrm {pot} }\rangle ={\frac {1}{2}}k_{B}T-\sum _{n=3}^{\infty }\left({\frac {n-2}{2}}\right)C_{n}\langle q^{n}\rangle $

並_不_能把平均位能以已知常數來表示。

### 布朗運動

均分定理可用於從朗之萬方程式中導出一粒子的布朗運動。根據該方程式，一質量為 _m_，速度為 _v_ 的粒子的運動是由牛頓第二定律支配的

　
: ${\frac {d\mathbf {v} }{dt}}={\frac {1}{m}}\mathbf {F} =-{\frac {\mathbf {v} }{\tau }}+{\frac {1}{m}}\mathbf {F} ^{\mathrm {rnd} },$

其中 **F**rnd 為一代表著粒子及其周圍分子的隨機碰撞的隨機力，而時間常數 τ 反映阻礙粒子在溶液中運動的阻力。阻力很多時候被寫成 **F**drag =-γ**v**，因此時間常數 τ 等於 _m_/γ。

取方程式與位置向量 _r_ 的點積，平均後，得布朗運動用的方程式

　
: ${\Bigl \langle }\mathbf {r} \cdot {\frac {d\mathbf {v} }{dt}}{\Bigr \rangle }+{\frac {1}{\tau }}\langle \mathbf {r} \cdot \mathbf {v} \rangle =0$

（由於 **F**rnd 跟位置向量 _r_ 不相關）。使用數學恆等式

　
: ${\frac {d}{dt}}\left(\mathbf {r} \cdot \mathbf {r} \right)={\frac {d}{dt}}\left(r^{2}\right)=2\left(\mathbf {r} \cdot \mathbf {v} \right)$

及

　
: ${\frac {d}{dt}}\left(\mathbf {r} \cdot \mathbf {v} \right)=v^{2}+\mathbf {r} \cdot {\frac {d\mathbf {v} }{dt}},$

布朗運動的基本方程式可被轉變成

　
: ${\frac {d^{2}}{dt^{2}}}\langle r^{2}\rangle +{\frac {1}{\tau }}{\frac {d}{dt}}\langle r^{2}\rangle =2\langle v^{2}\rangle ={\frac {6}{m}}k_{B}T,$

其中從最後的等式可由均分定理得平移動能：

　
: $\langle H^{\mathrm {kin} }\rangle ={\Bigl \langle }{\frac {p^{2}}{2m}}{\Bigr \rangle }=\langle {\tfrac {1}{2}}mv^{2}\rangle ={\tfrac {3}{2}}k_{B}T$ 。

以上 $\langle r^{2}\rangle $ 的微分方程式（有適當的初始條件）可被整解得

　
: $\langle r^{2}\rangle ={\frac {6k_{B}T\tau ^{2}}{m}}\left(e^{-t/\tau }-1+{\frac {t}{\tau }}\right)$ 。

在短時標上，_t_ <<_τ_，粒子的行動就像自由粒子一樣：由指數函數的泰勒級數得，距離的平方大約以_二次_增長：

　
: $\langle r^{2}\rangle \approx {\frac {3k_{B}T}{m}}t^{2}=\langle v^{2}\rangle t^{2}$ 。

但是在長時標上，_t_ >>_τ_，對數及常數項可被忽略，距離平方只以_一次_增長：

　
: $\langle r^{2}\rangle \approx {\frac {6k_{B}T\tau }{m}}t=6\gamma k_{B}Tt$ 。

這描述了粒子隨著時間的滲透。一條剛性子旋轉滲透用的類似方程式可用跟這個近似的方法導出。

### 恆星物理學

均分定理及相關的均功定理在很早以前就已被用於天體物理學。例如，均功定理可被用於由白矮星質量去估計恆星溫度或其錢德拉塞卡極限。

一顆恆星的平均溫度可由均分定理估計。由於大部份恆星都是球狀對稱，總重力位能可由積分法估算

　
: $H_{\mathrm {tot} }^{\mathrm {grav} }=-\int _{0}^{R}{\frac {4\pi r^{2}G}{r}}M(r)\,\rho (r)\,dr,$

其中 _M(r)_ 為半徑 _r_ 以內的質量而 _ρ(r)_ 則是於半徑 _r_ 時的恆星密度；_G_ 代表萬有引力常數，_R_ 為恆星總半徑。假定整個恆星內的密度一致，該積分可得公式

　
: $H_{\mathrm {tot} }^{\mathrm {grav} }=-{\frac {3GM^{2}}{5R}},$

其中 _M_ 為恆星總質量。因此一個單粒子平均位能為

　
: $H_{\mathrm {tot} }^{\mathrm {grav} }=-{\frac {3GM^{2}}{5R}},$

其中 _N_ 為恆星內的粒子數目。由於大部分恆星都是由離子化的氫所組成的，_N_ 大概等於 _(M/mp)_，其中 _mp_ 為一個質子的質量。應用均分定理可得一個恆星溫度的約值

　
: ${\Bigl \langle }r{\frac {\partial H^{\mathrm {grav} }}{\partial r}}{\Bigr \rangle }=\langle -H^{\mathrm {grav} }\rangle =k_{B}T={\frac {3GM^{2}}{5RN}}.$

將太陽的質量及半徑代入得太陽溫度的大約值 _T_ 等於一千四百萬克耳文，跟其核心溫度一千五百萬克耳文非常接近。但是，太陽被這個模型假設的要複雜得多 —— 它的溫度及質量都隨半徑而變動甚大 —— 而且如此極佳的切合度（≈7% 相對誤差）在一定情況下是偶然的。

### 恆星的形成

同一組方程式可被用於斷定巨型分子雲中的恆星形成狀態。這樣一個雲的本地密度起伏可以引致一個很壞的情況，當中分子雲會在自己的重力下向內倒塌。這樣一個倒塌在當均分定理 —— 或同樣地均功定理 —— 不再有效的時候發生，亦即當重力位能超越動能兩倍的時候

　
: ${\frac {3GM^{2}}{5R}}>3Nk_{B}T$

設雲的密度固定為 ρ

　
: $M={\frac {4}{3}}\pi R^{3}\rho $

得恆星收縮的最小質量，金斯質量 _MJ_

　
: $M_{J}^{2}=\left({\frac {5k_{B}T}{Gm_{p}}}\right)^{3}\left({\frac {3}{4\pi \rho }}\right)$

代入這樣的雲的典型觀測值（_T_=150 K，ρ=2×10−16g/cm3）可得一為 17 個太陽質量的最小質量估值，跟恆星形成的觀測一致。這個效應亦被稱為金斯不穩定性，以在 1902 年發表這理論的英國物理學家詹姆士・金斯命名。

## 推導

### 動能與馬克士威 - 波茲曼分布

均分定理的原公式化表述明確指出，在任何處熱平衡狀態的物理系統中，每一種粒子有著完全一樣的平均動能 _(3/2)kBT_。但事實上，這個結果僅適用於理想氣體，並且可使用馬克士威 - 波茲曼分布導出，方法如下。首先，我們先暫時只考慮粒子在 z 方向上的速度分布

　
: $f(v_{z})={\sqrt {\dfrac {m}{2\pi k_{B}T}}}\;e^{\frac {-m{v_{z}}^{2}}{2k_{B}T}}$

有了以上式子，便可以計算粒子在 z 方向上的方均速度

　
: $\langle {v_{z}}^{2}\rangle =\int _{-\infty }^{\infty }f(v_{z}){v_{z}}^{2}dv_{z}={\dfrac {k_{B}T}{m}}$

由於粒子在 x、y、z 各方向上的速度互相獨立且不互相影響，因此粒子的平均動能為

　
: $\langle E_{k}\rangle ={\dfrac {3}{2}}m\langle {v_{z}}^{2}\rangle ={\dfrac {3}{2}}k_{B}T$

注意，馬克士威 - 波茲曼分布不應與波茲曼分布混淆，因為前者必須假設粒子的總能等於其平移動能，才能由後者導出。

### 二次能量與配分函數

更一般地說，能量均分定理明確指出任何只在總能量 $H$ 中只以一簡單二次項 $Ax^{2}$ （此處 $A$ 為一常數）出現的自由度 $x$ ，於熱平衡時有平均能量 ${\frac {1}{2}}k_{B}T$ 。在這個案中均分定理可從配分函數 $Z(\beta )$ 中導出，其中 $\beta ={\frac {1}{k_{B}T}}$ ，為正則溫度倒數。從 $Z$ 的式子中，經變量 $x$ 取積分得因子

　
: $Z_{x}=\int _{-\infty }^{\infty }dx\ e^{-\beta Ax^{2}}={\sqrt {\frac {\pi }{\beta A}}}$ ，

而跟這個因子有關的平均能量為

　
: $\langle H_{x}\rangle =-{\frac {\partial \log Z_{x}}{\partial \beta }}={\frac {1}{2\beta }}={\frac {1}{2}}k_{B}T$

跟均分定理指出的一致。

### 一般證明

均分定理的一般證明可於很多統計力學的教科書中找到，包括微正則系綜用 及正則系綜用 的。它們都需要取系統相空間中的平均，而它是一個辛流形。

要解釋這些推導，需要引入以下的記法。首先，相空間是由廣義位置坐標 _q_j 跟它們的共軛動量 _p_j 一起描述的。量 _q_j 完整描述系統位形，而量 (_q_j,_p_j) 一起則描述完整描述系統的狀態。

第二，相空間的無限小體積

　
: $d\Gamma =\prod _{i}dq_{i}dp_{i}$

需被引入以用於定義

　
: $\Gamma (E,\Delta E)=\int _{H\in \left[E,E+\Delta E\right]}d\Gamma .$

在這式子中，_ΔE_ 被假定為非常小的，_ΔE<<E_。同樣地，Σ(_E_) 被定義為相空間中能量被 _E_ 低的總體積：

　
: $\Sigma (E)=\int _{H<E}d\Gamma .$

由於 _ΔE_ 非常小，以下這兩條積分式是等價的

　
: $\int _{H\in \left[E,E+\Delta E\right]}\ldots d\Gamma =\Delta E{\frac {\partial }{\partial E}}\int _{H<E}\ldots d\Gamma ,$

其中省略號代表被積函數。由此可得 Γ 與 ΔE 成正比

　
: $\Gamma =\Delta E\ {\frac {\partial \Sigma }{\partial E}}=\Delta E\ \rho (E),$

此處 _ρ(E)_ 為狀態密度。由統計力學常用的定義得，熵 _S_ 等於 _kB_log_Σ(E)_，而溫度被定義為

　
: ${\frac {1}{T}}={\frac {\partial S}{\partial E}}=k_{b}{\frac {\partial \log \Sigma }{\partial E}}=k_{b}{\frac {1}{\Sigma }}\,{\frac {\partial \Sigma }{\partial E}}$ 。

#### 正則系綜

在正則系綜中，系統正與一溫度為 _T_（以克耳文計算）的無限熱庫處與熱平衡。在相空間中每一態的機率由波茲曼因子乘上重整化因子 ${\mathcal {N}}$ 得出，而重整化因子的值是為使總機率成一而設的。

　
: ${\mathcal {N}}\int e^{-\beta H(p,q)}d\Gamma =1,$

其中 _β=1/kBT_。為一相空間變量 _xk_（可為 _qk_ 或 _pk_）使用分部積分法得方程式

　
: ${\mathcal {N}}\int \left[e^{-\beta H(p,q)}x_{k}\right]_{x_{k}=a}^{x_{k}=b}d\Gamma _{k}+{\mathcal {N}}\int e^{-\beta H(p,q)}x_{k}\beta {\frac {\partial H}{\partial x_{k}}}d\Gamma =1,$

其中 _dΓk=dΓ/dxk_ ，即第一個積分並沒有取 _xk_ 積分的。第一項通常為零，這是因為在極限下 _xk_ 為零，或因為能量在那極限下趨向無限。此時，可馬上得正則系綜的均分定理

　
: ${\mathcal {N}}\int e^{-\beta H(p,q)}x_{k}{\frac {\partial H}{\partial x_{k}}}\,d\Gamma ={\Bigl \langle }x_{k}{\frac {\partial H}{\partial x_{k}}}{\Bigr \rangle }={\frac {1}{\beta }}=k_{B}T.$

這裏被 $\langle \ldots \rangle $ 符號化的平均是正則系綜中的正則平均。

#### 微正則系綜

在微正則系綜中，系統與外界隔絕，或只是跟外界很弱的耦合起來。因此，其總能量實際上是恆定的；要明確的話，則需指出總能量 _H_ 局限於 _E_ 及 _E+ΔE_ 之間。對某固定能量 _E_ 及其範圍 _ΔE_ 而言，有一相空間區域 Γ，在內系統有該能量且每一態在該相空間區域的機率均等。已知該等定義，相空間變量 _xm_（可為 _qk_ 或 _pk_）的均分平均為

　
: ${\begin{aligned}{\Bigl \langle }x_{m}{\frac {\partial H}{\partial x_{n}}}{\Bigr \rangle }&={\frac {1}{\Gamma }}\,\int _{H\in \left[E,E+\Delta E\right]}x_{m}{\frac {\partial H}{\partial x_{n}}}\,d\Gamma \\&={\frac {\Delta E}{\Gamma }}\,{\frac {\partial }{\partial E}}\int _{H<E}x_{m}{\frac {\partial H}{\partial x_{n}}}\,d\Gamma \\&={\frac {1}{\rho }}\,{\frac {\partial }{\partial E}}\int _{H<E}x_{m}{\frac {\partial \left(H-E\right)}{\partial x_{n}}}\,d\Gamma ,\end{aligned}}$

此處能得出最後一等式是因為 _E_ 乃一常數，並不受 _x_n 影響。使用分部積分法得方程式

　
: ${\begin{aligned}\int _{H<E}x_{m}{\frac {\partial (H-E)}{\partial x_{n}}}\,d\Gamma &=\int _{H<E}{\frac {\partial }{\partial x_{n}}}{\bigl (}x_{m}(H-E){\bigr )}\,d\Gamma -\int _{H<E}\delta _{mn}(H-E)d\Gamma \\&=\delta _{mn}\int _{H<E}(E-H)\,d\Gamma ,\end{aligned}}$

這是由於第一行右方第一項為零的緣故（該項能被寫成在 _H_=_E_ 的超曲面上 _H_-_E_ 的一個積分）。

將該結果代入前面的方程式可得

　
: ${\Bigl \langle }x_{m}{\frac {\partial H}{\partial x_{n}}}{\Bigr \rangle }=\delta _{mn}{\frac {1}{\rho }}\,{\frac {\partial }{\partial E}}\int _{H<E}\left(E-H\right)\,d\Gamma =\delta _{mn}{\frac {1}{\rho }}\,\int _{H<E}\,d\Gamma =\delta _{mn}{\frac {\Sigma }{\rho }}$ 。

由於 $\rho ={\frac {\partial \Sigma }{\partial E}}$ ，可得均分定理

　
: ${\Bigl \langle }x_{m}{\frac {\partial H}{\partial x_{n}}}{\Bigr \rangle }=\delta _{mn}{\Bigl (}{\frac {1}{\Sigma }}{\frac {\partial \Sigma }{\partial E}}{\Bigr )}^{-1}=\delta _{mn}{\Bigl (}{\frac {\partial \log \Sigma }{\partial E}}{\Bigr )}^{-1}=\delta _{mn}k_{B}T$ 。

如是者就導出了均分定理的通用公式

也就是在之前應用中如此有用的那條公式。

## 限制

### 遍歷性需求

均分定律只對處於熱平衡的遍歷系統有效，這意味著同一能量的態被遷移的可能性必然一樣。故此，系統一定要可以讓它所有各形態的能量能夠互相交換，或在正則系綜中跟一熱庫一起。已被證明為遍歷的系統數量不多；雅科夫・西奈的硬球系統是一個有名的例子 。讓隔離系統保證其遍歷性 —— 因而，均分定理 —— 的需求已被研究過，同時研究還推動了動力系統混沌理論的發展。一混沌哈密頓系統不一定是遍歷系統，儘管假定它是通常也足夠準確。

有時候能量並_不_由它的各種形式所攤分，且此時均分定理在微正則系綜_不_成立，耦合諧波振蕩器系統就是在這狀況下常被引用的一個例子。如果系統跟外界隔絕，那每一個正常模態的能量是恆定的；能量並不由一個模態傳遞到另一模態的。因此在這樣一個系統中均分定理無效；每一個模態能量的量都被它的起始值所固定。如果能量函數中有著足夠強的非線性量的時候，能量可能可以在正常模態中傳遞，使系統走向遍歷並使均分定律有效。然而，柯爾莫哥洛夫 - 阿諾德 - 莫澤定理明確指出除非擾動夠強，否則能量不會交換；如擾動小的話，最低限度能量會繼續受困於一些模態中。

### 量子效應引起的失敗

當熱能 _kBT_ 比能階間的差要小得多的時候，均分法則就會失效。均分此時不再成立，是因為能階組成平滑連續能譜的這個假設跟實際情況不近似，而這假設在上面均分定理推導中有用到。歷史上，古典均分定理在解釋比熱及黑體輻射時的失敗，對表明需要一套物質及輻射的新理論（即量子力學及量子場論）起了關鍵性的作用。

要說明均分的失效，可考慮一單（量子）諧波振蕩子的平均能量，古典個案在上文已討論過。它的量子能階為 _En=nhν_，其中 _h_ 為普朗克常數，_ν_ 為振蕩子的基本頻率，而 _n_ 則為一整數。某指定能階正被置於正則系綜的機率可由其波茲曼因子得出

　
: $P(E_{n})={\frac {e^{-n\beta h\nu }}{Z}}$ ，

其中 _β_=1/_k_B_T_，而分母中的 _Z_ 為配分函數，此處為一等比級數

　
: $Z=\sum _{n=0}^{\infty }e^{-n\beta h\nu }={\frac {1}{1-e^{-\beta h\nu }}}$ 。

則其平均能量為

　
: $\langle H\rangle =\sum _{n=0}^{\infty }E_{n}P(E_{n})={\frac {1}{Z}}\sum _{n=0}^{\infty }nh\nu \ e^{-n\beta h\nu }=-{\frac {1}{Z}}{\frac {\partial Z}{\partial \beta }}=-{\frac {\partial \log Z}{\partial \beta }}$ 。

將 _Z_ 的式子代入得最後結果：

　
: $\langle H\rangle =h\nu {\frac {e^{-\beta h\nu }}{1-e^{-\beta h\nu }}}$ 。

於高溫時，當熱能 _kBT_ 比能階差 _hν_ 大得多的時候，指數變量 _βhν_ 比一要小得多，所以平均能量成了 _kBT_，跟均分定理一致（見圖十）。然而於低溫時，當 _hν_>>_kBT_ 的時候，平均能量走向零 —— 高頻能階被「凍結」了（見圖十）。作為另一例子，氫原子內部的受激電子態在室溫下並不提供任何比熱，這是由於熱能 _kBT_（大概是 0.025 eV）比最低及下一高能階之間的差（大概是 10 eV）要小得多的緣故。

相近的考量可用於任何能階差比熱能大得多的狀況下。例如，艾爾伯特・愛因斯坦就是用這套論證解決黑體輻射的紫外災變。由於在一封閉容器下的電磁場有無限個獨立模態，每一個都能被當作諧波振蕩器看待，因而就形成了悖論。如果每一個電磁模態皆有平均能量 _kBT_，容器內的能量將為無限大。然而，根據以上的論證，高 ω 模態的平均值當 ω 趨向無限時趨向零；而且描術模態實驗中能量分佈的普朗克黑體輻射定律，也是根據同一組論證所中得出的。

此外，更微妙的量子效應可引起均分定理的修正，例如全同粒子及連續對稱。全同粒子效應可在非常高密度且低溫時有著顯著的效果。比方說金屬的價電子可以有幾個電子伏的平均能量，正常情況一般對應數萬克耳文的溫度。如此的狀態，密度高得讓包立不相容原理使得古典門徑無效化，被稱為簡並態費米子氣體。這種氣體對白矮星及中子星的結構很重要。在低溫時，玻色 - 愛因斯坦凝聚（此凝聚中大量全同粒子佔據了低能量態）的費米子類比能夠形成；這種超流體電子是引起超導現象的成因。

## 另見

* 均功定理
* 分子運動論
* 統計力學
* 量子統計力學

## 註釋及參考資料

## 延伸閱讀

## 外部連結

* 一單原子及雙原子混合氣體的均分定理即時 Applet 展示 （頁面存檔備份，存於網際網路檔案館）（英文）
* 恆星物理中的能量均分定理（頁面存檔備份，存於網際網路檔案館），由耶路撒冷希伯來大學 Racah 物理研究所副教授 Nir J. Shaviv 所撰寫。（英文）