tecktom.com

何を削ったか

2026年10月15日8分で読めますAI生成教科書『虎の巻』第2部 第8章 ポリアセチレン:Peierls 転移と SSH 模型
目次
  1. 格子の量子的な揺らぎ
  2. 電子だけ・隣だけ・1次だけ
  3. 電子どうしの相互作用
  4. 鎖どうしの結合と温度
  5. 木下の層

AI 生成・木下が未勉強

ハミルトニアンの選び方で固体を考える 第2部 第8章 §2.8.5 ★★★ 格子を古典的に止め、π\pi 電子どうしの相互作用と鎖どうしの結合を捨てた。2量体化そのものはこの模型の中で厳密に証明されている

格子の量子的な揺らぎ

式 (2.8.1) では unu_n を数として扱い、原子核の運動エネルギーを捨てた。原子核が電子より十分重いという Born–Oppenheimer の極限である。実際の CH の組は零点振動をもち、その振幅は u0=0.04 A˚u_0=0.04\,\mathrm{\AA} と同じ桁になりうる。零点振動を入れると、2量体化の大きさは小さくなるが消えはしないと考えられている。この章の模型の中で、静止した格子・相互作用しない電子・半分詰まった輪について、全エネルギーを最小にする配置がちょうど2量体化(周期 22)になることはT. Kennedy and E. H. Lieb (1987), Phys. Rev. Lett. 59, 1309 が証明した。これは「周期 22 より複雑な歪みが、さらにエネルギーを下げることはない」という主張まで含む。サイト数が 4m+24m+2 の輪では、ホッピングと歪みのエネルギーの関数形によらず、最小は等間隔か2量体化のどちらかになることも示されている(E. H. Lieb and B. Nachtergaele (1995), Phys. Rev. B 51, 4777)。

近似の裏付け

型 B(極限) | 小さいもの (m/M)1/4→0(m/M)^{1/4}\to0(格子を古典的に止める。格子の零点振動を捨てる) | 証明 条件付き(静止した格子・相互作用しない電子・ホッピングが結合の長さの1次・調和なばねの範囲で、半分詰まった輪の最小が2量体化であることは Kennedy–Lieb 1987 が証明。格子の量子揺らぎを入れたときの2量体化の存続は要確認) | 行き先 近似の数理 第7章・近似の数理 第13章

π\pi 電子だけ・隣だけ・1次だけ

Hückel 法と同じく π\pi 電子だけを残し、σ\sigma 電子の役割はばね定数 KK にまとめた。ホッピングは隣どうしだけで、結合の長さの1次で打ち切った。σ\sigma 骨格のギャップ(10 eV10\,\mathrm{eV} の桁)が π\pi 電子と格子のエネルギー(1 eV1\,\mathrm{eV} 以下)より大きいことが、σ\sigma 電子を凍らせる根拠である。

近似の裏付け

型 A(射影) | 小さいもの 各炭素の pzp_z 軌道1つが張る部分空間(隣どうしのホッピングだけを残し、結合の長さの1次で打ち切る) | 証明 条件付き(深い井戸の極限でトンネル積分の形まで証明済み。実在の炭素の間隔での誤差は要確認) | 行き先 近似の数理 第4章・近似の数理 第12章

電子どうしの相互作用

式 (2.8.1) には π\pi 電子どうしの Coulomb 相互作用が無い。同じ炭素の上の斥力 UU(Hubbard 模型。第1部第1章)を足すと、半分詰まった鎖は格子が歪まなくても Mott 絶縁体になりうるので、ギャップの起源が格子の歪みか電子の相関かを分けて考える必要が出る。ポリアセチレンの UU は tt と同じ桁と見積もられており、この章の数値は相互作用を平均に押し込めた有効なパラメータとして読むべきものである。UU を入れた Peierls–Hubbard 模型でも、4m+24m+2 個の輪の最小が等間隔か2量体化のどちらかであることは Lieb–Nachtergaele 1995 が示しているが、どちらになるかは UU と λ\lambda の大きさによる。

近似の裏付け

型 B(極限) | 小さいもの U/t→0U/t\to0(π\pi 電子どうしの相互作用を捨てる) | 証明 要確認(U/tU/t は 1 の桁で小さくない。相互作用を入れた模型で2量体化が残るかは数値計算に頼る部分が多く、この章では前提を原論文で確かめていない) | 行き先 近似の数理 第13章

鎖どうしの結合と温度

ポリアセチレンの試料は多数の鎖の束で、隣の鎖との間にも弱い飛び移りがある。この章ではそれを捨て、1本の鎖だけを扱った。有限の温度では、1本の鎖の2量体化の型(u>0u>0 か u<0u<0)は熱で生じたソリトンによって区切られ、長い距離では乱れる。2つの基底状態のどちらかを選ぶ型の秩序は、短い距離の相互作用だけをもつ1次元の系では有限の温度で長い距離まで続かない。平均場の扱い(すべての鎖を1つの uu で表す)では2量体化が消える温度 TMFT_{\mathrm{MF}} が求まるが、これは本当の相転移の温度ではなく、ソリトンを作る費用と比べて意味を読む必要がある(例題 2.8.5)。

近似の裏付け

型 B(極限) | 小さいもの 鎖どうしの飛び移り →0\to0(1本の鎖を扱い、平均場で温度を入れる) | 証明 要確認(1次元の系で有限の温度の長距離秩序が無いことと、鎖どうしの結合で3次元の秩序が回復する条件は、この章では原論文で確かめていない) | 行き先 近似の数理 第6章

例題 2.8.5 平均場の転移温度と熱で生じるソリトンの間隔 ★★

オリジナル有限の温度では、ギャップの中央を化学ポテンシャルとし、サイト1個あたりの自由エネルギーを

F(y,T)N=−2kBT⋅1π∫−π/2π/2ln⁡[2cosh⁡ε+(ϕ)2kBT]dϕ+Ky22α2\frac{F(y,T)}{N}=-2 k_{\mathrm B} T\cdot\frac1\pi\int_{-\pi/2}^{\pi/2}\ln\Bigl[2\cosh\frac{\varepsilon _+(\phi)}{2 k_{\mathrm B} T}\Bigr]\mathrm{d}\phi+\frac{K y^2}{2\alpha^2}

とする(ε+\varepsilon _+ は式 (2.8.6)。準位 ±ε\pm\varepsilon の対の大きな正準ポテンシャル −kBTln⁡[(1+eβε)(1+e−βε)]=−2kBTln⁡[2cosh⁡(βε/2)]- k_{\mathrm B} T\ln\bigl[(1+\mathrm{e}^{\beta\varepsilon })(1+\mathrm{e}^{-\beta\varepsilon })\bigr]=-2 k_{\mathrm B} T\ln[2\cosh(\beta\varepsilon /2)] にスピンの2を掛けた)。

  1. (1) T→0T\to0 で F/NF/N が式 (2.8.8) に戻ることを確かめよ。
  2. (2) 式 (2.8.2) の値でこの FF を数値で最小にすると、yy が 00 になる温度は TMF≈4180 KT_{\mathrm{MF}}\approx4180\,\mathrm K である。Δ0/(kBTMF)\Delta_0/( k_{\mathrm B} T_{\mathrm{MF}}) を求め、弱い結合の極限の値 πe−γ=1.764\pi\mathrm{e}^{-\gamma}=1.764(γ\gamma は Euler の定数。BCS と同じ値)と比べよ。
  3. (3) ソリトン1個のエネルギーを Es=0.42 eVE_{\mathrm s}=0.42\,\mathrm{eV} とする。300 K300\,\mathrm K で熱で生じるソリトンの間隔は、おおよそ a eEs/(kBT)a\,\mathrm{e}^{E_{\mathrm s}/( k_{\mathrm B} T)} である。これを長さで求めよ(a=1.22 A˚a=1.22\,\mathrm{\AA})。
解答

(1) βε≫1\beta\varepsilon \gg1 で ln⁡[2cosh⁡(βε/2)]→βε/2\ln[2\cosh(\beta\varepsilon /2)]\to\beta\varepsilon /2 なので、第1項は −2kBT⋅1π∫ε+2kBTdϕ=−1π∫ε+ dϕ=−4tπE(1−z2)-2 k_{\mathrm B} T\cdot\frac1\pi\int\frac{\varepsilon _+}{2 k_{\mathrm B} T}\mathrm{d}\phi=-\frac1\pi\int\varepsilon _+\,\mathrm{d}\phi=-\frac{4t}{\pi}\mathrm E(1-z^2) になる。第2項は Ky2/(2α2)=K(2αu)2/(2α2)=2Ku2Ky^2/(2\alpha^2)=K(2\alpha u)^2/(2\alpha^2)=2Ku^2 で、式 (2.8.8) と一致する。

(2) kBTMF=8.617×10−5×4183=0.360 eVk_{\mathrm B} T_{\mathrm{MF}}=8.617\times10^{-5}\times4183=0.360\,\mathrm{eV} なので、Δ0/(kBTMF)=0.650/0.360=1.80\Delta_0/( k_{\mathrm B} T_{\mathrm{MF}})=0.650/0.360=1.80 である。弱い結合の値 1.7641.764 より約 2%2\% 大きく、λ=0.204\lambda=0.204 がまだ弱い結合の側にあることを示している。TMF≈4000 KT_{\mathrm{MF}}\approx4000\,\mathrm K は高分子が分解する温度より桁で高い。平均場で見る限り、ポリアセチレンは存在しうるどの温度でも2量体化している。

(3) kBT=0.02585 eVk_{\mathrm B} T=0.02585\,\mathrm{eV} で Es/(kBT)=16.2E_{\mathrm s}/( k_{\mathrm B} T)=16.2、e16.2≈1.1×107\mathrm{e}^{16.2}\approx1.1\times10^{7} なので、間隔は 1.1×107×1.22×10−10 m≈1.4 mm1.1\times10^7\times1.22\times10^{-10}\,\mathrm m\approx1.4\,\mathrm{mm} である。1本の鎖が途切れずに続く長さはこれよりはるかに短いので、室温の鎖の上に熱で生じたソリトンはほとんど無く、2量体化の型は鎖ごとにそろっている。1次元では無限に長い鎖の秩序は有限の温度で必ず壊れるが、壊れる長さが試料の鎖の長さを桁で超えるので、実際の鎖ではその効果は現れない。

木下の勉強ステータス

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

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

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