tecktom.com

Einstein の A,BA,B 係数:熱力学が自然放出を強制する

2026年10月15日4分で読めますAI生成教科書『基礎理論』第12部 第1章 輻射遷移:原子はなぜ・どのように光るか

AI 生成・木下が未勉強

量子力学入門 第12部 第1章 §12.1.2 ★★ 自然放出・AA/BB 係数(光学遷移の下支え)

自然放出は「光子場の零点振動が引き起こす誘導放出」であり、本当は電磁場の量子化(QED)が要る。しかし Einstein(1917)は、場を量子化せずに熱力学だけで自然放出の存在とそのレートを導いてみせた。論法を再現する。

2準位(下 11、上 22、縮退なし、E2−E1=ℏωE_2-E_1=\hbar\omega)の原子集団が、温度 TT の熱輻射(スペクトルエネルギー密度 u(ω)u(\omega))と平衡にあるとする。レートを

吸収:B12 u(ω)N1,誘導放出:B21 u(ω)N2,自然放出:A N2\text{吸収:}B_{12}\,u(\omega)N_1,\qquad \text{誘導放出:}B_{21}\,u(\omega)N_2,\qquad \text{自然放出:}A\,N_2

と置く(AA の項が無いと仮定してみると、平衡条件 B12uN1=B21uN2B_{12}uN_1=B_{21}uN_2 から N2/N1N_2/N_1 が温度によらなくなり、Boltzmann 分布 N2/N1=e−ℏω/kBTN_2/N_1=\mathrm{e}^{-\hbar\omega/k_BT} と矛盾する——自然放出は熱力学的に存在しなければならない)。AA 込みの平衡条件を uu について解くと

u(ω)=A/B21(B12/B21) eℏω/kBT−1.u(\omega)=\frac{A/B_{21}}{(B_{12}/B_{21})\,\mathrm{e}^{\hbar\omega/k_BT}-1}.

これが任意の温度で Planck 分布 u=ℏω3π2c31eℏω/kBT−1u=\frac{\hbar\omega^3}{\pi^2c^3}\frac{1}{\mathrm{e}^{\hbar\omega/k_BT}-1} に一致する条件は

 B12=B21,A=ℏω3π2c3 B21 \boxed{\ B_{12}=B_{21},\qquad A=\frac{\hbar\omega^3}{\pi^2c^3}\,B_{21}\ }

——吸収と誘導放出の係数は等しく、自然放出は誘導放出から普遍係数で決まる。B21B_{21} は通常の摂動論(第8部『時間依存摂動論とダイナミクス』)で計算できるから、これで AA が手に入る。結果(導出は BB の計算と等方平均 ∣dz∣2‾=∣d∣2/3\overline{|d_z|^2}=|\bm{d}|^2/3 を組むだけ):

 A=ω3 ∣dfi∣23πε0ℏc3 (SI 単位系)(12.1.1)\boxed{\ A=\frac{\omega^3\,|\bm{d}_{fi}|^2}{3\pi\varepsilon_0\hbar c^3}\ } \qquad(\text{SI 単位系})\tag{12.1.1}

励起状態の寿命は τ=1/∑fAf\tau=1/\sum_fA_{f}(全ての下向き遷移の和)。A∝ω3A\propto\omega^3 に注意——同じ双極子でも高い振動数の遷移ほど圧倒的に速く光る(ω2\omega^2 は終状態の光子の状態密度から、もう1つの ω\omega は光子1個のエネルギーからくる)。

例題 1.1(水素 2p2p 状態の寿命を計算する ★★★)

オリジナル水素の 2p→1s2p\to1s 遷移について双極子行列要素を計算し、寿命 τ=1/A\tau=1/A を求めよ。

解答

m=0m=0 の成分で代表させる(他の mm も等方性から同じ AA を与える)。⟨1s∣z^∣2p0⟩\langle 1s | \hat z | 2p_0 \rangle は動径積分と角度積分の積で、動径部分は R10=2a0−3/2e−r/a0R_{10}=2a_0^{-3/2}\mathrm{e}^{-r/a_0}、R21=126a0−3/2ra0e−r/2a0R_{21}=\frac{1}{2\sqrt6}a_0^{-3/2}\frac{r}{a_0}\mathrm{e}^{-r/2a_0} を使って

∫0∞R10R21 r3dr=1286243 a0≃1.290 a0,\int_0^\infty R_{10}R_{21}\,r^3\mathrm{d} r=\frac{128\sqrt6}{243}\,a_0 \simeq1.290\,a_0,

角度部分は ∫Y00∗cos⁡θ Y10dΩ=1/3\int Y_{00}^*\cos\theta\,Y_{10}\mathrm{d}\Omega=1/\sqrt3。よって ∣⟨1s∣z^∣2p0⟩∣=1282243a0≃0.745 a0|\langle 1s | \hat z | 2p_0 \rangle |=\frac{128\sqrt2}{243}a_0\simeq0.745\,a_0(動径積分の閉形式と数値積分の一致は確認済み)。ℏω=13.6 (1−14)=10.2 eV\hbar\omega=13.6\,(1-\tfrac14)=10.2\,\mathrm{eV} とともに (12.1.1) に代入すると

A≃6.3×108 s−1, τ≃1.60 ns A\simeq6.3\times10^8\ \mathrm{s^{-1}},\qquad \boxed{\ \tau\simeq1.60\ \mathrm{ns}\ }

——実測の水素 Lyman-α\alpha 寿命 1.6 ns1.6\,\mathrm{ns} と一致する(数値は SI 定数で独立計算して検証済み)。オーダーの覚え方:τ−1∼α3ω×(無次元数)\tau^{-1}\sim\alpha^3\omega\times(\text{無次元数})——電磁結合の微細構造定数 α≃1/137\alpha\simeq1/137 の3乗が「原子は振動周期(∼10−16\sim10^{-16} s)に比べてずっとゆっくり(10610^6 周期以上かけて)光る」ことを説明する。準位の自然幅は Γ=ℏ/τ∼4×10−7 eV\Gamma=\hbar/\tau\sim4\times10^{-7}\,\mathrm{eV}——準位間隔の 10−810^{-8} 倍で、「離散スペクトル」という描像が実際上成り立つ理由でもある。■\blacksquare

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『基礎理論』の記事一覧へ