光学と量子エレクトロニクス 第13章 §13.5 ★★★ 自由電子だけなら貴金属はどれも銀色になる。満ちた d バンドからのバンド間遷移を ε 2 \varepsilon_2 ε 2 の段として足し、ε 1 \varepsilon_1 ε 1 を Kramers–Kronig で作ると、遷移の始まりの位置で金・銀・銅の色が分かれる
自由電子だけでは色が出ない
金・銀・銅は伝導電子が原子1個あたり1個で、自由電子の数密度から出す ℏ ω p \hbar\omega_{\mathrm p} ℏ ω p はどれも 9 ∼ 11 e V 9\sim11\,\mathrm{eV} 9 ∼ 11 eV の紫外にある。可視光(1.65 ∼ 3.26 e V 1.65\sim3.26\,\mathrm{eV} 1.65 ∼ 3.26 eV )は緩和の帯に入り、(13.13) の反射率は周波数によらないので、自由電子だけの模型では3つとも白く(銀色に)見えるはずである。金が黄色く銅が赤いのは、伝導電子以外の電子が可視光に応答するからである。
これらの金属では、満ちた d バンドが Fermi 準位の数 eV 下にあり、そこから Fermi 準位の上の空いた状態へ電子を上げるバンド間遷移が、ある光子のエネルギー ℏ ω t \hbar\omega_{\mathrm t} ℏ ω t から始まる。H. Ehrenreich と H. R. Philipp は銀と銅の光学定数の測定を、自由電子の応答・バンド間遷移・集団振動の3つで解釈した1 1 H. Ehrenreich and H. R. Philipp (1962), Physical Review 128, 1622. 。遷移の始まる位置が、銀では紫外に、金と銅では可視の中にある。結晶のバンドからこの遷移を数える話(結合状態密度と吸収端の形)は ハミルトニアンの選び方で固体を考える 第4部第6章 で扱う。この節では、その結果をもっとも簡単な形の ε 2 \varepsilon_2 ε 2 として模型に入れる。
バンド間遷移を段で入れる
バンド間遷移の ε 2 \varepsilon_2 ε 2 を、ω t < ω < ω c \omega_{\mathrm t}<\omega<\omega_{\mathrm c} ω t < ω < ω c で一定値 A A A 、それ以外で 0 0 0 の段(箱)で近似する。実部は (13.9) で
ε 1 i n t e r ( ω ) = 2 π P ∫ ω t ω c A ω ′ ω ′ 2 − ω 2 d ω ′ = A π ln ∣ ω c 2 − ω 2 ω t 2 − ω 2 ∣ (13.17) \varepsilon_1^{\mathrm{inter}}(\omega)=\frac{2}{\pi}\,\mathrm P\!\int_{\omega_{\mathrm t}}^{\omega_{\mathrm c}}\frac{A\omega'}{\omega'^2-\omega^2}\,\mathrm{d}\omega'
=\frac{A}{\pi}\ln\lvert \frac{\omega_{\mathrm c}^2-\omega^2}{\omega_{\mathrm t}^2-\omega^2}\rvert\tag{13.17} ε 1 inter ( ω ) = π 2 P ∫ ω t ω c ω ′2 − ω 2 A ω ′ d ω ′ = π A ln ∣ ω t 2 − ω 2 ω c 2 − ω 2 ∣ ( 13.17 )
になる。ω < ω t \omega<\omega_{\mathrm t} ω < ω t で正で、ω → ω t \omega\to\omega_{\mathrm t} ω → ω t に向かって大きくなる。吸収の段の手前では、その段を作る電子が束縛された振動子のように分極し、正の誘電率を足す。(13.17) は吸収の端で対数的に発散する。端をなめらかにするには、ω \omega ω を上半面の z = ω + i η z=\omega+\mathrm{i}\eta z = ω + i η に置き換えた
ε r ( ω ) = 1 − ω p 2 ω ( ω + i γ ) + A π [ ln ( ω c 2 − z 2 ) − ln ( ω t 2 − z 2 ) ] , z = ω + i η (13.18) \varepsilon_{\mathrm r}(\omega)=1-\frac{\omega_{\mathrm p}^2}{\omega(\omega+\mathrm{i}\gamma)}+\frac{A}{\pi}\Bigl[\ln(\omega_{\mathrm c}^2-z^2)-\ln(\omega_{\mathrm t}^2-z^2)\Bigr],\qquad z=\omega+\mathrm{i}\eta\tag{13.18} ε r ( ω ) = 1 − ω ( ω + i γ ) ω p 2 + π A [ ln ( ω c 2 − z 2 ) − ln ( ω t 2 − z 2 ) ] , z = ω + i η ( 13.18 )
を使う(対数は主値の枝)。角括弧の関数は上半面で正則で、∣ z ∣ → ∞ \lvert z\rvert \to\infty ∣ z ∣ → ∞ で 0 0 0 に近づくので、η > 0 \eta>0 η > 0 でも Kramers–Kronig の関係を満たす。η → + 0 \eta\to+0 η → + 0 では虚部が ω t < ω < ω c \omega_{\mathrm t}<\omega<\omega_{\mathrm c} ω t < ω < ω c で A A A 、外で 0 0 0 に戻り、実部が (13.17) に戻る。先頭の 1 1 1 は、より深い内殻の寄与を無視した値である。
段が総和則に足す面積は ∫ ω ε 2 i n t e r d ω = A ( ω c 2 − ω t 2 ) / 2 \int\omega\varepsilon_2^{\mathrm{inter}}\mathrm{d}\omega=A(\omega_{\mathrm c}^2-\omega_{\mathrm t}^2)/2 ∫ ω ε 2 inter d ω = A ( ω c 2 − ω t 2 ) /2 で、これを (13.10) の形 π 2 ω p , i n t e r 2 \frac{\pi}{2}\omega_{\mathrm{p,inter}}^2 2 π ω p , inter 2 に書けば、段の電子の数を伝導電子の ( ω p , i n t e r / ω p ) 2 (\omega_{\mathrm{p,inter}}/\omega_{\mathrm p})^2 ( ω p , inter / ω p ) 2 倍と数えたことになる。模型の値は表 13.1 のように選ぶ。3つの金属で ℏ ω p \hbar\omega_{\mathrm p} ℏ ω p を共通の 9.0 e V 9.0\,\mathrm{eV} 9.0 eV にそろえ、遷移の始まり ℏ ω t \hbar\omega_{\mathrm t} ℏ ω t と段の高さ A A A だけを変えた。これは色の分かれ方を見るための模型の値で、測定値を再現するように合わせたものではない。
表 13.1 (13.18) の模型の値(単位は eV。A A A は無次元) ℏ ω p \hbar\omega_{\mathrm p} ℏ ω p ℏ γ \hbar\gamma ℏ γ A A A ℏ ω t \hbar\omega_{\mathrm t} ℏ ω t ℏ ω c \hbar\omega_{\mathrm c} ℏ ω c ℏ η \hbar\eta ℏ η 銀の模型 9.0 9.0 9.0 0.02 0.02 0.02 2.5 2.5 2.5 3.9 3.9 3.9 20 20 20 0.1 0.1 0.1 金の模型 9.0 9.0 9.0 0.07 0.07 0.07 4.0 4.0 4.0 2.4 2.4 2.4 20 20 20 0.1 0.1 0.1 銅の模型 9.0 9.0 9.0 0.09 0.09 0.09 5.0 5.0 5.0 2.1 2.1 2.1 20 20 20 0.1 0.1 0.1
図 13.4 表 13.1 の模型 (13.18) で計算した垂直入射の反射率。帯は可視光(1.65 ∼ 3.26 e V 1.65\sim3.26\,\mathrm{eV} 1.65 ∼ 3.26 eV )。銀の模型は可視で R ≈ 0.99 R\approx0.99 R ≈ 0.99 が平らに続き、3.85 e V 3.85\,\mathrm{eV} 3.85 eV に鋭い谷がある。金と銅の模型は遷移の始まり(2.4 2.4 2.4 と 2.1 e V 2.1\,\mathrm{eV} 2.1 eV )から上で反射率が落ち、青の側を反射しにくい
遮蔽されたプラズマ周波数と銀の谷
ω < ω t \omega<\omega_{\mathrm t} ω < ω t では、段の寄与 (13.17) は吸収を持たない背景の誘電率 ε b ( ω ) = 1 + ε 1 i n t e r ( ω ) \varepsilon_{\mathrm b}(\omega)=1+\varepsilon_1^{\mathrm{inter}}(\omega) ε b ( ω ) = 1 + ε 1 inter ( ω ) として振る舞う。ε 1 = 0 \varepsilon_1=0 ε 1 = 0 になる周波数は
ω ∗ = ω p ε b ( ω ∗ ) \omega^*=\frac{\omega_{\mathrm p}}{\sqrt{\varepsilon_{\mathrm b}(\omega^*)}} ω ∗ = ε b ( ω ∗ ) ω p
で、自由電子の ω p \omega_{\mathrm p} ω p より下がる。d 電子の分極が伝導電子の振動の電場を弱めるので、これを遮蔽されたプラズマ周波数と呼ぶ。銀の模型では ω ∗ \omega^* ω ∗ が遷移の始まり ω t \omega_{\mathrm t} ω t のすぐ下に来るので、ε 1 = 0 \varepsilon_1=0 ε 1 = 0 の点で ε 2 \varepsilon_2 ε 2 がまだ小さい。ε r ≈ 0 \varepsilon_{\mathrm r}\approx0 ε r ≈ 0 なら n ~ ≈ 0 \tilde n\approx0 n ~ ≈ 0 で、(13.8) は R ≈ 1 R\approx1 R ≈ 1 でなく 0 0 0 に近づく。反射率は ω ∗ \omega^* ω ∗ の上で鋭く落ち、谷を作る。ω ∗ \omega^* ω ∗ が吸収の中に入る金と銅では、ε 2 \varepsilon_2 ε 2 が大きいので谷は浅くぼやける。Ehrenreich と Philipp が「伝導電子と d 電子の両方が関わるプラズマ共鳴」と呼んだものが、この ω ∗ \omega^* ω ∗ の振動である。
例題 13.7 銀の模型の遮蔽されたプラズマ周波数 ★★
オリジナル 表 13.1 の銀の模型について次を求めよ。
(1) ε 1 ( ω ) = 0 \varepsilon_1(\omega)=0 ε 1 ( ω ) = 0 になる ℏ ω ∗ \hbar\omega^* ℏ ω ∗ (3 ∼ 3.9 e V 3\sim3.9\,\mathrm{eV} 3 ∼ 3.9 eV の範囲)と、そこでの背景の誘電率 Re ε b ( ω ∗ ) \operatorname{Re}\,\varepsilon_{\mathrm b}(\omega^*) Re ε b ( ω ∗ ) 。
(2) 静電場の極限の背景 ε b ( 0 ) = 1 + ( A / π ) ln ( ω c 2 / ω t 2 ) \varepsilon_{\mathrm b}(0)=1+(A/\pi)\ln(\omega_{\mathrm c}^2/\omega_{\mathrm t}^2) ε b ( 0 ) = 1 + ( A / π ) ln ( ω c 2 / ω t 2 ) を使った見積もり ω p / ε b ( 0 ) \omega_{\mathrm p}/\sqrt{\varepsilon_{\mathrm b}(0)} ω p / ε b ( 0 ) との違い。
(3) 段が総和則に足す面積を ℏ 2 ω p , i n t e r 2 = A ( ℏ 2 ω c 2 − ℏ 2 ω t 2 ) / π \hbar^2\omega_{\mathrm{p,inter}}^2=A(\hbar^2\omega_{\mathrm c}^2-\hbar^2\omega_{\mathrm t}^2)/\pi ℏ 2 ω p , inter 2 = A ( ℏ 2 ω c 2 − ℏ 2 ω t 2 ) / π で表し、伝導電子の ℏ 2 ω p 2 \hbar^2\omega_{\mathrm p}^2 ℏ 2 ω p 2 と比べよ。
解答 (1) (13.18) の実部を数値で 0 0 0 にすると ℏ ω ∗ = 3.79 e V \hbar\omega^*=3.79\,\mathrm{eV} ℏ ω ∗ = 3.79 eV で、そこでの Re ε b = 5.64 \operatorname{Re}\,\varepsilon_{\mathrm b}=5.64 Re ε b = 5.64 、9.0 / 5.64 = 3.79 9.0/\sqrt{5.64}=3.79 9.0/ 5.64 = 3.79 と矛盾なくそろう。ε 2 ( ω ∗ ) = 0.61 \varepsilon_2(\omega^*)=0.61 ε 2 ( ω ∗ ) = 0.61 と小さいので反射率の谷が深い(図 13.4 の最小は 3.85 e V 3.85\,\mathrm{eV} 3.85 eV で R = 0.09 R=0.09 R = 0.09 )。
(2) ε b ( 0 ) = 1 + ( 2.5 / π ) ln ( 400 / 15.21 ) = 3.60 \varepsilon_{\mathrm b}(0)=1+(2.5/\pi)\ln(400/15.21)=3.60 ε b ( 0 ) = 1 + ( 2.5/ π ) ln ( 400/15.21 ) = 3.60 で、9.0 / 3.60 = 4.74 e V 9.0/\sqrt{3.60}=4.74\,\mathrm{eV} 9.0/ 3.60 = 4.74 eV になる。段の手前で ε b ( ω ) \varepsilon_{\mathrm b}(\omega) ε b ( ω ) が大きくなるので、静電場の値で見積もると ω ∗ \omega^* ω ∗ を 1 eV ほど高く出してしまう。
(3) ℏ 2 ω p , i n t e r 2 = 2.5 × ( 400 − 15.21 ) / π = 306 e V 2 \hbar^2\omega_{\mathrm{p,inter}}^2=2.5\times(400-15.21)/\pi=306\,\mathrm{eV^2} ℏ 2 ω p , inter 2 = 2.5 × ( 400 − 15.21 ) / π = 306 e V 2 で、伝導電子の 81 e V 2 81\,\mathrm{eV^2} 81 e V 2 の 3.8 3.8 3.8 倍である。20 e V 20\,\mathrm{eV} 20 eV までの吸収に、伝導電子1個に対して d 電子 3.8 3.8 3.8 個分の重みを割り当てた模型になっている。
例題 13.8 金と銅の模型の反射の色 ★★
オリジナル 表 13.1 の3つの模型で、赤(1.9 e V 1.9\,\mathrm{eV} 1.9 eV , 653 n m 653\,\mathrm{nm} 653 nm )・緑(2.3 e V 2.3\,\mathrm{eV} 2.3 eV , 539 n m 539\,\mathrm{nm} 539 nm )・青(2.75 e V 2.75\,\mathrm{eV} 2.75 eV , 451 n m 451\,\mathrm{nm} 451 nm )の垂直入射の反射率を求め、白色光を反射したときの色を述べよ。
解答 (13.18) と (13.8) から次のようになる。
赤 1.9 e V 1.9\,\mathrm{eV} 1.9 eV 緑 2.3 e V 2.3\,\mathrm{eV} 2.3 eV 青 2.75 e V 2.75\,\mathrm{eV} 2.75 eV 銀の模型 0.994 0.994 0.994 0.992 0.992 0.992 0.986 0.986 0.986 金の模型 0.966 0.966 0.966 0.855 0.855 0.855 0.495 0.495 0.495 銅の模型 0.923 0.923 0.923 0.589 0.589 0.589 0.424 0.424 0.424
銀の模型は3色をほぼ同じに反射するので、反射光は白い。金の模型は青を半分しか反射せず、赤と緑を多く返すので、反射光は黄色に寄る。銅の模型は緑も6割に落ち、赤だけが強く残るので赤みが強い。3つの違いは ℏ ω t \hbar\omega_{\mathrm t} ℏ ω t を 3.9 3.9 3.9 , 2.4 2.4 2.4 , 2.1 e V 2.1\,\mathrm{eV} 2.1 eV と変えたことだけから出ている。遷移の始まりが可視の中に入ると、それより短い波長の光は電子をバンド間で上げることに使われ、反射されにくくなる。
1 H. Ehrenreich and H. R. Philipp (1962), Physical Review 128, 1622.
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com