光学と量子エレクトロニクス 第11章 §11.1 ★★★ 平面波1つと、振動子強度で重みをつけた Lorentz 振動子の集まり。薄い気体では原子1個の分極率を足すだけで感受率が決まる
第10章(第10章 )では、物質中の Maxwell 方程式を分極 P \bm{P} P で書き、線形で等方的な媒質の性質が比誘電率 ε r ( ω ) = 1 + χ ( ω ) \varepsilon_{\mathrm r}(\omega)=1+\chi(\omega) ε r ( ω ) = 1 + χ ( ω ) と複素屈折率 n ~ = ε r = n + i κ \tilde n=\sqrt{\varepsilon_{\mathrm r}}=n+\mathrm{i}\kappa n ~ = ε r = n + i κ にまとまることを示した。そこでは χ ( ω ) \chi(\omega) χ ( ω ) を与えられた関数として扱った。この章では、χ ( ω ) \chi(\omega) χ ( ω ) そのものを物質のモデルから計算する。
この章で取り出すモードは、無限に広い一様な媒質の中を進む単色の平面波である。物質のモデルは、単位体積あたり N N N 個の原子のそれぞれを Lorentz 振動子(第1章の式で、第1章 )とみなすもので、1つの原子が複数の共鳴を持つときは振動子強度で重みをつけて足す。結合は電子が感じる電場との相互作用 − e E -e\bm{E} − e E だが、ここで1つ問題が起こる。原子が密に並んだ固体や液体では、1つの原子が感じる電場は、巨視的な平均の電場 E \bm{E} E に、周りの原子の双極子が作る場を加えたものになる。これを局所場 と呼ぶ。
11.1節で、局所場を無視できる薄い気体について χ ( ω ) \chi(\omega) χ ( ω ) を書き下す。11.2 節で Lorentz の球を使って局所場を求め、Clausius–Mossotti の関係を導く。11.3 節と11.4 節が「出てくるもの」で、共鳴の近くの正常分散と異常分散、局所場による共鳴のずれと全反射の帯(残留線の帯)、位相速度と群速度、2つの振動子を結合させた遅い光を順に導く。11.5 節で置いた近似と実験で確かめる量を述べ、11.6 節で学部の光学で習う Cauchy の式・Sellmeier の式と並べ、11.7 節でまとめる。
モード・物質・結合の3つの選び方
場のモードとして、角周波数 ω \omega ω の平面波 E ( z , t ) = Re [ E 0 e i ( n ~ ω z / c − ω t ) ] \bm{E}(z,t)=\operatorname{Re}[\bm{E}_0\mathrm{e}^{\mathrm{i}(\tilde n\omega z/c-\omega t)}] E ( z , t ) = Re [ E 0 e i ( n ~ ω z / c − ω t ) ] を1つだけ取り出す。波長は原子の間隔よりはるかに長いとし、原子の間隔より大きく波長より小さい範囲で平均した巨視的な場を E \bm{E} E と書く。第10章で示したとおり、この平面波の強度は e − α a b s z \mathrm{e}^{-\alpha_{\mathrm{abs}}z} e − α abs z で減り、吸収係数は α a b s = 2 ω κ / c \alpha_{\mathrm{abs}}=2\omega\kappa/c α abs = 2 ω κ / c である。
物質のモデルは、第1章の Lorentz 振動子
m ( d 2 x d t 2 + γ d x d t + ω 0 2 x ) = − e E l o c ( t ) (11.1) m\Bigl(\frac{\mathrm{d}^2\bm{x}}{\mathrm{d} t^2}+\gamma\frac{\mathrm{d}\bm{x}}{\mathrm{d} t}+\omega_0^2\bm{x}\Bigr)=-e\bm{E}_{\mathrm{loc}}(t)\tag{11.1} m ( d t 2 d 2 x + γ d t d x + ω 0 2 x ) = − e E loc ( t ) ( 11.1 )
である。右辺の E l o c \bm{E}_{\mathrm{loc}} E loc は電子の位置にある電場で、局所場と呼ぶ。原子1個の双極子 d = − e x \bm{d}=-e\bm{x} d = − e x の複素振幅は α ( ω ) E l o c \alpha(\omega)\bm{E}_{\mathrm{loc}} α ( ω ) E loc になり、分極率は第1章の式と同じ α ( ω ) = ( e 2 / m ) / ( ω 0 2 − ω 2 − i γ ω ) \alpha(\omega)=(e^2/m)/(\omega_0^2-\omega^2-\mathrm{i}\gamma\omega) α ( ω ) = ( e 2 / m ) / ( ω 0 2 − ω 2 − i γ ω ) である。
原子の遷移は1つではない。第1章で見たとおり、量子力学の2準位系は弱い光に対して「電子 f f f 個分の Lorentz 振動子」として応答する(f f f は振動子強度)。原子の基底状態から出る遷移 j j j の角周波数を ω j \omega_j ω j 、振動子強度を f j f_j f j 、線幅を γ j \gamma_j γ j とすると、原子1個の分極率は遷移ごとの寄与の和
α ( ω ) = e 2 m ∑ j f j ω j 2 − ω 2 − i γ j ω (11.2) \alpha(\omega)=\frac{e^2}{m}\sum_j\frac{f_j}{\omega_j^2-\omega^2-\mathrm{i}\gamma_j\omega}\tag{11.2} α ( ω ) = m e 2 j ∑ ω j 2 − ω 2 − i γ j ω f j ( 11.2 )
で書ける。減衰を落とした形は、2次の摂動論で計算した量子力学の分極率と一致する(第1章の例題で、遷移が1つの場合を確かめた)。γ j \gamma_j γ j は現象論的に入れた減衰で、その量子論の起源は第24章(第24章 )と第34章(第34章 )で扱う。
薄い気体の感受率
原子の間隔が大きい気体では、周りの原子の双極子が作る場は巨視的な場に比べて小さく、E l o c ≈ E \bm{E}_{\mathrm{loc}}\approx\bm{E} E loc ≈ E とおける(どれだけ小さいかは 11.2 節で確かめる)。このとき分極は P = N α ( ω ) E \bm{P}=N\alpha(\omega)\bm{E} P = N α ( ω ) E で、P = ε 0 χ E \bm{P}=\varepsilon_0\chi\bm{E} P = ε 0 χ E と比べて
χ ( ω ) = N α ( ω ) ε 0 = ω p 2 ∑ j f j ω j 2 − ω 2 − i γ j ω , ω p 2 = N e 2 ε 0 m (11.3) \chi(\omega)=\frac{N\alpha(\omega)}{\varepsilon_0}=\omega_{\mathrm p}^2\sum_j\frac{f_j}{\omega_j^2-\omega^2-\mathrm{i}\gamma_j\omega},\qquad\omega_{\mathrm p}^2=\frac{Ne^2}{\varepsilon_0m}\tag{11.3} χ ( ω ) = ε 0 N α ( ω ) = ω p 2 j ∑ ω j 2 − ω 2 − i γ j ω f j , ω p 2 = ε 0 m N e 2 ( 11.3 )
を得る。ω p \omega_{\mathrm p} ω p は単位体積あたり N N N 個の自由電子の気体のプラズマ角周波数と同じ式で、第13章(第13章 )で金属の光学応答を決める量として再び現れる。∣ χ ∣ ≪ 1 \lvert \chi\rvert \ll1 ∣ χ ∣ ≪ 1 なら n ~ = 1 + χ ≈ 1 + χ / 2 \tilde n=\sqrt{1+\chi}\approx1+\chi/2 n ~ = 1 + χ ≈ 1 + χ /2 なので、
n − 1 ≈ 1 2 Re χ , κ ≈ 1 2 Im χ (11.4) n-1\approx\frac12\operatorname{Re}\chi,\qquad\kappa\approx\frac12\operatorname{Im}\chi\tag{11.4} n − 1 ≈ 2 1 Re χ , κ ≈ 2 1 Im χ ( 11.4 )
である。屈折率の 1 1 1 からのずれと吸収は、それぞれ感受率の実部と虚部の半分になる。時間の因子を e − i ω t \mathrm{e}^{-\mathrm{i}\omega t} e − i ω t に取ったので、Im χ > 0 \operatorname{Im}\chi>0 Im χ > 0 が吸収、Im χ < 0 \operatorname{Im}\chi<0 Im χ < 0 が増幅に対応する。
式 (11.3) は、光の周波数がどの遷移からも十分離れていれば γ j \gamma_j γ j を落として実数になり、媒質は透明になる。可視光に対して透明なガラスや空気は、共鳴が紫外と赤外にあって可視域がその間の谷にある物質である。その谷で屈折率がどう変わるかは 11.3 節と11.6 節で見る。
例題 11.1 原子状水素の気体の屈折率の見積もり ★★
オリジナル 0 ∘ ^\circ ∘ C・1 気圧の気体の数密度(Loschmidt 数)N = 2.687 × 10 25 m − 3 N=2.687\times10^{25}\,\mathrm{m^{-3}} N = 2.687 × 1 0 25 m − 3 の原子状水素を考える。第1章の例題で求めた 1 s → 2 p 1s\to2p 1 s → 2 p 遷移(ℏ ω 0 = 3 8 E h \hbar\omega_0=\frac38E_{\mathrm h} ℏ ω 0 = 8 3 E h 、E h = 4.360 × 10 − 18 J E_{\mathrm h}=4.360\times10^{-18}\,\mathrm J E h = 4.360 × 1 0 − 18 J 、振動子強度 f = 0.416 f=0.416 f = 0.416 )だけを残す。
(1) ω p 2 \omega_{\mathrm p}^2 ω p 2 と ω 0 \omega_0 ω 0 を求めよ。
(2) 波長 500 n m 500\,\mathrm{nm} 500 nm の光に対するこの遷移の寄与 n − 1 n-1 n − 1 を、減衰を落として求めよ。
(3) 残りの振動子強度 1 − f 1-f 1 − f はすべて ω 0 \omega_0 ω 0 より高い周波数の遷移にある。それを足すと n − 1 n-1 n − 1 は (2) より大きくなるか小さくなるか。
解答 (1) ω p 2 = N e 2 / ( ε 0 m ) = 2.687 × 10 25 × ( 1.602 × 10 − 19 ) 2 / ( 8.854 × 10 − 12 × 9.109 × 10 − 31 ) = 8.55 × 10 28 s − 2 \omega_{\mathrm p}^2=Ne^2/(\varepsilon_0m)=2.687\times10^{25}\times(1.602\times10^{-19})^2/(8.854\times10^{-12}\times9.109\times10^{-31})=8.55\times10^{28}\,\mathrm{s^{-2}} ω p 2 = N e 2 / ( ε 0 m ) = 2.687 × 1 0 25 × ( 1.602 × 1 0 − 19 ) 2 / ( 8.854 × 1 0 − 12 × 9.109 × 1 0 − 31 ) = 8.55 × 1 0 28 s − 2 。ω 0 = 3 8 E h / ℏ = 1.550 × 10 16 s − 1 \omega_0=\frac38E_{\mathrm h}/\hbar=1.550\times10^{16}\,\mathrm{s^{-1}} ω 0 = 8 3 E h /ℏ = 1.550 × 1 0 16 s − 1 (波長 121.5 n m 121.5\,\mathrm{nm} 121.5 nm の Lyman α \alpha α 線)。
(2) 光の角周波数は ω = 2 π c / λ = 3.767 × 10 15 s − 1 \omega=2\pi c/\lambda=3.767\times10^{15}\,\mathrm{s^{-1}} ω = 2 π c / λ = 3.767 × 1 0 15 s − 1 で、ω 0 2 − ω 2 = 2.262 × 10 32 s − 2 \omega_0^2-\omega^2=2.262\times10^{32}\,\mathrm{s^{-2}} ω 0 2 − ω 2 = 2.262 × 1 0 32 s − 2 。式 (11.3) と (11.4) から
n − 1 ≈ f ω p 2 2 ( ω 0 2 − ω 2 ) = 0.416 × 8.55 × 10 28 2 × 2.262 × 10 32 = 7.86 × 10 − 5 . n-1\approx\frac{f\omega_{\mathrm p}^2}{2(\omega_0^2-\omega^2)}=\frac{0.416\times8.55\times10^{28}}{2\times2.262\times10^{32}}=7.86\times10^{-5}. n − 1 ≈ 2 ( ω 0 2 − ω 2 ) f ω p 2 = 2 × 2.262 × 1 0 32 0.416 × 8.55 × 1 0 28 = 7.86 × 1 0 − 5 .
可視光は Lyman α \alpha α 線よりずっと低い周波数にあるので、分母の ω 2 \omega^2 ω 2 は ω 0 2 \omega_0^2 ω 0 2 の 6 % 6\% 6% にすぎず、n − 1 n-1 n − 1 は静電場に対する値とほとんど変わらない。
(3) 高い周波数の遷移 ω j > ω 0 \omega_j>\omega_0 ω j > ω 0 はどれも ω j 2 − ω 2 > 0 \omega_j^2-\omega^2>0 ω j 2 − ω 2 > 0 なので、正の寄与を足す。したがって n − 1 n-1 n − 1 は (2) より大きくなる。ただし分母が大きいので、同じ振動子強度でも寄与は小さい。たとえば残りの 0.584 0.584 0.584 がすべて ω 0 \omega_0 ω 0 の 1.5 1.5 1.5 倍の所にあれば、寄与は 0.584 / ( 2.25 − 0.059 ) 0.584/(2.25-0.059) 0.584/ ( 2.25 − 0.059 ) を 0.416 / ( 1 − 0.059 ) 0.416/(1-0.059) 0.416/ ( 1 − 0.059 ) と比べた比、つまり (2) の約 0.60 0.60 0.60 倍になる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com