プランクの法則

プランクの法則の最新ニュースをまとめて検索!

黒体放射スペクトル

より一般的な導入として、黒体も参照

プランクの法則(―のほうそく)とは物理学における黒体放射の公式であり、温度 T における黒体からの電磁輻射分光放射輝度を全波長領域において正しく説明することができる。プランクはこの法則の導出を考える中で、空洞の壁の振動子のエネルギーがエネルギー素量(energy element、現在はエネルギー量子とよばれる)ε = hν の整数倍になっていると仮定した。このエネルギーの量子仮説(量子化)はその後の量子力学の幕開けに大きな影響を与えた。

プランクの法則は周波数 ν と温度 T関数として

I(\nu,T) =\frac{2 h\nu^{3}}{c^2}\frac{1}{ e^{\frac{h\nu}{kT}}-1}

と表すことができる[1]

この関数は、放射面の単位面積立体角、周波数あたりの放射強度を表す。プランクの法則は放射強度の全立体角についての積分として u(ν,T) = πI(ν,T) という形で表すこともできる。さらに単位体積あたりのエネルギーとして u(ν,T) = 4πI(ν,T) / c と書くこともできる。

関数 I(ν,T)hν = 2.82kBT の位置にピークをもつ[2]。関数の値は高周波数においては指数関数的に、低周波数においては多項式的に減少する。

波長 λ の関数として、プランクの法則は(単位立体角あたり)

I'(\lambda,T) =\frac{2 hc^2}{\lambda^5}\frac{1}{ e^{\frac{hc}{\lambda kT}}-1}

という形であらわすことができる[3]

この関数は hc = 4.97λkBT の位置にピークをもつ。これはヴィーンの変位則でより一般的に用いられるピークである。

周波数範囲 12] または波長範囲 21] = [c / ν2,c / ν1] において放射される放射輝度はそれぞれの関数の積分として求められる。

\int_{\nu_1}^{\nu_2}I(\nu,T)\,d\nu=\int_{\lambda_2}^{\lambda_1}I'(\lambda,T)\,d\lambda

周波数が増加するとき波長は減少するため、2つの積分では上限・下限が入れ替わっている。

次の表に、数式の中に表れるそれぞれの記号の定義とSI単位cgs単位を示す。

記号 意味 SI単位系 cgs単位系
I, I' \, 分光放射輝度 または エネルギー(単位時間表面積立体角、周波数(波長)あたり) J⋅s-1⋅m-2⋅sr-1⋅Hz-1, または J⋅s-1⋅m-2⋅sr-1⋅m-1 erg⋅s-1⋅cm-2⋅Hz-1⋅sr-1, または erg⋅s-1⋅cm-2⋅sr-1⋅cm-1
\nu \, 周波数 ヘルツ (Hz) ヘルツ
\lambda \, 波長 メートル (m) センチメートル (cm)
T \, 黒体の温度 ケルビン (K) ケルビン
h \, プランク定数 ジュール⋅秒 (J⋅s) エルグ⋅秒 (erg⋅s)
c \, 光速 メートル毎秒 (m/s) センチメートル毎秒 (cm/s)
e \, 自然対数の底, 2.718281... 無次元量 無次元量
k \, ボルツマン定数 ジュール毎ケルビン (J/K) エルグ毎ケルビン (erg/K)

目次

[編集] 概要

波長と周波数は

\lambda = { c \over \nu }

という関係式によって結びついている[4]

プランクの法則は分光エネルギー密度に関して

u(\nu,T)    =   {  4 \pi \over c }   I(\nu,T)  =   \frac{8\pi h\nu^3 }{c^3}~\frac{1}{e^{\frac{h\nu}{kT}}-1}

と表すこともできる[5]。ここでエネルギー密度 u は単位体積、単位周波数あたりのエネルギーの単位(ジュール毎立方センチメートル毎ヘルツ)をもつ。この式を周波数について積分をすれば、全エネルギー密度を得る。黒体の輻射場は光子気体と考えることができ、その場合エネルギー密度は光子気体の熱力学変数の一つとなる。

また分光エネルギー密度を波長の関数として

u(\lambda,T) = {8\pi h c\over \lambda^5}{1\over e^{\frac{h c}{\lambda kT}}-1}

と表すこともできる。導出は後述する。

マックス・プランクは1900年(発行は1901年[6])にプランクの法則を生み出した。プランクの法則は最初、ヴィーンの放射法則より良い公式を得ようとする過程でみつけた公式だった。ヴィーンの放射法則はヴィルヘルム・ヴィーンが1896年に発表した公式であり、短波長(高周波数)においては実験データと一致するものの、長波長(低周波数)では一致しない。一方レイリー・ジーンズの法則(1900年に不完全な形でレイリーが発表した)は反対に長波長(低周波数)でよい一致を示す。プランクのみつけた前述の公式は、全ての波長領域において非常によく実験データと一致した。次にこの法則の導出を構築する過程で、プランクは物質中の荷電振動子の異なるモードについての電磁エネルギー分布を考えた。これらの振動子のエネルギーが離散的になっていると仮定したところ、プランクの法則を導出することができた。具体的には、エネルギーは振動数 ν に比例する単位エネルギー E、即ち

E = hν

の整数倍になっていなくてはならない。

プランクはこの量子化の仮定を、アルベルト・アインシュタイン光電効果の説明のために光子の存在を仮定するよりも5年に行っていた。この時点では、プランクは量子化は空洞壁面にあるであろう微小の共鳴子resonator、現在でいう原子)にのみ適用されるものであり、光それ自身が離散的なエネルギーの束や塊を伝播する性質を有しているとは仮定しなかった。更には、プランクはこの仮定にはなんら物理的重要性はなく、公式を導くための単なる数学的な道具に過ぎないと考えていた。

プランクの公式では黒体は全ての周波数の電磁波を放出するとしているが、これは非常に多数の光子が測定される実験でのみ実際に適用できる。例えば室温 (300 K) における表面積が1平方メートルの黒体は、1000年に一度程度しか可視領域の光子を放出せず、よって通常の実験などにおいては黒体は室温では可視光線を放出されないといっても差し支えない。実験データからプランクの法則を導出する際などのこの事実の重要性については[7]で議論されている。

この法則による影響として、プランクのエネルギー量子化の仮説とアインシュタインの光量子仮説は、ともに量子力学の発展における基礎となった。

[編集] 他の輻射法則との関係

プランクの法則は、以下にあげるように黒体輻射の法則と密接に関連している。

プランクの法則は高周波数(短波長)においては

I(\nu, T) \sim \frac{2h\nu^3}{c^2} e^{-\frac{h\nu}{kT}}
I'(\lambda, T) \sim \frac{2hc^2} {\lambda^5} e^{-\frac{hc}{\lambda kT}}

となり、ヴィーンの放射法則に漸近する。

また、低周波数(長波長)においては

I(\nu, T) \sim \frac{2\nu^2}{c^2}kT
I'(\lambda, T) \sim \frac{2c}{\lambda^4}kT

となり、レイリー・ジーンズの法則に漸近する。

また、プランクの法則の周波数(波長)についての積分

j^\star(T) = \int^\infty_0 I(\nu,T)d\nu \int d\Omega = \pi  \int^\infty_0 I(\nu,T)d\nu

= \sigma T^4, \quad \sigma = \frac{2\pi^5 k^4}{15 c^2 h^3}T^4

より、シュテファン=ボルツマンの法則が得られる。

さらに、\partial I'(\lambda, T) / \partial \lambda = 0 よりプランクの法則の分光放射輝度 I'(λ,T) が最大となる波長 λ を求めることにより、ヴィーンの変位則が得られる。

[編集] 導出

以下のプランクの法則の導出は[5]などでみられる。より一般的な導出は箱の中の気体 (en:gas in a box) の項目を参照。

伝導壁をもち電磁波で満たされた一辺の長さ L の立方体を考える。立方体の壁では、電場の平行成分と磁場の直行成分はあってはならない。箱の中の粒子 (particle in a box) の波動関数との類似により、場は周期的な関数の重ね合わせとして表される。壁に直行する3つの方向についての3つの波長 λ123

\lambda_{i} = \frac{2L}{n_{i}}

となる。ここで ni は整数である。ni のそれぞれの組について、2つの線型独立な解(モード)がある。量子論にしたがい、一つのモードのエネルギー準位は

E_{n_{1},n_{2},n_{3}}\left(r\right)=\left(r+\frac{1}{2}\right)\frac{hc}{2L}\sqrt{n_{1}^{2}+n_{2}^{2}+n_{3}^{2}} \qquad \mbox{(1)}

によって与えられる。

量子数 r はモードの中の光子数に対応している。ni のそれぞれの組の2つのモードはスピン1をもつ光子の2つの偏光状態に対応している。ここで注意すべきは、r = 0 においてもモードのエネルギーは零ではないことである。この電磁場の真空エネルギーはカシミール効果によるものである。これ以降、温度 T の箱の内部エネルギーを真空エネルギーとの相対値で計算してゆく。

統計力学に従い、特定のモードのエネルギー準位についての確率分布は

P_{r}=\frac{\exp\left(-\beta E\left(r\right)\right)}{Z\left(\beta\right)}

で与えられる。ここで

\beta\ \stackrel{\mathrm{def}}{=}\   1/\left(kT\right)

である。

分母 Z(β) は単モードの分配関数であり、Pr を正しく正規化する。

Z\left(\beta\right)=\sum_{r=0}^{\infty}\exp\left[-\beta E\left(r\right)\right]=\frac{1}{1-\exp\left[-\beta\varepsilon\right]}

ここで

\varepsilon\ \stackrel{\mathrm{def}}{=}\  \frac{hc}{2L}\sqrt{n_{1}^{2}+n_{2}^{2}+n_{3}^{2}}

を暗に定義した。これは単一光子のエネルギーである。あるモードにおける平均エネルギーは分配関数によって

\left\langle E\right\rangle=-\frac{d\log\left(Z\right)}{d\beta}=\frac{\varepsilon}{\exp\left(\beta\varepsilon\right)-1}

のように表される。

これはボース=アインシュタイン統計に従う粒子の場合の公式である。全光子数に制限がないため、化学ポテンシャルは零である。

箱の中の全エネルギーはあり得る全単一光子状態についての総和 \left\langle E\right\rangle に従う。これは L が無限大となる熱力学的極限において厳密に成り立つ。この極限では \varepsilon は連続となり、よって \left\langle E\right\rangle\varepsilon について積分することができる。この方法により箱の中の全エネルギーを計算するには、与えられたエネルギー範囲にどの程度の光子状態があるのかを評価する必要がある。今エネルギー \varepsilon\varepsilon + d\varepsilon の間にある単一光子状態の総数を g(\varepsilon)\,d\varepsilon と表すとする。ここで g(\varepsilon) は評価しようとする状態密度である。この場合には

U = \int_{0}^{\infty}\frac{\varepsilon}{\exp\left(\beta\varepsilon\right)-1}g(\varepsilon)\,d\varepsilon \qquad \mbox{(2)}

と書くことができる。

状態密度を計算するためには、等式(1)を

\varepsilon\ \stackrel{\mathrm{def}}{=}\  \frac{hc}{2L}n

と書き換える。ここで n はベクトルのノルム(長さ)

n=\sqrt{n_{1}^{2}+n_{2}^{2}+n_{3}^{2}}

である。

零以上の整数成分のベクトル \vec{n} について、それぞれ2つの光子状態がある。言い換えると、ある n-空間領域での光子状態の数は、その領域の体積の2倍である。d\varepsilon のエネルギー範囲は n-空間では dn= (2L/hc) d\varepsilon の厚さの殻に対応する。\vec{n} の要素は符号が正でなくてはならないため、この殻は丁度球の八分の一領域にわたる。よってエネルギー範囲 d\varepsilon にある光子状態の数 g(\varepsilon)\,d\varepsilon

g(\varepsilon)\,d\varepsilon=2\frac{1}{8}4\pi n^{2}\,dn=\frac{8\pi L^{3}}{h^{3}c^{3}}\varepsilon^{2}\,d\varepsilon

で与えられる。

この式を方程式 (2) に代入して

U =L^{3}\frac{8\pi}{h^{3}c^{3}}\int_{0}^{\infty}\frac{\varepsilon^{3}}{\exp\left(\beta\varepsilon\right)-1}\,d\varepsilon \qquad \mbox{(3)}

を得る。

この方程式から、周波数の関数 u(ν,T) または波長の関数 u(λ,T) として分光エネルギー密度を容易に導出することができる。

\frac{U}{L^3} = \int_0^{\infty}u(\nu,T)\, d\nu.

ここで

u(\nu,T) = {8\pi h\nu^3\over c^3}{1\over e^{h\nu/kT}-1}

である。この u(ν,T)黒体スペクトルとして知られる。これが単位周波数、単位体積あたりの分光エネルギー密度の関数である。

更に

\frac{U}{L^3} = \int_0^\infty u(\lambda,T)\, d\lambda

も導くことができる。ここで

u(\lambda,T) = {8\pi h c\over \lambda^5}{1\over e^{h c/\lambda kT}-1}

である。

これは同様に、単位波長、単位体積あたりの分光エネルギー密度の関数である。ボース気体とフェルミ気体の計算に現れるこの形の積分は多重対数関数によって表される。しかし今回の場合には閉形式の積分を初等関数を用いて表すことができる。方程式 (3) において

\varepsilon = kTx

と置換すると、積分変数を無次元量の割り算にすることができ

u(T) =\frac{8\pi (kT)^{4}}{(hc)^{3}} J

である。ここで J

J=\int_{0}^{\infty}\frac{x^{3}}{\exp\left(x\right)-1}\,dx = \frac{\pi^{4}}{15}

によって与えられる。この計算の証明は#附則において後述する。

よって箱の中の全電磁エネルギーは

{U\over V} = \frac{8\pi^5(kT)^4}{15 (hc)^3}

によって与えられる。ここで V = L3 は箱の体積である。(註: これはシュテファン=ボルツマンの法則ではない。即ち、黒体によって放射される全エネルギーではない。シュテファン=ボルツマンの法則の項目を参照。)放射は全方向 (4π) に等しく起き、またその伝播速度は光速 (c)であるため、分光放射輝度は

I(\nu,T) = \frac{u(\nu,T)\,c}{4\pi}

である。よって

I(\nu,T) = \frac{2 h\nu^3 }{c^2}~\frac{1}{e^{h\nu/kT}-1}

を得る。

この式を波長についての I'(λ,T) の形式へと変換するためには、νc / λ で置き換え、

I'(\lambda,T) = I(\nu,T)\left|\frac{d\nu}{d\lambda}\right|

の式を計算する。

[編集] 百分率

プランクの法則のグラフの形状は温度に依存しない。よって波長または周波数に温度をかけた値を基準として、全放射輝度の百分位点を示すことができる。以下の表では、1行目は放射輝度の百分位点、2行目には対応する波長を温度で割った値 x = \lambda \cdot T [μm⋅K] を示した。例えば、20% の 2676 というのは、0 - 2676 [μm⋅K] が全放射輝度の 20% を占めるということを意味している。

10% 20% 25.0% 30% 40% 41.8% 50% 60% 64.6% 70% 80% 90%
2195 2676 2898 3119 3582 3670 4107 4745 5099 5590 6864 9376

波長と周波数のピークはそれぞれ 25.0% と 64.6% にあり、表中に太字で示した。41.8% の点は波長と周波数の中間ピークである。これらはそれぞれ、プランクの式のうち 1 / λ5, ν3, (ν / λ)2 が最大となる点である。

どのピークを用いるかは応用する場合による。便利な選択は、ヴィーンの変位則による25.0%の波長ピークである。いくつかの目的には、全放射輝度を1/2に分ける中央値(50%の点)がより適している。放射輝度は短波長では指数的に、長波長では多項式的に減少するため、後者は波長ピークよりも周波数ピークにより近い。同じ理由により、中間ピークは中央値よりも短い波長に位置する。

太陽はT = 5778 K の黒体放射体とする近似が非常によく成り立ち、10% - 90% の百分位点を以下のように表にすることができる。2行目はナノメートル単位の波長である。

10% 20% 25.0% 30% 40% 41.8% 50% 60% 64.6% 70% 80% 90%
380 463 502 540 620 635 711 821 882 967 1188 1623

これは大気の上部に到達する放射輝度である。400 nm以下(紫外線)の放射輝度はおよそ12%であり、一方700 nm以上(赤外線)は全体の51%である。大気はこの分布を大きく変化させる。具体的には大部分の紫外線とかなりの赤外線を吸収し、可視光線の比率を上昇させる。

[編集] 歴史

一部の物理学の教科書を含む、量子理論の説明の多くは、プランクの法則の説明において重大な間違いを犯している。この間違いは1960年代よりも前に物理学史の研究者によって指摘されたものの、現状が示しているようにこの間違いを根絶するのは難しい。Helge Kraghの論文により、実際には何が起きたのかについてのはっきりした説明が与えられた[8]

広く知られている通説に反し、プランクは光を量子化しなかった。この根拠としては、プランクの1901年のオリジナルの論文[6]と、その中に参考文献としてあげられている彼の1901年以前の論文があげられる。また、著書 "Theory of Heat Radiation" においてプランク定数はヘルツ振動子 (Hertzian oscillator) を示していると説明している。量子化の概念は別の第三者によって、現在では量子力学として知られるものの中に開発された。この流れの中で次に重要な段階を踏んだのはアルベルト・アインシュタインであった。アインシュタインは光電効果を研究し、光は塊や光子として放出されるだけでなく、吸収もされるという模型と方程式を提出した。そして1924年、サティエンドラ・ボースがプランクの法則を理論的に導出することができる光子の統計力学を考え出した。

また別の通説に反し、プランクは「紫外発散」の問題を解決しようとしてこの法則を導いたわけではなかった。紫外発散とはポール・エーレンフェストによって与えられた用語であり、黒体放射に古典統計力学のエネルギー等配分の法則を適用すると、空洞の全エネルギーが無限大になってしまうという矛盾である。プランクは等配分則が普遍的に成り立っているとは考えておらず、よって「発散」の問題にも気づいていなかった。紫外発散は5年ほど後の1905年にアインシュタイン、レイリー卿ジェームズ・ジーンズによって独立に発見された。

[編集] 附則

積分

J=\int_{0}^{\infty}\frac{x^3}{e^x-1}\,dx

を計算する簡単な方法は、まず一般化した場合について計算し、最後に比較して答えを出すことである。積分

\int_{0}^{\infty}\frac{x^n}{e^x-1}\,dx = \int_{0}^{\infty}\frac{x^n e^{-x}}{1 - e^{-x}}\,dx

を考える。右辺の分母は常に1よりも小さいため、等比数列の総和の公式を用いて収束する数列へと式変形する。

 \frac{1}{1-e^{-x}} = \sum_{k=0}^{\infty} e^{-kx}

つまり、左辺の分数は無限等比級数の総和 1 + e^{-x} + e^{-2x} + e^{-3x} + \cdots に等しい。初項は 1、公比は e x である。

よって

\int_{0}^{\infty}x^{n} e^{-x} \sum_{k=0}^{\infty} e^{-kx}\,dx

とすることができた。

総和の左側の e x を総和の中に入れることを考える。すると 1 + e^{-x} + e^{-2x} + \cdotse^{-x} + e^{-2x} + e^{-3x} + \cdots へと変化するため、e x を取り除く代わりに総和の下限を 0 から 1 にずらし

\int_{0}^{\infty}x^{n}  \sum_{k=1}^{\infty} e^{-kx}\,dx

となる。

変数を u = kx と変えると、xndx はそれぞれ x^n = \frac{u^n}{k^n}dx = \frac{du}{k} となる。よって

\int_{0}^{\infty}\frac{u^n}{k^n} \sum_{k=1}^{\infty} e^{-u}\frac{du}{k}

または

\int_{0}^{\infty}u^n \sum_{k=1}^{\infty}\frac{1}{k^{n + 1}} e^{-u}du

となる。

総和のそれぞれの項は収束する積分であるため、総和を積分の外へ出すことができ

\sum_{k=1}^{\infty} \frac{1}{k^{n+1}} \int_{0}^{\infty}u^{n} e^{-u}\,du

と式変形することができる。

左の総和はゼータ関数 ζ(n + 1) であり、右の積分はガンマ関数 Γ(n + 1) である。以上より一般的な結果

\int_{0}^{\infty}\frac{x^{n}}{e^x-1}\,dx = \zeta(n+1) \Gamma{\left(n+1\right)}

または

\int_{0}^{\infty}\frac{x^{n-1}}{e^x - 1}\,dx = \zeta{\left(n\right)} \Gamma{\left(n\right)}

を得ることができた。

もともと考えていた問題は x3 の場合であったが、この場合の答えは

J=\zeta{\left(4\right)} \Gamma{\left(4\right)} = \frac{\pi^{4}}{90} \times 6 = \frac{\pi^4}{15}

である。以上より、問題の積分を計算することができた。

ここで

\sum_{n=1}^{\infty}\frac{1}{n^{4}} = \frac{\pi^4}{90}

という等式を用いた。この式はゼータ関数の特殊値である。この値はウォリス積 (Wallis product#Finding Zeta(4)) や閉路積分留数定理から導くことができる。また、ガンマ関数の計算は公式 Γ(n + 1) = n! により求めることができる。

今回扱った積分の異なる計算法はen:Stefan-Boltzmann law#Appendixを参照。また関連した積分として、ボース=アインシュタイン統計における積分計算で現れる多重対数関数も参照。

[編集] 関連項目

[編集] 参考文献と脚注

  1. ^ (Rybicki & Lightman 1979, p. 22)
  2. ^ Kittel, Thermal Physics p98
  3. ^ (Rybicki & Lightman 1979, p. 22)
  4. ^ (Rybicki & Lightman 1979, p. 1)
  5. ^ Brehm, J.J. and Mullin, W.J., "Introduction to the Structure of Matter: A Course in Modern Physics," (Wiley, New York, 1989) ISBN 047160531X.
  6. ^ Planck, Max, "On the Law of Distribution of Energy in the Normal Spectrum". Annalen der Physik, vol. 4, p. 553 ff (1901)
  7. ^ Ribaric, M. and Sustersic, L., arxiv:0810.0905.
  8. ^ Kragh, Helge Max Planck: The reluctant revolutionary Physics World, December 2000.

[編集] 参考文献

  • Rybicki, G. B.; Lightman, A. P. (1979), Radiative Processes in Astrophysics, New York: John Wiley & Sons, ISBN 0-471-82759-2 
  • Thornton, Stephen T., Andrew Rex (2002). Modern Physics. USA: Thomson Learning. ISBN 0-03-006049-4. 

[編集] より詳しく

  • Peter C. Milonni (1994). The Quantum Vacuum. Academic Press. 

[編集] 外部リンク

最終更新 2009年11月5日 (木) 13:29 (日時は個人設定で未設定ならばUTC)。
【プランクの法則】変更履歴

ご利用上の注意