ハミルトニアンの選び方で固体を考える 第4部 第8章 §4.8.6 ★★ 静的な乱れ・振動の量子性・ホッピングの揺らぎ・キャリアどうしの相互作用を削った。実験は移動度の温度変化と逆転領域で確かめる
何を削ったか
この章の結果は、(4.8.1) の H H H を2分子に切り、振動を古典的に扱って黄金律を当てたことに頼っている。そこで削ったものを、輸送に大きな順に並べる。
1つめは静的な乱れ である。薄膜の中の分子は向きも間隔もそろっておらず、分子ごとのサイトのエネルギー ε a \varepsilon _{\mathrm a} ε a が幅 σ d i s \sigma_{\mathrm{dis}} σ dis ほどばらつく。キャリアは低いエネルギーの分子に捕まり、そこから抜け出すには余分な活性化が要る。サイトのエネルギーを Gauss 分布とした数値計算から、移動度が exp [ − ( 2 σ d i s / 3 k B T ) 2 ] \exp[-(2\sigma_{\mathrm{dis}}/3 k_{\mathrm B} T)^2] exp [ − ( 2 σ dis /3 k B T ) 2 ] の形で温度に依存することが提案されている(H. Bässler (1993), Phys. Status Solidi B 175 , 15)。この依存性は (4.8.14) の e − λ / 4 k B T \mathrm{e}^{-\lambda/4 k_{\mathrm B} T} e − λ /4 k B T と形が違い、両者は温度変化の曲がり方で区別される。
2つめは振動の量子性 である。例題 4.8.2 で見たとおり、分子の再配置の一部は ℏ ω v ≫ k B T \hbar\omega_{\mathrm v}\gg k_{\mathrm B} T ℏ ω v ≫ k B T の骨格の振動による。この量子的な振動を1本残し、残りを古典的に扱うと、(4.8.11) は振動の量子数 n n n についての和
k E T = 2 π ℏ t h 2 4 π λ s k B T ∑ n = 0 ∞ e − S S n n ! exp [ − ( Δ G + λ s + n ℏ ω q ) 2 4 λ s k B T ] (4.8.20) k_{\mathrm{ET}}=\frac{2\pi}{\hbar}\frac{t_{\mathrm h}^2}{\sqrt{4\pi\lambda_{\mathrm s} k_{\mathrm B} T}}\sum_{n=0}^\infty\mathrm{e}^{-S}\frac{S^n}{n!}\exp\Bigl[-\frac{(\Delta G+\lambda_{\mathrm s}+n\hbar\omega_{\mathrm q})^2}{4\lambda_{\mathrm s} k_{\mathrm B} T}\Bigr]\tag{4.8.20} k ET = ℏ 2 π 4 π λ s k B T t h 2 n = 0 ∑ ∞ e − S n ! S n exp [ − 4 λ s k B T ( Δ G + λ s + n ℏ ω q ) 2 ] ( 4.8.20 )
に置き換わる(λ s \lambda_{\mathrm s} λ s は古典的な振動の再配置エネルギー、S = λ q / ℏ ω q S=\lambda_{\mathrm q}/\hbar\omega_{\mathrm q} S = λ q /ℏ ω q は量子的な振動の Huang–Rhys 因子)。この式は Marcus–Levich–Jortner の式と呼ばれる(J. Jortner (1976), J. Chem. Phys. 64 , 4860)。逆転領域では、余ったエネルギーを量子的な振動の励起(n ≥ 1 n\ge1 n ≥ 1 )として受け取る経路が開くので、速さの落ち方は (4.8.11) よりずっと緩やかになる(図 4.8.4 )。ℏ ω q ≪ k B T \hbar\omega_{\mathrm q}\ll k_{\mathrm B} T ℏ ω q ≪ k B T の極限で量子的な振動の熱励起まで入れた式に戻すと、(4.8.11) の λ = λ s + λ q \lambda=\lambda_{\mathrm s}+\lambda_{\mathrm q} λ = λ s + λ q に一致する(verify で 3 % 3\% 3% 以内に一致することを確かめた)。
3つめはホッピングの揺らぎ である。分子どうしの低い振動(ℏ ω ∼ 5 \hbar\omega\sim5 ℏ ω ∼ 5 〜10 m e V 10\,\mathrm{meV} 10 meV )は室温で大きく揺れ、t h t_{\mathrm h} t h そのものがその平均と同じ程度に揺らぐ。この非対角の乱れが電子を局在させ、結晶性のよい有機半導体の室温の輸送を決めているという見方が提案されている(A. Troisi and G. Orlandi (2006), Phys. Rev. Lett. 96 , 086601)。これは Holstein 模型の対角の結合では表せず、SSH 型の結合(第2部第8章 )を分子どうしの間に戻す必要がある。
4つめはキャリアどうしの相互作用と電極 である。例題 4.8.5 の密度は分子あたり 0.005 0.005 0.005 個で、キャリアどうしは遠い。しかし FET の電流には、電極から半導体への注入の抵抗と、界面の罠が加わる。(4.8.16) から読んだ移動度は、これらを含んだ見かけの値である。
図 4.8.4 跳び移りの速さの駆動力依存(t h = 5 m e V t_{\mathrm h}=5\,\mathrm{meV} t h = 5 meV 、T = 300 K T=300\,\mathrm K T = 300 K 、縦軸は Marcus の最大値 k max k_{\max} k m a x で割った値)。実線は古典の Marcus の式 (4.8.11) (λ = 0.2 e V \lambda=0.2\,\mathrm{eV} λ = 0.2 eV )。点線は再配置の半分を ℏ ω q = 0.18 e V \hbar\omega_{\mathrm q}=0.18\,\mathrm{eV} ℏ ω q = 0.18 eV の量子的な振動の分とした (4.8.20) (λ s = λ q = 0.1 e V \lambda_{\mathrm s}=\lambda_{\mathrm q}=0.1\,\mathrm{eV} λ s = λ q = 0.1 eV )。逆転領域で点線の落ち方は実線よりはるかに緩い。
裏付けの札
近似の裏付け
型 B(極限) | 小さいもの t h / λ → 0 t_{\mathrm h}/\lambda\to0 t h / λ → 0 と γ L Z → 0 \gamma_{\mathrm{LZ}}\to0 γ LZ → 0 (ホッピングの2次の黄金律。非断熱の移動) | 証明 条件付き(有限準位の系が浴に弱く結合する極限で黄金律の速さが出ることは Davies 1974。γ L Z ≪ 1 \gamma_{\mathrm{LZ}}\ll1 γ LZ ≪ 1 の非断熱の領域に限る。振動が毎回の跳び移りの間に平衡へ戻ることは次の札の前提) | 行き先 近似の数理 第14章
近似の裏付け
型 B(極限) | 小さいもの ℏ ω v / k B T → 0 \hbar\omega_{\mathrm v}/ k_{\mathrm B} T\to0 ℏ ω v / k B T → 0 (原子核の振動を古典的に扱う) | 証明 条件付き(振動が調和的で、電子が乗っても振動数が変わらず原点だけがずれる場合は、量子の式の極限として (4.8.11) が出ることを verify で数値的に確かめた。非調和性と振動数の変化は含まない) | 行き先 近似の数理 第7章
近似の裏付け
型 B(極限) | 小さいもの k E T / ω v → 0 k_{\mathrm{ET}}/\omega_{\mathrm v}\to0 k ET / ω v → 0 (跳び移りを記憶の無い乱歩として足し、Einstein の関係で移動度に直す) | 証明 要確認(振動の緩和が跳び移りより十分速いという仮定で、一般の分子性固体での誤差の評価は原論文で確かめていない) | 行き先 近似の数理 第14章
近似の裏付け
型 B(極限) | 小さいもの 分子間の相互作用の到達距離 → ∞ \to\infty → ∞ (Maier–Saupe・Landau の平均場) | 証明 条件付き(Kac 型の長距離の極限で平均場の自由エネルギーが厳密になることは Lebowitz–Penrose 1966 が示した。液晶の向きの模型への適用の前提と、短距離の相互作用での揺らぎの補正は要確認) | 行き先 近似の数理 第6章 ・近似の数理 第13章
何の実験の量で確かめるか
最も直接の量は移動度の温度変化 である。純度を極めたナフタレンの単結晶では、飛行時間法(膜の一方の面で光で作ったキャリアが反対の面に届くまでの時間を測る方法)で測った正孔の移動度が、冷やすほど T − n T^{-n} T − n で増え続け、10 K 10\,\mathrm K 10 K 付近で数百 c m 2 V − 1 s − 1 \mathrm{cm^2\,V^{-1}s^{-1}} c m 2 V − 1 s − 1 に達する(W. Warta and N. Karl (1985), Phys. Rev. B 32 , 1172)。これは (4.8.2) の下限を大きく上回るバンド的な振る舞いで、乱れた薄膜の活性化型の移動度とは逆の向きである。同じ分子の物質でも、結晶の質と温度によって4.8.1 節の2つの描像の間を行き来することが、この比較から読み取れる。
跳び移りの側では、(4.8.11) の駆動力依存、とくに逆転領域が、分子の中の電子移動の速さを駆動力を変えて測った実験で確かめられた(4.8.3 節の Miller–Calcaterra–Closs)。再配置エネルギー λ \lambda λ は、分子の中性と陽イオンの構造を量子化学で計算するか、光電子分光のスペクトルの振動構造の幅から見積もることができる。強誘電体では分極と電場の履歴曲線(P P P –E E E 曲線)から P s P_{\mathrm s} P s と抗電場を、液晶では複屈折や NMR の分裂から秩序変数 s s s を読む。
例題 4.8.8 量子的な振動で逆転領域が緩む・飛行時間法 ★★
オリジナル t h = 5 m e V t_{\mathrm h}=5\,\mathrm{meV} t h = 5 meV 、T = 300 K T=300\,\mathrm K T = 300 K で、再配置を λ s = 0.1 e V \lambda_{\mathrm s}=0.1\,\mathrm{eV} λ s = 0.1 eV (古典)と λ q = 0.1 e V \lambda_{\mathrm q}=0.1\,\mathrm{eV} λ q = 0.1 eV (ℏ ω q = 0.18 e V \hbar\omega_{\mathrm q}=0.18\,\mathrm{eV} ℏ ω q = 0.18 eV の量子的な振動)に分ける。
(1) Δ G = 0 \Delta G=0 Δ G = 0 で、(4.8.20) の n = 0 n=0 n = 0 の項の重みと障壁を求めよ。(4.8.20) 全体の k E T k_{\mathrm{ET}} k ET と、λ = 0.2 e V \lambda=0.2\,\mathrm{eV} λ = 0.2 eV の Marcus の式との比を求め、比が 1 より大きい理由を述べよ。
(2) − Δ G = 0.6 e V -\Delta G=0.6\,\mathrm{eV} − Δ G = 0.6 eV と 0.8 e V 0.8\,\mathrm{eV} 0.8 eV で、(4.8.20) と Marcus の式の比を求めよ。
(3) 例題 4.8.4 の移動度の膜(厚さ 10 μ m 10\,\mu\mathrm m 10 μ m )に 10 V 10\,\mathrm V 10 V をかけたとき、飛行時間法でキャリアが膜を横切る時間を求めよ。
解答 (1) S = λ q / ℏ ω q = 0.556 S=\lambda_{\mathrm q}/\hbar\omega_{\mathrm q}=0.556 S = λ q /ℏ ω q = 0.556 で、n = 0 n=0 n = 0 の重みは e − S = 0.574 \mathrm{e}^{-S}=0.574 e − S = 0.574 、障壁は λ s / 4 = 25 m e V \lambda_{\mathrm s}/4=25\,\mathrm{meV} λ s /4 = 25 meV である。和の全体は k E T = 2.89 × 10 11 s − 1 k_{\mathrm{ET}}=2.89\times10^{11}\,\mathrm{s^{-1}} k ET = 2.89 × 1 0 11 s − 1 で、Marcus の式(1.35 × 10 11 s − 1 1.35\times10^{11}\,\mathrm{s^{-1}} 1.35 × 1 0 11 s − 1 )の 2.14 2.14 2.14 倍になる。量子的な振動は室温で励起されないので、その分の再配置は熱の揺らぎで障壁を越えるのでなく、零点振動の波動関数の重なり(Franck–Condon 因子 e − S \mathrm{e}^{-S} e − S )として支払われる。障壁の高さが λ / 4 = 50 m e V \lambda/4=50\,\mathrm{meV} λ /4 = 50 meV から λ s / 4 = 25 m e V \lambda_{\mathrm s}/4=25\,\mathrm{meV} λ s /4 = 25 meV に下がる得が、重なりの損 0.574 0.574 0.574 を上回る。
(2) − Δ G = 0.6 e V -\Delta G=0.6\,\mathrm{eV} − Δ G = 0.6 eV で比は 89 89 89 ((4.8.20) は 3.63 × 10 10 s − 1 3.63\times10^{10}\,\mathrm{s^{-1}} 3.63 × 1 0 10 s − 1 、Marcus は 4.09 × 10 8 s − 1 4.09\times10^{8}\,\mathrm{s^{-1}} 4.09 × 1 0 8 s − 1 )、0.8 e V 0.8\,\mathrm{eV} 0.8 eV では 1.8 × 10 5 1.8\times10^{5} 1.8 × 1 0 5 である。余ったエネルギーを n = 2 , 3 , 4 n=2,3,4 n = 2 , 3 , 4 の振動の励起で受け取ると、古典の部分の障壁 ( Δ G + λ s + n ℏ ω q ) 2 / 4 λ s (\Delta G+\lambda_{\mathrm s}+n\hbar\omega_{\mathrm q})^2/4\lambda_{\mathrm s} ( Δ G + λ s + n ℏ ω q ) 2 /4 λ s を小さく保てるからである。
(3) 平均の速さは μ e V / L \mu_{\mathrm e}V/L μ e V / L なので、横切る時間は t T = L 2 / ( μ e V ) = ( 10 − 5 ) 2 / ( 1.31 × 10 − 6 × 10 ) = 7.6 μ s t_{\mathrm T}=L^2/(\mu_{\mathrm e}V)=(10^{-5})^2/(1.31\times10^{-6}\times10)=7.6\,\mu\mathrm s t T = L 2 / ( μ e V ) = ( 1 0 − 5 ) 2 / ( 1.31 × 1 0 − 6 × 10 ) = 7.6 μ s 。この時間を電流の波形から読み、μ e = L 2 / ( V t T ) \mu_{\mathrm e}=L^2/(Vt_{\mathrm T}) μ e = L 2 / ( V t T ) として移動度を出す。FET と違ってゲートで電荷をためないので、膜の内部(界面でない所)の移動度を測ることができる。
授業で習った形との対応
学部の固体物理の授業では、電気伝導を σ = n e μ e \sigma=ne\mu_{\mathrm e} σ = n e μ e 、μ e = e τ / m ∗ \mu_{\mathrm e}=e\tau/m^* μ e = e τ / m ∗ と書く。この章の跳び移りでも σ = n e μ e \sigma=ne\mu_{\mathrm e} σ = n e μ e はそのまま使えるが、移動度の中身が置き換わる:
バンド: μ e = e τ m ∗ , 跳び移り: μ e = e k B T k E T a 2 . \text{バンド:}\ \mu_{\mathrm e}=\frac{e\tau}{m^*},\qquad\text{跳び移り:}\ \mu_{\mathrm e}=\frac{e}{ k_{\mathrm B} T}\,k_{\mathrm{ET}}a^2 . バンド: μ e = m ∗ e τ , 跳び移り: μ e = k B T e k ET a 2 .
左は「加速と散乱のつり合い」、右は「乱歩の拡散を Einstein の関係で読んだもの」である。右の式に τ \tau τ も m ∗ m^* m ∗ も現れないのは、キャリアが加速される間もなく次の分子へ移るからである。2つの式が同じ大きさを与えるのは、およそ μ e ∼ e a 2 / ℏ \mu_{\mathrm e}\sim ea^2/\hbar μ e ∼ e a 2 /ℏ のあたり(例題 4.8.1 )で、そこが2つの描像の境目になる。
化学の授業では、反応の速さを Arrhenius の式 k = A e − E a / k B T k=A\mathrm{e}^{-E_{\mathrm a}/ k_{\mathrm B} T} k = A e − E a / k B T や遷移状態理論の k = k B T h e − Δ G ‡ / k B T k=\frac{ k_{\mathrm B} T}{h}\mathrm{e}^{-\Delta G^\ddagger/ k_{\mathrm B} T} k = h k B T e − Δ G ‡ / k B T で習う。Marcus の式 (4.8.11) はその特別な場合で、活性化エネルギー E a = ( λ + Δ G ) 2 / 4 λ E_{\mathrm a}=(\lambda+\Delta G)^2/4\lambda E a = ( λ + Δ G ) 2 /4 λ を2つの放物線の交点として駆動力の関数に書き、前の因子を電子の結合 t h 2 t_{\mathrm h}^2 t h 2 で決めたものである。遷移状態理論の前の因子 k B T / h k_{\mathrm B} T/h k B T / h は断熱の極限(4.8.3 節の γ L Z ≳ 1 \gamma_{\mathrm{LZ}}\gtrsim1 γ LZ ≳ 1 )に、Marcus の t h 2 t_{\mathrm h}^2 t h 2 の因子は非断熱の極限に対応する。
有機 FET の式 (4.8.16) は、電子回路の授業で習う MOSFET の2乗則と同じ形である。授業では μ e \mu_{\mathrm e} μ e を物質の定数として与えるが、有機 FET ではそれが (4.8.14) の温度と構造の関数になっている。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com