ハミルトニアンの選び方で固体を考える 第2部 第8章 §2.8.4 ★★ 弱い結合で終わる端と、2つの Kekulé 構造の境目には、ギャップの中央 ε=0 に孤立した準位が1つずつ現れる。その数は h(k) の巻き数で決まる
端の準位
2量体化した鎖を有限の長さで切ると、端の結合が強い(t1)か弱い(t2)かで違いが出る。サイト 1,2,…,N(N は偶数)の鎖の結合を、弱い・強い・弱い・…・弱いの順に並べる(N 個のサイトに N−1 本の結合があり、両端が弱い結合になる)。このとき ε=0 の固有状態を探す。Schrödinger 方程式の n 番目の成分は −tn−1ψn−1−tnψn+1=εψn で、ε=0 では奇数番目のサイトの振幅と偶数番目のサイトの振幅が別々の式に分かれる。偶数番目をすべて 0 とし、奇数番目について解くと、サイト2の式 −twψ1−tsψ3=0 から
ψ2j+1=(−tstw)jψ1,ψ2j=0(2.8.13)
を得る(ts=t+∣y∣, tw=t−∣y∣ は強い結合と弱い結合のホッピング)。tw<ts なので振幅は左端から指数関数的に減り、左端に局在した状態になる。減る長さは、単位胞(長さ 2a)ごとに tw/ts 倍になることから
ξ=ln(ts/tw)2a(2.8.14)
である。右端にも同じ形の状態があり、鎖が長ければ2つの状態は e−Na/ξ の桁の小さな重なりで ± にわずかに分かれるだけで、どちらもギャップの中央に残る(図 2.8.3)。逆に端が強い結合で終わる鎖では、同じ式で振幅が端から増えてしまうので規格化できず、端の準位は無い。
この違いは直観的にも分かる。弱い結合で終わる鎖では、端のサイトは隣と弱くしかつながっておらず、強い結合の対を作れずに1個だけ余っている。余った1個の pz 軌道が、ギャップの中央 εa=0 の準位として残る。

図 2.8.3 両端が弱い結合で終わる N=40 の鎖(ts=2.83eV, tw=2.17eV)を数値で対角化し、ε≈0 の2つの固有状態から左端に局在した組み合わせを作ったもの(左半分だけ示す)。振幅は奇数番目のサイトにだけあり、式 (2.8.13) のとおり単位胞ごとに −tw/ts=−0.770 倍になる。2つの固有値は ±6.1meV で、バルクの準位の最小 ±0.79eV に比べて桁違いに 0 に近い。
ソリトン:2つの Kekulé 構造の境目
長い鎖の途中で、左側は u>0、右側は u<0 の2量体化になっている配置を考える。境目の近くでは、強い結合の対を作れないサイトが1個余る。左側から見るとそのサイトは弱い結合で終わる端、右側から見ても弱い結合で終わる端になっており、上と同じ理由で境目に局在した ε=0 の準位が1つ現れる。実際の鎖では un は境目で急に符号を変えず、ある幅にわたって滑らかに +u0 から −u0 へ移る。この境目をソリトンと呼ぶ。SSH は数値計算で、ソリトンの幅が約 7a、ソリトン1個を作るエネルギーが 0.42eV になることを示した(Su–Schrieffer–Heeger 1980)。連続極限の模型での解析では、ソリトン1個のエネルギーは 2Δ0/π になる(H. Takayama, Y. R. Lin-Liu and K. Maki (1980), Phys. Rev. B 21, 2388)。Δ0=0.650eV なら 0.414eV で、ギャップの半分 Δ0 より小さい。電子を1個ギャップの上へ上げるより、ソリトンを作る方がエネルギーが安い。
ギャップの中央の準位は、スピンを含めて席が2つある。電子が1個入るとソリトンは電荷 0・スピン 21、2個入ると電荷 −e・スピン 0、空なら電荷 +e・スピン 0 になる(Su–Schrieffer–Heeger 1980)。電荷とスピンの組み合わせが普通の電子(電荷 −e・スピン 21)と逆になっている。ポリアセチレンにドーピングで電荷を入れると、電子はバンドへ入らず、電荷をもちスピンをもたないソリトンになると考えられている。
巻き数:端の準位の数を h(k) から数える
端の準位があるかどうかは、無限に長い鎖の h(k) だけから決まる。式 (2.8.5) の非対角成分を
q(k):=−(t1+t2e−2ika)
と書くと、k が Brillouin 域を1周するとき、q(k) は複素平面の上で中心 −t1、半径 t2 の円を1周する。t2<t1 なら円は原点を囲まず、t2>t1 なら原点を1回囲む。原点のまわりを回る回数
W:=2πi1∫−π/(2a)π/(2a)dkdlnq(k)dk={0−1(t1>t2)(t1<t2)(2.8.15)
が、h(k) の巻き数である(e−2ika は k が増えると時計回りに回るので、原点を囲むときの符号は −1 になる)。t1=t2(ギャップが閉じる所)を通らない限り、t1, t2 を連続に変えても W は整数のまま変わらない。
単位胞を (A, B) で区切った鎖を左端で切ると、端の結合は単位胞の中の結合 t1 の反対側、すなわち左端の A サイトは単位胞の中の t1 だけで隣とつながる。t1<t2(∣W∣=1)なら端の結合が弱く、端の準位が1つある。t1>t2(W=0)なら端の結合が強く、端の準位は無い。端の準位の数が無限の鎖の巻き数 ∣W∣ に等しいというこの対応は、第5部で扱うトポロジカル絶縁体の「バルクと端の対応」の最も簡単な例である(第5部第6章)。巻き数を数えるとき、h(k) の対角成分が 0 であること(A どうし・B どうしの飛び移りが無いこと)を使った。この性質(カイラル対称性)が崩れると、h(k) は2次元の円でなく3次元の球面の上を動くので、原点を囲むという言い方が使えなくなる。
例題 2.8.4 端の準位の局在の長さとソリトンの幅を比べる ★★
オリジナル式 (2.8.2) の値と 2.8.3節の y0=0.326eV を使う。
- (1) ts, tw と、式 (2.8.14) の局在の長さ ξ を a の単位で求めよ。
- (2) SSH が数値計算で得たソリトンの幅 7a と比べよ。
- (3) (1) の値で、ξ を連続極限の式 ξ≈ℏvF/Δ0(ℏvF=2ta は等間隔の鎖の Fermi 準位での dε/dk)と比べよ。
解答
(1) ts=2.826eV、tw=2.174eV、ln(ts/tw)=ln1.2999=0.2623 なので ξ=2a/0.2623=7.63a である。
(2) 端の準位の局在の長さ 7.6a は、ソリトンの幅 7a と同じ大きさである。どちらも「ギャップの中央の状態が、ギャップ Δ0 の絶縁体の中にどこまでしみ出すか」で決まる長さだからである。ソリトンの幅が格子間隔の数倍あるので、格子の離散性を無視した連続極限の模型(Takayama–Lin-Liu–Maki)がよい近似になる。
(3) ℏvF=2×2.5×a=5.0eV⋅a、Δ0=0.652eV として ξ≈5.0/0.652a=7.67a で、(1) と 0.5% で一致する。ts/tw=(1+z)/(1−z) で z=y/t が小さいとき ln1−z1+z≈2z なので、ξ≈2a/(2z)=at/y=2ta/(2y)=ℏvF/Δ0 となり、2つの式は z→0 で一致する。
ほかの巻では
- 近似の数理 第12章:端の準位の数と Wannier 関数の中心のずれ(Zak 位相)が同じ巻き数から来ることが読める。記号の対応:向こうの v,w はこの巻の t1,t2
- 相対エントロピーから情報と物理を組み立てる 第5部第1章:2量体化した鎖の基底状態のエンタングルメント・エントロピーが長さによらず飽和することが読める。記号の対応:向こうの δ はこの巻の y/t
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com