ビリアル定理

テンプレート:出典の明記 ビリアル定理（ビリアルていり、テンプレート:Lang-en-short）とは、多粒子系において、粒子が動き得る範囲が有限である場合に、古典力学、量子力学系のいずれにおいても成立する以下の関係式のことである。

$⟨K⟩=⟨\sum _{i=1}^N\frac{{𝐩}_i^2}{2m_i}⟩=\sum _{i=1}^N⟨\frac{{𝐩}_i^2}{2m_i}⟩=-\frac{1}{2}\sum _{i=1}^N⟨{𝐅}_i\cdot {𝐫}_i⟩$

テンプレート:Mvar は系の粒子数、テンプレート:Mvar は系全体の運動エネルギー

$K=\sum _{i=1}^N\frac{{𝐩}_i^2}{2m_i}$

で、テンプレート:Math は粒子 テンプレート:Mvar の運動量、テンプレート:Math は粒子 テンプレート:Mvar の位置座標、テンプレート:Math は粒子 テンプレート:Mvar に働く力、テンプレート:Mvar は粒子 テンプレート:Mvar の質量である。テンプレート:Math は物理量の平均操作（ここでは長時間平均）を意味する。

粒子 テンプレート:Mvar に働く力 テンプレート:Math が、系全体のポテンシャルエネルギー テンプレート:Math を用いて テンプレート:Math と表せるならば、ビリアル定理は、

$⟨K⟩=\frac{1}{2}\sum _{i=1}^N⟨{\mathrm{\nabla }}_{{𝐫}_i}V\cdot {𝐫}_i⟩$

という形で表せる。

ポテンシャルエネルギー テンプレート:Mvar が中心力ポテンシャルで、粒子間の距離のテンプレート:Math乗に比例する形

$V(𝐫)=a{r}^{n+1}(r=|𝐫|)$

で表せる(ここでべき指数は力の法則が$r^n$になるように選んだ)ならば、

$⟨K⟩=\frac{n+1}{2}⟨V⟩$

となる。中心力が電磁気力や重力の場合を考えると、テンプレート:Math であるから、

$⟨K⟩=-\frac{1}{2}⟨V⟩$

となる。ビリアル定理から次のことが言える。

• 系全体の運動エネルギー テンプレート:Mvar の時間平均は、系全体のポテンシャルエネルギー テンプレート:Mvar の時間平均の テンプレート:Math に等しい。

また、同等のこととして、

• 系全体のポテンシャルエネルギー テンプレート:Mvar の時間平均は、系全体の全エネルギーの時間平均に等しい。

• 系全体の運動エネルギー テンプレート:Mvar の時間平均と系全体の全エネルギーの時間平均を加えた物は テンプレート:Math。

ということが示される。

ビリアルとはラテン語で「力」という意味であり、ビリアル定理の名はそれに因む。ビリアル定理におけるビリアルとは、1870年にルドルフ・クラウジウスが導入した量で、各粒子の位置と運動量のドット積の総和 テンプレート:Math によって定義される テンプレート:Mvar を指す。

証明

古典力学系の場合のビリアル定理の証明。ビリアル

を時間で微分すると、

$\frac{dG}{dt}=\sum _i\frac{d{𝐫}_i}{dt}\cdot {𝐩}_i+\sum _i{𝐫}_i\cdot \frac{d{𝐩}_i}{dt}=\sum _i\frac{{𝐩}_i}{m_i}\cdot {𝐩}_i+\sum _i{𝐫}_i\cdot {𝐅}_i=2K+\sum _i{𝐫}_i\cdot {𝐅}_i$

より以下の関係が得られる。 テンプレート:NumBlk

この式の両辺を テンプレート:Math から時間 テンプレート:Mvar の範囲で積分して テンプレート:Mvar で割り、テンプレート:Math の極限をとって長時間平均する。すると、粒子が動き得る範囲は有限なのでビリアル テンプレート:Mvar も有限だから、左辺は 0 に収束する。

したがって、

$0=2⟨K⟩+⟨\sum _i{𝐫}_i\cdot {𝐅}_i⟩$

つまり、ビリアル定理

を得る。 次に、ポテンシャルエネルギー テンプレート:Mvar が中心力ポテンシャルで、粒子間の距離の テンプレート:Math 乗 (テンプレート:Math) に比例する形、すなわち、系のポテンシャル テンプレート:Mvar が各粒子対の相互作用の和

によって書き表される場合、粒子 テンプレート:Mvar に働く力 テンプレート:Math は、以下のように書ける。

ここで、

は、粒子 テンプレート:Mvar から粒子 テンプレート:Mvar に働く力である。これを、テンプレート:EquationNoteに代入すると、以下のようになる。

和は テンプレート:Math、テンプレート:Math の2重の和である。 この和を テンプレート:Math と テンプレート:Math に分け、

$\sum _{i\ne j}{𝐫}_i\cdot {𝐅}_{ji}=\sum _{i>j}{𝐫}_i\cdot {𝐅}_{ji}+\sum _{i<j}{𝐫}_i\cdot {𝐅}_{ji}$

第 2 項で添え字の入れ替えに対する反対称性 テンプレート:Math に注意すると、以下の様な形になる。

$=\sum _{i>j}{𝐫}_i\cdot {𝐅}_{ji}+\sum _{j<i}{𝐫}_j\cdot {𝐅}_{ij}=\sum _{i>j}{𝐫}_i\cdot {𝐅}_{ji}-\sum _{j<i}{𝐫}_j\cdot {𝐅}_{ji}=\sum _{i>j}({𝐫}_i-{𝐫}_j)\cdot {𝐅}_{ji}=-(n+1)V({𝐫}_1,\cdots ,{𝐫}_N).$

したがって、中心力ポテンシャルに関するビリアル定理は以下のようになる。

応用

ビリアル定理を太陽系や銀河を始めとする、非常に複雑な物理体系（重力多体系）に適用することにより、計算結果を簡素化することができるので非常に便利である。

また、ビリアル定理が成り立つ場合、次式から系の圧力を求めることができる。

$P=\frac{1}{3V}(2⟨K⟩+\sum _{i=1}^N⟨{𝐅}_i\cdot {𝐫}_i⟩)$

ここで、テンプレート:Mvar は圧力、テンプレート:Mvar は系の体積である。気体分子運動論では上式から圧力を求める。

一般化

一般化されたビリアル定理を、超ビリアル定理 テンプレート:En と言う。座標 テンプレート:Math と共役運動量 テンプレート:Math を考え、この 2 つの量を変数とした関数 テンプレート:Math を考える。この関数は、冒頭での粒子系と同様な境界条件の基で任意に選べるとする。ハミルトニアンを テンプレート:Mvar として、ポアソン括弧（詳細はハミルトン力学を参照）の時間平均、

$⟨[W,H]⟩=0$

となるのが古典的な超ビリアル定理である。量子力学では、上記交換子の基底状態における期待値がゼロとなる。

$⟨0|[W,H]|0⟩=0$

これが量子力学的な超ビリアル定理である。ここで、テンプレート:Mvar として上記のビリアルをとる。すなわち、

$W=𝐫\cdot 𝐏$

とすれば、通常のビリアル定理が導かれる。

関連項目

• 多体問題
• 統計力学
• 物性物理学
• ビリアル応力