ハミルトニアンの選び方で固体を考える 第4部 第8章 §4.8.3 ★★★ 跳び移りの速さは t h 2 t_{\mathrm h}^2 t h 2 と、2つの放物線の交点を熱の揺らぎで通る確率の積。駆動力 − Δ G -\Delta G − Δ G が λ \lambda λ を超えると遅くなる
黄金律で跳び移りを数える
t h t_{\mathrm h} t h が小さいとき、電子が分子1から分子2へ移る速さは、t h t_{\mathrm h} t h を摂動とする Fermi の黄金律(解析力学から様々な理論分野を組み立てる 第4部第12章 )で求まる。始めの状態は「電子が分子1にいて、原子核の座標が U 1 U_1 U 1 の熱平衡にある」状態、終わりの状態は「電子が分子2にいる」状態である。原子核を古典的に扱うと、電子が移る瞬間に原子核の座標 q q q は変わらないので、エネルギーの保存は U 1 ( q ) = U 2 ( q ) U_1(q)=U_2(q) U 1 ( q ) = U 2 ( q ) 、すなわち q q q が交点にあることを要求する。黄金律は
k E T = 2 π ℏ t h 2 ⟨ δ ( U 2 ( q ) − U 1 ( q ) ) ⟩ 1 (4.8.10) k_{\mathrm{ET}}=\frac{2\pi}{\hbar}t_{\mathrm h}^2\,\Bigl\langle\delta\bigl(U_2(q)-U_1(q)\bigr)\Bigr\rangle_1\tag{4.8.10} k ET = ℏ 2 π t h 2 ⟨ δ ( U 2 ( q ) − U 1 ( q ) ) ⟩ 1 ( 4.8.10 )
となる。⟨ ⋯ ⟩ 1 \langle\cdots\rangle_1 ⟨ ⋯ ⟩ 1 は U 1 U_1 U 1 の Boltzmann 分布 ∝ e − λ q 2 / k B T \propto\mathrm{e}^{-\lambda q^2/ k_{\mathrm B} T} ∝ e − λ q 2 / k B T での平均である。
エネルギーの差 Δ U ( q ) = U 2 − U 1 = λ ( 1 − 2 q ) + Δ G \Delta U(q)=U_2-U_1=\lambda(1-2q)+\Delta G Δ U ( q ) = U 2 − U 1 = λ ( 1 − 2 q ) + Δ G は q q q の1次式で、q q q は平均 0 0 0 ・分散 k B T / ( 2 λ ) k_{\mathrm B} T/(2\lambda) k B T / ( 2 λ ) の Gauss 分布に従う。したがって Δ U \Delta U Δ U も Gauss 分布に従い、平均は λ + Δ G \lambda+\Delta G λ + Δ G 、分散は 4 λ 2 ⋅ k B T / ( 2 λ ) = 2 λ k B T 4\lambda^2\cdot k_{\mathrm B} T/(2\lambda)=2\lambda k_{\mathrm B} T 4 λ 2 ⋅ k B T / ( 2 λ ) = 2 λ k B T である。Gauss 分布の Δ U = 0 \Delta U=0 Δ U = 0 での密度を書くと
k E T = 2 π ℏ t h 2 1 4 π λ k B T exp [ − ( λ + Δ G ) 2 4 λ k B T ] (4.8.11) k_{\mathrm{ET}}=\frac{2\pi}{\hbar}\,t_{\mathrm h}^2\,\frac{1}{\sqrt{4\pi\lambda k_{\mathrm B} T}}\,\exp\Bigl[-\frac{(\lambda+\Delta G)^2}{4\lambda k_{\mathrm B} T}\Bigr]\tag{4.8.11} k ET = ℏ 2 π t h 2 4 π λ k B T 1 exp [ − 4 λ k B T ( λ + Δ G ) 2 ] ( 4.8.11 )
を得る。これが Marcus の速度式 である(R. A. Marcus (1956), J. Chem. Phys. 24 , 966。溶液中の酸化還元反応のために作られ、溶媒の分極の揺らぎを反応座標にとった)。指数の肩は交点の高さ (4.8.9) を k B T k_{\mathrm B} T k B T で割ったものに等しく、前の因子 ( 4 π λ k B T ) − 1 / 2 (4\pi\lambda k_{\mathrm B} T)^{-1/2} ( 4 π λ k B T ) − 1/2 は Gauss 分布の幅から来る規格化である。
逆転領域
(4.8.11) は Δ G \Delta G Δ G の Gauss 関数で、Δ G = − λ \Delta G=-\lambda Δ G = − λ で最大になる。駆動力 − Δ G -\Delta G − Δ G を 0 0 0 から増やすと、はじめは障壁が下がって速くなるが、− Δ G > λ -\Delta G>\lambda − Δ G > λ では再び遅くなる。これを Marcus の逆転領域 と呼ぶ。図 4.8.1 で U 2 U_2 U 2 を下げていくと、交点は q = 0 q=0 q = 0 を通り過ぎて左側へ移り、U 1 U_1 U 1 の上で再び高くなる。「エネルギーを多く解放する反応ほど速い」という直観に反するこの予言は、剛直なスペーサーでつないだ分子の中の電子移動で確かめられた(J. R. Miller, L. T. Calcaterra and G. L. Closs (1984), J. Am. Chem. Soc. 106 , 3047)。
非断熱と断熱
(4.8.10) は、原子核が交点を1回通るたびに電子が移る確率が小さいことを前提にしている。原子核が交点を速さ q ˙ \dot q q ˙ で通るとき、2つの準位の差は ∣ d Δ U / d t ∣ = 2 λ ∣ q ˙ ∣ \lvert \mathrm{d}\Delta U/\mathrm{d} t\rvert =2\lambda\lvert \dot q\rvert ∣ d Δ U / d t ∣ = 2 λ ∣ q ˙ ∣ の速さで変わる。このときの乗り移りの確率は Landau–Zener の式(解析力学から様々な理論分野を組み立てる 第4部第9章 )で 1 − exp [ − 2 π t h 2 / ( ℏ ∣ d Δ U / d t ∣ ) ] 1-\exp[-2\pi t_{\mathrm h}^2/(\hbar\lvert \mathrm{d}\Delta U/\mathrm{d} t\rvert )] 1 − exp [ − 2 π t h 2 / ( ℏ ∣ d Δ U / d t ∣)] である。熱的な速さ ∣ q ˙ ∣ = ω v k B T / 2 λ \lvert \dot q\rvert =\omega_{\mathrm v}\sqrt{ k_{\mathrm B} T/2\lambda} ∣ q ˙ ∣ = ω v k B T /2 λ (λ q 2 \lambda q^2 λ q 2 を 1 2 M ω v 2 x 2 \frac12M\omega_{\mathrm v}^2x^2 2 1 M ω v 2 x 2 と読み、運動エネルギーの平均を 1 2 k B T \frac12 k_{\mathrm B} T 2 1 k B T とした)を入れると、指数の中は
γ L Z = 2 π t h 2 ℏ ω v 2 λ k B T (4.8.12) \gamma_{\mathrm{LZ}}=\frac{2\pi t_{\mathrm h}^2}{\hbar\omega_{\mathrm v}\sqrt{2\lambda k_{\mathrm B} T}}\tag{4.8.12} γ LZ = ℏ ω v 2 λ k B T 2 π t h 2 ( 4.8.12 )
となる。γ L Z ≪ 1 \gamma_{\mathrm{LZ}}\ll1 γ LZ ≪ 1 なら1回の通過での乗り移りは γ L Z \gamma_{\mathrm{LZ}} γ LZ の確率でしか起きず、速さは t h 2 t_{\mathrm h}^2 t h 2 に比例する (4.8.11) の形になる(非断熱 の移動)。γ L Z ≳ 1 \gamma_{\mathrm{LZ}}\gtrsim1 γ LZ ≳ 1 では交点を通るたびにほぼ確実に乗り移り、速さは t h t_{\mathrm h} t h によらず、原子核が交点へ来る頻度 ( ω v / 2 π ) e − E a / k B T (\omega_{\mathrm v}/2\pi)\mathrm{e}^{-E_{\mathrm a}/ k_{\mathrm B} T} ( ω v /2 π ) e − E a / k B T で頭打ちになる(断熱 の移動)。
例題 4.8.3 Marcus の速さと逆転領域の数値 ★★★
オリジナル t h = 5 m e V t_{\mathrm h}=5\,\mathrm{meV} t h = 5 meV 、λ = 0.2 e V \lambda=0.2\,\mathrm{eV} λ = 0.2 eV 、T = 300 K T=300\,\mathrm K T = 300 K とする。
(1) Δ G = 0 \Delta G=0 Δ G = 0 の k E T k_{\mathrm{ET}} k ET を求めよ。
(2) Δ G = − λ \Delta G=-\lambda Δ G = − λ の k E T k_{\mathrm{ET}} k ET と、(1) との比を求めよ。
(3) Δ G = − 0.4 e V \Delta G=-0.4\,\mathrm{eV} Δ G = − 0.4 eV と − 0.6 e V -0.6\,\mathrm{eV} − 0.6 eV の k E T k_{\mathrm{ET}} k ET を求め、逆転領域の意味を述べよ。
(4) 低い振動の ℏ ω v = 10 m e V \hbar\omega_{\mathrm v}=10\,\mathrm{meV} ℏ ω v = 10 meV を使って γ L Z \gamma_{\mathrm{LZ}} γ LZ を求め、非断熱の式が使えるかを判定せよ。t h = 20 m e V t_{\mathrm h}=20\,\mathrm{meV} t h = 20 meV ならどうか。
解答 (1) k B T = 25.85 m e V k_{\mathrm B} T=25.85\,\mathrm{meV} k B T = 25.85 meV 、4 π λ k B T = 0.255 e V \sqrt{4\pi\lambda k_{\mathrm B} T}=0.255\,\mathrm{eV} 4 π λ k B T = 0.255 eV 、e − λ / 4 k B T = e − 1.934 = 0.1446 \mathrm{e}^{-\lambda/4 k_{\mathrm B} T}=\mathrm{e}^{-1.934}=0.1446 e − λ /4 k B T = e − 1.934 = 0.1446 、2 π / ℏ = 9.55 × 10 15 e V − 1 s − 1 2\pi/\hbar=9.55\times10^{15}\,\mathrm{eV^{-1}s^{-1}} 2 π /ℏ = 9.55 × 1 0 15 e V − 1 s − 1 なのでk E T = 9.55 × 10 15 × ( 5 × 10 − 3 ) 2 / 0.255 × 0.1446 = 1.35 × 10 11 s − 1 k_{\mathrm{ET}}=9.55\times10^{15}\times(5\times10^{-3})^2/0.255\times0.1446=1.35\times10^{11}\,\mathrm{s^{-1}} k ET = 9.55 × 1 0 15 × ( 5 × 1 0 − 3 ) 2 /0.255 × 0.1446 = 1.35 × 1 0 11 s − 1 。1回の跳び移りに平均 7.4 p s 7.4\,\mathrm{ps} 7.4 ps かかる。
(2) 指数が 1 1 1 になるので k E T = 9.36 × 10 11 s − 1 k_{\mathrm{ET}}=9.36\times10^{11}\,\mathrm{s^{-1}} k ET = 9.36 × 1 0 11 s − 1 、(1) の e 1.934 = 6.92 \mathrm{e}^{1.934}=6.92 e 1.934 = 6.92 倍である。
(3) Δ G = − 0.4 e V \Delta G=-0.4\,\mathrm{eV} Δ G = − 0.4 eV では ( λ + Δ G ) 2 = ( 0.2 ) 2 (\lambda+\Delta G)^2=(0.2)^2 ( λ + Δ G ) 2 = ( 0.2 ) 2 が Δ G = 0 \Delta G=0 Δ G = 0 と同じなので、k E T = 1.35 × 10 11 s − 1 k_{\mathrm{ET}}=1.35\times10^{11}\,\mathrm{s^{-1}} k ET = 1.35 × 1 0 11 s − 1 に戻る。Δ G = − 0.6 e V \Delta G=-0.6\,\mathrm{eV} Δ G = − 0.6 eV では障壁が ( 0.4 ) 2 / ( 0.8 ) = 0.2 e V (0.4)^2/(0.8)=0.2\,\mathrm{eV} ( 0.4 ) 2 / ( 0.8 ) = 0.2 eV に上がり、k E T = 4.09 × 10 8 s − 1 k_{\mathrm{ET}}=4.09\times10^{8}\,\mathrm{s^{-1}} k ET = 4.09 × 1 0 8 s − 1 と、最大値の 2300 分の1になる。解放するエネルギーが λ \lambda λ より大きいと、そのエネルギーを原子核の振動として受け取る配置が平衡から遠くなり、そこへ揺らぎで行く確率が下がる。
(4) γ L Z = 2 π ( 0.005 ) 2 / ( 0.010 × 2 × 0.2 × 0.02585 ) = 0.154 \gamma_{\mathrm{LZ}}=2\pi(0.005)^2/(0.010\times\sqrt{2\times0.2\times0.02585})=0.154 γ LZ = 2 π ( 0.005 ) 2 / ( 0.010 × 2 × 0.2 × 0.02585 ) = 0.154 で、1 より小さいので非断熱の式は使える(誤差は1割程度)。t h = 20 m e V t_{\mathrm h}=20\,\mathrm{meV} t h = 20 meV では γ L Z \gamma_{\mathrm{LZ}} γ LZ は 16 倍の 2.47 2.47 2.47 になり、断熱の側に入る。このとき (4.8.11) は t h 2 t_{\mathrm h}^2 t h 2 に比例して速さを過大に見積もるので、上限の目安 ( ω v / 2 π ) e − λ / 4 k B T = 3.5 × 10 11 s − 1 (\omega_{\mathrm v}/2\pi)\mathrm{e}^{-\lambda/4 k_{\mathrm B} T}=3.5\times10^{11}\,\mathrm{s^{-1}} ( ω v /2 π ) e − λ /4 k B T = 3.5 × 1 0 11 s − 1 を使う(ここでは障壁の t h t_{\mathrm h} t h だけの低下も無視した)。結晶性のよい有機半導体の多くは、この2つの領域の境目にある。
ほかの巻では
近似の数理 第7章 :断熱定理とそのギャップの前提が読め、Landau–Zener の断熱と非断熱の分かれ目の数学の側が分かる
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com