レイリー・ジーンズの法則

レイリー・ジーンズの法則（レイリー・ジーンズのほうそく、英: Rayleigh–Jeans Law）は、黒体輻射におけるエネルギー密度の理論式の1つである。レイリー・ジーンズの公式とも呼ばれ、イギリスの物理学者であるレイリー卿とジェームズ・ジーンズに因む。レイリー・ジーンズの公式は黒体から熱放射される電磁波、すなわち輻射場のある温度におけるエネルギー密度のスペクトル分布を与える。輻射場を一次元調和振動子の集まりとして扱い、古典統計力学を適用することで導かれる。この公式が適用できるのは長波長（低振動数）の限られたスペクトル領域のみであり、不完全な理論式である。全スペクトル領域で成り立つ完全な理論式は、量子論に基づくプランクの公式によって与えられる。しかしながら、レイリー・ジーンズの公式は古典物理学の限界を浮き彫りにし、輻射の理論や現代物理学の発展に重要な役割を果たした。 溶鉱炉で熱した鉄が光を放射するように、物質を高温にしていくと物質原子から電磁波の熱輻射が生じる。逆に輻射された電磁波（輻射場）を受けた物質は、電磁波を吸収し、エネルギーは物質原子の熱運動に変わる。キルヒホッフの法則によれば、輻射場と熱平衡状態にある物質の輻射能と吸収能の比は、物質によらず、振動数と温度のみに依存する普遍的な関数となる。これは特に黒体輻射の場合と一致する。また、この関数は分光放射輝度であり、輻射場のエネルギー密度に比例する。黒体はすべての振動数の電磁波を完全に吸収する理想的な物体であるが、黒体輻射は空洞炉内での熱輻射（空洞輻射）の形で実現できる。空洞炉内で輻射場は壁の物質からの放射と吸収を通じて、熱平衡状態になる。 空洞炉内での輻射場が温度 T の熱平衡状態にあるとする。振動数が ν から ν+dνの間にある輻射場の単位体積あたりのエネルギー密度を u(ν,T)dν とすると、レイリー・ジーンズの公式は、 u ( ν , T ) = 8 π ν 2 k T c 3 {\displaystyle u(\nu ,T)={\frac {8\pi \nu ^{2}kT}{c^{3}}}} と表される。ここで、k は ボルツマン定数、c は 光速である。振動数 ν と波長 λ の関係 λ=c/ν から、 波長が λ から λ+dλ の間にある輻射場の単位体積あたりのエネルギー密度を u(λ,T)dλ とすると、レイリー・ジーンズの公式は u ( λ , T ) = 8 π k T λ 4 {\displaystyle u(\lambda ,T)={\frac {8\pi kT}{\lambda ^{4}}}} と表すこともできる。 この式は 輻射場は空洞炉の壁の物質と熱平衡にある電磁波（光）であり、波として振る舞う 電磁波の全てのモードに対してエネルギー等配分の法則が成り立つ という2つの古典物理学的な仮定から導出される。しかし、上式は振動数が低い（波長が長い）領域では実験と良く一致するが、振動数が高く（波長が短く）なればなるほど実験結果とズレが大きくなる。また、放射の全エネルギー密度 ∫ 0 ∞ u ( ν , T ) d ν ( = ∫ 0 ∞ u ( λ , T ) d λ ) {\displaystyle \int _{0}^{\infty }u(\nu ,T)d\nu \left(=\int _{0}^{\infty }u(\lambda ,T)d\lambda \right)} を計算しようとすると発散して無限大になってしまう。このことは、黒体放射の問題に対して古典物理学が破綻することを端的に示している。 レイリー・ジーンズの法則は、レイリー卿が1900年に最初に発表した。 レイリーは、1900年に『完全輻射の法則についての注意』と題する2ページの短い論文の中で、レイリー・ジーンズの公式の原型となる式を提案した。この中で、ヴィーンの公式が、長波長、高温の領域では温度が上昇しても u(λ,T) が一定の極限値に近づき、それ以上は増加しない問題を指摘した。レイリーは空気の振動である音との類似性から、振動数 ν と ν+dν の間にある輻射場の振動子のモード数が ν2dν に比例するとした。さらに振動子の各モードにマクスウェル＝ボルツマンらによるエネルギー等分配則を適用することで u ( ν , T ) d ν ∝ T ν 2 d ν {\displaystyle u(\nu ,T)d\nu \varpropto T\nu ^{2}d\nu } もしくは波長で表した u ( λ , T ) d λ ∝ T λ 4 d λ {\displaystyle u(\lambda ,T)d\lambda \varpropto {\frac {T}{\lambda ^{4}}}d\lambda } という結果を得た。また、明確な根拠を示さずに指数因子を加えた u ( λ , T ) d λ ∝ T λ 4 e − c λ T d λ {\displaystyle u(\lambda ,T)d\lambda \varpropto {\frac {T}{\lambda ^{4}}}e^{-{\frac {c}{\lambda T}}}d\lambda } を提案している。これはヴィーンの公式で λ−5 の項をλ−4T に置き換えた修正に相当する。但し、これらの結果については、具体的な比例係数を求めておらず、議論も不完全なものであった。その後、1905年にレイリーは係数まで含めた形で導出を行ったが、係数が正しい結果と8倍違っていた。同年、ジェームズ・ジーンズが係数に誤りがあることを指摘した。1905年、レイリーは論文『気体とエーテルの力学理論』において、レイリー・ジーンズの公式の係数の具体的な導出を行った。彼は電磁波を伝播させる媒質として考えられていたエーテルの振動について考察し、弦の振動との類似性から波長 λ の持つ振動のモード数を計算した。そして、これらのモードにエネルギー等分配則を適用することで、結果として係数に8倍の誤りを含んでいたものの、レイリー・ジーンズの公式を得た。 レイリー・ジーンズの法則は、波長が短いときに実験結果と合わない。逆に、ヴィーンの公式（1896年）は波長が短いときに実験と一致するが長波長領域では実験と合わない。ヴィーンの公式を改良したプランクの公式（1900年） u ( λ , T ) = 8 π h c λ 5   1 e h c / λ k T − 1 {\displaystyle u(\lambda ,T)={\frac {8\pi hc}{\lambda ^{5}}}~{\frac {1}{e^{hc/\lambda kT}-1}}} は、波長が短いときも長いときも実験結果とよく合う。ここで h は プランク定数、c は光速度である。また、放射の全エネルギー密度も有限の値になる。 プランクの公式は、波長の長い時、または高温の時レイリー・ジーンズの式に近づく。 1905年、アルベルト・アインシュタインは『光の生成と変換に関する、ひとつの発見法的観点について』と題した論文で、プランクの公式についての議論を行い、別の方法でレイリー・ジーンズの公式を導いた。プランクは輻射場と熱平衡状態にある荷電共鳴子（ヘルツの共鳴子）についてのエントロピーの議論からプランクの公式を導いていたが、共鳴子の平均エネルギーEνと輻射場のエネルギー密度の間には、 E ν ¯ = c 3 8 π ν 2 u ( ν , T ) {\displaystyle {\overline {E_{\nu }}}={\frac {c^{3}}{8\pi \nu ^{2}}}u(\nu ,T)} の関係が成りたつ。アインシュタインはこの共鳴子について、エネルギー等分配則を適用することでレイリー・ジーンズの公式が得られることを示した。また、レイリー・ジーンズの公式では、振動数について積分した全エネルギー密度が無限大になることを指摘した。 電荷や電流が分布しない真空中では、電磁場は無限個の一次元調和振動子の集まりとして扱える。この調和振動子の集団が温度 T の熱平衡状態にあるとき、古典統計力学を適用するとレイリー・ジーンズの法則が導かれる。古典統計力学のエネルギー等分配則によれば、振動数 ν の一次元調和振動子のエネルギー εν の熱平均は、 ⟨ ϵ ν ⟩ = k T {\displaystyle \langle \epsilon _{\nu }\rangle =kT} で与えられる。したがって、振動数が ν ∼ν+dν にある電磁波のモード数を求め、それに kT を乗じれば、考えている空間全体での振動数が ν ∼ν+dν にある輻射場のエネルギーとなる。これを空間の単位体積当たりに直せば、エネルギー密度となる。 振動数が ν から ν+dν の間にある電磁波の単位体積あたりのモード数は、モード密度 m(ν) によって、m(ν)dν で与えられる。したがって、振動数が ν から ν+dν の間にある輻射場のエネルギー密度は、振動数 ν の振動子の平均エネルギー〈εν〉=kT と単位体積あたりのモード数 m(ν)dν の積により u ( ν , T ) d ν = m ( ν ) ⟨ ϵ ν ⟩ d ν = m ( ν ) k T d ν {\displaystyle u(\nu ,T)d\nu =m(\nu )\langle \epsilon _{\nu }\rangle d\nu =m(\nu )kTd\nu } と表される。電磁波の波数の取りうる条件と、電磁波が横波で偏光の2自由度が存在することを考慮すると、モード密度は m ( ν ) = 8 π ν 2 c 3 {\displaystyle m(\nu )={\frac {8\pi \nu ^{2}}{c^{3}}}} となる。その結果、 u ( ν , T ) d ν = 8 π ν 2 c 3 k T d ν {\displaystyle u(\nu ,T)d\nu ={\frac {8\pi \nu ^{2}}{c^{3}}}kTd\nu } が得られる。 しかしながら、古典統計力学を適用するアプローチは正しくない。無限個の振動子が kT のエネルギーを持つため、振動数について、全振動数の領域 ν=0∼∞ について積分すると発散する。すなわち、紫外破綻の問題が生じる原因となっている。 u(ν,T) を全振動数（または u(λ,T) を全波長）の領域に亘って積分すれば、輻射場の全エネルギー密度を求めることができる。レイリー・ジーンズの法則に基づけば、 U V = ∫ 0 ∞ u ( ν , T ) d ν = 8 π k T c 3 ∫ 0 ∞ ν 2 d ν = ∞ ( T ≠ 0 ) {\displaystyle {\begin{aligned}{\frac {U}{V}}&=\int _{0}^{\infty }u(\nu ,T)d\nu ={\frac {8\pi kT}{c^{3}}}\int _{0}^{\infty }\nu ^{2}d\nu \\&=\infty \quad (T\neq 0)\end{aligned}}} となり、温度が T=0 以外の場合には、無限大となる。すなわち、真空中に無限大のエネルギーが存在すると誤った結果を導く。このことは、アルベルト・アインシュタインによって1905年の論文で指摘された。この問題をオランダの物理学者ポール・エーレンフェストは、紫外破綻の問題と呼んだ。 プランクの公式では、エネルギー密度は u ( ν , T ) = 8 π h ν 3 c 3 1 e h ν / k T − 1 {\displaystyle u(\nu ,T)={\frac {8\pi h\nu ^{3}}{c^{3}}}{\frac {1}{e^{h\nu /kT}-1}}} で与えられる。プランクの公式は高温、または低振動度数で hν<<kT が満たされる場合については、 u ( ν , T ) ∼ 8 π ν 2 k T c 3 ( h ν k T << 1 ) {\displaystyle u(\nu ,T)\sim {\frac {8\pi \nu ^{2}kT}{c^{3}}}\quad {\biggl (}{\frac {h\nu }{kT}}<<1{\biggr )}} とレイリー・ジーンズの公式で近似される。また、プランクの公式では、u(ν,T) を全振動数の領域に亘って積分した全エネルギー密度は U V = ∫ 0 ∞ u ( ν , T ) d ν = 8 π 5 15 ( k T ) 4 ( h c ) 3 ∝ T 4 {\displaystyle {\frac {U}{V}}=\int _{0}^{\infty }u(\nu ,T)d\nu ={\frac {8\pi ^{5}}{15}}{\frac {(kT)^{4}}{(hc)^{3}}}\varpropto T^{4}} と有限の値をとる。 レイリー・ジーンズの公式とプランクの公式のエネルギー密度は両者ともに u ( ν , T ) = m ( ν ) ⟨ ϵ ν ⟩ {\displaystyle u(\nu ,T)=m(\nu )\langle \epsilon _{\nu }\rangle } と、電磁波の体積当たりのモード数であるモード密度 m(ν) と振動数 ν の調和振動子の温度 T での平均エネルギー〈εν〉 の積となっている。モード密度は m ( ν ) = 8 π ν 2 c 3 {\displaystyle m(\nu )={\frac {8\pi \nu ^{2}}{c^{3}}}} と共通であるが、レイリー・ジーンズの公式が古典統計力学に基づく ⟨ ϵ ν ⟩ = k T {\displaystyle \langle \epsilon _{\nu }\rangle =kT} を適用しているのに対し、プランクの公式は光子の量子統計性に基づく ⟨ ϵ ν ⟩ = h ν e h ν / k T − 1 {\displaystyle \langle \epsilon _{\nu }\rangle ={\frac {h\nu }{e^{h\nu /kT}-1}}} を適用している。 量子論の枠組みでは輻射場は量子化された場で表される。電荷や電流が分布しない真空中では、量子化された電磁場の各モードが量子力学的な調和振動子に対応する。この量子力学的な調和振動子は光子を表しており、振動数 ν の光子のエネルギーは ϵ ν ( n ) = ( n + 1 2 ) h ν ( n = 0 , 1 , 2 , ⋯ ) {\displaystyle \epsilon _{\nu }(n)={\biggl (}n+{\frac {1}{2}}{\biggr )}h\nu \quad (n=0,1,2,\cdots )} という離散的な値をとる。各 n に対応するエネルギー状態は調和振動子の第 n 励起状態であるが、これは輻射場の振動数 ν のモードに n 個の光子が励起した状態である。 n=0 での基底状態でエネルギーは有限の値 hν/2 をとるが、これは零点エネルギーと呼ばれる。観測に掛かるのはエネルギーの基底状態からの差であり、零点エネルギーの効果は以降の議論で無視できる。零点エネルギーを無視すると、振動数 ν のモードの光子は、hν, 2hν, 3hν… という hν の整数倍のエネルギーのみをとる。 温度 T の平衡状態でエネルギーεν(n)=nhν を持つ状態にある確率はボルツマン因子を用いて、 P ( n ) = exp ⁡ ( − n h ν / k T ) ∑ n = 1 ∞ exp ⁡ ( − n h ν / k T ) = exp ⁡ ( − n h ν / k T ) Z {\displaystyle P(n)={\frac {\exp(-nh\nu /kT)}{\sum _{n=1}^{\infty }\exp {(-nh\nu /kT)}}}={\frac {\exp {(-nh\nu /kT)}}{Z}}} で与えられることから、その期待値は ⟨ ϵ ν ⟩ = ∑ n = 1 ∞ ϵ ν ( n ) P ( n ) = h ν e h ν / k T − 1 {\displaystyle \langle \epsilon _{\nu }\rangle =\sum _{n=1}^{\infty }\epsilon _{\nu }(n)P(n)={\frac {h\nu }{e^{h\nu /kT}-1}}} となり、プランクの公式を与える結果が得られる。離散的なエネルギーの間隔をゼロとし、エネルギーが連続的であるとする極限 hν → 0 では ⟨ ϵ ν ⟩ → k T {\displaystyle \langle \epsilon _{\nu }\rangle \rightarrow kT} となり、レイリー・ジーンズの公式を与える結果になる。 輻射場で満たされた空洞炉内が電荷や電流が存在しない真空であるとする。このとき、電場と磁場は波動方程式を満たし、これは電磁波の伝播を記述する。この電磁波は、波数ベクトル k と偏光の2つの自由度 γ=1,2 で特徴付けられるモードに展開できる。このとき、電磁波のエネルギーは無限個の正準座標{Qk,γ}と正準運動量{Pk,γ}の組を用いて、調和振動子の集まりとして表すことができる。 電場と磁場は、ベクトルポテンシャルによって、記述することができる。クーロンゲージの条件(divA=0)を適用するとベクトルポテンシャル A(r,t) は波動方程式を満たし、 A ( r , t ) = ∑ k ( A k e i ( k ⋅ r − ω k t ) + A k ∗ e − i ( k ⋅ r − ω k t ) ) {\displaystyle \mathbf {A} (\mathbf {r} ,t)=\sum _{\mathbf {k} }(\mathbf {A} _{\mathbf {k} }e^{i(\mathbf {k} \cdot \mathbf {r} -\omega _{k}t)}+\mathbf {A} _{\mathbf {k} }^{\ast }e^{-i(\mathbf {k} \cdot \mathbf {r} -\omega _{k}t)})} と展開できる。但し、A(r,t) には周期境界条件を課したほか、 ω k = c | k | = c k {\displaystyle \omega _{k}=c|\mathbf {k} |=ck} とした。一方、空洞炉内の電磁波のエネルギーは電場 E(r,t)、磁場 B(r,t) によって U = 1 2 ∫ { ϵ 0 E 2 ( r , t ) + 1 μ 0 B 2 ( r , t ) } d 3 r = 2 ϵ 0 V ∑ k ω k 2 ( A k ⋅ A k ∗ ) {\displaystyle {\begin{aligned}U&={\frac {1}{2}}\int \left\{\epsilon _{0}\mathbf {E} ^{2}(\mathbf {r} ,t)+{\frac {1}{\mu _{0}}}\mathbf {B} ^{2}(\mathbf {r} ,t)\right\}d^{3}\mathbf {r} \\&=2\epsilon _{0}V\sum _{\mathbf {k} }\omega _{k}^{\,2}(\mathbf {A} _{\mathbf {k} }\cdot \mathbf {A} _{\mathbf {k} }^{\ast })\end{aligned}}} で与えられる。ここで実の正準変数 Qk(t), Pk(t) の組を Q k ( t ) = ϵ 0 V ( A k e − i ω k t + A k ∗ e i ω k t ) P k ( t ) = Q ˙ k ( t ) = − i ω k ϵ 0 V ( A k e − i ω k t − A k ∗ e i ω k t ) {\displaystyle {\begin{aligned}\mathbf {Q} _{\mathbf {k} }(t)&={\sqrt {\epsilon _{0}V}}(\mathbf {A} _{\mathbf {k} }e^{-i\omega _{k}t}+\mathbf {A} _{\mathbf {k} }^{\ast }e^{i\omega _{k}t})\\\mathbf {P} _{\mathbf {k} }(t)&={\dot {\mathbf {Q} }}_{\mathbf {k} }(t)\\&=-i\omega _{\mathbf {k} }{\sqrt {\epsilon _{0}V}}(\mathbf {A} _{\mathbf {k} }e^{-i\omega _{k}t}-\mathbf {A} _{\mathbf {k} }^{\ast }e^{i\omega _{k}t})\end{aligned}}} で導入すると、U をこれらの正準変数で書き表したハミルトニアン H として H = ∑ k H k = ∑ k 1 2 ( P k 2 ( t ) + ω k 2 Q k 2 ( t ) ) {\displaystyle {\begin{aligned}H&=\sum _{\mathbf {k} }H_{\mathbf {k} }\\&=\sum _{\mathbf {k} }{\frac {1}{2}}(\mathbf {P} _{\mathbf {k} }^{\,2}(t)+\omega _{k}^{\,2}\mathbf {Q} _{\mathbf {k} }^{\,2}(t))\end{aligned}}} が得られる。Hk は調和振動子のハミルトニアンそのものである。 クーロンゲージの条件から波数ベクトル k は Ak と垂直である。 k に垂直で、互いに直交する2つの単位ベクトル ek,1, ek,2 を取れば、Akは A k = ∑ γ = 1 , 2 A k , γ e k , γ {\displaystyle \mathbf {A} _{\mathbf {k} }=\sum _{\gamma =1,2}A_{\mathbf {k} ,\gamma }\mathbf {e} _{\mathbf {k} ,\gamma }} と偏光の2自由度に対応する形に展開できる。Qk(t), Pk(t)についても、同様に ek,γで展開できる。その係数を Qk,γ(t), Pk,γ(t) とすれば、これらも正準変数である。ハミルトニアンは最終的に H = ∑ k , γ H k , γ = ∑ k , γ 1 2 ( P k , γ 2 ( t ) + ω k 2 Q k , γ 2 ( t ) ) {\displaystyle {\begin{aligned}H&=\sum _{\mathbf {k} ,\gamma }H_{\mathbf {k} ,\gamma }\\&=\sum _{\mathbf {k} ,\gamma }{\frac {1}{2}}(P_{\mathbf {k} ,\gamma }^{\,2}(t)+\omega _{k}^{\,2}Q_{\mathbf {k} ,\gamma }^{\,2}(t))\end{aligned}}} となるが、これは波数ベクトル k と偏光の2つの自由度 γ=1,2 で特徴付けられる調和振動子の集まりである。