tecktom.com

この章のまとめ

2026年10月15日10分で読めますAI生成教科書『虎の巻』第2部 第8章 ポリアセチレン:Peierls 転移と SSH 模型

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第2部 第8章 §2.8.7 ★★★ SSH 模型の全エネルギーを歪みの大きさについて最小にすると、半分詰まった1次元の鎖は必ず2量体化して絶縁体になり、端と境目にギャップの中央の準位が残る

この章では、原子の位置 unu_n を数として H^\hat H に戻した SSH 模型 (2.8.1) を選んだ。ホッピングは結合の長さの1次で変わり、σ\sigma 骨格はばねとして入る。1つおきに逆向きにずれる配置 (2.8.3) では単位胞が原子2個になり、Bloch の定理で kk ごとの 2×22\times2 の行列 (2.8.5) に分かれる。そのバンド (2.8.6) は、等間隔の鎖のバンドを半分の Brillouin 域へ折り返し、ゾーンの端(等間隔の鎖の Fermi 波数)に幅 4∣y∣4\lvert y\rvert のギャップを開いたものである。

全エネルギーは楕円積分 (2.8.8) で書け、小さな歪みでは電子の得が y2ln⁡(1/y)y^2\ln(1/y)、ばねの損が y2y^2 に比例する。1次元の Fermi 面が2点だけで、歪みの波数 2kF2k_{\mathrm F} が Fermi 面全体を結ぶことが対数の起源で、等間隔の鎖は必ず不安定になる(Peierls 転移)。ギャップは Δ0≈8t e−1−1/(2λ)\Delta_0\approx8t\,\mathrm{e}^{-1-1/(2\lambda)} で、ポリアセチレンの値では u0=0.0397 A˚u_0=0.0397\,\mathrm{\AA}、Eg=1.30 eVE_{\mathrm g}=1.30\,\mathrm{eV} になる。2量体化した鎖の端が弱い結合で終わるとき、また2つの Kekulé 構造の境目(ソリトン)には、ギャップの中央に局在した準位が1つ現れ、その有無は h(k)h(k) の巻き数 (2.8.15) で決まる。

削ったものは、格子の量子揺らぎ(静止した格子の範囲では2量体化が厳密に最小であることが証明されている)、π\pi 電子以外と遠くへの飛び移り、電子どうしの相互作用、鎖どうしの結合の4つである。模型は結合の長さの交替と光学ギャップを桁でとらえるが、結合の長さの差は実験の約 1.71.7 倍になる。

次の章(第2部第9章)では、2次元のグラフェンで、同じ 2×22\times2 の行列(単位胞に炭素2個)からゾーンの角でギャップの閉じた Dirac コーンが現れることを見る。グラフェンの h(k)h(\bm{k}) も対角成分が 00 の形をもち、この章の巻き数の考え方がそのまま Dirac 点のまわりの位相として現れる。

この章の要点

  • SSH 模型:H^=−∑[t−α(un+1−un)](c^n+1†c^n+h.c.)+K2∑(un+1−un)2\hat H=-\sum[t-\alpha(u_{n+1}-u_n)](\hat c^\dagger_{n+1}\hat c_n+\mathrm{h.c.})+\frac K2\sum(u_{n+1}-u_n)^2。2量体化 un=(−1)nuu_n=(-1)^nu で t1,2=t±yt_{1,2}=t\pm y、y=2αuy=2\alpha u。
  • バンド:ε±(k)=±2t2cos⁡2ka+y2sin⁡2ka\varepsilon _\pm(k)=\pm2\sqrt{t^2\cos^2ka+y^2\sin^2ka}、ギャップ Eg=4∣y∣=2Δ0E_{\mathrm g}=4\lvert y\rvert =2\Delta_0。
  • 全エネルギー:E/N=−4tπE(1−z2)+2Ku2E/N=-\frac{4t}\pi\mathrm E(1-z^2)+2Ku^2。電子の得 ∝z2ln⁡(4/z)\propto z^2\ln(4/z) がばねの損 ∝z2\propto z^2 に勝つ。Δ0≈8t e−1−1/(2λ)\Delta_0\approx8t\,\mathrm{e}^{-1-1/(2\lambda)}、λ=2α2/(πKt)\lambda=2\alpha^2/(\pi Kt)。
  • 端の準位:弱い結合で終わる端に ψ2j+1∝(−tw/ts)j\psi_{2j+1}\propto(-t_{\mathrm w}/t_{\mathrm s})^j、局在の長さ ξ=2a/ln⁡(ts/tw)≈ℏvF/Δ0\xi=2a/\ln(t_{\mathrm s}/t_{\mathrm w})\approx\hbar v_{\mathrm F}/\Delta_0。数は巻き数 ∣W∣\lvert W\rvert。
  • ソリトン:2つの Kekulé 構造の境目。電荷とスピンの組が電子と逆。エネルギー ≈2Δ0/π\approx2\Delta_0/\pi。

総合問題 炭素の輪はいつ2量体化するか:ベンゼンとポリアセチレンをつなぐ ★★★

オリジナルサイト NN 個(NN は偶数)の輪に π\pi 電子 NN 個を入れ、式 (2.8.1) の SSH 模型(周期境界)を考える。2量体化 un=(−1)nuu_n=(-1)^nu を仮定し、y=2αuy=2\alpha u, z=y/tz=y/t とする。ばねのエネルギーは輪全体で 2NKu2=NKy2/(2α2)2NKu^2=NKy^2/(2\alpha^2) である。

  1. (1) (使う道具:第2部第3章・第2部第2章)N=6N=6(ベンゼン)で y=0y=0 のときの1電子の準位を −2tcos⁡(2πν/6)-2t\cos(2\pi\nu/6) から求め、π\pi 電子のエネルギー −8t-8t を確かめよ。
  2. (2) (使う道具:2.8.2節の 2×22\times2 の行列)N=6N=6 で単位胞は3個、2ka∈{0,±2π/3}2ka\in\{0,\pm2\pi/3\} である。各 kk の固有値 ±∣t1+t2e−2ika∣\pm\lvert t_1+t_2\mathrm{e}^{-2\mathrm{i} ka}\rvert を求め、電子のエネルギーが −4t−4t2+3y2-4t-4\sqrt{t^2+3y^2} になることを示せ。
  3. (3) (使う道具:2.8.3節の得と損の比較)N=6N=6 の全エネルギーの変化を zz の2次まで展開し、ベンゼンが2量体化する条件を λ=2α2/(πKt)\lambda=2\alpha^2/(\pi Kt) で書け。式 (2.8.2) の値ではどうなるか。さらに、1+3z2−1≤32z2\sqrt{1+3z^2}-1\le\frac32z^2 を使って、この値ではどの大きさの zz でも2量体化しないことを示せ。
  4. (4) (使う道具:2.8.2節・例題 2.8.2)N=4N=4 では 2ka∈{0,π}2ka\in\{0,\pi\} である。電子のエネルギーが −4t−4∣y∣-4t-4\lvert y\rvert になることを示し、輪は λ\lambda によらず必ず歪むことと、歪みの大きさ ∣y∣=α2/K\lvert y\rvert =\alpha^2/K を求めよ。
  5. (5) (使う道具:2.8.3節・第2部第3章)N=4m+2N=4m+2 の等間隔の輪の HOMO と LUMO の差は 4tsin⁡(π/N)4t\sin(\pi/N) である。これが無限の鎖のギャップ 2Δ0=1.30 eV2\Delta_0=1.30\,\mathrm{eV} より小さくなると輪は2量体化する、と見積もって、2量体化する最小の NN の目安を求めよ。数値で対角化して全エネルギーを最小にすると、4m+24m+2 の輪で最初に2量体化するのは N=18N=18 である。見積もりと比べよ。
解答

(1) 準位は ν=0\nu=0 で −2t-2t、ν=1,5\nu=1,5 で −t-t、ν=2,4\nu=2,4 で tt、ν=3\nu=3 で 2t2t である。電子6個は下の3準位にスピンを逆にして2個ずつ入り、エネルギーは 2(−2t−t−t)=−8t2(-2t-t-t)=-8t である。

(2) 2ka=02ka=0 では ∣t1+t2∣=2t\lvert t_1+t_2\rvert =2t。2ka=±2π/32ka=\pm2\pi/3 では ∣t1+t2e∓2πi/3∣2=t12+t22+2t1t2cos⁡2π3=t12+t22−t1t2\lvert t_1+t_2\mathrm{e}^{\mp2\pi\mathrm{i}/3}\rvert ^2=t_1^2+t_2^2+2t_1t_2\cos\frac{2\pi}3=t_1^2+t_2^2-t_1t_2 で、t12+t22=2t2+2y2t_1^2+t_2^2=2t^2+2y^2, t1t2=t2−y2t_1t_2=t^2-y^2 を入れると t2+3y2t^2+3y^2 になる。下の3準位 −2t-2t, −t2+3y2-\sqrt{t^2+3y^2}(2つ)に2個ずつ入れて、電子のエネルギーは −4t−4t2+3y2-4t-4\sqrt{t^2+3y^2} である。y=0y=0 で −8t-8t に戻る。

(3) 全エネルギーの変化は

ΔE=−4t(1+3z2−1)+3Kt2α2z2≈(−6t+3Kt2α2)z2\Delta E=-4t\bigl(\sqrt{1+3z^2}-1\bigr)+\frac{3Kt^2}{\alpha^2}z^2\approx\Bigl(-6t+\frac{3Kt^2}{\alpha^2}\Bigr)z^2

である(ばねは 6Ky2/(2α2)=3Kt2z2/α26Ky^2/(2\alpha^2)=3Kt^2z^2/\alpha^2)。z2z^2 の係数が負になる条件は 2α2/(Kt)>12\alpha^2/(Kt)>1、すなわち λ>1/π=0.318\lambda>1/\pi=0.318 である。式 (2.8.2) では 2α2/(Kt)=33.62/52.5=0.640<12\alpha^2/(Kt)=33.62/52.5=0.640<1(λ=0.204\lambda=0.204)なので、小さな2量体化は全エネルギーを上げる。さらに 1+x≤1+x/2\sqrt{1+x}\le1+x/2 から 1+3z2−1≤32z2\sqrt{1+3z^2}-1\le\frac32z^2 なので、ΔE≥(−6t+3Kt2/α2)z2\Delta E\ge(-6t+3Kt^2/\alpha^2)z^2 がすべての zz で成り立つ。係数は −15+23.42=8.42 eV>0-15+23.42=8.42\,\mathrm{eV}>0 なので、z≠0z\neq0 ではつねに ΔE>0\Delta E>0 である。ポリアセチレンと同じ tt, α\alpha, KK を使っても、ベンゼンは6本の結合の長さがそろった形にとどまる。等間隔の輪ですでに HOMO(−t-t)と LUMO(tt)の間に 2t2t の差があり、対数の得が生じないからである。

(4) 2ka=02ka=0 で ±2t\pm2t、2ka=π2ka=\pi で ±∣t1−t2∣=±2∣y∣\pm\lvert t_1-t_2\rvert =\pm2\lvert y\rvert である。電子4個は −2t-2t と −2∣y∣-2\lvert y\rvert に2個ずつ入り、電子のエネルギーは −4t−4∣y∣-4t-4\lvert y\rvert になる。ばねは 4Ky2/(2α2)=2Ky2/α24Ky^2/(2\alpha^2)=2Ky^2/\alpha^2 なので ΔE=−4∣y∣+2Ky2/α2\Delta E=-4\lvert y\rvert +2Ky^2/\alpha^2。∣y∣\lvert y\rvert が小さいとき1次の得が2次の損に必ず勝つので、λ\lambda によらず歪む。dΔE/d∣y∣=−4+4K∣y∣/α2=0\mathrm{d}\Delta E/\mathrm{d}\lvert y\rvert =-4+4K\lvert y\rvert /\alpha^2=0 から ∣y∣=α2/K=16.81/21=0.800 eV\lvert y\rvert =\alpha^2/K=16.81/21=0.800\,\mathrm{eV}、∣u∣=∣y∣/(2α)=0.0976 A˚\lvert u\rvert =\lvert y\rvert /(2\alpha)=0.0976\,\mathrm{\AA} である。y=0y=0 で 00 の準位が2つ縮退して電子2個が半分だけ埋めていたのが、歪みで分かれる。Jahn–Teller 効果と同じ仕組みである。この歪みは 2.8.3節の u0u_0 の2倍以上あり、ホッピングを uu の1次で打ち切った近似の範囲を超えかけている。

(5) 4tsin⁡(π/N)<2Δ04t\sin(\pi/N)<2\Delta_0 から sin⁡(π/N)<1.30/10=0.130\sin(\pi/N)<1.30/10=0.130、π/N<0.1304\pi/N<0.1304、N>24.1N>24.1 である。4m+24m+2 の形で最小のものは N=26N=26 になる。数値の答え N=18N=18 はこの見積もりより小さい。見積もりは「対数の得が始まるのは、等間隔の輪の差が鎖のギャップより小さくなってから」という粗い判定で、実際には差がギャップと同じ桁(N=18N=18 で 4tsin⁡(π/18)=1.74 eV4t\sin(\pi/18)=1.74\,\mathrm{eV})でも、Fermi 準位から離れた多数の準位の2次の得の和が損を上回る。それでも、N=6N=6 の輪(ベンゼン)が等間隔で、数十個の炭素の輪から先が2量体化するという大きさの区分は見積もりでとらえられている。N→∞N\to\infty では輪のサイト1個あたりのエネルギーは無限の鎖の式 (2.8.8) に近づき、4m4m と 4m+24m+2 の区別は消える。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ