ハミルトニアンの選び方で固体を考える 第3部 第3章 §3.3.3 ★★★ 光は波数が小さいので Brillouin 域の中心 k ≈ 0 \bm{k}\approx\bm{0} k ≈ 0 だけを見る。音響の枝の傾き(音速)を Brillouin 散乱、光学の枝の Γ \Gamma Γ 点を Raman 散乱で測る
光の散乱を分極率で書く
光の電場 E \bm{E} E は電子を揺らし、結晶に双極子モーメント p = α E \bm{p}=\alpha\bm{E} p = α E を誘起する(α \alpha α は分極率で、一般には2階のテンソル)。原子がずれると電子の雲の形が変わるので、分極率は基準振動の座標 Q Q Q の関数で、α ( Q ) = α 0 + ( ∂ α / ∂ Q ) 0 Q + ⋯ \alpha(Q)=\alpha_0+(\partial\alpha/\partial Q)_0\,Q+\cdots α ( Q ) = α 0 + ( ∂ α / ∂ Q ) 0 Q + ⋯ と展開することができる。電場を E 0 cos ω i t E_0\cos\omega_{\mathrm i}t E 0 cos ω i t 、振動を Q 0 cos Ω t Q_0\cos\Omega t Q 0 cos Ω t とすると
p ( t ) = α 0 E 0 cos ω i t + 1 2 ( ∂ α ∂ Q ) 0 Q 0 E 0 [ cos ( ω i − Ω ) t + cos ( ω i + Ω ) t ] (3.3.10) p(t)=\alpha_0E_0\cos\omega_{\mathrm i}t+\frac12\Bigl(\frac{\partial\alpha}{\partial Q}\Bigr)_{\!0}Q_0E_0\bigl[\cos(\omega_{\mathrm i}-\Omega)t+\cos(\omega_{\mathrm i}+\Omega)t\bigr]\tag{3.3.10} p ( t ) = α 0 E 0 cos ω i t + 2 1 ( ∂ Q ∂ α ) 0 Q 0 E 0 [ cos ( ω i − Ω ) t + cos ( ω i + Ω ) t ] ( 3.3.10 )
となる。第1項は周波数の変わらない Rayleigh 散乱、第2項は振動の周波数だけずれた散乱で、低い側 ω i − Ω \omega_{\mathrm i}-\Omega ω i − Ω を Stokes 線 、高い側を anti-Stokes 線 と呼ぶ。C. V. Raman と K. S. Krishnan は液体の散乱光にこのずれた線を見いだした(C. V. Raman and K. S. Krishnan (1928), Nature 121 , 501)。
量子力学では Q Q Q を Q ^ = ℏ / ( 2 Ω ) ( a ^ + a ^ † ) \hat Q=\sqrt{\hbar/(2\Omega)}(\hat a+\hat a^\dagger) Q ^ = ℏ/ ( 2Ω ) ( a ^ + a ^ † ) (質量を含めた基準座標)に置き換える。Stokes 線はフォノンを1個作る過程で行列要素の2乗が n B + 1 n_{\mathrm B}+1 n B + 1 に、anti-Stokes 線は1個吸収する過程で n B n_{\mathrm B} n B に比例する。中性子の式 (3.3.9) と同じ温度の因子が現れ、光でも k \bm{k} k 一定・反転対称の結晶なら、Stokes と anti-Stokes の強度比から詳細つり合いによって温度が分かる。
運動量の保存と2種類の光散乱
結晶の中で光の波数は n ω / c n\omega/c nω / c (n n n は屈折率)である。入射光 k i \bm{k}_{\mathrm i} k i が散乱光 k s \bm{k}_{\mathrm s} k s になり、フォノン k p h \bm{k}_{\mathrm{ph}} k ph を作るとき、k i = k s + k p h \bm{k}_{\mathrm i}=\bm{k}_{\mathrm s}+\bm{k}_{\mathrm{ph}} k i = k s + k ph である(光の波数は G \bm{G} G よりはるかに小さいので、G = 0 \bm{G}=\bm{0} G = 0 の場合しか起きない)。フォノンのエネルギーは光のエネルギーの数 % 以下なので ∣ k s ∣ ≈ ∣ k i ∣ \lvert \bm{k}_{\mathrm s}\rvert \approx\lvert \bm{k}_{\mathrm i}\rvert ∣ k s ∣ ≈ ∣ k i ∣ とおけ、散乱角を θ \theta θ とすると
∣ k p h ∣ = 2 n ω i c sin θ 2 (3.3.11) \lvert \bm{k}_{\mathrm{ph}}\rvert =2\,\frac{n\omega_{\mathrm i}}{c}\,\sin\frac\theta2\tag{3.3.11} ∣ k ph ∣ = 2 c n ω i sin 2 θ ( 3.3.11 )
である。例題 3.3.1 で見たように、これは Brillouin 域の大きさの 10 − 3 10^{-3} 1 0 − 3 の程度しかない。光が見るのは k p h ≈ 0 \bm{k}_{\mathrm{ph}}\approx\bm{0} k ph ≈ 0 のフォノンだけである。
k ≈ 0 \bm{k}\approx\bm{0} k ≈ 0 には2種類のフォノンがある。音響の枝は ω = v ∣ k ∣ \omega=v\lvert \bm{k}\rvert ω = v ∣ k ∣ (v v v は音速)で 0 に近づくので、ずれは
Ω B = v ∣ k p h ∣ = 2 n v ω i c sin θ 2 (3.3.12) \Omega_{\mathrm B}=v\lvert \bm{k}_{\mathrm{ph}}\rvert =\frac{2nv\,\omega_{\mathrm i}}{c}\,\sin\frac\theta2\tag{3.3.12} Ω B = v ∣ k ph ∣ = c 2 n v ω i sin 2 θ ( 3.3.12 )
と GHz の程度で、散乱角とともに変わる。これを Brillouin 散乱 と呼ぶ(L. Brillouin (1922), Ann. Phys. (Paris) 17 , 88)。θ \theta θ を変えて Ω B \Omega_{\mathrm B} Ω B を測ると、音響の枝の原点の傾き=音速が分かる。光学の枝は k → 0 \bm{k}\to\bm{0} k → 0 でも有限の ω T O \omega_{\mathrm{TO}} ω TO や ω L O \omega_{\mathrm{LO}} ω LO に近づくので、ずれは THz の程度で散乱角によらない。これが結晶の Raman 散乱 で、これで光学の枝の Γ \Gamma Γ 点(Brillouin 域の中心)の値が分かる。分光の慣習では、ずれを波長の逆数の差 ν ˉ ≔ 1 / λ i − 1 / λ s \bar\nu\coloneqq1/\lambda_{\mathrm i}-1/\lambda_{\mathrm s} ν ˉ : = 1/ λ i − 1/ λ s (単位 c m − 1 \mathrm{cm^{-1}} c m − 1 )で書き、Raman シフトと呼ぶ。ℏ Ω = h c ν ˉ \hbar\Omega=hc\bar\nu ℏΩ = h c ν ˉ で、1 m e V 1\,\mathrm{meV} 1 meV は 8.07 c m − 1 8.07\,\mathrm{cm^{-1}} 8.07 c m − 1 に当たる。
選択則:反転対称と相互排他
分極率が振動で変わらなければ(∂ α / ∂ Q = 0 \partial\alpha/\partial Q=0 ∂ α / ∂ Q = 0 )、式 (3.3.10) の第2項は無く、その振動は Raman 散乱に現れない。一方、赤外線を直接吸収するには、振動で双極子モーメントが変わる必要がある(∂ p / ∂ Q ≠ 0 \partial\bm{p}/\partial Q\ne\bm{0} ∂ p / ∂ Q = 0 )。この2つの条件は、反転対称で整理することができる。
結晶(または分子)が反転中心を持つとき、k = 0 \bm{k}=\bm{0} k = 0 の基準振動は反転で Q → + Q Q\to+Q Q → + Q になるもの(偶)と Q → − Q Q\to-Q Q → − Q になるもの(奇)に分かれる。双極子モーメントはベクトルで反転で符号を変え、分極率は2階のテンソルで反転で符号を変えない。奇の振動では、反転した配置は − Q -Q − Q の配置と同じなのでα ( − Q ) = α ( Q ) \alpha(-Q)=\alpha(Q) α ( − Q ) = α ( Q ) 、したがって ∂ α / ∂ Q = 0 \partial\alpha/\partial Q=0 ∂ α / ∂ Q = 0 で Raman 散乱に現れない。偶の振動では、反転した配置が元の配置そのものなので p ( Q ) = − p ( Q ) \bm{p}(Q)=-\bm{p}(Q) p ( Q ) = − p ( Q ) 、すなわち振動に沿って双極子は 0 のままで ∂ p / ∂ Q = 0 \partial\bm{p}/\partial Q=\bm{0} ∂ p / ∂ Q = 0 となり、赤外吸収に現れない。反転中心のある系では、同じ振動が赤外と Raman の両方に現れることはない。これを相互排他律 と呼ぶ。
ダイヤモンド構造の Si では、2つの原子の中点が反転中心で、Γ \Gamma Γ 点の光学振動(2つの原子が逆向きに動く)は反転で原子が入れ替わって偶になる。したがって Si の Γ \Gamma Γ 点の光学フォノンは Raman 散乱に現れ、1個のフォノンによる赤外吸収は無い。室温の Si の Raman 線は 520.2 ± 0.5 c m − 1 520.2\pm0.5\,\mathrm{cm^{-1}} 520.2 ± 0.5 c m − 1 に観測されている(J. H. Parker, Jr., D. W. Feldman, and M. Ashkin (1967), Phys. Rev. 155 , 712)。塩化ナトリウム型の結晶では各原子の位置が反転中心で、Γ \Gamma Γ 点の光学振動は奇になる。こちらは赤外で強く吸収し、1個のフォノンの Raman 線は現れない。
例題 3.3.3 CO2 _2 2 の振動の赤外・Raman 活性 ★★
オリジナル 直線分子 O–C–O の4つの振動(対称伸縮、逆対称伸縮、2方向の曲げ)を考える(基準振動の求め方は 解析力学から様々な理論分野を組み立てる 第2部第1章 )。C の位置が反転中心である。
(1) 対称伸縮(2つの O が C から逆向きに同じだけ離れる)と逆対称伸縮(2つの O が同じ向きに動き、C が逆向きに動く)は、反転で偶か奇か。
(2) 各振動が赤外・Raman のどちらに現れるかを述べよ。
(3) O の電荷を − δ -\delta − δ 、C の電荷を + 2 δ +2\delta + 2 δ の点電荷とみなすと、逆対称伸縮で双極子が変わり、対称伸縮で変わらないことを式で確かめよ。
解答 (1) O の変位を u 1 u_1 u 1 、u 3 u_3 u 3 、C の変位を u 2 u_2 u 2 (軸方向)とする。反転は位置 x → − x x\to-x x → − x と変位 u → − u u\to-u u → − u を同時に行い、2つの O を入れ替える:( u 1 , u 2 , u 3 ) → ( − u 3 , − u 2 , − u 1 ) (u_1,u_2,u_3)\to(-u_3,-u_2,-u_1) ( u 1 , u 2 , u 3 ) → ( − u 3 , − u 2 , − u 1 ) 。対称伸縮 ( − d , 0 , d ) (-d,0,d) ( − d , 0 , d ) はこれで ( − d , 0 , d ) (-d,0,d) ( − d , 0 , d ) に戻るので偶。逆対称伸縮 ( d , − d ′ , d ) (d,-d',d) ( d , − d ′ , d ) は ( − d , d ′ , − d ) (-d,d',-d) ( − d , d ′ , − d ) になるので奇。曲げ(O が同じ向き、C が逆向きに軸と垂直に動く)も同じ形で奇である。
(2) 相互排他律から、偶の対称伸縮は Raman 散乱だけに、奇の逆対称伸縮と曲げは赤外吸収だけに現れる。
(3) 軸方向の双極子は p = ∑ j q j x j p=\sum_jq_jx_j p = ∑ j q j x j の変化 Δ p = − δ u 1 + 2 δ u 2 − δ u 3 \Delta p=-\delta u_1+2\delta u_2-\delta u_3 Δ p = − δ u 1 + 2 δ u 2 − δ u 3 。対称伸縮 ( − d , 0 , d ) (-d,0,d) ( − d , 0 , d ) では Δ p = δ d − δ d = 0 \Delta p=\delta d-\delta d=0 Δ p = δ d − δ d = 0 。逆対称伸縮 ( d , − d ′ , d ) (d,-d',d) ( d , − d ′ , d ) では Δ p = − 2 δ d − 2 δ d ′ ≠ 0 \Delta p=-2\delta d-2\delta d'\ne0 Δ p = − 2 δ d − 2 δ d ′ = 0 で、双極子が振動とともに振れるので赤外を吸収する。
例題 3.3.4 Brillouin シフトと Raman 線で温度を読む ★★
オリジナル
(1) 屈折率 n = 1.5 n=1.5 n = 1.5 、縦波の音速 v = 6000 m s − 1 v=6000\,\mathrm{m\,s^{-1}} v = 6000 m s − 1 の透明な結晶に、真空中の波長 532 n m 532\,\mathrm{nm} 532 nm の光を当て、真後ろ(θ = π \theta=\pi θ = π )に散乱された光を測る。Brillouin シフト Ω B / 2 π \Omega_{\mathrm B}/2\pi Ω B /2 π と、関わるフォノンの波数を求めよ。
(2) Si の Raman 線 520.2 c m − 1 520.2\,\mathrm{cm^{-1}} 520.2 c m − 1 のフォノンのエネルギー(meV)と周波数を求め、波長 532 n m 532\,\mathrm{nm} 532 nm の光で測ったときの Stokes 線と anti-Stokes 線の波長を求めよ。
(3) 300 K 300\,\mathrm K 300 K で anti-Stokes 線と Stokes 線の光子の数の比を求めよ。双極子放射の出す光子の数が周波数の3乗に比例すること(放射の強さが4乗、光子1個が1乗)を含めよ。
解答 (1) 式 (3.3.11) で sin ( θ / 2 ) = 1 \sin(\theta/2)=1 sin ( θ /2 ) = 1 、∣ k p h ∣ = 2 n ⋅ 2 π / λ = 3.54 × 10 7 m − 1 \lvert \bm{k}_{\mathrm{ph}}\rvert =2n\cdot2\pi/\lambda=3.54\times10^{7}\,\mathrm{m^{-1}} ∣ k ph ∣ = 2 n ⋅ 2 π / λ = 3.54 × 1 0 7 m − 1 。Ω B / 2 π = v ∣ k p h ∣ / 2 π = 2 n v / λ = 33.8 G H z \Omega_{\mathrm B}/2\pi=v\lvert \bm{k}_{\mathrm{ph}}\rvert /2\pi=2nv/\lambda=33.8\,\mathrm{GHz} Ω B /2 π = v ∣ k ph ∣ /2 π = 2 n v / λ = 33.8 GHz 。フォノンのエネルギーは 0.14 m e V 0.14\,\mathrm{meV} 0.14 meV で、光の 2.3 e V 2.3\,\mathrm{eV} 2.3 eV の 6 × 10 − 5 6\times10^{-5} 6 × 1 0 − 5 にすぎない。
(2) ℏ Ω = h c ν ˉ = 64.5 m e V \hbar\Omega=hc\bar\nu=64.5\,\mathrm{meV} ℏΩ = h c ν ˉ = 64.5 meV 、Ω / 2 π = c ν ˉ = 15.6 T H z \Omega/2\pi=c\bar\nu=15.6\,\mathrm{THz} Ω/2 π = c ν ˉ = 15.6 THz 。1 / λ s = 1 / λ i ∓ ν ˉ 1/\lambda_{\mathrm s}=1/\lambda_{\mathrm i}\mp\bar\nu 1/ λ s = 1/ λ i ∓ ν ˉ から、Stokes 線は 547.1 n m 547.1\,\mathrm{nm} 547.1 nm 、anti-Stokes 線は 517.7 n m 517.7\,\mathrm{nm} 517.7 nm 。
(3) β ℏ Ω = 2.495 \beta\hbar\Omega=2.495 β ℏΩ = 2.495 で e − β ℏ Ω = 0.0825 \mathrm{e}^{-\beta\hbar\Omega}=0.0825 e − β ℏΩ = 0.0825 。周波数の比 ( ω A S / ω S ) 3 = ( 547.1 / 517.7 ) 3 = 1.181 (\omega_{\mathrm{AS}}/\omega_{\mathrm S})^3=(547.1/517.7)^3=1.181 ( ω AS / ω S ) 3 = ( 547.1/517.7 ) 3 = 1.181 を掛けて、光子の数の比は 0.0974 0.0974 0.0974 。測った比から ( ω A S / ω S ) 3 (\omega_{\mathrm{AS}}/\omega_{\mathrm S})^3 ( ω AS / ω S ) 3 を割り戻せば e − β ℏ Ω \mathrm{e}^{-\beta\hbar\Omega} e − β ℏΩ が分かり、試料の温度を接触せずに読むことができる。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com