αヘリカル ペプチド の熱運動 。この不安定な運動はランダムかつ複雑で、特定の原子のエネルギーは大きく変動する可能性がある。しかしながら、等分配定理を用いることで、各原子の平均運動エネルギーと、多くの振動モードの平均位置エネルギーを計算することができる。灰色、赤色、青色の球はそれぞれ炭素 、 酸素 、窒素原子を表し、小さな白色 の球 は 水素原子を表す 。古典 統計力学 において、等分配定理は 系の温度とその平均 エネルギー を関連付けます。等分配定理は、等分配の法則 、エネルギーの等分配 、あるいは単に等分配 とも呼ばれます。等分配の本来の考え方は、熱平衡状態 において、エネルギーはその様々な形態の間で均等に分配されるというものでした。例えば、分子の並進運動 における自由度 あたりの平均運動エネルギーは、 回転運動 におけるそれと等しくなければなりません。
等分配定理は定量的な予測を行います。ビリアル定理 と同様に、与えられた温度における系の運動エネルギーと位置エネルギーの平均値を与え、そこから系の熱容量を計算できます。しかし、等分配定理は、特定の粒子の運動エネルギーや単一の バネ の位置エネルギーなど、エネルギーの個々の構成要素の平均値も与えます。例えば、単原子 理想気体 中のすべての原子の平均運動エネルギーはであると予測されます。3 / 2 熱平衡状態における k B T で、 k B はボルツマン定数 、 Tは (熱力学的)温度 です、どんなに複雑な熱平衡状態 にある古典系 にも適用できます。これは、理想気体の法則 や、固体の比熱容量 に関するデュロン・プティの法則 を導くために使用できます。 [ 1 ] 等分配定理は相対論的 効果を考慮し ても成り立つため、白色矮星 や中性子星などの 恒星 の特性を予測するためにも使用できます
等分配定理は特定の条件下では正確な予測を行いますが、低温など量子効果 が顕著な場合には不正確になります。特定の自由度における 熱エネルギー k B T が 量子エネルギー間隔よりも小さい場合、この自由度の平均エネルギーと熱容量は等分配定理によって予測される値よりも小さくなります。熱エネルギーがこの間隔よりもはるかに小さい場合、そのような自由度は「凍結されている」と言われます。例えば、固体の熱容量は低温では減少しますが、これは様々な種類の運動が凍結されるためであり、等分配定理によって予測されるように一定のままであるのに対し、固体は減少します。このような熱容量の減少は、19世紀の物理学者にとって、古典物理学が誤りであり、より繊細で新しい科学的モデルが必要であることを示す最初の兆候の一つでした。他の証拠とともに、等分配法では黒体放射( 紫外線大惨事 としても知られる)をモデル化できなかったことから、マックス・プランクは 、光を発する物体の振動子のエネルギーは量子化されていると提唱しました。これは量子力学 と量子場理論 の発展を促した革命的な仮説でした。
基本的な概念と簡単な例 図2. 298.15 K (25 °C )における4種の希ガス の分子速度の確率密度関数。4 種のガスは、ヘリウム ( 4 He )、ネオン ( 20 Ne )、アルゴン ( 40 Ar )、キセノン ( 132 Xe ) である。上付き文字はそれぞれの質量数を示す。これらの確率密度関数は、確率×速度の逆数の 次元 を持つ。確率は無次元であるため、秒/メートルの単位で表すことができる。 「等分配」という名称は「均等な分割」を意味し、先行詞æquus(「等しい、あるいは均等な」)から派生したラテン語の equiと、名詞 partitio (「分割、部分」)から派生したpartitionに由来する。[ 2 ] [ 3 ] 等分配の本来の概念は、系が熱平衡に達すると、系の全運動エネルギーは 平均して すべての独立した部分の間で均等に分配されるというものである。等分配はこれらのエネルギーを定量的に予測する。例えば、温度T で熱平衡状態にある不活性希ガス の各原子は、平均並進運動エネルギーが であると予測される。 3 / 2 k B T 、ここでk B はボルツマン定数 です。結果として、運動エネルギーは1 ⁄ 2 (質量)(速度) 2に等しいため、同じ温度では、 キセノン の重い原子はヘリウム の軽い原子よりも平均速度が低くなります、4つの希ガス中の原子の速度の マクスウェル・ボルツマン分布を 示しています
この例では、運動エネルギーが速度の2乗に比例するという点が重要です。等分配定理は、熱平衡状態において、エネルギーの2乗にしか現れない自由度 (例えば、粒子の位置や速度の成分)は平均エネルギーが1 ⁄ 2 k B T であり、したがって系の熱容量に 1 ⁄ 2 k B 寄与することを示しています。これは多くの応用があります。
並進エネルギーと理想気体 質量m 、速度v の粒子の(ニュートン力学)運動エネルギーは、v x 、v y 、v z が速度v の直交座標成分である として、次のように表されます 。ここで、Hは ハミルトニアン(Hamiltonian) の略で、ハミルトニアン形式は 等分割定理の 最も一般的な形式 において中心的な役割を果たすため、以降はエネルギーの記号として使用されます。H 親族 = 1 2 m | v | 2 = 1 2 m ( v × 2 + v y 2 + v z 2 ) 、 {\displaystyle H_{\text{kin}}={\tfrac {1}{2}}m|\mathbf {v} |^{2}={\tfrac {1}{2}}m\left(v_{x}^{2}+v_{y}^{2}+v_{z}^{2}\right),}
運動エネルギーは速度の各成分の2乗に比例するため、等分配により、これら3つの成分はそれぞれ熱平衡状態における平均運動エネルギーに1 ⁄ 2 k B T 寄与する。したがって、粒子の平均運動エネルギーは 3 / 2 k B T 、上記の希ガスの例と同じです
より一般的には、単原子理想気体では、全エネルギーは純粋に(並進)運動エネルギーから成ります。仮定によれば、粒子は内部自由度を持たず、互いに独立して運動します。したがって、等分配則から、N 個の粒子からなる理想気体の全エネルギーは、 3 / 2 N k B T 。
したがって、気体の熱容量は 3 / 2 N k B であり、したがって、特に、そのようなガス粒子1モル の熱容量は 3 / 2 N A k B = 3 / 2 R 、ここでN A はアボガドロ定数 、 Rは 気体定数 である。R ≈ 2 cal /( mol · K ) なので、 等分配 法則から 理想気体の モル熱容量 はおよそ 3 cal/(mol · K )と予測される。この予測は、単原子気体と比較した実験によって確認されている。 [ 4 ]
平均運動エネルギーを用いることで、気体粒子の 二乗平均速度 v rms を計算することも可能となる。 ここで、M = N A m は気体粒子1モルあたりの質量である。この結果は、ウラン 濃縮 法を規定するグラハムの 流出 の法則など、多くの応用に役立つ。[ 5 ] v 実効値 = ⟨ v 2 ⟩ = 3 k B T m = 3 R T M 、 {\displaystyle v_{\text{rms}}={\sqrt {\left\langle v^{2}\right\rangle}}={\sqrt {\frac {3k_{\text{B}}T}{m}}}={\sqrt {\frac {3RT}{M}}},}
溶液中の回転エネルギーと分子の回転 同様の例として、主慣性モーメントがそれぞれ I 1 、I 2 、I 3 である回転分子が挙げられます。古典力学によれば、このような分子の 回転エネルギー は次のように表されます。 ここで、 ω 1 、ω 2 、ω 3 は角速度 の主成分です。並進の場合と全く同じ理屈で、等分配は熱平衡状態において各粒子の平均回転エネルギーが H r o t = 1 2 ( I 1 ω 1 2 + I 2 ω 2 2 + I 3 ω 3 2 ) 、 {\displaystyle H_{\mathrm {rot} }={\tfrac {1}{2}}(I_{1}\omega _{1}^{2}+I_{2}\omega _{2}^{2}+I_{3}\omega _{3}^{2}),} 3 / 2 k B T . 同様に、等分配定理により、分子の平均(より正確には、二乗平均平方根)角速度を計算することができます。 [ 6 ]
剛体分子の回転、すなわち溶液中の分子のランダムな回転は、核磁気共鳴 、特にタンパク質NMR と残留双極子結合 によって観測される緩和 において重要な役割を果たしている。[ 7 ] 回転拡散は、蛍光異方性 、流動複屈折 、誘電分光法 などの他の生物物理学的プローブによっても観測することができる。[ 8 ]
位置エネルギーと調和振動子 等分配は運動エネルギー だけでなく位置エネルギーにも適用されます。重要な例としては、 バネ などの調和振動子 が挙げられます。バネの 弾性定数a はバネの剛性を表し、q は平衡状態からの偏差です。このような1次元系の質量がm の場合、その運動エネルギーH kin は 、 v とp = mv は振動子の速度と運動量を表します。これらの項を組み合わせると、全エネルギーが得られます[ 9 ] 。したがって、等 分配は、熱平衡状態において振動子が平均エネルギーを持つことを意味し、 角括弧は括弧内の量の平均を表します[ 10 ]。 H ポット = 1 2 a q 2 、 {\displaystyle H_{\text{pot}}={\tfrac {1}{2}}aq^{2},\,} H 親族 = 1 2 m v 2 = p 2 2 m 、 {\displaystyle H_{\text{kin}}={\frac {1}{2}}mv^{2}={\frac {p^{2}}{2m}},} H = H kin + H pot = p 2 2 m + 1 2 a q 2 . {\displaystyle H=H_{\text{kin}}+H_{\text{pot}}={\frac {p^{2}}{2m}}+{\frac {1}{2}}aq^{2}.} ⟨ H ⟩ = ⟨ H kin ⟩ + ⟨ H pot ⟩ = 1 2 k B T + 1 2 k B T = k B T , {\displaystyle \langle H\rangle =\langle H_{\text{kin}}\rangle +\langle H_{\text{pot}}\rangle ={\tfrac {1}{2}}k_{\text{B}}T+{\tfrac {1}{2}}k_{\text{B}}T=k_{\text{B}}T,} ⟨ … ⟩ {\displaystyle \left\langle \ldots \right\rangle }
この結果は、振り子 、振動する分子、受動的な電子振動子 など、あらゆるタイプの調和振動子に当てはまります。このような振動子のシステムは多くの状況で発生します。等分配により、各振動子は平均全エネルギーk B T を受け取り、したがってシステムの熱容量に k B を寄与します。これは、ジョンソン・ナイキスト雑音の式 [ 11 ] と固体の熱容量に関するデュロン・プティの法則 を導くために使用できます。特に後者の応用は、等分配の歴史において重要な意味を持ちました。
図3. 結晶中の原子は、格子 内の平衡位置を中心に振動することができる。このような振動は、結晶誘電体の 熱容量 の大部分を占める。金属 の場合、電子 も熱容量に寄与する。
固体の比熱容量 等分配定理の重要な応用の一つは、結晶固体の比熱容量である。このような固体中の各原子は3つの独立した方向に振動することができるため、固体は3 N 個の独立した単振動子の 系とみなすことができる。ここで、Nは 格子中の原子数を表す。各振動子の平均エネルギーはk B T であるため、固体の平均全エネルギーは3 N k B T であり、その熱容量は3 N k B である。
アボガドロ定数 N A をN とし、気体定数 R とボルツマン定数k B の関係式R = N A k B を用いることで、固体の比熱容量 に関するデュロン・プティの法則 を説明できます。この法則は、固体元素の比熱容量(単位質量あたり)は原子量 に反比例するというものです。現代版では、固体のモル熱容量は3R ≈ 6 cal/(mol·K)とされています。
しかし、この法則は量子効果のために低温では不正確であり、また、実験的に導かれた熱力学の第三法則 とも矛盾している。この法則によれば、いかなる物質のモル熱容量も温度が絶対零度に近づくとゼロになるはずである。[ 11 ] 量子効果を取り入れたより正確な理論は、アルバート・アインシュタイン (1907年)とピーター・デバイ (1911年)によって開発された。[ 12 ]
他の多くの物理システムは、結合した振動子 の集合としてモデル化できます。このような振動子の運動は、ピアノの弦 の振動モードやオルガンのパイプの 共鳴のように、 基準振動モード に分解できます。一方、このようなシステムでは、基準振動モード間でエネルギーの交換が行われないため、等分配則がしばしば破綻します。極端な状況では、各モードは独立しており、したがってそれらのエネルギーは独立に保存されます。これは、正式にはエルゴード性 と呼ばれる、ある種のエネルギー混合が、等分配則の成立に重要であることを示しています。
粒子の沈降 位置エネルギーは必ずしも位置の2乗の関数ではありません。しかし、等分配定理は、自由度xがエネルギーに x・ s (一定の実数s に対して)の倍数のみを寄与する場合、熱平衡状態においてその部分の平均エネルギーはk・ B・ T / s であることを示しています
この拡張は、重力 下における粒子の沈降 に簡単に応用できる。[ 13 ] 例えば、ビールに時々見られる濁りは、光 を散乱させる タンパク質の 塊によって発生する可能性がある。[ 14 ] 時間が経つにつれて、これらの塊は重力の影響を受けて下方に沈降し、ボトルの上部よりも底部に近い方に濁りが多くなる。しかし、反対方向に働くプロセスでは、粒子はボトルの上部に向かって再び拡散する 。平衡に達すると、等分配定理を使用して、浮力のある質量 m b の特定の塊の平均位置を決定できる。無限に高いビール瓶の場合、重力による位置エネルギーは 次のように与えられる。 ここで、 z はボトル内のタンパク質の塊の高さ、g は重力による加速度 である。 s = 1 なので、タンパク質の塊の平均位置エネルギーはk B T に等しい。したがって、浮遊質量10 MDa ( ウイルス とほぼ同じ大きさ)のタンパク質凝集体は、 平衡状態で平均高さ約2cmのヘイズを生成する。このような沈降過程はメイソン・ウィーバーの式 で記述される。[ 15 ] H g r a v = m b g z {\displaystyle H^{\mathrm {grav} }=m_{\text{b}}gz}
歴史 運動エネルギーの等分配は、1843年、より正確には1845年にジョン・ジェームズ・ウォーターストン によって最初に提案されました。[ 16 ] 1859年、ジェームズ・クラーク・マクスウェルは 、気体の運動熱エネルギーは直線エネルギーと回転エネルギーに均等に分割されると主張しました。[ 17 ] 1876年、ルートヴィヒ・ボルツマンは この原理を拡張し、平均エネルギーは系内のすべての独立した運動成分に均等に分割されることを示しました。[ 18 ] [ 19 ] ボルツマンは等分配定理を応用して、固体の 比熱容量 に関するデュロン・プティの法則 を理論的に説明しました
図4. 二原子気体のモル比熱 と温度の関係を理想化したグラフ。高温では等分配から予測される(7/2) R ( R は気体定数 )と一致するが、(5/2) R まで減少し、その後、 3 / 2 低温では、振動運動と回転 運動のモードが「凍結」されるため、 Rに等しい。等分配定理の破綻は 量子力学 によってのみ解決されたパラドックスにつながった。ほとんどの分子では、遷移温度T rot は 室温よりはるかに低いが、T vib は 10倍以上高くなることがある。典型的な例は一酸化炭素 (CO)で、T rot ≈ 2.8 K 、T vib ≈ 3103 K である。非常に大きな原子や弱く結合した原子を持つ分子では、T vib は 室温(約300 K)に近くなることがある。例えば、ヨウ素 ガス(I 2 )ではT vib ≈ 308 Kである。 [ 20 ] 等分配定理の歴史は、19世紀に研究された比熱容量の歴史と絡み合っています。1819年、フランスの物理学者 ピエール=ルイ・デュロン とアレクシ・テレーズ・プティは 、室温での固体元素の比熱容量はその元素の原子量に反比例することを発見しました。[ 21 ] 彼らの法則は、原子量を測定する技術として長年使用されました。[ 12 ] しかし、その後のジェームズ ・デュワーとハインリッヒ・フリードリヒ・ウェーバー による研究で、このデュロン=プティの法則は 高温 でのみ成立することが示されました。[ 22 ] より低い温度、またはダイヤモンド などの非常に硬い固体では、比熱容量は低くなります。[ 23 ]
気体の比熱容量の実験的観察によって、等分配定理の妥当性に対する懸念も生じた。定理によれば、単純な一原子気体のモル熱容量はおよそ 3 cal/(mol·K) であるのに対し、二原子気体のモル熱容量はおよそ 7 cal/(mol·K) である。実験によって前者の予測は確認されたが[ 4 ] 、二原子気体のモル熱容量は典型的には約 5 cal/(mol·K) であり[ 24 ] 、非常に低温では約 3 cal/(mol·K) まで低下することが判明した[ 25 ] 。マクスウェルは 1875 年に、実験と等分配定理の不一致はこれらの数字が示唆するよりもはるかに大きいと指摘した。[ 26 ] 原子には内部部品があるため、熱エネルギーはこれらの内部部品の運動に使われるはずで、単原子気体と二原子気体の比熱はそれぞれ3cal/(mol·K)と7cal/(mol·K)よりもはるかに高くなります。
3つ目の矛盾は金属の比熱に関するものであった。[ 27 ] 古典的なドルーデモデル によれば、金属電子はほぼ理想気体として振る舞うため、 3 / 2 等 分配定理により、熱容量は電子数N e k Bに比例します。 しかし 、 実験 的 には、電子は熱容量にほとんど寄与しません。多くの導体と絶縁体のモル熱容量はほぼ同じです。 [ 27 ]
等分配定理がモル熱容量を説明できない理由については、いくつかの説明が提案された。ボルツマンは 等分配定理の導出は正しいと主張したが、エーテル との相互作用により気体が熱平衡状態 にない可能性を示唆した。[ 28 ] ケルビン卿は 、等分配定理の導出は実験と矛盾するため誤りであると示唆したが、その理由を示すことはできなかった。[ 29 ] 1900年、レイリー卿は、等分配定理と熱平衡の実験的仮定は どちらも 正しいという、より急進的な見解を提唱した。彼は、両者を調和させるためには、等分配定理の「破壊的な単純さからの脱出」を可能にする新しい原理が必要であると指摘した。[ 30 ] アルバート・アインシュタインは 1906年に比熱の異常は量子効果、特に固体の弾性モードにおけるエネルギーの量子化によるものであることを示してその逃げ道を提供した。[ 31 ] アインシュタインは等分配の失敗を利用して、物質の新しい量子論の必要性を主張した。[ 12 ] ネルンストの 1910年の低温での比熱の測定[ 32 ] はアインシュタインの理論を支持し、物理学者の間で量子論 が広く受け入れられることになった。[ 33 ]
等分割定理の最も一般的な形式は、適切な仮定(後述)の下で、 ハミルトン エネルギー関数H と自由度x n を持つ物理系において、すべての指数m とn に対して熱平衡状態で次の等分割式が成り立つことを述べています。[ 6 ] [ 10 ] [ 13 ] ここでδ mn はクロネッカーのデルタ であり、 m = n の場合には1 、それ以外の場合には0です。平均括弧は位相空間全体のアンサンブル平均、または エルゴード性 の仮定の下では単一系の時間平均であると仮定 されます。⟨ x m ∂ H ∂ x n ⟩ = δ m n k B T . {\displaystyle \left\langle x_{m}{\frac {\partial H}{\partial x_{n}}}\right\rangle =\delta _{mn}k_{\text{B}}T.} ⟨ … ⟩ {\displaystyle \left\langle \ldots \right\rangle }
一般的な等分配定理は、系 の全エネルギーが一定の場合のミクロカノニカル集団 [ 10 ] と、系がエネルギー交換可能な熱浴に結合されている場合 の カノニカル集団 [ 6 ] [ 34 ] の両方において成立する。一般式の導出は本稿の後半で 示す。
一般的な式は次の 2 つと同等です。
⟨ x n ∂ H ∂ x n ⟩ = k B T for all n {\displaystyle \left\langle x_{n}{\frac {\partial H}{\partial x_{n}}}\right\rangle =k_{\text{B}}T\quad {\text{for all }}n} ⟨ x m ∂ H ∂ x n ⟩ = 0 for all m ≠ n . {\displaystyle \left\langle x_{m}{\frac {\partial H}{\partial x_{n}}}\right\rangle =0\quad {\text{for all }}m\neq n.} 自由度x n が ハミルトニアンH において二次項a n x n 2 としてのみ現れる場合、これらの式の最初の式は、 この自由度が平均エネルギーに与える寄与の2倍であるを意味する。したがって、二次エネルギーを持つ系の等分配定理は、この一般公式から容易に導かれる。同様の議論は、 2 をs に置き換えれば、 a n x n s という形式のエネルギーにも当てはまる。 k B T = ⟨ x n ∂ H ∂ x n ⟩ = 2 ⟨ a n x n 2 ⟩ , {\displaystyle k_{\text{B}}T=\left\langle x_{n}{\frac {\partial H}{\partial x_{n}}}\right\rangle =2\left\langle a_{n}x_{n}^{2}\right\rangle ,} ⟨ H ⟩ {\displaystyle \langle H\rangle }
自由度x n は 系の位相空間 上の座標であり、したがって、一般的には一般化位置 座標q k と一般化運動量 座標p k に細分される。ここで、p k はq k の共役運動量 である。この場合、式1はすべてのk に対して、 ⟨ p k ∂ H ∂ p k ⟩ = ⟨ q k ∂ H ∂ q k ⟩ = k B T . {\displaystyle \left\langle p_{k}{\frac {\partial H}{\partial p_{k}}}\right\rangle =\left\langle q_{k}{\frac {\partial H}{\partial q_{k}}}\right\rangle =k_{\text{B}}T.}
ハミルトン力学 の方程式を用いると[ 9 ] 、これらの式は次のようにも書ける。 ⟨ p k d q k d t ⟩ = − ⟨ q k d p k d t ⟩ = k B T . {\displaystyle \left\langle p_{k}{\frac {dq_{k}}{dt}}\right\rangle =-\left\langle q_{k}{\frac {dp_{k}}{dt}}\right\rangle =k_{\text{B}}T.}
同様に、式2を用いて 、 ⟨ p j ∂ H ∂ p k ⟩ = ⟨ q j ∂ H ∂ q k ⟩ = 0 for all j ≠ k . {\displaystyle \left\langle p_{j}{\frac {\partial H}{\partial p_{k}}}\right\rangle =\left\langle q_{j}{\frac {\partial H}{\partial q_{k}}}\right\rangle =0\quad {\text{ for all }}\,j\neq k.} ⟨ p j ∂ q k ∂ t ⟩ = − ⟨ q j ∂ p k ∂ t ⟩ = 0 for all j ≠ k . {\displaystyle \left\langle p_{j}{\frac {\partial q_{k}}{\partial t}}\right\rangle =-\left\langle q_{j}{\frac {\partial p_{k}}{\partial t}}\right\rangle =0\quad {\text{ for all }}\,j\neq k.}
ビリアル定理との関係 一般的な等分割定理は、ビリアル定理 (1870年提唱[ 35 ] ) の拡張であり、 t が 時間 を 表すと述べている。[ 9 ] 2つの重要な違いは、ビリアル定理は個々の 平均ではなく合計された 平均を互いに関連付け、それらを温度 T に結び付けないことである。もう1つの違いは、ビリアル定理の従来の導出では時間にわたる平均が使用されるのに対し、等分割定理の導出では位相空間 にわたる平均が使用されることである。 ⟨ ∑ k q k ∂ H ∂ q k ⟩ = ⟨ ∑ k p k ∂ H ∂ p k ⟩ = ⟨ ∑ k p k d q k d t ⟩ = − ⟨ ∑ k q k d p k d t ⟩ , {\displaystyle \left\langle \sum _{k}q_{k}{\frac {\partial H}{\partial q_{k}}}\right\rangle =\left\langle \sum _{k}p_{k}{\frac {\partial H}{\partial p_{k}}}\right\rangle =\left\langle \sum _{k}p_{k}{\frac {dq_{k}}{dt}}\right\rangle =-\left\langle \sum _{k}q_{k}{\frac {dp_{k}}{dt}}\right\rangle ,}
応用
理想気体の法則 理想気体は 等分配定理の重要な応用例である。 粒子あたりの平均運動エネルギーの式を与えるだけでなく 、等分配定理は古典力学から理想気体の法則 を導くためにも使用できる。[ 6 ] q = ( q x , q y , q z ) およびp = ( p x , p y , p z )を 気体中の粒子の位置ベクトルと運動量、 Fを その粒子に働く正味の力とすると、 最初の等式はニュートンの第二法則 、2行目はハミルトン方程式と等分配定理を用いる。N 個 の粒子 系について合計すると、⟨ H k i n ⟩ = 1 2 m ⟨ p x 2 + p y 2 + p z 2 ⟩ = 1 2 ( ⟨ p x ∂ H k i n ∂ p x ⟩ + ⟨ p y ∂ H k i n ∂ p y ⟩ + ⟨ p z ∂ H k i n ∂ p z ⟩ ) = 3 2 k B T {\displaystyle {\begin{aligned}\langle H^{\mathrm {kin} }\rangle &={\frac {1}{2m}}\langle p_{x}^{2}+p_{y}^{2}+p_{z}^{2}\rangle \\&={\frac {1}{2}}\left(\left\langle p_{x}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{x}}}\right\rangle +\left\langle p_{y}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{y}}}\right\rangle +\left\langle p_{z}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{z}}}\right\rangle \right)={\frac {3}{2}}k_{\text{B}}T\end{aligned}}} ⟨ q ⋅ F ⟩ = ⟨ q x d p x d t ⟩ + ⟨ q y d p y d t ⟩ + ⟨ q z d p z d t ⟩ = − ⟨ q x ∂ H ∂ q x ⟩ − ⟨ q y ∂ H ∂ q y ⟩ − ⟨ q z ∂ H ∂ q z ⟩ = − 3 k B T , {\displaystyle {\begin{aligned}\langle \mathbf {q} \cdot \mathbf {F} \rangle &=\left\langle q_{x}{\frac {dp_{x}}{dt}}\right\rangle +\left\langle q_{y}{\frac {dp_{y}}{dt}}\right\rangle +\left\langle q_{z}{\frac {dp_{z}}{dt}}\right\rangle \\&=-\left\langle q_{x}{\frac {\partial H}{\partial q_{x}}}\right\rangle -\left\langle q_{y}{\frac {\partial H}{\partial q_{y}}}\right\rangle -\left\langle q_{z}{\frac {\partial H}{\partial q_{z}}}\right\rangle =-3k_{\text{B}}T,\end{aligned}}} 3 N k B T = − ⟨ ∑ k = 1 N q k ⋅ F k ⟩ . {\displaystyle 3Nk_{\text{B}}T=-\left\langle \sum _{k=1}^{N}\mathbf {q} _{k}\cdot \mathbf {F} _{k}\right\rangle .}
図5. 特定の分子の運動エネルギーは大きく変動する 可能性がありますが、等分配定理を用いることで、任意の温度における平均エネルギーを計算できます。また、等分配定理は、気体の 圧力 、体積 、温度を 関連付ける式である理想気体の法則 を導出します。(この図では、5つの分子が運動を追跡するために赤色で表示されています。この色分けには他に意味はありません。) ニュートンの第三法則 と理想気体の仮定により、系に働く正味の力は容器の壁によって加えられる力であり、この力は気体の 圧力P によって与えられます。 したがって、 d S は容器の壁に沿った微小面積要素です。位置ベクトルq の発散 は なので 、 発散定理 から次の式が成り立ちます。 ここで、dV は容器内の微小体積、V は 容器の全体積です。 − ⟨ ∑ k = 1 N q k ⋅ F k ⟩ = P ∮ surface q ⋅ d S , {\displaystyle -\left\langle \sum _{k=1}^{N}\mathbf {q} _{k}\cdot \mathbf {F} _{k}\right\rangle =P\oint _{\text{surface}}\mathbf {q} \cdot d\mathbf {S} ,} ∇ ⋅ q = ∂ q x ∂ q x + ∂ q y ∂ q y + ∂ q z ∂ q z = 3 , {\displaystyle {\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 ∮ s u r f a c e q ⋅ d S = P ∫ v o l u m e ( ∇ ⋅ q ) d V = 3 P V , {\displaystyle P\oint _{\mathrm {surface} }\mathbf {q} \cdot \mathbf {dS} =P\int _{\mathrm {volume} }\left({\boldsymbol {\nabla }}\cdot \mathbf {q} \right)\,dV=3PV,}
これらの式をまとめると、 N 個の粒子に対する理想気体の法則 が直ちに導かれる。 ここで、n = N / N A は気体のモル数、R = N A k B は気体定数 である。等分配法は理想気体の法則と内部エネルギーの簡単な導出を提供するが、分配関数 を用いた別の方法でも同じ結果が得られる。[ 36 ] 3 N k B T = − ⟨ ∑ k = 1 N q k ⋅ F k ⟩ = 3 P V , {\displaystyle 3Nk_{\text{B}}T=-\left\langle \sum _{k=1}^{N}\mathbf {q} _{k}\cdot \mathbf {F} _{k}\right\rangle =3PV,} P V = N k B T = n R T , {\displaystyle PV=Nk_{\text{B}}T=nRT,}
二原子気体 二原子気体は、 剛性 a のバネ で結合された2つの 質量m⁻¹ とm⁻² としてモデル化することができ、これは剛体回転子-調和振動子近似 と呼ばれます。[ 20 ] この系の古典エネルギーは、 p⁻¹ と p⁻²が 2 つの原子の運動量、qが 原子間距離の平衡値からの偏差である 場合に、次式で表されます。エネルギーの自由度はすべて2乗であり、したがって、全平均エネルギーに1 ⁄ 2 k B T 、熱容量に 1 ⁄ 2 k B 寄与するはずです。したがって、 N 個の二原子分子の気体の熱容量は7 N · 1 ⁄ 2 k B と予測されます。運動量p⁻¹ とp⁻² は それぞれ3つの自由度を寄与し、拡張qは7番目 の 自由度を寄与します。したがって、他の自由度を持たない1モルの二原子分子の熱容量は、 H = | p 1 | 2 2 m 1 + | p 2 | 2 2 m 2 + 1 2 a q 2 , {\displaystyle H={\frac {\left|\mathbf {p} _{1}\right|^{2}}{2m_{1}}}+{\frac {\left|\mathbf {p} _{2}\right|^{2}}{2m_{2}}}+{\frac {1}{2}}aq^{2},} 7 / 2 N A k B = 7 / 2 R であり、したがって、予測されるモル熱容量はおよそ 7 cal/(mol·K) となるはずです。しかし、二原子気体のモル熱容量の実験値は典型的には約 5 cal/(mol·K) [ 24 ] であり、非常に低い温度では 3 cal/(mol·K) まで低下します。 [ 25 ] 等分配予測とモル熱容量の実験値との間のこの不一致は、より複雑な分子モデルを用いて説明することはできません。なぜなら、自由度を追加しても予測される比熱は増加するだけで、減少することはないからです。 [ 26 ] この不一致は、物質の量子論 の必要性を示す重要な証拠でした。
図6. かに星雲 のX線画像と可視光画像を組み合わせたもの。この星雲の中心には、太陽 の約1.5倍の質量を持ちながら、直径はわずか25kmしかない高速回転する中性子星 が存在する。等分配定理は、このような中性子星の特性を予測するのに役立つ。
極限相対論的理想気体 上記では、等分配法則を使用して、ニュートン力学 から古典的な理想気体の法則 を導出しました。しかし、白色矮星 や中性子星 などの一部のシステムでは相対論的効果が 支配的になり、[ 10 ] 理想気体の方程式を修正する必要があります。等分配定理は、極端に相対論的な理想気体 に対応する法則を導出する便利な方法を提供します。[ 6 ] このような場合、単一粒子 の運動エネルギーは次の式で与えられます。 H をp x 運動量成分について微分すると、次の式が得られます 。 p y およびp z 成分 についても同様です。 3 つの成分を加算すると、次の式が得られます 。最後の等式は、等分配法則から得られます。したがって、極端に相対論的な気体の平均全エネルギーは、非相対論的な場合の 2 倍、つまりN 粒子の場合は3 Nk B T になります。 H k i n ≈ c p = c p x 2 + p y 2 + p z 2 . {\displaystyle H_{\mathrm {kin} }\approx cp=c{\sqrt {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}}.} p x ∂ H k i n ∂ p x = c p x 2 p x 2 + p y 2 + p z 2 {\displaystyle 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}}}}} ⟨ H k i n ⟩ = ⟨ c p x 2 + p y 2 + p z 2 p x 2 + p y 2 + p z 2 ⟩ = ⟨ p x ∂ H k i n ∂ p x ⟩ + ⟨ p y ∂ H k i n ∂ p y ⟩ + ⟨ p z ∂ H k i n ∂ p z ⟩ = 3 k B T {\displaystyle {\begin{aligned}\langle H_{\mathrm {kin} }\rangle &=\left\langle c{\frac {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}{\sqrt {p_{x}^{2}+p_{y}^{2}+p_{z}^{2}}}}\right\rangle \\&=\left\langle p_{x}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{x}}}\right\rangle +\left\langle p_{y}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{y}}}\right\rangle +\left\langle p_{z}{\frac {\partial H^{\mathrm {kin} }}{\partial p_{z}}}\right\rangle \\&=3k_{\text{B}}T\end{aligned}}}
非理想気体 理想気体では、粒子は衝突のみで相互作用すると仮定される。等分配定理は、粒子が保存力 によっても相互作用する「非理想気体」のエネルギーと圧力を導くためにも用いられる。この場合、保存力のポテンシャルU ( r ) は粒子間の距離rのみに依存する。 [ 6 ] この状況は、まず単一の気体粒子に注目し、残りの気体を球対称 分布で近似することで記述できる。そして、ある粒子からr の距離に別の粒子が存在する確率密度が 4πr 2 ρg ( r ) となるような放射状分布関数 g ( r ) を導入するのが通例である。ここでρ = N / V は気体の平均密度 である。 [ 37 ] 与えられた粒子とガスの残りの部分との相互作用に関連する平均位置エネルギーは、 したがって、ガス全体の平均位置エネルギーは、N はガス中の粒子の数であり、すべての粒子についての合計は各相互作用を2回カウントするため、係数1 ⁄ 2が 必要である。運動エネルギーと位置エネルギーを加算し、等分配を適用すると、エネルギー方程式 が得られる。 同様の議論[ 6 ] は、圧力方程式 を導くのにも使用できる。⟨ h p o t ⟩ = ∫ 0 ∞ 4 π r 2 ρ U ( r ) g ( r ) d r . {\displaystyle \langle h_{\mathrm {pot} }\rangle =\int _{0}^{\infty }4\pi r^{2}\rho U(r)g(r)\,dr.} ⟨ H pot ⟩ = 1 2 N ⟨ h p o t ⟩ {\displaystyle \langle H_{\text{pot}}\rangle ={\tfrac {1}{2}}N\langle h_{\mathrm {pot} }\rangle } H = ⟨ H k i n ⟩ + ⟨ H p o t ⟩ = 3 2 N k B T + 2 π N ρ ∫ 0 ∞ r 2 U ( r ) g ( r ) d r . {\displaystyle H=\langle H_{\mathrm {kin} }\rangle +\langle H_{\mathrm {pot} }\rangle ={\frac {3}{2}}Nk_{\text{B}}T+2\pi N\rho \int _{0}^{\infty }r^{2}U(r)g(r)\,dr.} 3 N k B T = 3 P V + 2 π N ρ ∫ 0 ∞ r 3 U ′ ( r ) g ( r ) d r . {\displaystyle 3Nk_{\text{B}}T=3PV+2\pi N\rho \int _{0}^{\infty }r^{3}U'(r)g(r)\,dr.}
非調和振動子 非調和振動子(単純調和振動子とは対照的に)は、位置エネルギーが拡張q (平衡状態からのシステムの偏差を測定する一般化された位置 )において二次ではない振動子です。このような振動子は、等分配定理に補完的な観点を提供します。[ 38 ] [ 39 ] 簡単な例として、 C とs が 任意の実定数 で ある形の位置エネルギー関数が挙げられます 。これらの場合、等分配の法則によれば 、平均位置エネルギーはk B T / s となり、二次調和振動子(s = 2 )の場合のような k B T /2 にはなりません。 H p o t = C q s , {\displaystyle H_{\mathrm {pot} }=Cq^{s},\,} k B T = ⟨ q ∂ H p o t ∂ q ⟩ = ⟨ q ⋅ s C q s − 1 ⟩ = ⟨ s C q s ⟩ = s ⟨ H p o t ⟩ . {\displaystyle k_{\text{B}}T=\left\langle q{\frac {\partial H_{\mathrm {pot} }}{\partial q}}\right\rangle =\langle q\cdot sCq^{s-1}\rangle =\langle sCq^{s}\rangle =s\langle H_{\mathrm {pot} }\rangle .}
より一般的には、1 次元システムの典型的なエネルギー関数は、 非負の整数 n に対して、拡大q : のテイラー展開を 持ちます。平衡点 では、正味の力はなく、したがってエネルギーの 1 階微分は 0なので、 n = 1 の項はありません。平衡位置でのエネルギーは慣例により 0 に設定できるため、 n = 0 の項を含める必要はありません。この場合、等分配の法則により、[ 38 ] が予測されます。 ここで引用した他の例とは対照的に、等分配式では、 平均ポテンシャルエネルギーを既知の定数で表すこと はできません。 H p o t = ∑ n = 2 ∞ C n q n {\displaystyle H_{\mathrm {pot} }=\sum _{n=2}^{\infty }C_{n}q^{n}} k B T = ⟨ q ∂ H p o t ∂ q ⟩ = ∑ n = 2 ∞ ⟨ q ⋅ n C n q n − 1 ⟩ = ∑ n = 2 ∞ n C n ⟨ q n ⟩ . {\displaystyle k_{\text{B}}T=\left\langle q{\frac {\partial H_{\mathrm {pot} }}{\partial q}}\right\rangle =\sum _{n=2}^{\infty }\langle q\cdot nC_{n}q^{n-1}\rangle =\sum _{n=2}^{\infty }nC_{n}\langle q^{n}\rangle .} ⟨ H p o t ⟩ = 1 2 k B T − ∑ n = 3 ∞ ( n − 2 2 ) C n ⟨ q n ⟩ {\displaystyle \langle H_{\mathrm {pot} }\rangle ={\frac {1}{2}}k_{\text{B}}T-\sum _{n=3}^{\infty }\left({\frac {n-2}{2}}\right)C_{n}\langle q^{n}\rangle }
ブラウン運動 図7. 三次元における粒子のブラウン運動の例 等分配定理は、ランジュバン方程式 から粒子のブラウン運動を 導くために用いることができる。[ 6 ] この方程式によれば、質量m 、速度v の粒子の運動はニュートン力学第二法則 に従う。 ここで、F rnd は 粒子と周囲の分子とのランダム衝突を表すランダム力であり、時定数τは溶液中における粒子の運動に抗する 抗力 を表す。抗力はF drag = − γ v と表記されることが多いため、時定数τは m / γ に等しい。 d v d t = 1 m F = − v τ + 1 m F r n d , {\displaystyle {\frac {d\mathbf {v} }{dt}}={\frac {1}{m}}\mathbf {F} =-{\frac {\mathbf {v} }{\tau }}+{\frac {1}{m}}\mathbf {F} _{\mathrm {rnd} },}
この方程式と位置ベクトルr のドット積を平均化すると、ブラウン運動の方程式が得られます (ランダム力F rnd は 位置r と無相関であるため)。 数学的な恒等式 と ブラウン運動の基本方程式を使用すると、次のように変形できます。 ここで、最後の等式は、並進運動エネルギーの等分割定理から得られます。 の上記の微分方程式 は(適切な初期条件を使用すれば) 正確に解くことができます。 t ≪ τ の小さな時間スケールでは、粒子は自由に移動する粒子として動作します。指数関数 のテイラー級数 により、距離の二乗はほぼ2 乗で 増加します。 ただし、 t ≫ τ の長い時間スケールでは、指数項と定数項は無視でき、距離の二乗は線形に のみ増加します。 これは、時間の経過に伴う粒子の拡散 を表します。 剛体分子の回転拡散の類似の方程式も、同様の方法で導出できます。 ⟨ r ⋅ d v d t ⟩ + 1 τ ⟨ r ⋅ v ⟩ = 0 {\displaystyle \left\langle \mathbf {r} \cdot {\frac {d\mathbf {v} }{dt}}\right\rangle +{\frac {1}{\tau }}\langle \mathbf {r} \cdot \mathbf {v} \rangle =0} d d t ( r ⋅ r ) = d d t ( r 2 ) = 2 ( r ⋅ v ) {\displaystyle {\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)} d d t ( r ⋅ v ) = v 2 + r ⋅ d v d t , {\displaystyle {\frac {d}{dt}}\left(\mathbf {r} \cdot \mathbf {v} \right)=v^{2}+\mathbf {r} \cdot {\frac {d\mathbf {v} }{dt}},} d 2 d t 2 ⟨ r 2 ⟩ + 1 τ d d t ⟨ r 2 ⟩ = 2 ⟨ v 2 ⟩ = 6 m k B T , {\displaystyle {\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_{\text{B}}T,} ⟨ H k i n ⟩ = ⟨ p 2 2 m ⟩ = ⟨ 1 2 m v 2 ⟩ = 3 2 k B T . {\displaystyle \langle H_{\mathrm {kin} }\rangle =\left\langle {\frac {p^{2}}{2m}}\right\rangle =\langle {\tfrac {1}{2}}mv^{2}\rangle ={\tfrac {3}{2}}k_{\text{B}}T.} ⟨ r 2 ⟩ {\displaystyle \langle r^{2}\rangle } ⟨ r 2 ⟩ = 6 k B T τ 2 m ( e − t / τ − 1 + t τ ) . {\displaystyle \langle r^{2}\rangle ={\frac {6k_{\text{B}}T\tau ^{2}}{m}}\left(e^{-t/\tau }-1+{\frac {t}{\tau }}\right).} ⟨ r 2 ⟩ ≈ 3 k B T m t 2 = ⟨ v 2 ⟩ t 2 . {\displaystyle \langle r^{2}\rangle \approx {\frac {3k_{\text{B}}T}{m}}t^{2}=\langle v^{2}\rangle t^{2}.} ⟨ r 2 ⟩ ≈ 6 k B T τ m t = 6 k B T t γ . {\displaystyle \langle r^{2}\rangle \approx {\frac {6k_{\text{B}}T\tau }{m}}t={\frac {6k_{\text{B}}Tt}{\gamma }}.}
恒星物理学 等分配定理と関連するビリアル定理は、 天体物理学 のツールとして長い間使用されてきました。[ 40 ] 例として、ビリアル定理は恒星の温度や白色 矮星の質量に関するチャンドラセカール限界を 推定するために使用されることがあります。[ 41 ] [ 42 ]
恒星の平均温度は等分配定理から推定できる。[ 43 ] ほとんどの恒星は球対称なので、総重力 位置エネルギーは 次の積分で推定できる。 ここで、M ( r )は半径 r 内の質量、ρ ( r )は半径 r における恒星の密度、Gは 重力定数 、Rは 恒星の全半径である。恒星全体で密度が一定であると仮定すると、この積分から次の式が得られる。 ここで、 M は恒星の全質量である。したがって、単一粒子の平均位置エネルギーは 次の式で表される。 ここで、 N は恒星内の粒子の数である。ほとんどの恒星は主に イオン化した 水素 で構成されているため、Nは おおよそM / m p に等しくなり、m p は陽子1個の質量である。等分配定理を適用することで、恒星の温度を推定することができます。 太陽 の質量と半径を代入すると、太陽温度はT = 1400万ケルビンと推定され 、これは太陽核の温度である1500万ケルビンに非常に近い値となります。しかし、太陽はこのモデルで想定されているよりもはるかに複雑であり、温度と密度は半径によって大きく変化します。そのため、このような良好な一致(相対誤差 約7% )は、部分的には偶然の産物です。[ 44 ] H g r a v = − ∫ 0 R 4 π r 2 G r M ( r ) ρ ( r ) d r , {\displaystyle H_{\mathrm {grav} }=-\int _{0}^{R}{\frac {4\pi r^{2}G}{r}}M(r)\,\rho (r)\,dr,} H g r a v = − 3 G M 2 5 R , {\displaystyle H_{\mathrm {grav} }=-{\frac {3GM^{2}}{5R}},} ⟨ H g r a v ⟩ = H g r a v N = − 3 G M 2 5 R N , {\displaystyle \langle H_{\mathrm {grav} }\rangle ={\frac {H_{\mathrm {grav} }}{N}}=-{\frac {3GM^{2}}{5RN}},} ⟨ r ∂ H g r a v ∂ r ⟩ = ⟨ − H g r a v ⟩ = k B T = 3 G M 2 5 R N . {\displaystyle \left\langle r{\frac {\partial H_{\mathrm {grav} }}{\partial r}}\right\rangle =\langle -H_{\mathrm {grav} }\rangle =k_{\text{B}}T={\frac {3GM^{2}}{5RN}}.}
巨大分子雲における 星形成 の条件を決定する際にも、同じ公式を適用することができます。[ 45 ] このような雲の密度の局所的な変動は、雲が自身の重力によって内側に崩壊する暴走状態につながる可能性があります。このような崩壊は、等分配定理(または同等のビリアル定理 )がもはや有効でなくなったとき、つまり重力による位置エネルギーが運動エネルギーの2倍を超えたときに発生します。 雲の 密度ρが 一定であると仮定すると 、恒星の収縮の最小質量、つまりジーンズ質量M J が得られます。 このような雲で典型的に観測される値(T = 150 K 、ρ = 3 G M 2 5 R > 3 N k B T . {\displaystyle {\frac {3GM^{2}}{5R}}>3Nk_{\text{B}}T.} M = 4 3 π R 3 ρ {\displaystyle M={\frac {4}{3}}\pi R^{3}\rho } M J 2 = ( 5 k B T G m p ) 3 ( 3 4 π ρ ) . {\displaystyle M_{\text{J}}^{2}=\left({\frac {5k_{\text{B}}T}{Gm_{p}}}\right)^{3}\left({\frac {3}{4\pi \rho }}\right).} 2 × 10 −16 g/cm 3 )は、最小質量が太陽質量の17倍と推定され、観測された星形成と一致する。この効果は、1902年に発表したイギリスの物理学者ジェームズ・ホップウッド・ジーンズ にちなんで、ジーンズ不安定性としても知られている。 [ 46 ]
導出
運動エネルギーとマクスウェル・ボルツマン分布等 分配定理の元の定式化は、熱平衡状態 にあるあらゆる物理系において、すべての粒子が正確に同じ平均並進運動エネルギー を持つというものです3 / 2 k B T . [ 47 ] しかし、これは理想気体の 場合にのみ当てはまり、同じ結果はマクスウェル・ボルツマン分布 からも導き出せます。まず、z成分の速度のマクスウェル・ボルツマン分布のみを考慮することにしますf ( v z ) = m 2 π k B T e − m v z 2 2 k B T {\displaystyle f(v_{z})={\sqrt {\dfrac {m}{2\pi k_{\text{B}}T}}}\;e^{\frac {-m{v_{z}}^{2}}{2k_{\text{B}}T}}}
この式で、 Z 成分 の平均二乗速度を計算できる。⟨ v z 2 ⟩ = ∫ − ∞ ∞ f ( v z ) v z 2 d v z = k B T m {\displaystyle \langle {v_{z}}^{2}\rangle =\int _{-\infty }^{\infty }f(v_{z}){v_{z}}^{2}dv_{z}={\dfrac {k_{\text{B}}T}{m}}}
速度の異なる成分は互いに独立しており、同じ分布関数を共有するため、平均並進運動エネルギーは次のように表される。 ⟨ E k ⟩ = 3 2 m ⟨ v z 2 ⟩ = 3 2 k B T {\displaystyle \langle E_{k}\rangle ={\dfrac {3}{2}}m\langle {v_{z}}^{2}\rangle ={\dfrac {3}{2}}k_{\text{B}}T}
注意:マクスウェル-ボルツマン分布を ボルツマン分布 と混同しないでください。前者は、粒子のエネルギーがその並進運動エネルギーに等しいと仮定することで後者から導き出すことができます。
等分割定理が述べているように、粒子が特定の量子エネルギー状態にある確率を用いて粒子エネルギーを平均化することによっても同じ結果が得られます。[ 36 ]
二次エネルギーと分配関数 より一般的には、等分配定理は、全エネルギーH に単純な2次項Ax 2 ( A は定数)としてのみ現れる自由度 xは 、熱平衡状態で平均エネルギーが1 ⁄ 2 k B T を持つことを述べています。この場合、等分配定理は分配関数 Z ( β ) ( β = 1/( k B T ) は標準的な温度逆数 )から導くことができます。[ 48 ] 変数x についての積分は、Z の式における因子となります 。この因子に関連する平均エネルギーは、等分配定理によれば、次のように与えられます 。 Z x = ∫ − ∞ ∞ d x e − β A x 2 = π β A , {\displaystyle Z_{x}=\int _{-\infty }^{\infty }dx\ e^{-\beta Ax^{2}}={\sqrt {\frac {\pi }{\beta A}}},} ⟨ H x ⟩ = − ∂ log Z x ∂ β = 1 2 β = 1 2 k B T {\displaystyle \langle H_{x}\rangle =-{\frac {\partial \log Z_{x}}{\partial \beta }}={\frac {1}{2\beta }}={\frac {1}{2}}k_{\text{B}}T}
一般的な証明 等分割定理の一般的な導出は、ミクロカノニカル集団 [ 6 ] [ 10 ] とカノニカル集団 [ 6 ] [ 34 ] の両方について、多くの統計力学の 教科書に記載されています。それらは、 シンプレクティック多様体 である系の位相空間 にわたって平均を取ることを伴います
これらの導出を説明するために、以下の表記法を導入する。まず、位相空間は一般化位置座標 q j とその共役運動量 p j によって記述される。q j は系の構成を 完全に記述し、 ( q j , p j )は系の 状態を 完全 に記述する。
次に、位相空間の微小体積 を導入し、系のエネルギーHが2つの限界値 E とE + Δ E の間にある位相空間の部分の体積Σ( E , Δ E ) を定義する。 この式では、Δ E は非常に小さいと仮定し、Δ E ≪ E とする。同様に、Ω( E )は、エネルギーが E 未満である位相空間の全体積として定義される。 d Γ = ∏ i d q i d p i {\displaystyle d\Gamma =\prod _{i}dq_{i}\,dp_{i}\,} Σ ( E , Δ E ) = ∫ H ∈ [ E , E + Δ E ] d Γ . {\displaystyle \Sigma (E,\Delta E)=\int _{H\in \left[E,E+\Delta E\right]}d\Gamma .} Ω ( E ) = ∫ H < E d Γ . {\displaystyle \Omega (E)=\int _{H<E}d\Gamma .}
Δ E は非常に小さいので、 楕円が積分関数を表す以下の積分は同値である 。このことから、 Σは Δ E に比例することがわかる。 ここでρ ( E )は 状態密度 である。統計力学 の通常の定義によれば、エントロピー Sは k B log Ω( E ) に等しく、温度 T は次のように定義される 。∫ H ∈ [ E , E + Δ E ] … d Γ = Δ E ∂ ∂ E ∫ H < E … d Γ , {\displaystyle \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 , Δ E ) = Δ E ∂ Ω ∂ E = Δ E ρ ( E ) , {\displaystyle \Sigma (E,\Delta E)=\Delta E\ {\frac {\partial \Omega }{\partial E}}=\Delta E\ \rho (E),} 1 T = ∂ S ∂ E = k B ∂ log Ω ∂ E = k B 1 Ω ∂ Ω ∂ E . {\displaystyle {\frac {1}{T}}={\frac {\partial S}{\partial E}}=k_{\text{B}}{\frac {\partial \log \Omega }{\partial E}}=k_{\text{B}}{\frac {1}{\Omega }}\,{\frac {\partial \Omega }{\partial E}}.}
正統派アンサンブル 正準集団 では、システムは温度 T (ケルビン)の無限熱浴と熱平衡状態にある。 [ 6 ] [ 34 ] 位相空間 における各状態の確率は、ボルツマン因子と 正規化因子の 積で与えられ、 β = 1/( k B T ) となるように確率の合計が1になるように選択される 。位相空間変数x k の部分積分 を用いると、上記は次のように表すことができる 。 ここで、d Γ k = d Γ/ dx k 、すなわち、最初の積分はx k 上では実行されない。2つの限界a とb の間で最初の積分を実行し、2番目の積分を簡略化すると、次式が得られる。 N {\displaystyle {\mathcal {N}}} N ∫ e − β H ( p , q ) d Γ = 1 , {\displaystyle {\mathcal {N}}\int e^{-\beta H(p,q)}d\Gamma =1,} N ∫ e − β H ( p , q ) d Γ = N ∫ d [ x k e − β H ( p , q ) ] d Γ k − N ∫ x k ∂ e − β H ( p , q ) ∂ x k d Γ , {\displaystyle {\mathcal {N}}\int e^{-\beta H(p,q)}d\Gamma ={\mathcal {N}}\int d[x_{k}e^{-\beta H(p,q)}]d\Gamma _{k}-{\mathcal {N}}\int x_{k}{\frac {\partial e^{-\beta H(p,q)}}{\partial x_{k}}}d\Gamma ,} N ∫ [ e − β H ( p , q ) x k ] x k = a x k = b d Γ k + N ∫ e − β H ( p , q ) x k β ∂ H ∂ x k d Γ = 1 , {\displaystyle {\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,}
最初の項は通常ゼロになる。これは、x k が極限でゼロになるか、あるいはエネルギーがそれらの極限で無限大になるからである。その場合、正準集団の等分配定理は直ちに成立する。 N ∫ e − β H ( p , q ) x k ∂ H ∂ x k d Γ = ⟨ x k ∂ H ∂ x k ⟩ = 1 β = k B T . {\displaystyle {\mathcal {N}}\int e^{-\beta H(p,q)}x_{k}{\frac {\partial H}{\partial x_{k}}}\,d\Gamma =\left\langle x_{k}{\frac {\partial H}{\partial x_{k}}}\right\rangle ={\frac {1}{\beta }}=k_{\text{B}}T.}
ここで、 で表される平均は、標準アンサンブル に対して取られたアンサンブル平均 です。 ⟨ … ⟩ {\displaystyle \langle \ldots \rangle }
ミクロカノニカルアンサンブル ミクロカノニカルアンサンブルでは、系は他の世界から孤立しているか、少なくとも非常に弱く結合しています。[ 10 ] したがって、その全エネルギーは事実上一定です。明確にするために、全エネルギーHは E とE + dE の間に閉じ込められていると言います。与えられたエネルギーE と広がりdE に対して、系がそのエネルギーを持つ位相空間 Σ の領域があり、その位相空間 の領域内の各状態の確率は、ミクロカノニカルアンサンブルの定義により等しくなります。これらの定義から、位相空間変数x m ( q k またはp k のいずれか)とx n の等分配平均は次のように与えられます
⟨ x m ∂ H ∂ x n ⟩ = 1 Σ ∫ H ∈ [ E , E + Δ E ] x m ∂ H ∂ x n d Γ = Δ E Σ ∂ ∂ E ∫ H < E x m ∂ H ∂ x n d Γ = 1 ρ ∂ ∂ E ∫ H < E x m ∂ ( H − E ) ∂ x n d Γ , {\displaystyle {\begin{aligned}\left\langle x_{m}{\frac {\partial H}{\partial x_{n}}}\right\rangle &={\frac {1}{\Sigma }}\,\int _{H\in \left[E,E+\Delta E\right]}x_{m}{\frac {\partial H}{\partial x_{n}}}\,d\Gamma \\&={\frac {\Delta E}{\Sigma }}\,{\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 に依存しない定数であるため成り立ちます。部分積分を行うと、 最初の行の右側の最初の項がゼロであるため、関係式が得られます (これは、 H = E となる超曲面 上のH − E の積分として書き直すことができます)。 ∫ H < E x m ∂ ( H − E ) ∂ x n d Γ = ∫ H < E ∂ ∂ x n ( x m ( H − E ) ) d Γ − ∫ H < E δ m n ( H − E ) d Γ = δ m n ∫ H < E ( E − H ) d Γ , {\displaystyle {\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}}}
この結果を前の式に代入すると、 ⟨ x m ∂ H ∂ x n ⟩ = δ m n 1 ρ ∂ ∂ E ∫ H < E ( E − H ) d Γ = δ m n 1 ρ ∫ H < E d Γ = δ m n Ω ρ . {\displaystyle \left\langle x_{m}{\frac {\partial H}{\partial x_{n}}}\right\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 {\Omega }{\rho }}.}
等分割定理は次のようになります。 ρ = ∂ Ω ∂ E {\displaystyle \rho ={\frac {\partial \Omega }{\partial E}}} ⟨ x m ∂ H ∂ x n ⟩ = δ m n ( 1 Ω ∂ Ω ∂ E ) − 1 = δ m n ( ∂ log Ω ∂ E ) − 1 = δ m n k B T . {\displaystyle \left\langle x_{m}{\frac {\partial H}{\partial x_{n}}}\right\rangle =\delta _{mn}\left({\frac {1}{\Omega }}{\frac {\partial \Omega }{\partial E}}\right)^{-1}=\delta _{mn}\left({\frac {\partial \log \Omega }{\partial E}}\right)^{-1}=\delta _{mn}k_{\text{B}}T.}
このようにして、上で説明した応用 で非常に役立つ等分割定理の一般的な定式化を導き出しました 。 ⟨ x m ∂ H ∂ x n ⟩ = δ m n k B T , {\displaystyle \left\langle x_{m}{\frac {\partial H}{\partial x_{n}}}\right\rangle =\delta _{mn}k_{\text{B}}T,}
制限事項 図9.理想的に結合した振動子 の孤立系では、様々な基準モード間でエネルギーは共有され ません 。各モードのエネルギーは一定であり、他のモードのエネルギーとは独立しています。したがって、等分配定理は、ミクロカノニカル集団 (孤立している場合)ではこのような系には成立しませんが、 カノニカル集団 (熱浴に結合している場合)では成立します。しかし、モード間に十分に強い非線形結合を加えることで、エネルギーは共有され、両方の集団で等分配が成立します
エルゴード性の要件 等分配の法則は熱平衡状態にある エルゴード 系にのみ成立し、これは同じエネルギーを持つすべての状態が等しく存在する可能性がなければならないことを意味します。[ 10 ] したがって、系内の様々な形態の間で、または正準集団 内の外部熱浴 との間でエネルギーを交換できなければなりません。厳密にエルゴードであることが証明された物理系の数は少なく、有名な例はヤコフ・シナイ の剛体球系 です。[ 49 ] 孤立系がエルゴード性 、つまり等分配を保証するための要件は研究されており、現代の力学系 カオス理論 の動機となっています。カオスハミルトン系は エルゴードである必要はありませんが、それは通常良い仮定です。[ 50 ]
エネルギーが様々な形態間で共有されず 、ミクロカノニカル集団で等分配が成立しない 反例としてよく挙げられるのは、結合した調和振動子のシステムである。 [ 50 ] システムが他の部分から隔離されている場合、各通常モード のエネルギーは一定であり、エネルギーは1つのモードから別のモードに転送されない。したがって、このようなシステムでは等分配は成立せず、各通常モードのエネルギー量は初期値で固定される。エネルギー 関数に十分に強い非線形項が存在する場合、エネルギーは通常モード間で転送され、エルゴード性をもたらし、等分配の法則が有効になる可能性がある。しかし、コルモゴロフ・アーノルド・モーザーの定理に よれば、非線形摂動が十分に強くない限りエネルギーは交換されず、非線形摂動が小さすぎると、エネルギーは少なくともいくつかのモードに閉じ込められたままになる。
もう一つの簡単な例は、円形容器内で衝突する有限個の粒子からなる理想気体である。容器の対称性により、このような気体の角運動量は保存される。したがって、同じエネルギーを持つ状態がすべて存在するわけではない。その結果、平均粒子エネルギーはこの粒子の質量だけでなく、他のすべての粒子の質量にも依存する。[ 51 ]
エルゴード性が破れるもう一つの方法は、非線形ソリトン 対称性の存在である。1953年、フェルミ 、パスタ 、ウラム 、ツィンゴウは、非線形項(あるテストでは2次、別のテストでは3次、さらに3次への区分線形近似)を含む振動弦の コンピュータシミュレーション を行った。彼らは、システムの挙動が等分割に基づく直感から予想されたものとは全く異なることを発見した。モードのエネルギーが等しく共有されるのではなく、システムは非常に複雑な準周期的挙動を示した。この不可解な結果は、最終的に1965年にクラスカルとザブスキーによって論文で説明され、シミュレーションされたシステムをコルテヴェク・ド・フリース方程式 に結び付けることによってソリトン数学の発展につながった。
量子効果による失敗 熱エネルギーk B T がエネルギー準位間の間隔よりも著しく小さい場合、等分配の法則は破綻する。エネルギー準位が滑らかな連続体を形成すると仮定することは、 上記の等分配定理の導出 に必要な近似としては不十分であるため、等分配はもはや成立しない。[ 6 ] [ 10 ] 歴史的に、古典的な等分配定理が比熱 と黒体放射を説明できなかったことは、物質と放射に関する新しい理論、すなわち 量子力学 と量子場の理論 の必要性を示す上で決定的な役割を果たした。[ 12 ]
図10. 量子力学的振動子(赤で表示)の平均エネルギーを温度の関数として対数プロットしたlog-logプロット。比較のために、等分配定理によって予測される値を黒で示している。高温では両者はほぼ完全に一致するが、低温(k B T ≪ hν) では量子力学的値ははるかに急速に減少する。これにより、紫外線カタストロフィー の問題が解消される。つまり、与えられた温度において、高周波モード(hν ≫ k B T )のエネルギーはほぼゼロとなる。 等分配の崩壊を説明するために、古典的なケースで上で議論した単一の(量子)調和振動子の平均エネルギーを考えてみましょう。無関係な零点エネルギー 項は確率分布に含まれる指数関数から因数分解できるため無視すると、量子調和振動子のエネルギー準位はE n = nhν で与えられます。ここで、h はプランク定数 、ν は振動子の基本周波数、 nは整数です。与えられたエネルギー準位が 正準集団 に占められる確率は、そのボルツマン因子 で与えられます。 ここでβ = 1/ k B T であり、分母Z は分配関数 です。ここでは等比級数です。 P ( E n ) = e − n β h ν Z , {\displaystyle P(E_{n})={\frac {e^{-n\beta h\nu }}{Z}},} Z = ∑ n = 0 ∞ e − n β h ν = 1 1 − e − β h ν . {\displaystyle Z=\sum _{n=0}^{\infty }e^{-n\beta h\nu }={\frac {1}{1-e^{-\beta h\nu }}}.}
その平均エネルギーは次のように与えられる。 ⟨ H ⟩ = ∑ n = 0 ∞ E n P ( E n ) = 1 Z ∑ n = 0 ∞ n h ν e − n β h ν = − 1 Z ∂ Z ∂ β = − ∂ log Z ∂ β . {\displaystyle \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 に式を代入すると最終結果が得られる[ 10 ] ⟨ H ⟩ = h ν e − β h ν 1 − e − β h ν . {\displaystyle \langle H\rangle =h\nu {\frac {e^{-\beta h\nu }}{1-e^{-\beta h\nu }}}.}
高温では、熱エネルギーk B T が エネルギー準位間の間隔hν よりもはるかに大きいため、指数関数引数βhν は1よりはるかに小さくなり、平均エネルギーはk B T となり、等分配定理と一致します(図10)。しかし、低温ではhν ≫ k B T となり、平均エネルギーはゼロになります。つまり、高周波エネルギー準位は「凍結」されます(図10)。別の例として、水素原子の内部励起電子状態は、室温でのガスとしての比熱には寄与しません。これは、熱エネルギーk B T (およそ0.025 eV )が、最低電子エネルギー準位と次に高い電子エネルギー準位間の間隔(およそ10 eV)よりもはるかに小さいためです。
エネルギー準位間隔が熱エネルギーよりもはるかに大きい場合は、同様の考慮が当てはまります。この考え方は、マックス・プランク やアルベルト・アインシュタインなどが 黒体放射 の紫外線カタストロフィー を解決するために使用しました。[ 52 ] このパラドックスは、閉じた容器内に電磁場 の独立したモードが無数にあり、それぞれが調和振動子として扱うことができるために生じます。各電磁モードの平均エネルギーがk B T であれば、容器内のエネルギーは無限大になります。[ 52 ] [ 53 ] しかし、上記の考え方によれば、高周波モードの平均エネルギーはν が 無限大になるとゼロになります。さらに、モードにおけるエネルギーの実験的分布を記述するプランクの 黒体放射の法則も、同じ考え方から導かれます。[ 52 ]
その他のより微妙な量子効果は、同一粒子 や連続対称性 など、等分配の修正につながる可能性があります。同一粒子の効果は、非常に高い密度と低温で支配的になることがあります。たとえば、金属の価電子は、数 電子ボルト の平均運動エネルギーを持つことができますが、これは通常、数万ケルビンの温度に相当します。パウリの排他原理により古典的なアプローチが無効になるほど密度が高くなるこのような状態は、 縮退フェルミオン気体 と呼ばれます。このような気体は、白色矮星 と中性子星 の構造にとって重要です。低温では、ボーズ・アインシュタイン凝縮 (多数の同一粒子が最低エネルギー状態を占める)のフェルミオン類似体が形成される可能性があり、このような 超流動電子が 超伝導 の原因です。
参照
注記と参考文献 ^ ストーン、A.ダグラス、「アインシュタインと量子」、第13章「凍結した振動」、2013年。ISBN 978-0691139685 ^ 「equi-」 。オンライン語源辞典。 2008年12月20 日 閲覧 ^ "partition" . オンライン語源辞典 . 2008年12月20日 閲覧。 。^ a b Kundt, A ; Warburg E (1876). 「Über die specifische Wärme des Quecksilbergases (水銀ガスの比熱について)」 Annalen der Physik (ドイツ語). 157 (3): 353– 369. 書誌コード : 1876AnP...233..353K . doi : 10.1002 / andp.18762330302 ^ ウラン濃縮に関するファクトシート 米国原子力規制委員会。2007年4月30日アクセス^ a b c d e f g h i j k l 英国パスリア (1972)。 統計力学 。ペルガモンプレス。ページ 43–48、73–74。ISBN 0-08-016747-0 。^ Cavanagh J、Fairbrother WJ、Palmer AG 3rd、Skelton NJ、Rance M (2006). タンパク質NMR分光法:原理と実践 (第2版). Academic Press. ISBN 978-0-12-164491-8 。^ Cantor, CR; Schimmel PR (1980). 生物理化学 第2部 生物学的構造と機能の研究のための技術 WH Freeman. ISBN 978-0-7167-1189-6 。^ a b c Goldstein, H (1980). 古典力学 (第2版). Addison-Wesley. ISBN 0-201-02918-9 。^ a b c d e f g h i Huang, K (1987). 統計力学 (第2版). John Wiley and Sons. pp. 136– 138. ISBN 0-471-81518-7 。^ a b Mandl, F (1971). 統計物理学 . John Wiley and Sons. pp. 213–219 . ISBN 0-471-56658-6 。^ a b c d ペイス、A (1982). 主は微妙なり 。オックスフォード大学出版局. ISBN 0-19-853907-X 。^ a b Tolman, RC (1918). 「エネルギー分配の一般理論と量子論への応用」 (PDF) . Physical Review . 11 (4): 261– 275. Bibcode : 1918PhRv...11..261T . doi : 10.1103 /PhysRev.11.261 ^ Miedl M, Garcia M, Bamforth C (2005). 「モデルビールシステムにおけるヘイズ形成」. J. Agric. Food Chem . 53 (26): 10161–5 . Bibcode : 2005JAFC...5310161M . doi : 10.1021/jf0506941 . PMID 16366710 . ^ Mason, M; Weaver W (1924). 「流体中における微粒子の沈降」. Physical Review . 23 (3): 412– 426. Bibcode : 1924PhRv...23..412M . doi : 10.1103/PhysRev.23.412 . ^ Brush, SG (1976). 『熱と呼ぶべき運動』第1巻 . アムステルダム: 北ホラント. pp. 134– 159. ISBN 978-0-444-87009-4 。ブラッシュ、SG(1976年)。『熱と呼ぶべき運動』第2巻 。アムステルダム:北ホラント。336 ~ 339ページ。ISBN 978-0-444-87009-4 。 ウォーターストン, JJ (1846). 「運動状態にある自由弾性分子からなる媒質の物理学について」 Proc . R. Soc. Lond . 5 : 604. doi : 10.1098/rspl.1843.0077 (抄録のみ)。全文掲載:Waterston, JJ; Rayleigh, L. (1893). 「運動状態にある自由かつ完全弾性の分子からなる媒体の物理学について」 . Philosophical Transactions of the Royal Society . A183 : 1– 79. Bibcode : 1892RSPTA.183....1W . doi : 10.1098/rsta.1892.0001 . JSハルデーン編(1928年)の 再版。ジョン・ジェームズ・ウォーターストンの科学論文集 。エディンバラ:オリバー&ボイド。 ウォーターストン、JJ (1843)『精神機能についての考察 』(彼の論文集 、3、167、183 に再録)ウォーターストン判事 (1851年) 『英国協会報告書 』 21 :6。 ウォーターストンの重要な論文は1845年に執筆され、王立協会 に提出されました。協会は彼の論文の出版を拒否した後、原稿の返却も拒否し、それをファイルに保管しました。原稿は1891年にレイリー卿 によって発見され、彼は元の査読者がウォーターストンの研究の重要性を認識していなかったと批判しました。ウォーターストンは1851年に自身のアイデアを公表したため、等分割定理の最初のバージョンを発表した点でマクスウェルよりも優先権を有しています。^ マクスウェル、JC (2003). 「気体の動的理論の図解」. WDニーヴン編. ジェームズ・クラーク・マクスウェルの科学論文集 . ニューヨーク: ドーバー. 第1巻, pp. 377–409. ISBN 978-0-486-49560-6 。 1859年9月21日、アバディーンで開催された英国協会の会合でマクスウェル教授によって朗読されました^ ボルツマン、L (1871)。 「Einige allgemeine Sätze über Wärmegleichgewicht (熱平衡に関するいくつかの一般的な記述)」。 ウィーン・ベリヒテ (ドイツ語)。 63 : 679–711 . この予備研究において、ボルツマンは、システムが外部の調和力によって作用されるとき、平均総運動エネルギーが平均総位置エネルギーに等しいことを示しました。^ ボルツマン、L (1876)。 「Über die Natur der Gasmoleküle (気体分子の性質について)」。 ウィーン・ベリヒテ (ドイツ語)。 74 : 553–560 . ^ a b McQuarrie, DA (2000). 統計力学 (改訂第2版). University Science Books. pp. 91–128 . ISBN 978-1-891389-15-3 。^ プチ、AT ; デュロン PL (1819)。 「Recherches sur quelques point importants de la théorie de la Chaleur (熱理論の重要な点に関する研究)」 。 Annales de Chimie et de Physique (フランス語)。 10 : 395–413 . ^ Dewar, J (1872). 「高温における炭素の比熱」. Philosophical Magazine . 44 : 461. ウェーバー、HF (1872)。「Die specifische Wärme des Kohlenstoffs (炭素の比熱)」 。Annalen der Physik (ドイツ語)。147 (10): 311–319 。ビブコード : 1872AnP...223..311W 。土井 :10.1002/andp.18722231007 。ウェーバー、HF (1875)。「Die specifische Wärmen der Elemente Kohlenstoff, Bor und Silicium (元素状炭素、ホウ素、シリコンの比熱)」 。Annalen der Physik (ドイツ語)。154 ( 3): 367–423、553–582 。ビブコード : 1875AnP... 230..367W 。土井 :10.1002/andp.18752300307 。^ デ・ラ・リヴェ、A;マーセット F (1840)。 「Quelques recherches sur la chaleur spécifique (比熱に関するいくつかの研究)」 。 Annales de Chimie et de Physique (フランス語)。 75 。マソン: 113–144 . ハイバーサス州レグノー (1841 年)。「Recherches sur la Chaleur spécifique des corps simples et des corps composés (deuxième Mémoire) (単純体と複合体の比熱の研究)」 。アナール・ド・シミーとフィジーク 。 (3me Série) (フランス語)。1 :129~ 207 1841 年 1 月 11 日に科学アカデミーで朗読。ウィガンド、A (1907)。 「Über Temperaturabhängigkeit der spezifischen Wärme fester Elemente (固体の比熱の温度依存性について)」 。 Annalen der Physik (ドイツ語)。 22 (1): 99–106 。 ビブコード : 1906AnP...327...99W 。 土井 : 10.1002/andp.19063270105 。 ^ a b ヴュラー、A (1896)。 Lehrbuch der Experimentalphysik (実験物理学の教科書) (ドイツ語)。ライプツィヒ:トイブナー。 Vol. 2、507ff。 ^ a b オイケン、A (1912)。 「Die Molekularwärme des Wasserstoffs beitiefen Tempen (低温における水素の分子比熱)」。 Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften (ドイツ語)。 1912 : 141–151 . ^ a b マクスウェル、JC (1890). 「物体の分子構成の動的証拠について」. WDニーヴン編. 『ジェームズ・クラーク・マクスウェルの科学論文集』 . ケンブリッジ大学出版局. 第2巻, pp.418–438. ISBN 0-486-61534-0 . ASIN B000GW7DXY. 1875年2月18日、マクスウェル教授が化学協会で行った講演^ a b Kittel, C (1996). 固体物理学入門 . ニューヨーク: John Wiley and Sons. pp. 151– 156. ISBN 978-0-471-11181-8 。^ Boltzmann, L (1895). 「気体理論に関するいくつかの問題について」 . Nature . 51 ( 1322): 413–415 . Bibcode : 1895Natur..51..413B . doi : 10.1038/051413b0 . S2CID 4037658 ^ トムソン, W (1904). ボルチモア講演集 . ボルチモア: ジョンズ・ホプキンス大学出版局. Sec. 27. ISBN 0-8391-1022-7 。1987年にMITプレスから『ケルビンのボルチモア講義と現代理論物理学:歴史的および哲学的視点』 (ロバート・カーゴン、ピーター・アチンシュタイン編)として再出版されました。ISBN 978-0-262-11117-1 ^ レイリー、JWS (1900). 「 運動エネルギーの分配の法則」 . 哲学雑誌 . 49 (296): 98–118 . doi : 10.1080/14786440009463826 ^ アインシュタイン、A (1906)。 「Die Plancksche Theorie der Strahlung und die Theorie der spezfischen Wärme (プランク輻射理論と比熱理論)」 。 Annalen der Physik (ドイツ語)。 22 (1): 180–190 。 Bibcode : 1906AnP...327..180E 。 土井 : 10.1002/andp.19063270110 。 アインシュタイン、A (1907)。「Berichtigung zu meiner Arbeit: 'Die Plancksche Theorie der Strahlung und die Theorie der spezfischen Wärme' (前回の記事の訂正)」 。Annalen der Physik (ドイツ語)。22 (4): 800。Bibcode : 1907AnP...327..800E 。土井 :10.1002/andp.19073270415 。S2CID 122548821 。アインシュタイン、A (1911)。「Eine Beziehung zwischen dem elastischen Verhalten and der spezifischen Wärme bei festen Körpern mit einatomigem Molekül (単一原子分子による固体の弾性挙動と比熱との関係)」 。Annalen der Physik (ドイツ語)。34 (1): 170–174 。Bibcode : 1911AnP...339..170E 。土井 :10.1002/andp.19113390110 。S2CID 122512507 。アインシュタイン、A (1911)。「Bemerkung zu meiner Arbeit: 'Eine Beziehung zwischen dem elastischen Verhalten and der spezifischen Wärme bei festen Körpern mit einatomigem Molekül' (前回の記事へのコメント)」 。Annalen der Physik (ドイツ語)。34 (3): 590。Bibcode : 1911AnP...339..590E 。土井 :10.1002/andp.19113390312 。アインシュタイン、A (1911)。「Elementare Betrachtungen über die thermische Molekularbewegung in festen Körpern (固体中の分子の熱運動に関する初等観察)」 。Annalen der Physik (ドイツ語)。35 (9): 679–694 。Bibcode : 1911AnP...340..679E 。土井 :10.1002/andp.19113400903 。^ ネルンスト、W (1910)。 「Untersuchungen über die spezfische Wärme beitiefen Tempen. II. (低温における比熱の調査)」。 Sitzungsberichte der Königlich Preussischen Akademie der Wissenschaften (ドイツ語)。 1910 : 262–282 . ^ ヘルマン、アルミン (1971)。 量子論の創世記 (1899–1913) (原題: Frühgeschichte der Quantentheorie (1899–1913) 、クロード・W・ナッシュ編訳)。マサチューセッツ州ケンブリッジ:MIT Press。 124–145 ページ 。 ISBN 0-262-08047-8 LCCN 73151106 ^ a b c トルマン, RC (1938). 統計力学の原理 . ニューヨーク: ドーバー出版. pp. 93– 98. ISBN 0-486-63896-0 。^ Clausius, R (1870). 「温度変化の機械的衝撃に対する応答」 . Annalen der Physik (ドイツ語). 141 (9): 124– 130. 書誌コード : 1870AnP...217..124C . doi : 10.1002/andp.18702170911 クラウジウス, RJE (1870). 「熱に応用可能な力学的定理について」.哲学雑誌 . シリーズ4. 40 : 122–127 .^ a b Vu-Quoc, L., Configuration integral (statistical mechanics) , 2008. このウィキサイトはダウンしています。2012年4月28日のウェブアーカイブにあるこの記事を 参照してください。 ^ McQuarrie, DA (2000). 統計力学 (改訂第2版). University Science Books. pp. 254–264 . ISBN 978-1-891389-15-3 。^ a b Tolman , RC (1927). 統計力学とその物理学および化学への応用 . Chemical Catalog Company. pp. 76–77 ^ Terletskii, YP (1971). 統計物理学 (N. Fröman編訳)アムステルダム: 北ホラント. pp. 83– 84. ISBN 0-7204-0221-2 LCCN 70157006 ^ Collins, GW (1978). The Virial Theorem in Stellar Astrophysics . Pachart Press. Bibcode : 1978vtsa.book ..... C ^ チャンドラセカール, S. (1939). 『恒星構造研究入門 』シカゴ: シカゴ大学出版局. pp. 49– 53. ISBN 0-486-60413-6 。^ クルガノフ、V (1980). 『先端天体物理学入門 』 ドルドレヒト、オランダ:D. ライデル社. pp . 59–60 , 134–140 , 181–184 ^ Chiu, HY (1968). 恒星物理学, 第1巻 . マサチューセッツ州ウォルサム: Blaisdell Publishing. LCCN 67017990 . ^ ノイズ, RW (1982). 『太陽、われらの星 』 ケンブリッジ, MA: ハーバード大学出版局. ISBN 0-674-85435-7 。^ キャロル、ブラッドリー・W.、オストリー、デール・A. (1996). 『現代恒星天体物理学入門』 マサチューセッツ州レディング:アディソン・ウェスリー. ISBN 0-201-59880-9 。^ Jeans, JH (1902). 「球状星雲の安定性」. Philosophical Transactions of the Royal Society A . 199 ( 312– 320): 1– 53. Bibcode : 1902RSPTA.199....1J . doi : 10.1098/rsta.1902.0012 . ^ McQuarrie, DA (2000). 統計力学 (改訂第2版). University Science Books. pp. 121–128 . ISBN 978-1-891389-15-3 。^ Callen, HB (1985). 熱力学と熱統計学入門 . ニューヨーク: John Wiley and Sons. pp. 375– 377. ISBN 0-471-86256-8 。^ Arnold, VI ; Avez A (1957). Théorie ergodique des systèms dynamiques (フランス語). Gauthier-Villars, Paris. (英語版: Benjamin-Cummings, Reading, Mass. 1968). ^ a b ライヒル, LE (1998). 統計物理学の現代コース (第2版). ワイリー・インターサイエンス. pp. 326– 333. ISBN 978-0-471-59520-5 。^ Naplekov, Dmitry M.; Yanovsky, Vladimir V. (2023-02-28). 「等分配を欠く理想気体におけるエネルギー分布」 . Scientific Reports . 13 (1): 3427. Bibcode : 2023NatSR..13.3427N . doi : 10.1038/ s41598-023-30636-6 . ISSN 2045-2322 . PMC 9974969. PMID 36854979 ^ a b c アインシュタイン、A (1905)。 「Über einen die Erzeugung und Verwandlung des Lichtes betreffenden heuristischen Gesichtspunkt (光の生成と変換のヒューリスティック モデル)」 。 Annalen der Physik (ドイツ語)。 17 (6): 132–148 。 Bibcode : 1905AnP...322..132E 。 土井 : 10.1002/andp.19053220607 。 英語翻訳は Wikisourceから 入手できます。^ レイリー, JWS (1900). 「完全放射の法則に関する考察」 . 哲学雑誌 . 49 : 539–540 . Bibcode : 1900PMag...49..539R . doi : 10.1080/14786440009463878 .
さらに読む
外部リンク